Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_013198	985,365	A→G	I904I (ATT→ATC)	pacL	cation‑transporting ATPase

Read alignment evidence...
	seq id	position		change	freq	score	reads	annotation	genes	product
*	NC_013198	985,365	0	A→G	100.0%	241.4	62	I904I (ATT→ATC)	pacL	cation‑transporting ATPase
Reads supporting (aligned to +/- strand):  new base (20/42):  ref base (0/0):  total (20/42)

AAAGTGACACTAAAAAATGGCCTCGGAAAGTTCAAATTCCAAGGCCATTTTCTTAGTGTTGCTTATTCGACACGCTTTACCTTCGCATGGGCTAAGCGCCGACGATTGAGCCATTTCCGTCCTTCTTCGATTAAGAAGAGTGGAATCGGGATGATGGCTAGGAAGATCCAATCGGTGAGGGCTAGTGGGCCCGTGTGGAAGAGTTCCTGTAAGAACGGGGTGTAACTGAGTAGGGCCAGCAGGAAAATTTCAAAGATAATGCCGAACCAGACAAGTCGATTGGCAAACAA  .  NC_013198/985222‑985511                                                                                                                                                |                                                                                                                                                   aaaGTGACACTAAAAAATGGCCTCGGAAAGTTCAAATTCCAAGGCCATTTTCTTAGTGTTGCTTATTCGACACGCTTTACCTTCGCATGGGCTAAGCGCCGACGATTGAGCCATTTCCGTCCTTCTTCGATTAAGAAGAGTGGGATCGGga                                                                                                                                             <  1:567370/151‑1       acacTAAAAAATGGCCTCGGAAAGTTCAAATTCCAAGGCCATTTTCTTAGTGTTGCTTATTCGACACGCTTTACCTTCGCATGGGCTAAGCGCCGACGATTGAGCCATTTCCGTCCTTCTTCGATTAAGAAGAGTGGGATCGGGATGATgg                                                                                                                                       <  2:495971/151‑1             aaaacTGACCTGGGAAAGATCAAATCCCCAGGCCATTCTGTTAATGTTGCTTATTCGACACGTTTTGCCTTTTTATGGGCAACGCGCTGAAGATTGAGCCATTTCCGTCCGTCTTCGATTAAGACGAGTGCGATCGCGATGAGGGATAgca                                                                                                                                 <  2:96500/151‑3             aaaaaTGGCCTCGGAAAGTTCAAATTCCAAGGCCATTTTCTTAGTGTTGCTTATTCGACACGCTTTACCTTCGCATGGGCTAAGCGCCGACGATTGAGCCATTTCCGTCCTTCTTCGATTAAGAAGAGTGGGATCGGGATGATGGCTAGGa                                                                                                                                 <  2:96448/151‑1              aaaaTGGCCTCGGAAAGTTCAAATTCCAAGGCCATTTTCTTAGTGTTGCTTATTCGACACGCTTTACCTTCGCATGGGCTAAGCGCCGACGATTGAGCCATTTCCGTCCTTCTTCGATTAAGAAGAGTGGGATCGGGATGATGGCTAGGaa                                                                                                                                <  1:76612/151‑1                    gCCTCGGAAAGTTCAAATTCCAAGGCCATTTTCTTAGTGTTGCTTATTCGACACGCTTTACCTTCGCATGGGCTAAGCGCCGACGATTGAGCCATTTCCGTCCTTCTTCGATTAAGAAGAGTGGGATCGGGATGATGGCTAGGAAGATCCa                                                                                                                          >  1:89323/1‑151                       tCGGAAAGTTCAAATTCCAAGGCCATTTTCTTAGTGTTGCTTATTCGACACGCTTTTCCTTCGCATGGGCTAAGCGCCGACGATTGAGCCATTTCCGTCCTTCTTCGATTAAGAAGAGTGGGATCGGGATGATGGCTAGGAAGATCCAATc                                                                                                                       <  1:829703/151‑1                           aaaGTTCAAATTCCAAGGCCATTTTCTTAGTGTTGCTTATTCGACACGCTTTACCTTCGCATGGGCTAAGCGCCGACGATTGAGCCATTTCCGTCCTTCTTCGATTAAGAAGAGTGGGATCGGGATGATGGCTAGGAAGATCCAATCGGTg                                                                                                                   <  1:949958/151‑1                              attCAAATTCCAAGGCCATTTTCTTAGTGTTGCTTATTCCACACGCTTTACCTTCGAATGGGCTAAGCGCCGACGATTGAGCCATTTCCGTCCTTCTTCAATTAAGAAGCGTGGGATCAGGCTGCTGGCTAGGAAGATCCAATCggagcgg                                                                                                                >  2:121148/2‑146                               ttCAAATTCCAAGGCCATTTTCTTAGTGTTGCTTATTCGACACGCTTTACCTTCGCATGGGCTAAGCGCCGACGATTGAGCCATTTCCGTCCTTCTTCGATTAAGAAGAGTGGGATCGGGATGATGGCTAGGAAGATCCAATCGGTGAggg                                                                                                               <  1:299168/151‑1                                 cAAATTCCAAGGCCATTTTCTTAGTGTTGCTTATTCGACACGCTTTACCTTCGCATGGGCTAAGCGCCGACGATTGAGCCATTTCCGTCCTTCTTCGATTAAGAAGAGTGGGATCGGGATGATGGCTAGGAAGATCCAATCGGTGAGGGCt                                                                                                             <  1:319425/151‑1                                     ttCCAAGGCCATTTTCTTAGTGTTGGTTATTCGGCACGCTTTACCTTCGCATGGGCTAAGCGCCGACGATTGAGCCATTTCCGTCCTTCTTCGATTAAGAAGAGTGGGATCGGGATGATGGCTAGGAAGATCCAATCGGTGAGGGCTAGTg                                                                                                         <  2:878627/151‑1                                     ttCCAAGGCCATTTTCTTAGTGTTGCTTATTCGACACGCTTTACCTTCGCATGGGCTAAGCGCCGACGATTGAGCCATTTCCGTCCTTCTTCGATTAAGAAGAGTGGGATCGGGATGATGGCTAGGAAGATCCAATCGGTGAGGGCTAGTg                                                                                                         <  2:634082/151‑1                                     ttCCAAGGCCATTTTCTTAGTGTTGCTTATTCGACACGCTTTACCTTCGCATGGGCTAAGCGCCGACGATTGAGCCATTTCCGTCCTTCTTCGATTAAGAAGAGTGGGATCGGGATGATGGCTAGGAAGATCCAATCGGTGAGGGCTAGTg                                                                                                         <  1:723101/151‑1                                     ttCCAAGGCCATTTTCTTAGTGTTGCTTATTCGACACGCTTTACCTTCGCATGGGCTAAGCGCCGACGATTGAGCCATTTCCGTCCTTCTTCGATTAAGAAGAGTGGGATCGGGATGATGGCTAGGAAGATCCAATCGGTGAGGGCTAGTg                                                                                                         <  2:358037/151‑1                                      tCCAAGGCCATTTTCTTAGTGTTGCTTATTCGACACGCTTTACCTTCGCATGGGCTAAGCGCCGACGATTGAGCCATTTCCGTCCTTCTTCGATTAAGAAGAGTGGGATCGGGATGATGGCTAGGAAGATCCAATCGGTGAGGGCTAGTgg                                                                                                        <  1:725788/151‑1                                         aaGGCCATTTTCTTAGTGTTGCTTATTCGACACGCTTTACCTTCGCATGGGCTAAGCGCCGACGATTGAGCCATTTCCGTCCTTCTTCGATTAAGAAGAGGGGGATCGGGATGATGGCTAGGAAGATCCAATCGGTGAGGGCTAGTGGGcc                                                                                                     <  1:307617/151‑1                                          aGGCCATTTTCTGAGTGCTTGTTATTCGACACGCTTTGGCTTCGTGCGGGCTATGTGGCGCCCATTCAGCAATTTCCGTGCGGTTGCGATTGAGAAAGGAGGGATCTGGATAATGGCTAGGAAGATCCTATCGGTGATGGCTAGTgggtcc                                                                                                    <  2:281778/151‑4                                            gCCATTTTCTTAGTGTTGCTTATTCGACACGCTTTACCTTCGCATGGGCTAAGCGCCGACGATTGAGCCATTTCCGTCCTTCTTCGATTAAGAAGAGTGGGATCGGGATGATGGCTAGGAAGATCCAATCGGTGAGGGCTAGTGGGCCCgt                                                                                                  >  2:448643/1‑151                                            gCCATTTTCTTAGTGTTGCTTATTCGACACGCTTTACCTTCGCATGGGCTAAGCGCCGACGATTGAGCCATTTCCGTCCTTCTTCGATTAAGAAGAGTGGGATCGGGATGATGGCTAGGAAGATCCAATCGGTGAGGGCTAGTGGGCCCgt                                                                                                  >  1:898800/1‑151                                            gCCATTTTCTTAGTGTTGCTTATTCGACACGCTTTACCTTCGCAAGGGCTAAGCGCCGACGATTGAGCCATTTCCGTCCTTCTTCGATTAAGAAGAGTGGGATCGGGATGATGGCTAGGAAGATCCAATCGGTGAGGGCTAGTGGGCCCgt                                                                                                  >  1:714791/1‑151                                                   tCTTAGTGTTGCTTATTCGACACGCTTTACCTTCGCATGGGCTAAGCGCCGACGATTGAGCCATTTCCGTCCTTCTTCGATTAAGAAGAGTGGGATCGGGATGATGGCTAGGAAGATCCAATCGGTGAGGGCTAGTGGGCCCGTGTGGaan                                                                                           <  1:811909/151‑2                                                   tCTTAGTGTTGCTTATTCGACACGCTTTACCTTCGCATGGGCTAAGCGCCGACGATTGAGCCATTTCCGTCCTTCTTCGATTAAGAAGAGTGGGATCGGGATGATGGCTAGGAAGATCCAATCGGTGAGGGCTAGTGGGCCCGTGTGGAag                                                                                           <  1:636094/151‑1                                                   tCTTAGTGTTGCTTATTCGACACGCTTTACCTTCGCATGGGCTAAGCGCCGACGATTGAGCCATTTCCGTCCTTCTTCGATTAAGAAGAGTGGGATCGGGATGATGGCTAGGAAGATCCAATCGGTGAGGGCTAGTGGGCCCGTGTGGAag                                                                                           <  1:217296/151‑1                                                   tCTTAGTGTTGCTTATTCGACACGCTTTACCTTCGCATGGGCTAAGCGCCGACGATTGAGCCATTTCCGTCCTTCTTCGATTAAGAAGAGTGGGATCGGGATGATGGCTAGGAAGATCCAATCGGTGAGGGCTAGTGGGCCCGTGTGGAag                                                                                           <  1:762166/151‑1                                                   tCTTAGTGTTGCTTATTCGACACGCTTTACCTTCGCATGGGCTAAGCGCCGACGATTGAGCCATTTCCGTCCTTCTTCGATTAAGAAGAGTGGGATCGGGATGATGGCTAGGAAGATCCAATCGGTGAGGGCTAGTGGGCCCGTGTGGAag                                                                                           <  2:378681/151‑1                                                     ttAGTGTTGCTTATTCGACACGCTTTACCTTCGCATGGGCTAAGCGCCGACGATTGAGCCATTTCCGTCCTTCTTCGATTAAGAAGAGTGGGATCGGGATGATGGCTAGGAAGATCCAATCGGTGAGGGCTAGTGGGCCCGTGTGGAagag                                                                                         <  2:730834/151‑1                                                     ttAGTGTTGCTGATTCGACACGCTTTACCTTCGCATGGGCTAAGCGCCGACGATTGAGCCATTTCCGTCCTTCTTCGATTAAGAAGAGTGGGATCGGGATGATGGCTAGGAAGATCCAATCGGTGAGGGCTAGTGGGCCCGTGTGGAagag                                                                                         <  2:829298/151‑1                                                     gtAGTGTTGCTTATTCGACGCGCTGTTCCTTCGCATGGGCTAAGCTCCGACGATTGAGCCATTTCCGTCCTTCTTCGATTAAGAAGAGTGGGATCGGGATGTTGGCTAGGAAGAGCCAATCGGTGAGGGCTAGTGGGCCCGTGTGGAagag                                                                                         <  2:582290/150‑1                                                           ttGCTTATTCGACACGCTTTACCTTCGCATGGGCTAAGCGCCGACGATTGAGCCATTTCCGTCCTTCTTCGATTAAGAAGAGTGGGATCGGGATGATGGCTAGGAAGATCCAATCGGTGAGGGCTAGTGGGCCCGTGTGGAAGAGTTCCTg                                                                                   <  2:1010849/151‑1                                                                         cGCTTTACCTTCGCATGGGCTAAGCGCCGACGATTGAGCCATTTCCGTCCTTCTTCGATTAAGAAGAGTGGGATCGGGATGATGGCTAGGAAGATCCAATCGGTGAGGGCTAGTGGGCCCGTGTGGAAGAGTTCCTGTAAGAACGGGgtgt                                                                     >  2:38605/1‑151                                                                         cGCTTTACCTTCGCATGGGCTAAGCGCCGACGATTGAGCCATTTCCGTCCTTCTTCGATTAAGAAGAGTGGGATCGGGATGATGGCTAGGAAGATCCAATCGGTGAGGGCTAGTGGGCCCGTGTGGAAGAGTTCCTGTAAGAACGGGgtgt                                                                     >  2:35204/1‑151                                                                         cGCTTTACCTTCGCATGGGCTAAGCGCCGACGATTGAGCCATTTCCGTCCTTCTTCGATTAAGAAGAGTGGGATCGGGATGATGGCTAGGAAGATCCAATCGGTGAGGGCTAGTGGGCCCGTGTGGAAGAGTTCCTGTAAGAACGGGgtgt                                                                     >  1:601140/1‑151                                                                         cGCTTTACCTTCGCATGGGCTAAGCGCCGACGATTGAGCCATTTCCGTCCTTCTTCGATTAAGAAGAGTGGGATCGGGATGATGGCTAGGAAGATCCAATCGGTGAGGGCTAGTGGGCCCGTGTGGAAGAGTTCCTGTAAGAACGGGgtgt                                                                     >  2:207008/1‑151                                                                            tttACCTTCGCATGGGCTAAGCGCCGACGATTGAGCCATTTCCGTCCTTCTTCGATTAAGAAGAGTGGGATCGGGATGATGGCTAGGAAGATCCAATCGGTGAGGGCTAGTGGGCCCGTGTGGAAGAGTTCCTGTAAGAACGGGGTGTAAc                                                                  <  2:435905/151‑1                                                                            tttACCTTCGCATGGGCTAAGCGCCGACGATTGAGCCATTTCCGTCCTTCTTCGATTAAGAAGAGTGGGATCGGGATGATGGCTAGGAAGATCCAATCGGTGAGGGCTAGTGGGCCCGTGTGGAAGAGTTCCTGTAAGAACGGGGTGTAAc                                                                  <  1:986077/151‑1                                                                            tttACCTTCGCATGGGCTAAGCGCCGACGATTGAGCCATTTCCGTCCTTCTTCGATTAAGAAGAGTGGGATCGGGATGATGGCTAGGAAGATCCAATCGGTGAGGGCTAGTGGGCCCGTGTGGAAGAGTTCCTGTAAGAACGGGGTGTAAc                                                                  <  1:763163/151‑1                                                                            tttACCTTCGCATGGGCTAAGCGCCGACGATTGAGCCATTTCCGTCCTTCTTCGATTAAGAAGAGTGGGATCGGGATGATGGCTAGGAAGATCCAATCGGTGAGGGCTAGTGGGCCCGTGTGGAAGAGTTCCTGTAAGAACGGGGTGTAAc                                                                  <  1:762595/151‑1                                                                            tttACCTTCGCATGGGCTAAGCGCCGACGATTGAGCCATTTCCGTCCTTCTTCGATTAAGAAGAGTGGGATCGGGATGATGGCTAGGAAGATCCAATCGGTGAGGGCTAGTGGGCCCGTGTGGAAGAGTTCCTGTAAGAACGGGGTGTAAc                                                                  <  1:633946/151‑1                                                                            tttACCTTCGCATGGGCTAAGCGCCGACGATTGAGCCATTTCCGTCCTTCTTCGATTAAGAAGAGTGGGATCGGGATGATGGCTAGGAAGATCCAATCGGTGAGGGCTAGTGGGCCCGTGTGGAAGAGTTCCTGTAAGAACGGGGTGTAAc                                                                  <  1:318579/151‑1                                                                            tttACCTTCACATGGGCTAAGCGCCGACGATTGAGCCATTTCCGTCCTTCTTCGATTAAGAAGAGTGGGATCGGGATGATGGCTAGGAAGATCCAATCGGTGAGGGCTAGTGGGCCCGTGGGGAAGAGTTCCTGTAAGAACGGGGTGTAAc                                                                  <  2:687765/151‑1                                                                              tACCTTCGCATGGGCTAAGCGCCGACGATTGAGCCATTTCCGTCCTTCTTCGATTAAGAAGAGTGGGATCGGGATGATGGCTAGGAAGATCCAATCGGTGAGGGCTAGTGGGCCCGTGTGGAAGAGTTCCTGTAAGAACGGGGTGTAACTg                                                                <  2:452382/151‑1                                                                              tACCTTCGCATGGGCTAAGCGCCGACGATTGAGCCATTTCCGTCCTTCTTCGATTAAGAAGAGTGGGATCGGGATGATGGCTAGGAAGATCCAATCGGTGAGGGCTAGTGGGCCCGTGTGGAAGAGTTCCTGTAAGAACGGGGTGTAACTg                                                                <  2:959133/151‑1                                                                              tACCTTCGCATGGGCTAAGCGCCGACGATTGAGCCATTTCCGTCCTTCTTCGATTAAGAAGAGTGGGATCGGGATGATGGCTAGGAAGATCCAATCGGTGAGGGCTAGTGGGCCCGTGTGGAAGAGTTCCTGTAAGAACGGGGTGTAACTg                                                                <  1:1030353/151‑1                                                                              tACCTTCGCATGGGCTAAGCGCCGACGATTGAGCCATTTCCGTCCTTCTTCGATTAAGAAGAGTGGGATCGGGATGATGGCTAGGAAGATCCAATCGGTGAGGGCTAGTGGGCCCGTGTGGAAGAGTTCCTGTAAGAACGGGGTGTAACTg                                                                <  1:1040254/151‑1                                                                              tACCTTCGCATGGGCTAAGCGCCGACGATTGAGCCATTTCCGTCCTTCTTCGATTAAGAAGAGTGGGATCGGGATGATGGCTAGGAAGATCCAATCGGTGAGGGCTAGTGGGCCCGTGTGGAAGAGTTCCTGTAAGAACGGGGTGTAACTg                                                                <  2:113471/151‑1                                                                                    cGCATGGGCTAAGCGCCGACGATTGAGCCATTTCCGTCCTTCTTCGATTAAGAAGAGTGGGATCGGGATGATGGCTAGGAAGATCCAATCGGTGAGGGCTAGTGGGCCCGTGTGGAAGAGTTCCTGTAAGAACGGGGTGTAACTGAGTAgg                                                          >  2:1075100/1‑151                                                                                      cATGGGCTAAGCGCCGACGATTGAGCCATTTCCGTCCTTCTTCGATTAAGAAGAGTGGGATCGGGATGATGGCTAGGAAGATCCAATCGGTGAGGGCTAGTGGGCCCGTGTGGAAGAGTTCCTGTAAGAACGGGGTGTAACTGAGTAGGGc                                                        <  2:278882/151‑1                                                                                          ggCTAAGCGCCGACGATTGAGCCATTTCCGTCCTTCTTCGATTAAGAAGAGTGGGATCGGGATGATGGCTAGGAAGATCCAATCGGTGAGGGCTAGTGGGCCCGTGTGGAAGAGTTCCTGTAAGAACGGGGTGTAACTGAGTAGGGCcagc                                                    >  2:892350/1‑151                                                                                          ggCTAAGCGCCGACGATTGAGCCATTTCCGTCCTTCTTCGATTAAGAAGAGTGGGATCGGGATGATGGCTAGGAAGATCCAATCGGTGAGGGCTAGTGGGCCCGTGTGGAAGAGTTCCTGTAAGAACGGGGGGTAACTGAGTAGGGCcagc                                                    >  2:841186/1‑151                                                                                          ggCTAAGCGCCGACGATTGAGCCATTTCCGTCCTTCTTCGATTAAGAAGAGCGGGATCGGGATGATGGCTAGGAAGATCCAATCGGCGAGGGCTAGTGGGCCCGTGTGGAAGAGTTCCTGTAAGAACGGGGTGTAACTGAGTCGGGCcagc                                                    >  2:566023/1‑151                                                                                           gCTAAGCGCCGACGATTGAGCCATTTCCGTCCTTCTTCGATTAAGAAGAGTGGGATCGGGATGATGGCTAGGAAGATCCAATCGGTGAGGGCTAGTGGGCCCGTGTGGAAGAGTTCCTGTAAGAACGGGGTGTAACTGAGTAGGGCcagca                                                   <  1:1071097/151‑1                                                                                                                 aTTTCCGTCCTTCTTCGATTAAGAAGAGTGGGATCGGGATGATGGCTAGGAAGATCCAATCGGTGAGGGCTAGTGGGCCCGTGTGGAAGAGTTCCTGTAAGAACGGGGTGTAACTGAGTAGGGCCAGCAGGAAAATTTCAAAGATAATGcc                             <  1:328942/151‑1                                                                                                                          cttcttCGATTAAGAAGAGTGGGATCGGGATGATGGCTAGGAAGATCCAATCGGTGAGGGCTAGTGGGCCCGTGTGGAAGAGTTCCTGTAAGAACGGGGTGTAACTGAGTAGGGCCAGCAGGAAAATTTCAAAGATAATGCCGAACCAGAc                    <  2:933118/151‑1                                                                                                                          cttcttCGATTAAGAAGAGTGGGATCGGGATGATGGCTAGGAAGATCCAATCGGTGAGGGCTAGTGGGCCCGTGTGGAAGAGTTCCTGTAAGAACGGGGTGTAACTGAGTAGGGCCAGCAGGAAAATTTCAAAGATAATGCCGAACCAGAc                    <  1:1063542/151‑1                                                                                                                                  aTTAAGAAGAGTGGGATCGGGATGATGGCTAGGAAGATCCAATCGGTGAGGGCTAGTGGGCCCGTGTGGAAGAGTTCCTGTAAGAACGGGGTGTAACTGAGTAGGGGCAGCAGGAAAAATTCAAAGATAATGCAGAACCAgccgaggcgct            >  2:691174/1‑141                                                                                                                                  aTTAAGAAGAGTGGGATCGGGATGATGGCTAGGAAGATCCAATCGGTGAGGGCTAGTGGGCCCGTGTGGAAGAGTTCCTGTAAGAACGGGGTGTAACTGAGTAGGGCCAGCAGGAAAATTTCAAAGATAATGCCGAACCAGACAAGTCGAt            >  1:995577/1‑151                                                                                                                                  aTTAAGAAGAGTGGGATCGGGATGATGGCTAGGAAGATCCAATCGGTGAGGGCTAGTGGGCCCGTGTGGAAGAGTTCCTGTAAGAACGGGGTGTAACTGAGTAGGGCCAGCAGGAAAATTTCAAAGATAATGCCGAACCAGACAAGTCGAt            >  2:645522/1‑151                                                                                                                                  aTTAAGAAGAGTGGGATCGGGATGATGGCTAGGAAGATCCAATCGGTGAGGGCTAGTGGGCCCGTGTGGAAGAGTTCCTGTAAGAACGGGGTGTAACTGAGTAGGGCCAGCAGGAAAATTTCAAAGATAATGCCGAACCAGACAAGTCGAt            >  1:173216/1‑151                                                                                                                                  aTTAAGAAGAGTGGGATCGGGATGATGGCTAGGAAGATCCAATCGGTGAGGGCTAGTGGGCCCGTGTGGAAGAGTTCCTGTAAGAACGGGGTGTAACTGAGTAGGGCCAGCAGGAAAATTTCAAAGATAATGCCGAACCAGACAAGTCGAt            >  2:140958/1‑151                                                                                                                                          gagTGGGATCGGGATGATGGCTAGGAAGATCCAATCGGTGAGGGCTAGTGGGCCCGTGTGGAAGAGTTCCTGTAAGAACGGGGTGTAACTGAGTAGGGCCAGCAGGAAAATTTCAAAGATAATGCCGAACCAGACAAGTCGATTGGcaaaa    >  2:93432/1‑150                                                                                                                                            gTGGGATCGGGATGATGGCTAGGAAGATCCAATCGGTGAGGGCTAGTGGGCCCGTGTGGAAGAGTTCCTGTAAGAACGGGGTGTAACTGAGTAGGGCCAGCAGGAAAATTTCAAAGATAATGCCGAACCAGACAAGTCGATTGGcaaacaa  >  1:448981/1‑151                                                                                                                                                |                                                                                                                                                   AAAGTGACACTAAAAAATGGCCTCGGAAAGTTCAAATTCCAAGGCCATTTTCTTAGTGTTGCTTATTCGACACGCTTTACCTTCGCATGGGCTAAGCGCCGACGATTGAGCCATTTCCGTCCTTCTTCGATTAAGAAGAGTGGAATCGGGATGATGGCTAGGAAGATCCAATCGGTGAGGGCTAGTGGGCCCGTGTGGAAGAGTTCCTGTAAGAACGGGGTGTAACTGAGTAGGGCCAGCAGGAAAATTTCAAAGATAATGCCGAACCAGACAAGTCGATTGGCAAACAA  .  NC_013198/985222‑985511 Base quality scores:  ATCG  < 3 ≤  ATCG  < 7 ≤  ATCG  < 22 ≤  ATCG  < 33 ≤  ATCG  < 41 ≤  ATCG