Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,640,9810CT48.1% 55.1 / 47.5 81F40F (TTC→TTT) dtpBdipeptide and tripeptide permease B
Reads supporting (aligned to +/- strand):  ref base C (10/32);  new base T (39/0);  total (49/32)
Fisher's exact test for biased strand distribution p-value = 4.11e-14
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GGTATTGGAAGGTTAATAAATATGAATACAACAACACCCATGGGGATGCTGCAGCAACCTCGCCCATTTTTCATGATCTTTTTTGTCGAGTTATGGGAGCGATTCGGCTACTACGGCGTGCAGGGCGTACTGGC‑GGTTTTCTTC‑GTTAAACAGCTTGGATTCTCGCAAGAGCAGGC‑TTTTGTC‑ACTTTTGGTGCTTTTGCTGCGCTGGTCTAT‑GGCCTCATTTCCA‑TTGGCGGCTATGTCGGCGACCACCTGCTGGGGACCAAACGCACCATTGTTCTTGGAGCACTTGTGC  >  NC_000913/3640841‑3641132
                                                                                                                                             |                                                                                                                                                            
ggTATTGGAAGGTTAATAAATATGAATACAACAACACCCATGGGGATGCTGCAGCAACCTCGCCCATTTTTCATGATCTTTTTTGTCGAGTTATGGGAGCGATTCGGCTACTACGGCGTGCAGGGCGTACTGGC‑GGGtttttttttttaaa                                                                                                                                                    >  1:443304/1‑143 (MQ=255)
ggTATTGGAAGGTTAATAAATATGAATACAACAACACCCATGGGGATGCTGCAGCAACCTCGCCCATTTTTCATGATCTTTTTTGTCGAGTTATGGGAGCGATTCGGCTACTACGGCGTGCAGGGCGTACTGGC‑GGGTTTTTTT‑GTTtaaa                                                                                                                                                   >  1:937561/1‑150 (MQ=255)
 gTATTGGAAGGTTAATAAATATGAATACAACAACACCCATGGGGATGCTGCAGCAACCTCGCCCATTTTTCATGATCTTTTTTGTCGAGTTATGGGAGCGATTCGGCTACTACGGCGTGCAGGGCGTACTGGCGGGTTTTCTTC‑GTTaaca                                                                                                                                                   >  2:55078/1‑149 (MQ=255)
 gTATTGGAAGGTTAATAAATATGAATACAACAACACCCATGGGGATGCTGCAGCAACCTCGCCCATTTTTCATGATCTTTTTTGTCGAGTTATGGGAGCGATTCGGCTACTACGGCGTGCAGGGCGTACTGGC‑GGTTTTTTTT‑GTTtaaaa                                                                                                                                                  >  1:10721/1‑149 (MQ=255)
 gTATTGGAAGGTTAATAAATATGAATACAACAACACCCATGGGGATGCTGCAGCAACCTCGCCCATTTTTCATGATCTTTTTTGTCGAGTTATGGGAGCGATTCGGCTACTACGGCGTGCAGGGCGTACTGGC‑GGTTTTTTTCTGTTaaaa                                                                                                                                                   >  2:1022441/1‑150 (MQ=255)
 gTATTGGAAGGTTAATAAATATGAATACAACAACACCCATGGGGATGCTGCAGCAACCTCGCCCATTTTTCATGATCTTTTTTGTCGAGTTATGGGAGCGATTCGGCTACTACGGCGTGCAGGGCGTACTGGC‑GGTTTTCTTC‑GTTAAAc                                                                                                                                                   <  1:826364/150‑1 (MQ=255)
 gTATTGGAAGGTTAATAAATATGAATACAACAACACCCATGGGGATGCTGCAGCAACCTCGCCCATTTTTCATGATCTTTTTTGTCGAGTTATGGGAGCGATTCGGCTACTACGGCGTGCAGGGCGTACTGGC‑GGTTTTCTTC‑GTTAAAc                                                                                                                                                   <  1:1073804/150‑1 (MQ=255)
 gTATTGGAAGGTTAATAAATATGAATACAACAACACCCATGGGGATGCTGCAGCAACCTCGCCCATTTTTCATGATCTTTTTTGTCGAGTTATGGGAGCGATTCGGCTACTACGGCGTGCAGGGCGTACTGGC‑GGTTTTCTTC‑GTTAAAc                                                                                                                                                   <  2:413739/150‑1 (MQ=255)
 gTATTGGAAGGTTAATAAATATGAATACAACAACACCCATGGGGATGCTGCAGCAACCTCGCCCATTTTTCATGATCTTTTTTGTCGAGTTATGGGAGCGATTCGGCTACTACGGCGTGCAGGGCGTACTGGC‑GGGTttctttttttaaaa                                                                                                                                                   >  1:852791/1‑142 (MQ=255)
 gTATTGGAAGGTTAATAAATATGAATACAACAACACCCATGGGGATGCTGCAGCAACCTCGCCCATTTTTCATGATCTTTTTTGTCGAGTTATGGGAGCGATTCGGCTACTACGGCGTGCAGGGCGTACTGGC‑GGGTTTTTTT‑GTTtaaaa                                                                                                                                                  >  1:263171/1‑149 (MQ=255)
 gTATTGGAAGGTTAATAAATATGAATACAACAACACCCATGGGGATGCTGCAGCAACCTCGCCCATTTTTCATGATCTTTTTTGTCGAGTTATGGGAGCGATTCGGCTACTACGGCGTGCAGGGCGTACTGGC‑GGGTTTTTTC‑TTTtaaaa                                                                                                                                                  >  1:145028/1‑149 (MQ=255)
 gTATTGGAAGGCTAATAAATATGAATACAACAACACCCATGGGGATGCTGCAGCAACCTCGCCCATTTTTCATGATCTTTTTTGTCGAGTTATGGGAGCGATTCGGCTACTACTGCGTGCAGGGCGTACTGGC‑GGGTttgtttctttataa                                                                                                                                                   >  2:758472/1‑142 (MQ=255)
          ggTTAATAAATATGAATACAACAACACCCATGGGGATGCTGCAGCAACCTCGCCCATTTTTCATGATCTTTTTTGTCGAGTTATGGGAGCGATTCGGCTACTACGTCGTGCAGGGCGTACTGGC‑GGGTTTTTTT‑GTTTAACAACtctgcgt                                                                                                                                         >  2:185783/1‑145 (MQ=255)
          ggTTAATAAATATGAATACAACAACACCCATGGGGATGCTGCAGCAACCTCGCCCATTTTTCATGATCTTTTTTGTCGAGTTATGGGAGCGATTCGGCTACTACGGCGTGCAGGGCGTACTGGC‑GGTTTTTTTT‑GTTTAACAGCTTGGGtt                                                                                                                                         >  1:109290/1‑151 (MQ=255)
          ggTTAATAAATATGAATACAACAACACCCATGGGGATGCTGCAGCAACCTCGCCCATTTTTCATGATCTTTTTTGTCGAGTTATGGGAGCGATTCGGCTACTACGGCGTGCAGGGCGTACTGGC‑GGTTTTTTTT‑GTTTAAAAACtttggtt                                                                                                                                         >  1:1077949/1‑146 (MQ=255)
          ggTTAATAAATATGAATACAACAACACCCATGGGGATGCTGCAGCAACCTCGCCCATTTTTCATGATCTTTTTTGTCGAGTTATGGGAGCGATTCGGCTACTACGGCGTGCAGGGCGTACTGGC‑GGTTTTTTTT‑GTTAAAAAGCtttggtt                                                                                                                                         >  2:531032/1‑146 (MQ=255)
               ataaataTGAATACAACAACACCCATGGGGATGCTGCAGCAACCTCGCCCATTTTTCATGATCTTTTTTGTCGAGTTATGGGAGCGATTCGGCTACTACGGCGTGCAGGGCGTACTGGC‑GGTTTTCTTT‑GTTAAACACCTTTGGTtctcc                                                                                                                                     >  2:111532/1‑149 (MQ=255)
                   atatGAATACAACAACACCCATGGGGATGCTGCAGCAACCTCGCCCATTTTTCATGATCTTTTTTGTCGAGTTATGGGAGCGATTCGGCTACTACGGCGTGCAGGGCGTACTGGG‑GGTTTTCTTT‑GTTAAACAACATGGGTTTTCTCaaaa                                                                                                                                >  1:942939/1‑149 (MQ=255)
                   atatGAATACAACAACACCCATGGGGATGCTGCAGCAACCTCGCCCATTTTTCATGATCTTTTTTGTCGAGTTATGGGAGCGATTCGGCTACTACGGCGTGCAGGGCGTACTGGC‑GGGTTTTTTT‑GTTAAAAAGCTTGGGGTCTCGCAaga                                                                                                                                >  1:633810/1‑151 (MQ=255)
                    tatGAATACAACAACACCCATGGGGATGCTGCAGCAACCTCGCCCATTTTTCATGATCTTTTTTGTCGAGTTATGGGAGCGATTCGGCTACTACGGCGTGCAGGGCGTACTGGC‑GGTTTTTTTT‑TTTTAAAAACTTGGGTTTTTGCAagag                                                                                                                               >  1:1178094/1‑151 (MQ=255)
                    tatGAATACAACAACACCCATGGGGATGCTGCAGCAACCTCGCCCATTTTTCATGATCTTTTTTGTCGAGTTATGGGAGCGATTCGGCTACTACGGCGTGCAGGGCGTACTGGC‑GGTTTTTTTT‑CTTTAAAAACTTTGGTTTTCCCAagag                                                                                                                               >  1:824000/1‑151 (MQ=255)
                    tatGAATACAACAACACCCATGGGGATGCTGCAGCAACCTCGCCCATTTTTCATGATCTTTTTTGTCGAGTTATGGGAGCGATTCGGCTACTACGGCGTGCAGGGCGTACTGGC‑GGGTTTTTTC‑TTTTAAAAACTTggggtttcgccaaag                                                                                                                               >  1:503264/1‑138 (MQ=255)
                    tatGAATACAACAACACCCATGGGGATGCTGCAGCAACCTCGCCCATTTTTCATGATCTTTTTTGTCGAGTTATGGGAGCGATTCGGCTACTACGGCGCGCAGGGCGGGCTGGC‑GGCTTTTTTC‑TTTAAAAAGCCTgggatttcccaaaag                                                                                                                               >  1:269083/1‑138 (MQ=255)
                       gAATACAACAACACCCATGGGGATGCTGCAGCAACCTCGCCCATTTTTCATGATCTTTTTTGTCGAGTTATGGGAGCGATTCGGCTACTACGGCGTGCAGGGCGTACTGGCGGGTTTTTTTG‑TTTAAAAACCTTGGGtttctgaagagccg                                                                                                                             >  2:1102057/1‑139 (MQ=255)
                       gAATACAACAACACCCATGGGGATGCTGCAGCAACCTCGCCCATTTTTCATGATCTTTTTTGTCGAGTTATGGGAGCGATTCGGCTACTACGGCGTGCAGGGCGTACTGGC‑GGTTTTTTTC‑GTTTAAAAACTTTGATTTTCGCCAGAGCAg                                                                                                                            >  2:992874/1‑151 (MQ=255)
                       gAATACAACAACACCCATGGGGATGCTGCAGCAACCTCGCCCATTTTTCATGATCTTTTTTGTCGAGTTATGGGAGCGATTCGGCTACTACGGCGTGCAGGGCGTACTGGC‑GGTTTTCTTT‑GTTTAACAGGTTGGGTTTTTGCAAAAGCAg                                                                                                                            >  1:632058/1‑151 (MQ=255)
                       gAATACAACAACACCCATGGGGATGCTGCAGCAACCTCGCCCATTTTTCATGATCTTTTTTGTCGAGTTATGGGAGCGATTCGGCTACTACGGCGTGCAGGGCGTACTGGC‑GGGTTTTTTT‑GTTTAAAAGCTTGGGTTTTCGCAAGAGCAg                                                                                                                            >  1:878801/1‑151 (MQ=255)
                       gAATACAACAACACCCATGGGGATGCTGCAGCAACCTCGCCCATTTTTCATGATCTTTTTTGTCGAGTTATGGGAGCGATTCGGCTACTACGGCGTGCAGGGCGTACTGGC‑GGGTTTTTTT‑GTTAAAAACCCTGGAATTCCTCAAgaggag                                                                                                                            >  1:456763/1‑151 (MQ=255)
                       gAATACAACAACACCCATGGGGATGCTGCAGCAACCTCGCCCATTTTTCATGATATTTTTTGGCGAGTTATGGGAGCGATTCGGCTACTACGGCGCGCAGGGCGTACCGGC‑GGGTTTCCTT‑CTTTAACAGCGTTGGTGTTCGCAAGAGCAg                                                                                                                            >  1:559816/1‑151 (MQ=255)
                            caccaacaCCCATGGGGATGCTGCAGCAACCTCGCCCATTTTTCATGATCTTTTTTGTCGAGTTATGGGAGCGATTCGGCTACTACGGCGTGCAGGGCGTACTGGC‑GGTTTTCTTC‑GTTTAACAGCTTGGGTTCTCGCACGAGAAGGC‑ttt                                                                                                                      >  2:5220/4‑151 (MQ=255)
                            caacaacaCCCATGGGGATGCTGCAGCAACCTCGCCCATTTTTCATGATCTTTTTTGTCGAGTTATGGGAGCGATTCGGCTACTACGGCGTGCAGGGCGTACTGGC‑GGTTTTTTTT‑TTTTAAAAAGTTGGGTTTTCCCCAAAGCAGGG‑ttt                                                                                                                      >  1:88931/1‑151 (MQ=255)
                            caacaacaCCCATGGGGATGCTGCAGCAACCTCGCCCATTTTTCATGATCTTTTTTGTCGAGTTATGGGAGCGATTCGGCTACTACGGCGTGCAGGGCGTACTGGC‑GGGTTTTTTT‑GTTAAAAAACCTGGGTTTTCGGAAAAGGAGGG‑ttt                                                                                                                      >  1:294578/1‑151 (MQ=255)
                            caacaacaCCCATGGGGATGCTGCAGCAACCTCGCCCATTTTTCATGATCTTTTTTGTCGAGTTATGGGAGCGATTCGGCTACTACGGCGTGCAGGGCGTACTGGC‑GGGTTTTTTT‑CTTAAACAGCCTTGGTTTTCGCCAAAGCAGGC‑gtt                                                                                                                      >  2:373224/1‑151 (MQ=255)
                              acaacaCCCATGGGGATGCTGCAGCAACCTCGCCCATTTTTCATGATCTTTTTTGTCGAGTTATGGGAGCGATTCGGCTACTACGGCGTGCAGGGCGTACTGGC‑GGGTTTTTTT‑GTTTAAAAACTTTGGTTTTCGCAAGAGCAGGG‑ttttt                                                                                                                    >  1:843442/1‑150 (MQ=255)
                               caacaCCCATGGGGATGCTGCAGCAACCTCGCCCATTTTTCATGATCTTTTTTGTCGAGTTATGGGAGCGATTCGGCTACTACGGCGTGCAGGGCGTACTGGC‑GGTTTTCTTC‑GTTAAACAGCTTGGATTCTCGCAAGAGCAGGC‑TTTTGt                                                                                                                   <  1:830632/151‑1 (MQ=255)
                                aacaCCCATGGGGATGCTGCAGCAACCTCGCCCATTTTTCATGATCTTTTTTGTCGAGTTATGGGAGCGATTCGGCTACTACGGCGTGCAGGGCGTACTGGC‑GGTTTTCTTC‑GTTAAACAGCTTGGATTCTCGCAAGAGCAGGC‑TTTTGTc                                                                                                                  <  2:1051207/151‑1 (MQ=255)
                                  cacCCATGGGGATGCTGCAGCAACCTCGCCCATTTTTCATGATCTTTTTTGTCGAGTTATGGGAGCGATTCGGCTACTACGGCGTGCAGGGCGTACTGGC‑GGTTTTCTTC‑GTTAAACAGCTTGGATTCTCGCAAGAGCAGGC‑TTTCGTC‑Ac                                                                                                               <  2:280126/151‑1 (MQ=255)
                                    cccATGGGGATGCTGCAGCAACCTCGCCCATTTTTCATGATCTTTTTTGTCGAGTTATGGGAGCGATTCGGCTACTACGGCGTGCAGGGCGTACTGGC‑GGTTTTTTTT‑GTTAAAAAGCTTGGGTTTTCGCAAAAGCAGGG‑TTTTTTC‑CCtt                                                                                                             >  1:203412/1‑151 (MQ=255)
                                    cccATGGGGATGCTGCAGCAACCTCGCCCATTTTTCATGATCTTTTTTGTCGAGTTATGGGAGCGATTCGGCTACTACGGCGTGCAGGGCGTACTGGC‑GGTTTTTTTC‑GTTTAAAAACTTTGGTTCTCGCAAAAGCCGGC‑gtttttttctt                                                                                                              >  2:143940/1‑144 (MQ=255)
                                    cccATGGGGATGCTGCAGCAACCTCGCCCATTTTTCATGATCTTTTTTGTCGAGTTATGGGAGCGATTCGGCTACTACGGCGTGCAGGGCGTACTGGC‑GGGTTTTTTT‑TTTTTACAACTTTGATTGTCCGCAGAGCAggggttttttcctt                                                                                                               >  2:798162/1‑139 (MQ=255)
                                     ccATGGGGATGCTGCAGCAACCTCGCCCATTTTTCATGATCTTTTTTGTCGAGTTATGGGAGCGATTCGGCTACTACGGCGTGCAGGGCGTACTGGC‑GGTTTTTTTC‑TTTTAAAAACTTTGGTTTTCTCCAGAGCGGGG‑tttttttccttt                                                                                                             >  1:674485/1‑143 (MQ=255)
                                      cATGGGGATGCTGCAGCAACCTCGCCCATTTTTCATGATCTTTTTTGTCGAGTTATGGGAGCGATTCGGCTACTACGGCGTGCAGGGCGTACTGGC‑GGTTTTTTTC‑GTTCAACAGCTTTGGTTCTCCCCAGAGCCGGG‑CTTTGTC‑TCcttt                                                                                                           >  2:584509/1‑151 (MQ=255)
                                      cATGGGGATGCTGCAGCAACCTCGCCCATTTTTCATGATCTTTTTTGTCGAGTTATGGGAGCGATTCGGCTACTACGGCGTGCAGGGCGTACTGGC‑GGTTTTCTTC‑GTTAAAc                                                                                                                                                   <  1:764686/113‑1 (MQ=255)
                                                  gcagcaACCTCGCCCATTTTTCATGATCTTTTTTGTCGAGTTATGGGAGCGATTCGGCTACTACGGCTTGCAGGGCGTACTGGC‑GGTTTTTTTT‑GTTTAACAGCCTGGGTTCTCTCAACACCAGGC‑TTTTTTT‑CCTTTTGGTTGttttttt                                                                                               >  1:477150/1‑148 (MQ=255)
                                                       aaCCTCGCCCATTTTTCATGATCTTTTTTGTCGAGTTATGGGAGCGATTCGGCTACTACGGCGTGCAGGGCGTACTGGCGGGTTTTTTTC‑TTTAAA‑AACTTGGGTTTTCGCAAGAGAAGGC‑TTTTTTC‑ACTTTTTGGGCTTTttctggg                                                                                            >  1:1013913/1‑147 (MQ=255)
                                                           tCGCCCATTTTTCATGATCTTTTTTGTCGAGTTATGGGAGCGATTCGGCTACTACGGCGTGCAGGGCGTACTGGC‑GGTTTTCTTC‑GTTAAACAGCTTGGATTCTCGCAAGAGCAGGC‑TTTTGTC‑ACTTTTGGTGCTTTTGCTGCGCTGGTc                                                                                      <  1:993305/151‑1 (MQ=255)
                                                                cATTTTTCATGATCTTTTTTGTCGAGTTATGGGAGCGATTCGGCTACTACGGCGTGCAGGGCGTACTGGC‑GGTTTTCTTC‑GTTAAACAGCTTGGATTCTCGCAAGAGCAGGC‑TTTTGTC‑ACTTTTGGTGCTTTTGCTGCGCTGGTCTAT‑gg                                                                                <  2:277293/151‑1 (MQ=255)
                                                                cATTTTTCATGATCTTTTTTGTCGAGTTATGGGAGCGATTCGGCTACTACGGCGTGCAGGGCGTACTGGC‑GGTTTTCTTC‑GTTAAACAGCTTGGATTCTCGCAAGAGCAGGC‑TTTTGTC‑ACTTTTGGTGCTTTTGCTGCGCTGGTCTAT‑gg                                                                                <  1:906675/151‑1 (MQ=255)
                                                                  tttttCATGATCTTTTTTGTCGAGTTATGGGAGCGATTCGGCTACTACGGCGTGCAGGGCGTACTGGC‑GGTTTTTTTT‑GTTTAACCGCTTTGGTTTTCCCCAGAGCAGGC‑CTTTTTT‑ACTTTTTTTGCTTTTTCTTCGCTGGTCTTT‑GGc                                                                               >  1:193769/1‑150 (MQ=255)
                                                                  tttttCATGATCTTTTTTGTCGAGTTATGGGAGCGATTCGGCTACTACGGCGTGCAGGGCGTACTGGC‑GGTTTTCTTC‑GTTTAAAAACTTGGGTTTTCGCAAGAGCAGGG‑TTTTGTT‑ACTTTTGGTGGTTTTTCTGCGCTGGGCTAT‑TGc                                                                               >  1:694017/1‑150 (MQ=255)
                                                                  tttttCATGATCTTTTTTGTCGAGTTATGGGAGCGATTCGGCTACTACGGCGTGCAGGGCGTACTGGC‑GGTTTTCTTC‑GTTAAACAGCTTGGATTCTCGCAAGAGCAGGC‑TTTTGTC‑ACTTTTGGTGCTTTTGCTGCGCTGGTCTAT‑GGcc                                                                              <  1:340904/151‑1 (MQ=255)
                                                                        aTGATCTTTTTTGTCGAGTTATGGGAGCGATTCGGCTACTACGGCGTGCAGGGCGTACTGGC‑GGTTTTCTTC‑GTCAAACAGCTTGGATTCTCGCAAGAGCAGGC‑TTTTGTC‑ACTTTTGGTGCTTTTGCTGCGCTGGTCTAT‑GGCCTCAttt                                                                        <  1:735698/151‑1 (MQ=255)
                                                                          gATCTTTTTTGTCGAGTTATGGGAGCGATTCGGCTACTACGGCGTGCAGGGCGTACTGGC‑GGTTTTTTTT‑GTTTAACAACTTGGGTTTTCGCAAGAGCAGGG‑TTTTTTT‑ACTTTTTGGGGTTTTTCTGCGCTGGGCTTT‑TGCCCCAtttt                                                                       >  1:795347/1‑149 (MQ=255)
                                                                           aTCTTTTTTGTCGAGTTATGGGAGCGATTCGGCTACTACGGCGTGCAGGGCGTACTGGC‑GGTTTTTTTT‑GTTTAAAAGCTTGGGTTTTCGCCAAAGGAGGC‑TTTTTTT‑ACTTTTTGGGGTTTTTCTGGGCTGGGCTTT‑GGGCTCCTTTTCa                                                                     >  1:463430/1‑151 (MQ=255)
                                                                                 tttGTCGAGTTATGGGAGCGATTCGGCTACTACGGCGTGCAGGGCGTACTGGC‑GGTTTTCTTC‑GTTAAACAGCTTGGATTCTCGCAAGAGCAGGC‑TTTTGTC‑ACTTTTGGTGCTTTTGCTGCGCTGGTCTAT‑GGCCTCATTTCCA‑TTggcg                                                              <  2:284746/151‑1 (MQ=255)
                                                                                 tttGTCGAGTTATGGGAGCGATTCGGCTACTACGGCGTGCAGGGCGTACTGGC‑GGTTTTCTTC‑GTTAAACAGCTTGGATTCTCGCAAGAGCAGGC‑TTTTGTC‑ACTTTTGGTGCTTTTGCTGCGCTGGTCTAT‑GGCCTCATTTCCA‑TTggcg                                                              <  2:1084736/151‑1 (MQ=255)
                                                                                  ttGTCGAGTTATGGGAGCGATTCGGCTACTACGGCGTGCAGGGCGTACTGGC‑GGTTTTCTTC‑GTTAAACAGCTTGGATTCTCGCAAGAGCAGGC‑TTTTGTC‑ACTTTTGGTGCTTTTGCTGCGCTGGTCTAT‑GGCCTCATTTCCA‑TTggcgg                                                             <  1:949250/151‑1 (MQ=255)
                                                                                        aGTTATGGGAGCGATTCGGCTACTACGGCGTGCAGGGCGTACTGGCGGGTTTTTTTC‑GTT‑AACAGGTTTGATTTTCGCAAGAGCAGGCGTTTTTTC‑CCTTTT‑GCGCTTTTTCTGCGCTTGTCTAT‑GGGCTCATTTTCA‑TTtgcgggggttt                                                       >  2:880363/1‑145 (MQ=255)
                                                                                        aGTTATGGGAGCGATTCGGCTACTACGGCGTGCAGGGCGTACTGGC‑GGTTTTTTTT‑GTTCAACAGCCTGGGATTTCGCAAAAGCAGGC‑TTTTTTCTTCTTTTGGT‑CTTTTTCTGCGCTTGGGTAT‑GGCGTCATTTTCA‑TTggcggggtttt                                                       >  2:780341/1‑145 (MQ=255)
                                                                                          ttATGGGAGCGATTCGGCTACTACGGCGTGCAGGGCGTACTGGC‑GGTTTTCTTC‑GTTAAACAGCTTGGATTCTCGCAAGAGCAGGC‑TTTTGTC‑ACTTTTGGTGCTTTTGCTGCGCTGGTCTAT‑GGCCTCATTTCCA‑TTGGCGGCTATGt                                                       <  1:409766/149‑1 (MQ=255)
                                                                                            aTGGGAGCGATTCGGCTACTACGGCGTGCAGGGCGTACTGGC‑GGTTTTCTTC‑GTTAAACAGCTTGGATTCTCGCAAGAGCAGGC‑TTTTGTC‑ACTTTTGGTGCTTTTGCTGCGCTGGTCTAT‑GGCCTCATTTCCA‑TTGGCGGCTATGTCGGc                                                   <  1:1057363/151‑1 (MQ=255)
                                                                                            aTGGGAGCGATTCGGCTACTACGGCGTGCAGGGCGTACTGGC‑GGGTTTTTTT‑CTTAAACAACTTGGGTTTTCGCCAGAGCAGGC‑TTTTTTC‑CCTTTTTGGGCTTTTTCTTCGCTGGTCTCT‑GGCCTCCTTTTCATTTGGGGGCTTTTTCgg                                                    >  1:1037122/1‑151 (MQ=255)
                                                                                                gAGCGATTCGGCTACTACGGCGTGCAGGGCGTACTGGC‑GGTTTTCTTC‑GTTAAACAGCTTGGATTCTCGCAAGAGCAGGC‑TTTTGTC‑ACTTTTGGTGCTTTTGCTGCGCTGGTCTAT‑GGCCTCATTTCCA‑TTGGCGGCTATGTCGGCGac                                                <  2:953235/150‑1 (MQ=255)
                                                                                                         ggCTACTACGGCGTGCAGGGCGTACTGGC‑GGTTTTTTTT‑GTTAAACAGCTTGGATTCTTGCCAGGGCAGGG‑TTTTGTT‑ACTTTTTGGGGTTTTTGTGGGCTGGGCTATGGGCCCCTTTTCCT‑TTGG‑GGGTATTTTGGGGACCACCCTCTgg                                      >  1:901401/1‑151 (MQ=255)
                                                                                                         ggCTACTACGGCGTGCAGGGCGTACTGGC‑GGTTTTCTTC‑GTTTAACCACTTTGGTTTTTGCAAGAGGAGGC‑TTTTGTT‑ACTTTTTGTGCTTTTGTTGCGCTGGTCTAT‑GGCCTCATTTCCA‑TTGGCGGCTATGTTGGTGagcacccta                                         >  2:906148/1‑145 (MQ=255)
                                                                                                             actacGGCGTGCAGGGCGTACTGGC‑GGTTTTCTTC‑GTTAAACAGCTTGGATTCTCGCAAGAGCAGGC‑TTTTGTC‑ACTTTTGGTGCTTTTGCTGCGCTGGTCTAT‑GGCCTCATTTCCA‑TTGGCGGCTATGTCGGCGACCACCTGCTGGGGAc                                  <  2:937167/151‑1 (MQ=255)
                                                                                                             actacGGCGTGCAGGGCGTACTGGC‑GGGTTTTTTC‑GTTAAACAGCTTGGATTCTCGCAAGAGCAGGG‑TTTTTTC‑ACTTTTGGGGCTTTTTCTGCGCTGGGCTAT‑GGCCTCCTTTCCA‑TTTGGGGGTATTTTGGGGgccaccccccgggggg                                  >  1:993964/1‑141 (MQ=255)
                                                                                                              ctgcGGCGTGCAGGGCGTACTGGC‑GGTTTTCTTC‑GTTAAACAGCTTGGATTCTCGCAAGAGCAGGC‑TTTTGTC‑ACTTTTGGTGCTTTTGCTGCGCTGGTCTAT‑GGCCTCATTTCCA‑TTGGCGGCTATGTCGGCGACCACCTGCTGGGGAcc                                 <  2:108891/148‑1 (MQ=255)
                                                                                                              ctacGGCGTGCAGGGCGTACTGGC‑GGTTTTCTTC‑GTTAAACAGCTTGGATTCTCGCAAGAGCAGGC‑TTTTGTC‑ACTTTTGGTGCTTTTGCTGCGCTGGTCTATGGGCCTCATTTCCA‑TTGGCGGCTATGTCGGCGACCACCTGCTGGGGAc                                  <  1:584905/151‑1 (MQ=255)
                                                                                                              ctacGGCGTGCAGGGCGTACTGGC‑GGTTTTCTTC‑GTTAAACAGCTTGGATTCTCGCAAGAGCAGGC‑TTTTGTC‑ACTTTTGGTGCTTTTGCTGCGCTGGTCTAT‑GGCCTCATTTCCA‑TTGGCGGCTATGTCGGCGACCACCTGCTGGGGAcc                                 <  2:1058747/151‑1 (MQ=255)
                                                                                                                  ggCGTGCAGGGCGTACTGGC‑GGTTTTTTTC‑GTTAAAAAGCTTGGGTTTTTGCAAAGGCAGGC‑GTTTTTC‑CCTTTTGGTGTTTTTTCTGCGCTGGGCTAT‑GGGCTCCTTTTCA‑TTGGCGGCTATTTCTGCGGCCCCCCCCTggggggccaaa                             >  1:34919/1‑144 (MQ=255)
                                                                                                                   gCGTGCAGGGCGTACTGGC‑GGTTTTCTTC‑GTTAAACAGCTTGGATTCTCGCAAGAGCAGGC‑TTTTGTC‑ACTTTTGGTGCTTTTGCTGCGCTGGTCTAT‑GGCCTCATTTCCA‑TTGGCGGCTATGTCGGCGACCACCTGCTGGGGACCAAAc                             <  1:313239/150‑1 (MQ=255)
                                                                                                                     gTGCAGGGCGTACTGGC‑GGTTTTTTTT‑GGTAAACAGCTTGGGTTTTCCCAAGAGCAGGG‑TTTTTTC‑CCTTTTTGTGCTTTTTTTGCGCCGGTCTAT‑GGCCCCCTTTCCA‑TTTGCGGGTTTTTCGGCGGCCCCCCCCCTGGGGCCAAACGc                           >  1:557192/1‑150 (MQ=255)
                                                                                                                             cGTACTGGC‑GGTTTTCTTC‑GTTAAACAGCTTGGATTCTCGCAAGAGCAGGC‑TTTTGTC‑ACTTTTGGTGCTTTTGCTGCGCTGGTCTAT‑GGCCTCATTTCCA‑TTGGCGGCTATGTCGGCGACCACCTGCTGGGGACCAAACGCACCATTGt                   <  1:579805/150‑1 (MQ=255)
                                                                                                                               tACTGGC‑GGTTTTCTTC‑GTTAAACAGCTTGGATTCTCGCAAGAGCAGGC‑TTTTGTC‑ACTTTTGGTGCTTTTGCTGCGCTGGTCTAT‑GGCCTCATTTCCA‑TTGGCGGCTATGTCGGCGACCACCTGCTGGGGACCAAACGCACCATTGTTc                 <  2:294488/150‑1 (MQ=255)
                                                                                                                               tACTGGC‑GGTTTTCTTC‑GTTAAACAGCTTGGATTCTCGCAAGAGCAGGC‑TTTTGTC‑ACTTTTGGTGCTTTTGCTGCGCTGGTCTAT‑GGCCTCATTTCCA‑TTGGCGGCTATGTCGGCGACCACCTGCTGGGGACCAAACGCACCATTGTTc                 <  1:410889/150‑1 (MQ=255)
                                                                                                                               tACTGGC‑GGTTTTCTTC‑GTTAAACAGCTTGGATTCTCGCAAGAGCAGGC‑TTTTGTC‑ACTTTTGGTGCTTTTGCTGCGCTGGTCTAT‑GGCCTCATTTCCA‑TTGGCGGCTATGTCGGCGACCACCTGCTGGGGACCAAACGCACCATTGTTc                 <  2:96569/150‑1 (MQ=255)
                                                                                                                               tACTGGC‑GGTTTTCTTC‑GTTAAACAGCTTGGATTCTCGCAAGAGCAGGC‑TTTTGTC‑ACTTTTGGTGCTTTTGCTGCGCTGGTCTAT‑GGCCTCATTTCCA‑TTGGCGGCTATGTCGGCGACCACCTGCTGGGGACCAAACGCACCATTGTTc                 <  1:845267/150‑1 (MQ=255)
                                                                                                                               tACTGGC‑GGTTTTCTTC‑GTTAAACAGCTTGGATTCTCGCAAGAGCAGGC‑TTTTGTC‑ACTTTTGGTGCTTTTGCTGCGCTGGTCTAT‑GGCCTAATTTCCG‑TTGGCGGCTATGTCGGCGACCACCTGCTGGGGACCAAACGCACCATTGTTc                 <  1:268745/150‑1 (MQ=255)
                                                                                                                                   ggC‑GGTTTTCTTC‑GTTAAACAGCTTGGATTCTCGCAAGAGCAGGC‑TTTTGTC‑ACTTTTGGTGCTTTTGCTGCGCTGGTCTAT‑GGCCTCATTTCCA‑TTGGCGGCTATGTCGGCGACCACCTGCTGGGGACCAAACGCACCATTGTTCTTGGa            >  1:1136406/1‑151 (MQ=255)
                                                                                                                                   ggC‑GGTTTTCTTC‑GTTAAACAGCTTGGATTCTCGCAAGAGCAGGC‑TTTTGTC‑ACTTTTGGTGCTTTTGCTGCGCTGGTCTAT‑GGCCTCATTTCCA‑TTGGCGGCTATGTCGGCGAACACCTGCTGGGGACCAAAAGCAACATTGTTCTTggg            >  1:813919/1‑150 (MQ=255)
                                                                                                                                    gC‑GGTTTTCTTC‑GTTAAACAGCTTGGATTCTCGCAAGAGCAGGC‑TTTTGTC‑ACTTTTGGTGCTTTTGCTGCGCTGGTCTAT‑GGCCTCATTTCCA‑TTGGCGGCTATGTCGGCGACCACCTGCTGGGGACCAAACGCACCATTGTTCTTGGa            <  2:535972/150‑1 (MQ=255)
                                                                                                                                    gC‑GGTTTTCTTC‑GTTAAACAGCTTGGATTCTCGCAAGAGCAGGC‑TTTTGTC‑ACTTTTGGTGCTTTTGCTGCGCTGGTCTAT‑GGCCTCATTTCCA‑TTGGCGGCTATGTCGGCGACCACCTGCTGGGGACCAAACGCACCATTGTTCTTGGa            <  2:787247/150‑1 (MQ=255)
                                                                                                                                         ttttCTTC‑GTTAAACAGCTTGGATTCTCGCAAGAGCAGGC‑TTTTGTC‑ACTTTTGGTGCTTTTGCTGCGCTGGTCTAT‑GGCCTCATTTCCA‑TTGGCGGCTATGTCGGCGACCACCTGCTGGGGACCAAACGCACCATTGTTCTTGGAGCAc        <  2:1146107/150‑1 (MQ=255)
                                                                                                                                             cttc‑GTTAAACAGCTTGGATTCTCGCAAGAGCAGGC‑TTTTGTC‑ACTTTTGGTGCTTTTGCTGCGCTGGTCTAT‑GG‑CTCATTTCCA‑TTGGCGGCTATGTCGGCGACCACCTGCTGGGGACCAAACGCACCATTGTTCTTGGAGCACTTGTGc  <  2:109261/151‑1 (MQ=255)
                                                                                                                                             |                                                                                                                                                            
GGTATTGGAAGGTTAATAAATATGAATACAACAACACCCATGGGGATGCTGCAGCAACCTCGCCCATTTTTCATGATCTTTTTTGTCGAGTTATGGGAGCGATTCGGCTACTACGGCGTGCAGGGCGTACTGGC‑GGTTTTCTTC‑GTTAAACAGCTTGGATTCTCGCAAGAGCAGGC‑TTTTGTC‑ACTTTTGGTGCTTTTGCTGCGCTGGTCTAT‑GGCCTCATTTCCA‑TTGGCGGCTATGTCGGCGACCACCTGCTGGGGACCAAACGCACCATTGTTCTTGGAGCACTTGTGC  >  NC_000913/3640841‑3641132

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