Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,556,200	0	G	A	44.3%	80.2 / 37.4	70	S217F (TCC→TTC)	rtcA	RNA 3'‑terminal phosphate cyclase
Reads supporting (aligned to +/- strand):  ref base G (32/7);  new base A (0/31);  total (32/38)
Fisher's exact test for biased strand distribution p-value = 2.26e-13
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ACCAACTGTGCCGCAACCACCTCGGCACTGACGCGCTTTTCACCGACGACAAAAAAGCG‑TTCGGTGATATTTTCACTTTCGACTTCAAGCGAGACGGTATTACCCGGCCCCTGGTCGCGCGGCAGGTTATGAATATTCTGTTCAT‑GCAGGG‑AAAAACTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGTTAAACGATGCCACGGGCGAGACTTCC  >  NC_000913/3556050‑3556343                                                                                                                                                         |                                                                                                                                                 aCCAACTGTGCCGCAACCACCTCGGCACTGACGCGCTTTTCACCGACGACAAAAAAGCG‑TTCGGTGATATTTTCACTTTCGACTTCAAGCGAGACGGTATTACCCGGCCCCTGGTCGCGCGGCAGGTTATGAATATTCTGTTCAT‑GCAggg                                                                                                                                                  >  2:490083/1‑151 (MQ=255)   cAACTGTGCCGCAACCACCTCGGCACTGACGCGCTTTTCACCGACGACAAAAAAGCG‑TTCGGTGATATTTTCACTTTCGACTTCAAGCGAGACGGTATTACCCGGCCCCTGGTCGCGCGGCAGGTTATGAATATTCTGTTCAT‑GCAGGG‑a                                                                                                                                                >  1:155739/1‑150 (MQ=255)     aCTGTGCCGCAACCACCTCGGCACTTACGCGCTTTTCACCGACTACAAAAAAGCG‑TTCGGTGATATTTGCACTTTCGACTGCAAGCGAGACGGTATTACCCGGCCCCTGGTCGCGCGGCAGGTTATGAATATTCTGTTCAT‑GCAGCG‑aaaa                                                                                                                                             <  1:599098/151‑1 (MQ=255)        gtgCCGCAACCACCTCGGCACTGACGCGCTTTTCACCGACGACAAAAAAGCG‑TTCGGTGATATTTTCACTTTCGACTTCAAGCGAGACGGTATTACCCGGCCCCTGGTCGCGCGGCAGGTTATGAATATTCTGTTCAT‑GCAGGG‑AAAAAc                                                                                                                                           <  2:766409/150‑1 (MQ=255)          gCCGCAACCACCTCGGCACTGACGCGCTTTTCACCGACGACAAAAAAGCG‑TTCGGTGATATTTTCACTTTCGACTTCAAGCGAGACGGTATTACCCGGCCCCTGGTCGCGCGGCAGGTTATGAATATTCTGTTCAT‑GCAGGG‑AAAAACTcc                                                                                                                                        <  2:274079/151‑1 (MQ=255)           cctcAACCACCTCGGCACTGACGCGCTTTTCACCTACGACAAAAAAGCG‑TTCGGTGATATTTTCACTTTCGACTTCAAGCGAGACGGTATTACCCGGCCCCTTGTCGCGCGGCAGGTTATGAATATTCTGTTCAT‑GCAGGG‑AAAAACTc                                                                                                                                         <  2:961766/146‑1 (MQ=255)             gCAACCACCTCGGCACTGACGCGCTTTTCACCGACGACAAAAAAGCG‑TTCGGTGATATTTTCACTTTCGACTTCAAGCGAGACGGTATTACCCGGCCCCTGGTCGCGCGGCAGGTTATGAATATTCTGTTCAT‑GCAGGG‑AAAAACTCCCCg                                                                                                                                     >  2:1038273/1‑151 (MQ=255)                 ccaccTCGGCACTGACGCGCTTTCCACCGACGACAAAAAAGCG‑TTCGGTGATATTTTCACTTTCGACTTCAAGCGAGACGGTATTCCCCGGCCCCTGGTCGCGCGGCGGGTTATGAATATTCTGTTCAT‑GCAGGG‑AAAAACTCCCCGCCAg                                                                                                                                 <  1:34113/151‑1 (MQ=255)                   gccTCGGCACTGACGCGCTTTTCACCGACGACAAAAAAGCG‑TTCGGTGATATTTTCACTTTCGACTTCAAGCGAGACGGTATTACCCGGCCCCTGGTCGCGCGGCAGGTTATGAATATTCTGTTCAT‑GCAGGG‑AAAAACTCCCCGCCAgtg                                                                                                                               >  1:800155/2‑151 (MQ=255)                     cTCGGCACTGACGCGCTTTTCACCGACGACAAAAAAGCG‑TTCGGTGATATTTTCACTTTCGACTTCAAGCGAGACGGTATTACCCGGCCCCTGGTCGCGCGGCAGGTTATGAATATTCTGTTCAT‑GCAGGG‑AAAAACTCCCCGCCAgtgt                                                                                                                              >  1:1087748/1‑150 (MQ=255)                     cTCGGCACTGACGCGCTTTTCACCGACGACAAAAAAGCG‑TTCGGTGATATTTTCACTTTCGACTTCAAGCGAGACGGTATTACCCGGCCCCTGGTCGCGCGGCAGGTTATGAATATTCTGTTCAT‑GCAGGG‑AAAAACTCCCCGCCAgtgt                                                                                                                              >  1:1176217/1‑150 (MQ=255)                        ggCACTGACGCGCTTTTCACCGACGACAAAAAAGCG‑TTCGGTGATATTTTCACTTTCGACTTCAAGCGAGACGGTATTACCCGGCCCCTGGTCGCGCGGCAGGTTATGAATATTCTGTTCAT‑GCAGGG‑AAAAACTCCCCGCCAGTGTAGCg                                                                                                                          >  2:479702/1‑151 (MQ=255)                             tGACGCGCTTTTCACCGACGACAAAAAAGCG‑TTCGGTGATATTTTCACTTTCGACTTCAAGCGAGACGGTATTACCCGGCCCCTGGTCGCGCGGCAGGTTATGAATATTCTGTTCAT‑GCAGGG‑AAAAACTCCCCGCCAGTGTAGCGATTTc                                                                                                                     >  1:927770/1‑151 (MQ=255)                                  cgcTTTTCACCGACGACAAAAAAGCG‑TTCGGTGATATTTTCACTTTCGACTTCAAGCGAGACGGTATTACCCGGCCCCTGGTCGCGCGGCAGGTTATGAATATTCTGTTCAT‑GCAGGG‑AAAAACTCCCCGCCAGTGTAGCGATTTCACGc                                                                                                                 >  1:570830/1‑150 (MQ=255)                                     ttttCCCCAACGAAAAAAAAGCG‑TTGGGTGATTTTTTCCCTTTCGCCTTCAAGCGAGACGGTATTACCCGGCCCCGGGTCGCGCGGCAGTTTATGAATATTCTGTTCAT‑GCAGGA‑AAAACCTCCCCGCCAGTGTAGCGATTTCACGCTCa                                                                                                              <  2:62349/150‑1 (MQ=255)                                     tttcccccGACAAAAAAAAAGCG‑TTCGGTGATATTTTCCCTTTCGACTTCAAGCGAGACGGTATTCCCCGCCCCCTGGTCGCGCGGCAGTTTATGAATATTCTGTTCAT‑GCAGGA‑AAAACCTCCCCGCCAGTGTAGCGATTTCACg                                                                                                                  <  2:339871/140‑1 (MQ=255)                                     ttcccACCGCCGACAAAAACGGT‑TGGGGTGTATTTTTCCCTTCCGCCTACAAGCGAGAGGTTATTCCCCGGCCCCGGGTCGCGCGGCGGGTTATGATTATTCTTTTCAT‑GAGGGA‑AAAAACCCCCCGCCAGTGTAGCGATTTCACGCTc                                                                                                               <  1:1139377/145‑1 (MQ=255)                                      ttgaccccgccaaaaaaaaGGGTTTCGGTGATATTTTCACTTTCGACTCAAAGCGAGCGGTTATTCCCCGCCCCCTGGTCGCGCGGCGGTTTATAAATTTTCTTTTCGG‑GCAGGA‑AAAACCTCCCCGCCAGTGTAGCGAtt                                                                                                                       <  1:803026/128‑1 (MQ=255)                                         ccccGACAAAAAAAAGGCG‑TTCGGTGATTTTTTCCCTTTCGACTCCAGGCGAGACGGTATTCCCCGGCCCCGGGTCGCGCGGCGGGTTATAAATTTTCTTTTCAT‑GCAGGA‑AAAACCCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAAc                                                                                                         <  2:326030/149‑1 (MQ=255)                                              acgacAAAAAAGCG‑TTCGGTGATATTTTCACTTTCGACTTCAAGCGAGACGGTATTACCCGGCCCCTGGTCGCGCGGCAGGTTATGAATATTCTGTTCAT‑GCAGGG‑AAAAACTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGgc                                                                                                    >  2:9461/1‑151 (MQ=255)                                               caaaaaaagaGCG‑TTGGTTAATTTTTTCTCTTTCGTCTTAAAGCGAGGCGTTATCCCCCGGCCCCGGGTCGCGCGGGAGTTTATAAATATTTTGTGCAG‑GCAGGA‑AAAACCCCCCCGCCAGTGTAGCGATTTCACGCTCAGc                                                                                                            <  1:249646/138‑1 (MQ=255)                                                 aaaaaaaaGGG‑TTCGGTGATATTTTCCCTTTCGACTTAAAGCGAGACGGTATTCCCGGGCCCGGGGTCGCGCGGCAGGTTATGAATATTTTGTTCAG‑GCGGGA‑AAAACCCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGcgcg                                                                                                 <  2:186734/149‑1 (MQ=255)                                                   aaaaaaGCG‑TTCGGTGATATTTTCACTTTCGACTTCAAGCGAGACGGTATTACCCGGCCCCTGGTCGCGCGGCAGGTTATGAATATTCTGTTCAT‑GCAGGG‑AAAAACTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGc                                                                                               >  2:782306/1‑151 (MQ=255)                                                    aaaaaGCG‑TTGGGTAATATTTTCACTTTCGACTACAAGCGAGACGGTATTCCCCGCCCCCTGGTCGCGCGGCAGGTTATAAATATTCTGTTCAT‑GCAGGA‑AAAACCCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGca                                                                                              <  1:120704/151‑1 (MQ=255)                                                    aaaaaGCG‑TTCGGTGATATTTTCACTTTCGACTTCAAGCGAGACGGTATTACCCGGCCCCTGGTCGCGCGGCAGGTTATGAATATTCTGTTCAT‑GCAGGG‑AAAAACTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGc                                                                                               >  2:352432/1‑150 (MQ=255)                                                    aaaaaCCG‑TGCGGTTATTTTTCCCTTTTCGGTTTCAAAGGAGACGGTTTTCCCCGGCCCCGGGCCGCGCGGGGGTTTATAATTATTC‑GTTTATCGGAGGA‑AAAACCTCCCCGCCAGTGTAGCGATTTCACGCTCAGCACCATGgc                                                                                                    <  1:424542/145‑1 (MQ=255)                                                           g‑TTCGGTGATATTTTCACTTTCGACTTCAAGCGAGACGGTATTACCCGGCCCCTGGTCGCGCGGCAGGTTATGAATATTCTGTTCAT‑GCAGGG‑AAAAACTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGc                                                                                       >  2:1016134/1‑151 (MQ=255)                                                             ttCGGTGATTTTTTCACTTCCGCCTTCAAGCGAGACGTTATTCCCCGCCCCCTGGTCGCGCGGCGGTTTATGAATATTCTGTTCAT‑GCAGGA‑AAAACCCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCa                                                                                         <  2:81976/148‑1 (MQ=255)                                                             ttCGGTGATATTTTCACTTTCGACTTCAAGCGAGACGGTATTACCCGGCCCCTGGTCGCGCGGCAGGTTATGAATATTCTGTTCAT‑GCAGGG‑AAAAACTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCt                                                                                      >  2:229717/1‑151 (MQ=255)                                                                ggTGATTTTTTCCCTTCCGACTGCAAGCGGGAGTGTATCCCCCGCCCCCTGGCCGGGCGGCGGTTTATGAATATTCTGTTCAG‑GGGGGA‑AAAACCCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTaa                                                                                    <  1:924818/150‑1 (MQ=255)                                                                       ttttCACTTTCGACTTCAAGCGAGACGGTATTACCCGGCCCCTGGTCGCGCGGCAGGTTATGAATATTCTGTTCAT‑GCAGGG‑AAAAACTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATgg                                                                                                     >  2:972917/1‑126 (MQ=255)                                                                           cccTTTCGCCTTCAAGCGAGACGGTATCCCCCGCCCCCTGGCCGCGCGGCGGTTTATGAATATTCTGTTCAT‑GCAGGA‑AAAACCCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTcc                                                                        <  2:757729/149‑1 (MQ=255)                                                                            ccgttctctttCAGGCGAGCGGGTATCCCGCGGCCCCGGTCCGGGCGGCTGTTTAGGAAATTTCTTTCCAG‑GCGAGA‑AAAACCCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCa                                                                       <  1:430661/142‑1 (MQ=255)                                                                                tcatctccaggCGAGACGGTATCCCCCGGCCCCGGCTCGCGCGGCATGTTATAAATATTCTGTCCAT‑GCGGGA‑AAAACCCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGc                                                                    <  1:494535/140‑1 (MQ=255)                                                                                 cGACTTCAAGCGAGACGTTTTTCCCCGGCCCCGGGTCGCGCGGCAGGTTATAAATATTCTGTTCAT‑GCAGGA‑AAAACCCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCAc                                                                      <  2:439723/147‑1 (MQ=255)                                                                                 cGACTCCAGGCGATCGGGTATCCCCCGGCCCCGGGTCGCGGGGCAGGTTATGATTTTTCTGTTCAG‑CCAGGA‑AAAACCCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCa                                                                   <  1:858336/150‑1 (MQ=255)                                                                                  gACTTCAAGCGAGACGGTATTACCCGGCCCCTGGTCGCGCGGCAGGTTATGAATATTCTGTTCAT‑GCAGGG‑AAAAACTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATc                                                                 >  1:1075898/1‑151 (MQ=255)                                                                                          gCGGGAGGGTATTCCCCGGCCCCGGGCCCGCCGGAAGTTTATGATTATTCTTTCAAT‑GCAGGA‑AAAACCCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGca                                                             <  1:1111322/147‑1 (MQ=255)                                                                                            gagaCGGTATTACCCGGCCCCTGGTCGCGCGGCAGGTTATGAATATTCTGTTCAT‑GCAGGG‑AAAAACTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGt                                                       >  2:1138074/1‑151 (MQ=255)                                                                                                 ggTATTACCCGGCCCCTGGTCGCGCGGCAGGTTATGAATATTCTGTTCAT‑GCAGGG‑AAAAACTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCcg                                                  >  1:1097161/1‑151 (MQ=255)                                                                                                 ggTATTACCCGGCCCCTGGTCGCGCGGCAGGTTATGAATATTCTGTTCAT‑GCAGGG‑AAAAACTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCcg                                                  >  1:180954/1‑151 (MQ=255)                                                                                                     ttcccgggCCCCTGGTCGCGCGGCGGGTTATAAATATTCTGTTCAT‑GCAGGA‑AAAACCCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCgcgc                                               <  1:248099/144‑1 (MQ=255)                                                                                                      tACCCGGCCCCTGGTCGCGCGGCAGGTTATGAATATTCTGTTCAT‑GCAGGG‑AAAAACTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCgcgcgc                                             >  1:505695/1‑149 (MQ=255)                                                                                                        cccGGCCCCTGGTCGCGCGGCAGTTTAGGAATATCCTGTTCAT‑GAAGGA‑AAAACCCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTc                                             <  1:245745/149‑1 (MQ=255)                                                                                                           ggCCCCTGGTCGCGCGGCAGGTTATGAATATTCTGTTCAT‑GCAGGG‑AAAAACTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCa                                         >  1:680446/1‑150 (MQ=255)                                                                                                           gccccgTGGCCGGGCGGAGGTTTTTGAATTTTCTTTCCAG‑GCAGAA‑AAAACCCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGcc                                          <  2:120523/147‑1 (MQ=255)                                                                                                             ccccTGGTCGCGCGGCAGGTTATGAATATTCTGTTCAT‑GCAGGG‑AAAAACTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGt                                      >  2:1155763/1‑151 (MQ=255)                                                                                                             ccccTGGTCGCGCGGCAGGTTATGAATATTCTGTTCAG‑GCAGGA‑AAAACCCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCaa                                        <  1:370729/149‑1 (MQ=255)                                                                                                             ccccGGGCCGCCGGGCAGGTTATGAATATTTTGTGCA‑‑GAAGGGAAAAAACCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGt                                      <  1:1138358/151‑1 (MQ=255)                                                                                                              cccGGGCCGCGGGGGAGTTTAGAAATATTCTTTTAAG‑GAAGGA‑AAAACCCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGtt                                     <  1:1037977/151‑1 (MQ=255)                                                                                                                    gcgcgagcggttatATGAAAATTTTTCTCAG‑GAAGGA‑AAAAACTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAg                                <  1:1133240/137‑1 (MQ=255)                                                                                                                     cgcgcgGCAGGTTATGAATATTCTGTTCAT‑GCAGGG‑AAAAACTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGgt                              >  1:483319/1‑151 (MQ=255)                                                                                                                     cgcgcgGCAGGTTATGAATATTCTGTTCAT‑GCAGGG‑AAAAACTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGgt                              >  2:512941/1‑151 (MQ=255)                                                                                                                      gcgcgGGAGTTTATGATTATCCTGTCAAT‑GCAGGG‑AAACCCTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGgt                              <  1:97771/150‑1 (MQ=255)                                                                                                                       cggggCAGTTTATAAATATTCTGTTCAT‑GCAGGA‑AAAACCCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGgt                              <  2:386833/146‑1 (MQ=255)                                                                                                                        gcgggggtttATAAATATTTTTTTCAG‑GCAGGA‑AAAACCTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAg                                <  2:1130503/138‑1 (MQ=255)                                                                                                                         cgGCAGGTTATGAATATTCTGTTCAT‑GCAGGG‑AAAAACTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGtt                           >  1:242388/1‑150 (MQ=255)                                                                                                                          ggCGGGTTATAAATATTTTGTTCAT‑GCAGGA‑AAAACCTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGTTa                          <  2:1038028/150‑1 (MQ=255)                                                                                                                          ggCAGGTTATGAATATTCTGTTCAT‑GCAGGG‑AAAAACTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGTTa                          >  2:1124924/1‑150 (MQ=255)                                                                                                                           gaggattATAAATTTTCTTTCCAG‑GCAGAA‑AAAAACCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGTTaaa                        <  1:901277/146‑1 (MQ=255)                                                                                                                                ttacaaaTATTCTTCTCAG‑GCAGGA‑AAAACCCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGTTAAACGAt                    <  1:209610/145‑1 (MQ=255)                                                                                                                                 tATGAATATTCTGTTCAT‑GCAGGG‑AAAAACTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGTTAAACGATGc                  >  2:443486/1‑151 (MQ=255)                                                                                                                                     aaTATTCTGTTCAT‑GCAGGG‑AAAAACTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGTTAAACGATGCCACg              >  1:615830/1‑151 (MQ=255)                                                                                                                                         tccTGTTCAT‑GCAGGA‑AAAACCCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGTTAAACGATGCCACGGGc           <  1:152244/148‑1 (MQ=255)                                                                                                                                         gttcttgtagggaaGGA‑AAAACCCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGTTAAACGATGcc                 <  1:1045929/132‑1 (MQ=255)                                                                                                                                          tCTGTTCAT‑GCAGGG‑AAAAACTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGTTAAACGATGCCACGGGc           >  1:551972/1‑149 (MQ=255)                                                                                                                                            ttttCAT‑GGGGGA‑AAAACCCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGTTAAACGATGCCACGGGCgag        <  2:1041534/148‑1 (MQ=255)                                                                                                                                            tGTTCAT‑GCAGGG‑AAAAACTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGTTAAACGATGCCACGGGCg          >  1:1070607/1‑148 (MQ=255)                                                                                                                                                cAT‑GCAGGG‑AAAAACTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGTTAAACGATGCCACGGGCGAGACTTc   >  1:822859/1‑151 (MQ=255)                                                                                                                                                cAT‑GCAGGG‑AAAAACTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGTTAAACGATGCCACGGGCGAGACTTc   >  1:528364/1‑151 (MQ=255)                                                                                                                                                 aT‑GCAGGG‑AAAAACTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGTTAAACGATGCCACGGGCGAGACTTcc  >  2:1156108/1‑151 (MQ=255)                                                                                                                                                       aag‑AAAATCTCCCCGATAGTGTAGCGTATCGACGCTCAGCAACAGGGCGCGGCCCACCAGCGAATAGAACTTCTCCTCGCATCTGCACAATGTTCCCGCGCTAGACCAGTTGCAAGATGTTATACGATGCTACGGGTGAGACtt    <  1:212966/142‑1 (MQ=255)                                                                                                                                                         |                                                                                                                                                 ACCAACTGTGCCGCAACCACCTCGGCACTGACGCGCTTTTCACCGACGACAAAAAAGCG‑TTCGGTGATATTTTCACTTTCGACTTCAAGCGAGACGGTATTACCCGGCCCCTGGTCGCGCGGCAGGTTATGAATATTCTGTTCAT‑GCAGGG‑AAAAACTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGTTAAACGATGCCACGGGCGAGACTTCC  >  NC_000913/3556050‑3556343

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑