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

ACCAACTGTGCCGCAACCACCTCGGCACTGACGCGCTTTTCACCGACGACAAAAAAGC‑GTTCGGTGATATTTTCACTTTCGACTTCAAGCGAGACGGTATTACCCGGCCCCT‑GGTCGCGCGGCAGGTTATGAATATTCTGTTCATGCAGGGAAAAA‑CTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGTTAAACGATGCCACGGGCGAGACTT  >  NC_000913/3556050‑3556341                                                                                                                                                         |                                                                                                                                               aCCAACTGTGCCGCAACCACCTCGGCACTGACGCGCTTTTCACCGACGACAAAAAAGC‑GTTCGGTGATATTTTCACTTTCGACTTCAAGCGAGACGGTATTACCCGGCCCCT‑GGTCGCGCGGCAGGTTATGAATATTCTGTTCATGCAggg                                                                                                                                                >  2:425756/1‑151 (MQ=255) aCCAACTGTGCCGCAACCACCTCGGCACTGACGCGCTTTTCACCGACGACAAAAAAGC‑GTTCGGTGATATTTTCACTTTCGACTTCAAGCGAGACGGTATTACCCGGCCCCT‑GGTCGCGCGGCAGGTTATGAATATTCTGTTCATGCAggg                                                                                                                                                >  1:535837/1‑151 (MQ=255)  ccAACTGTGCCGCAACCACCTCGGCACTGACGCGCTTTTCACCGACGACAAAAAAGC‑GTTCGGTGATATTTTCACTTTCGACTTCAAGCGAGACGGTATTACCCGGCCCCT‑GGTCGCGCGGCAGGTTATGAATATTCTGTTCATGCAggg                                                                                                                                                <  2:788823/150‑1 (MQ=255)  ccAACTGTGCCGCAACCACCTCGGCACTGACGCGCTTTTCACCGACGACAAAAAAGC‑GTTCGGTGATATTTTCACTTTCGACTTCAAGCGAGACGGTATTACCCGGCCCCT‑GGTCGCGCGGCAGGTTATGAATATTCTGTTCATGCAGGGa                                                                                                                                               <  2:18604/151‑1 (MQ=255)      cTGTGCCGCAACCACCTCGGCACTGACGCGCTTTTCACCGACGACAAAAAAGC‑GTTCGGTGATATTTTCACTTTCGACTTCAAGCGAGACGGTATTACCCGGCCCCT‑GGTCGCGCGGCAGGTTATGAATATTCTGTTCATGCAGGGaaaa                                                                                                                                            >  2:455108/1‑150 (MQ=255)        gtgCCGCAACCACCTCGGCACTGACGCGCTTTTCACCGACGACAAAAAAGC‑GTTCGGTGATATTTTCACTTTCGACTTCAAGCGAGACGGTATTACCCGGCCCCT‑GGTCGCGCGGCAGGTTATGAATATTCTGTTCATGCAGGGAAAAA‑c                                                                                                                                         <  2:539827/150‑1 (MQ=255)          gCCGCAACCAACTCGGCACTGACGCGCTTTTCACCGACGACAAAAAAGC‑GTTCGGTGATATTTTCACTTTCGACTACAAGCGAGACGGTATTACCCGGCCCCT‑GGTCGCGCGGCAGGTTATGAATATTCTGTTCATGCAGGGAAAAA‑CTcc                                                                                                                                      >  1:351624/1‑151 (MQ=255)              cAACCACCTCGGCACTGACGCGCTTTTCACCGACGACAAAAAAGC‑GTTCGGTGATATTTTCCCTTTCGACTTCAAGCGAGACGGTATTACCCGGCCCCT‑GGTCGCGCGGCAGGTTATGAATATTCTGTTCATGCAGGGAAAAA‑CTCCCCGc                                                                                                                                  <  2:429781/151‑1 (MQ=255)              cAACCACCTCGGCACTGACGCGCTTTTCACCGACGACAAAAAAGC‑GTTCGGTGATATTTTCACTTTCGACTTCAAGCGAGACGGTATTACCCGGCCCCT‑GGTCGCGCGGCAGGTTATGAATATTCTGTTCATGCAGGGAAAAA‑CTCCCCGc                                                                                                                                  >  2:94678/1‑151 (MQ=255)               aaCCACCTCGGCACTGACGCGCTTTTCACCGACGACAAAAAAGC‑GTTCGGTGATATTTTCACTTTCGACTTCAAGCGAGACGGTATTACCCGGCCCCT‑GGTCGCGCGGCAGGTTATGAATATTCTGTTCATGCAGGGAAAAA‑CTCCCCGc                                                                                                                                  >  2:729743/1‑150 (MQ=255)                accaccTCGGCACTGACGCGCTTTTCACCGACGACAAAAAAGC‑GTTCGGTGATATTTTCACTTTCGACTTCAAGCGAGACGGTATTACCCGGCCCCT‑GGTCGCGCGGCAGGTTATGAATATTCTGTTCATGCAGGGAAAAA‑CTCCCCGcc                                                                                                                                 >  2:415618/1‑150 (MQ=255)                   cccccgccACTGACGCTTTTTTCCCCGCCAAAAAAAAAGC‑GTGGGGTGATTTTTTCCCTTGCGACTTCAAGCGAGAGGGTATTCCCCGCCCCCT‑GGCCGCGCGGCAGTTTATGAATATTCTGTTCATGCGGGAAAAAA‑CCCCCCGCCAgtg                                                                                                                             <  1:1005347/144‑1 (MQ=255)                        gacccTGACGCGCTTTCCCCCGACGACAAAAAAGG‑GTTCGGTGATATTTTCCCTTTCGACTTCAAGCGAGCCGGTATTCCCCGCCCCCT‑GGCCGCGCGGCAGGTTATGAATATTCTTTTCTTGCAGGAAAAAC‑CTCCCCGCCAGTGTAGCg                                                                                                                        <  2:676413/147‑1 (MQ=255)                              gACGCGCTTTTCACCGACGACAAAAAAGC‑GTTCGGTGATATTTTCACTTTCGACTTCAAGCGAGACGGTATTACCCGGCCCCT‑GGTCGCGCGGCAGGTTATGAATATTCTGTTCATGCAGGGAAAAA‑CTCCCCGCCAGTGTAGCGATTTc                                                                                                                   >  2:943766/1‑150 (MQ=255)                                     ttttCACCGACGACAAAAAAGC‑TTTCGGTAATATTTTCACTTTCGACTTCAAGCAAGACGGTATTCCCCGCCCCTT‑GGTCGCGCGGAAGTTTATAAATATTCTGTTCAGGCAGGAAAAAC‑CTCCCCGCCAGTGTAGCGATTTCACGCTCa                                                                                                            <  2:384528/150‑1 (MQ=255)                                        ccccccaccaaaaaaaaGC‑GTGCGGTGTTTTTTCCTTTTCCGACTTCAAGCGAGACGTTTTTCCCCGCCCCTT‑GGCCGCGCGGCAGTTTATGAATTTTCTTTACATGAAGGAAAAAC‑CTCCCCGCCAGTGTAGCGATTTCACGCTc                                                                                                             <  1:916555/135‑1 (MQ=255)                                               cgacAAAAAGGC‑GTTCGGTGATTTTTTCCCTTTCGACTCCAAGCGAGACGTTTTCCCCCGCCCCCG‑GGTCGCGCGGCGGGTTATAAATATTCTTTTCATGAAGGAAAAAA‑CCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACAt                                                                                                     <  1:124349/147‑1 (MQ=255)                                                gacAAAAAAGC‑GTTCGGTGATATTTTCACTTTCGACTTCAAGCGAGACGGTATTACCCGGCCCCT‑GGTCGCGCGGCAGGTTATGAATATTCTGTTCATGCAGGGAAAAA‑CTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGgcgc                                                                                                >  1:650007/1‑151 (MQ=255)                                                  aaaaaaaGC‑GTTGGGTGATTTTTTCACTTGCGACTCCAAGCGGGACGTTATCCCCCGCCCCCT‑GTTCGCGCGGCGGTTTATAAATATTTTTTTCTTGCAGGAAAAAC‑CCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGgcgc                                                                                                <  2:125384/148‑1 (MQ=255)                                                   aaaaaaGC‑GTTCGGTGATATTTTCACTTTCGACTTCAAGCGAGACGGTATTACCCGGCCCCT‑GGTCGCGCGGCAGGTTATGAATATTCTGTTCATGCAGGGAAAAA‑CTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGc                                                                                             >  2:875517/1‑151 (MQ=255)                                                     aaaaGCGGTTGGTTGATTTTTCCCCTTTCGCCTCAAAGCGAGACGGTATCCCCCGGCCCCT‑GGTCGCGCGGCGGGTTATAAATATTCTGTTCAGCCAGGAAAAAA‑CCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGca                                                                                            <  2:813823/151‑1 (MQ=255)                                                      aaggC‑GTTGGGTTATTTTTTCCCTTCCGCCTCCAGGGGGGAGGTTTTTCCCCGCCCCCG‑GGCCGCGCGGCAGTTTTTGAATTTTCTGTTCAGGGAGGAAAAAC‑CCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGcaca                                                                                          <  1:977626/148‑1 (MQ=255)                                                      aaaGC‑GTTCGGTGATATTTTCACTTTCGACTTCAAGCGAGACGGTATTACCCGGCCCCT‑GGTCGCGCGGCAGGTTATGAATATTCTGTTCATGCAGGGAAAAA‑CTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGcac                                                                                           >  2:946110/1‑150 (MQ=255)                                                          c‑GTTCGGTGATATTTTCACTTTCGACTTCAAGCGAGACGGTATTACCCGGCCCCT‑GGTCGCGCGGCAGGTTATGAATATTCTGTTCATGCAGGGAAAAA‑CTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACAcc                                                                                        >  1:753425/1‑149 (MQ=255)                                                             ttCGGTGATATTTTCACTTTCGACTTCAAGCGAGACGGTATTACCCGGCCCCT‑GGTCGCGCGGCAGGTTATGAATATTCTGTTCATGCAGGGAAAAA‑CTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCt                                                                                    >  2:160005/1‑151 (MQ=255)                                                                 gTGATATTTTCACTTTCGACTTCAAGCGAGACGGTATTACCCGGCCCCT‑GGTCGCGCGGCAGGTTATGAATATTCTGTTCATGCAGGGAAAAA‑CTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAAt                                                                                 >  1:564048/1‑150 (MQ=255)                                                                     ttttttCCCTTTCGACTTAAAGCGAGCCGTTATTCCCCGCCCCCT‑GGTCGCGCGGCGGTTTATAAATTTTCTGTTCAGGCAGGAAAAAC‑CCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGaa                                                                             <  1:518763/148‑1 (MQ=255)                                                                     tatTTCCACTTTGGACTTAAAGCGAGAGGGTATTCCCCGCCCCCT‑GGCCGCGGGGCGGGTTATAAATATTCTGTTCATGCAGGAAAAAC‑CCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGaa                                                                             <  2:514119/150‑1 (MQ=255)                                                                        tttCCCTTCCGACTTCAAGCGAGACGTTATTCCCCGCCCCCT‑GGTCGCGCGGCAGTTTATGAATTTTCTGTTCATGCAGGAAAAAA‑CCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACtt                                                                          <  2:706373/150‑1 (MQ=255)                                                                         ttCACTTTCGACTTCAAGCGAGACGGTATTACCCGGCCCCT‑GGTCGCGCGGCAGGTTATGAATATTCTGTTCATGCAGGGAAAAA‑CTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTtct                                                                        >  1:129804/1‑151 (MQ=255)                                                                         ttCACTTTCGACTTCAAGCGAGACGGTATTACCCGGCCCCT‑GGTCGCGCGGCAGGTTATGAATATTCTGTTCATGCAGGGAAAAA‑CTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTtct                                                                        >  2:165238/1‑151 (MQ=255)                                                                         ttCACTTTCGACTTCAAGCGAGACGGTATTACCCGGCCCCT‑GGTCGCGCGGCAGGTTATGAATATTCTGTTCATGCAGGGAAAAA‑CTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTtct                                                                        >  1:809075/1‑151 (MQ=255)                                                                              tttCGACTTCAAGCGAGACGGTATTACCCGGCCCCT‑GGTCGCGCGGCAGGTTATGAATATTCTGTTCATGCAGGGAAAAA‑CTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACg                                                                   >  1:202781/1‑151 (MQ=255)                                                                               ttcacctccaagcaagACGGTTTTCCCCGGCCCCG‑GGCCGCCGGGGAGTTTAGAAATTTTCTTTTCATGCGGGAAAAAC‑CCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCa                                                                     <  1:981444/140‑1 (MQ=255)                                                                               tccGACTTAAAGGGAGACGGTATTACCCGCCCCCT‑GGCCGCGCGGCAGGTTATAAATATTCTGTTCAGGAAGGAAAAAC‑CCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCAc                                                                    <  2:709128/147‑1 (MQ=255)                                                                                tCGACTTCAAGCGAGACGGTATTACCCGGCCCCT‑GGTCGCGCGGCAGGTTATGAATATTCTGTTCATGCAGGGAAAAA‑CTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGc                                                                  >  1:855687/1‑150 (MQ=255)                                                                                  gACTTCAAGCGAGACGGTATTCCCCGGCCCCT‑GGTCGCGCGGCAGTTTATGAATATTCTGTTCATGCAGGAAAAAC‑CTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCa                                                                 <  2:116059/149‑1 (MQ=255)                                                                                   aCTTCAAGCGAGACGGTATTACCCGGCCCCT‑GGTCGCGCGGCAGGTTATGAATATTCTGTTCATGCAGGGAAAAA‑CTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCt                                                              >  2:709395/1‑151 (MQ=255)                                                                                   aCTTCAAGCAAAACGTTTTTCCCCGGCCCCT‑GGCCGCGCGGCAGGTTATAAATTTTCTGTCCA‑GCAAGGAAAAACCCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCa                                                                 <  2:746128/148‑1 (MQ=255)                                                                                    cTTAAAGGGAGAGGGTTTTCCCCGGCCCCT‑GGTCCCGGGGCAGTTTATAAATTTTTTTTTCAGGAGGGAAAAAA‑CCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCa                                                                 <  1:992957/147‑1 (MQ=255)                                                                                     ttaaaGCGAGACGGTATTCCCCGGCCCCT‑GGTCGCGCGGAAGGTTATAAATATTCTTTTCAGGCAGGAAAAAA‑CTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCt                                                              <  2:938931/146‑1 (MQ=255)                                                                                      tCAAGCGAGACGGTATTACCCGGCCCCT‑GGTCGCGCGGCAGGTTATGAATATTCTGTTCATGCAGGGAAAAA‑CTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGc                                                            >  2:495745/1‑150 (MQ=255)                                                                                      caaaGCGAGAGGTTATCCCCCGGCCCCG‑GGTCGCGCGGCAGGTTATGAATTTTCTGTTCAGGCAGGAAAAAC‑CCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGca                                                           <  2:534397/149‑1 (MQ=255)                                                                                       cAAGCGAGACGGTATTACCCGGCCCCT‑GGTCGCGCGGCAGGTTATGAATATTCTGTTCATGCAGGGAAAAA‑CTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGcac                                                          >  2:123550/1‑151 (MQ=255)                                                                                        aaGCGAGACGGTATTCCCCGGCCCCT‑GGTCGCGCGGCAGTTTATAAATTTTCTGTTCATGCAGGAAAAAC‑CTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTg                                                             <  2:436244/147‑1 (MQ=255)                                                                                        aaGCGAGACGGTATTACCCGGCCCCT‑GGTCGCGCGGCAGGTTATGAATATTCTGTTCATGCAGGGAAAAA‑CTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGcac                                                          >  1:157153/1‑150 (MQ=255)                                                                                            gagaCGGTATTACCCGGCCCCT‑GGTCGCGCGGCAGGTTATGAATATTCTGTTCATGCAGGGAAAAA‑CTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGt                                                     >  1:672370/1‑151 (MQ=255)                                                                                            gagaCGGTATTACCCGGCCCCT‑GGTCGCGCGGCAGGTTATGAATATTCTGTTCATGCAGGGAAAAA‑CTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGt                                                     >  2:775527/1‑151 (MQ=255)                                                                                            gagaCGGTATCCCCCGCCCCCT‑GGTCGCGCGGCGGTTTATAATTTTTCTTTCCATGAAGGAAAAAC‑CTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTg                                                             <  1:607727/143‑1 (MQ=255)                                                                                                 ggTATTACCCGGCCCCT‑GGTCGCGCGGCAGGTTATGAATATTCTGTTCATGCAGGGAAAAA‑CTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCcg                                                >  1:297188/1‑151 (MQ=255)                                                                                                  gTATTACCCGGCCCCT‑GGTCGCGCGGCAGGTTATGAATATTCTGTTCATGCAGGGAAAAA‑CTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCcgc                                               >  1:72177/1‑151 (MQ=255)                                                                                                      tACCCCGCCCCTGGGTCGCGCGGGAGTTTAATAATTTTCTTTTCGGGAAGGTAAAAA‑CTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCcgc                                               <  1:189361/148‑1 (MQ=255)                                                                                                               ccT‑GGTCGCGCGGCAGGTTATGAATATTCTGTTCATGCAGGGAAAAA‑CTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTg                                  >  2:646425/1‑151 (MQ=255)                                                                                                                cT‑GGTCGCGCGGCAGGTTATGAATATTCTGTTCATGCAGGGAAAAA‑CTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGc                                 >  1:704841/1‑151 (MQ=255)                                                                                                                 t‑GGTCGCGCGGCAGGTTATGAATATTCTGTTCATGCAGGGAAAAA‑CTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGc                                 >  1:587018/1‑150 (MQ=255)                                                                                                                       gcgcgGCCGGTTTTTAATTTGCTCTACATGGCGGAAAAAA‑CCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGGTGCAAGgt                            <  1:5430/150‑1 (MQ=255)                                                                                                                              aGGTTATGAATATTCTGTTCATGCAGGGAAAAA‑CTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGTTAAACg                    >  1:907667/1‑151 (MQ=255)                                                                                                                               ggTTATGAATATTCTTTTCATGCAGGAAAAAA‑CTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGTTaaa                      <  2:14778/148‑1 (MQ=255)                                                                                                                                 ttATAAATTTTCTGTTCATGCAGGAAAAAC‑CCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGTTAAACGa                   <  2:500928/149‑1 (MQ=255)                                                                                                                                          ttgcgtgtttggggggAAAAC‑CCCCCCGCCAGTGTAGCGATGTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAGCTTCTCCACGCATCCGCACAATGTGCCCGCGCTCGCCAAGTTGCAAGGTGTTAAACGATGCCACggg          <  1:742285/136‑1 (MQ=255)                                                                                                                                          ttCTTTTCAGGCAGGAAAAAC‑CTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGTTAAACGATGCCAc             <  2:862192/146‑1 (MQ=255)                                                                                                                                          ttCTTTTAATGCAGGAAAAAC‑CCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGTTAAACGATGCCACgg           <  2:788035/148‑1 (MQ=255)                                                                                                                                          ttCTGTTCATGCAGGAAAAAA‑CCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGTTAAACGATGCCACGGGc         <  2:592077/150‑1 (MQ=255)                                                                                                                                              gTTCATGCAGGGAAAAA‑CTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGTTAAACGATGCCACGGGCGAGAc    >  1:694302/1‑151 (MQ=255)                                                                                                                                                tCATGCAGGGAAAAA‑CTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGTTAAACGATGCCACGGGCGAGACtt  >  2:664957/1‑151 (MQ=255)                                                                                                                                                 caggCAGGAAAACC‑CCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGTTAAACGATGCCACGGGCga       <  1:679836/142‑1 (MQ=255)                                                                                                                                                         |                                                                                                                                               ACCAACTGTGCCGCAACCACCTCGGCACTGACGCGCTTTTCACCGACGACAAAAAAGC‑GTTCGGTGATATTTTCACTTTCGACTTCAAGCGAGACGGTATTACCCGGCCCCT‑GGTCGCGCGGCAGGTTATGAATATTCTGTTCATGCAGGGAAAAA‑CTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGTTAAACGATGCCACGGGCGAGACTT  >  NC_000913/3556050‑3556341

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