Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,640,981	0	C	T	46.7%	47.6 / 37.6	60	F40F (TTC→TTT)	dtpB	dipeptide and tripeptide permease B
Reads supporting (aligned to +/- strand):  ref base C (8/24);  new base T (28/0);  total (36/24)
Fisher's exact test for biased strand distribution p-value = 3.80e-10
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GTATTGGAAGGTTAATAAATATGAATACAACAACACCCATGGGGATGCTGCAGCAACCTCGCCCATTTTTCATGATCTTTTTTGTCGAGTTATGGGAGCGATTCGGCTACTACGGCGTGCAGGGCGTACTGGC‑GGTTTTCTTCG‑TTAAACAGCTTGGATTCTCGCAAGA‑GCA‑GGCTTTTGTCACTTTTGGTGCTTTTGCTGCGCTGGTCTA‑TGGCCTCATTTCCATTGGCGGCTATGTCGGCGACCACCTGCTGGGGACCAAACGCACCATTGTTCTTGGAGCACTTG  >  NC_000913/3640842‑3641129                                                                                                                                             |                                                                                                                                                         gTATTGGAAGGTTAATAAATATGAATACAACAACACCCATGGGGATGCTGCAGCAACCTCGCCCATTTTTTATGATCTTTTTTGTCGAGTTATGGGAGCGATTCGGCTACTACGGCGTGCAGGGCGTACTGGT‑GGGttttttttcttcaca                                                                                                                                               >  2:1037335/1‑142 (MQ=255) gTATTGGAAGGTTAATAAATATGAATACAACAACACCCATGGGGATGCTGCAGCAACCTCGCCCATTTTTCATGATCTTTTTTGTCGAGTTATGGGAGCGATTCGGCTACTACGGCGTGCAGGGCGTACTGGC‑GGTttttttttgttaaaa                                                                                                                                               >  2:732731/1‑142 (MQ=255)          ggTTAATAAATATGAATACAACAACACCCATGGGGATGCTGCAGCAACCTCGCCCATTTTTCATGATCTTTTTTGTCGAGTTATGGGAGCGATTCGGCTACTACGGCGTGCAGGGCGTACTGGC‑GGGTTTTTTTT‑TTTAACaaccttggtt                                                                                                                                     >  1:484045/1‑142 (MQ=255)           gTTAATAAATATGAATACAACAACACCCATGGGGATGCTGCAGCAACCTCGCCCATTTTTCATGATCTTTTTTGTCGAGTTATGGGAGCGATTCGGCTACTACGGCGTGCATGGCGTACTGGC‑GGGGTTTTTTG‑TTTAAAAGCTTGGGttt                                                                                                                                    >  2:556970/1‑150 (MQ=255)           gTTAATAAATATGAATACAACAACACCCATGGGGATGCTGCAGCAACCTCGCCCATTTTTCATGATCTTTTTTGTCGAGTTATGGGAGCGATTCGGCTACTACGGCGTGCAGGGCGTACTGGC‑GGTTTTCTTCG‑TTAAACAGCTTGGAtt                                                                                                                                     <  1:931415/150‑1 (MQ=255)           gTTAATAAATATGAATACAACAACACCCATGGGGATGCTGCAGCAACCTCGCCCATTTTTCATGATCTTTTTTGTCGAGTTATGGGAGCGATTCGGCTACTACGGCGTGCAGGGCGTACTGGC‑GGTTTTCTTCG‑TTAAACAGCTTGGAtt                                                                                                                                     <  1:338966/150‑1 (MQ=255)           gTTAATAAATATGAATACAACAACACCCATGGGGATGCTGCAGCAACCTCGCCCATTTTTCATGATCTTTTTTGTCGAGTTATGGGAGCGATTCGGCTACTACGGCGTGCAGGGCGTACTGGC‑GGTGTTTTTCG‑TTTAACAGCTTGGCttt                                                                                                                                    >  1:847187/1‑150 (MQ=255)           gTTAATAAATATGAATACAACAACACCCATGGGGATGCTGCAGCAACCTCGCCCATTTTTCATGATCTTTTTTGTCGAGTTATGGGAGCGATTCGGCTACTACGGCGTGCAGGGCGTACTGGC‑GGGTTTTTTTT‑TTTAAAAACTTGGGttt                                                                                                                                    >  2:849636/1‑150 (MQ=255)           gTTAATAAATATGAATACAACAACACCCATGGGGATGCTGCAGCAACCTCGCCCATTTTTCATGATCTTTTTTGTCGAGTTATGGGAGCGATTCGGCTACTACGGCGTGCAGGGCGTACTGGC‑GGGTTTTTTTG‑TTTAAAAACtttgggtt                                                                                                                                    >  1:364555/1‑145 (MQ=255)                   atatGAATACAACAACACCCATGGGGATGCTGCAGCAACCTCGCCCATTTTTCATGATCTTTTTTGTCGAGTTATGGGAGCGATTCGGCTACTACGGCGTGCAGGGCGTACTGGC‑GGTTTTTTTTT‑TTTAAAAGCtttggttttctcgaga                                                                                                                            >  2:1053558/1‑137 (MQ=255)                    tatGAATACAACAACACCCATGGGGATGCTGCAGCAACCTCGCCCATTTTTCATGATCTTTTTTGTCGAGTTATGGGAGCGATTCGGCTACTACGGCGTGCAGGGCGTACTGGC‑GGTTTTTTTTT‑TTTAAACGCTTTGATTTTCGCAaga‑g                                                                                                                          >  2:954660/1‑151 (MQ=255)                    tatGAATACAACAACACCCATGGGGATGCTGCAGCAACCTCGCCCATTTTTCATGATCTTTTTTGTCGAGTTATGGGAGCGATTCGGCTACTACGGCGTGCAGGGCGTACTGGC‑GGTTTTTTTTGTTTAAACA‑CTTGGATTTTCGCAggaa                                                                                                                           >  1:124335/1‑150 (MQ=255)                    tatGAATACAACAACACCCATGGGGATGCTGCAGCAACCTCGCCCATTTTTCATGATCTTTTTTGTCGAGTTATGGGAGCGATTCGGCTACTACGGCGTGCAGGGCGTACTGGC‑GGTTTTCTTCG‑TTAAACAGCTTGGATTCTCGCAaga‑g                                                                                                                          <  1:28813/151‑1 (MQ=255)                       gAATACAACAACACCCATGGGGATGCTGCAGCAACCTCGCCCATTTTTCATGATCTTTTTTGTCGAGTTATGGGAGCGATTCGGCTACTACGGCGTGCAGGGCGTACTGGC‑GGTTTTTTTTT‑TTAAAAAGCTTGGGTTTTTGCCAGA‑GCA‑g                                                                                                                      >  2:171685/1‑151 (MQ=255)                       gAATACAACAACACCCATGGGGATGCTGCAGCAACCTCGCCCATTTTTCATGATCTTTTTTGTCGAGTTATGGGAGCGATTCGGCTACTACGGCGTGCAGGGCGTACTGGC‑GGTTTTTTTCG‑TTTAACAGCCTGGGTTCTCGCCAGA‑GCA‑g                                                                                                                      >  1:733450/1‑151 (MQ=255)                       gAATACAACAACACCCATGGGGATGCTGCAGCAACCTCGCCCATTTTTCATGATCTTTTTTGTCGAGTTATGGGAGCGATTCGGCTACTACGGCGTGCAGGGCGTACTGGC‑GGGTTTTTTTC‑TTAAAAAAATTGGGGTCTCGCCAGA‑GCA‑g                                                                                                                      >  2:321436/1‑151 (MQ=255)                         aTACAACAACACCCATGGGGATGCTGCAGCAACCTCGCCCATTTTTCATGATCTTTTTTGTCGAGTTATGGGAGCGATTCGGCTACTACGGCGTGCAGGGCGTACTGGC‑GGTTTTCTTCG‑TTAAACAGCTTGGATTCTCGCAAGA‑GCA‑GGc                                                                                                                    <  2:90154/151‑1 (MQ=255)                            caacaacaCCCATGGGGATGCTGCAGCAACCTCGCCCATTTTTCATGATCTTTTTTGTCGAGTTATGGGAGCGATTCGGCTACTACGGCGTGCAGGGCGTACTGGC‑GGTTTTTTTTG‑TTTAAAAGACTGGGGTTTCGCAAAA‑GCC‑GGGttt                                                                                                                 >  1:445438/1‑151 (MQ=255)                            caacaacaCCCATGGGGATGCTGCAGCAACCTCGCCCATTTTTCATGATCTTTTTTGTCGAGTTATGGGAGCGATTCGGCTACTACGGCGTGCAGGGCGTACTGGC‑GGTTTTTTTCG‑TTAAACAGCTTGGATTCTCGCCAGA‑GCA‑GGGttt                                                                                                                 >  1:585087/1‑151 (MQ=255)                              acaacaCCCATGGGGATGCTGCAGCAACCTCGCCCATTTTTCATGATCTTTTTTGTCGAGTTATGGGAGCGATTCGGCTACTACGGCGTGCAGGGCGTACTGGC‑GGTTTTTTTTG‑TTTAAAAGCTTTGGTTTTCTCCAGA‑GCA‑GGGttttt                                                                                                               >  2:291126/1‑150 (MQ=255)                              acaacaCCCATGGGGATGCTGCAGCAACCTCGCCCATTTTTCATGATCTTTTTTGTCGAGTTATGGGAGCGATTCGGCTACTACGGCGTGCAGGGCGTACTGGC‑GGTTTTTTTTC‑TTTAACAACCTGGGGTTTTGCAAGAGGCG‑GGCtttt                                                                                                                >  1:596625/1‑151 (MQ=255)                               caacaCCCATGGGGATGCTGCAGCAACCTCGCCCATTTTTCATGATCTTTTTTGTCGAGTTATGGGAGCGATTCGGCTACTACGGCGTGCAGGGCGTACTGGC‑GGTTTTCTTCG‑TTAAACAGCTTGGATTCTCGCAAGA‑GCA‑GGCTTTTGt                                                                                                              <  1:997583/151‑1 (MQ=255)                                 acacCCATGGGGATGCTGCAGCAACCTCGCCCATTTTTCATGATCTTTTTTGTCGAGTTATGGGAGCGATTCGGCTACTACGGCGTGCAGGGCGTACTGGC‑GGTTTTTTTTG‑TTTAAAAGCTTGGGTTTTCGCCAGA‑ACC‑GGGTTTTTTc                                                                                                             >  2:246670/1‑150 (MQ=255)                                 acacCCATGGGGATGCTGCAGCAACCTCGCCCATTTTTCATGATCTTTTTTGTCGAGTTATGGGAGCGATTCGGCTACTACGGCGTGCAGGGCGTACTGGC‑GGTTTTTTTCG‑TTTAAAAGCTTTGGTTCTCGCCAAA‑GCC‑GGCgtttttt                                                                                                             >  1:846363/1‑147 (MQ=255)                                 acacCCATGGGGATGCTGCAGCAACCTCGCCCATTTTTCATGATCTTTTTTGTCGAGTTATGGGAGCGATTCGGCTACTACGGCGTGCAGGGCGTACTGGC‑GGTTTTCTTCG‑TTAAACAGCTTGGATTCTCGCAAGA‑GCA‑GGCTTTTGTCa                                                                                                            <  2:267279/151‑1 (MQ=255)                                    cccATGGGGATGCTGCAGCAACCTCGCCCATTTTTCATGATCTTTTTTGTCGAGTTATGGGAGCGATTCGGCTACTACGGCGTGCAGGGCGTACTGGC‑GGTTTTTTTTG‑TTTAACCACTTGGGTTTTTGCCAAA‑GGA‑GGGtttttttcctt                                                                                                         >  2:395649/1‑144 (MQ=255)                                    cccATGGGGATGCTGCAGCAACCTCGCCCATTTTTCATGATCTTTTTTGTCGAGTTATGGGAGCGATTCGGCTACTACGGCGTGCAGGGCGTACTGGC‑GGTTTTTTTTG‑TTAAAAAGCTTGGGTTTTTGCAAGA‑GCC‑GGGCTTTTTCCCtt                                                                                                         >  2:144593/1‑151 (MQ=255)                                    cccATGGGGATGCTGCAGCAACCTCGCCCATTTTTCATGATCTTTTTTGTCGAGTTATGGGAGCGATTCGGCTACTACGGCGTGCAGGGCGTACTGGC‑GGTTTTCTTCT‑TTAAACAACCTGGCTTTTTCCCAGG‑GCAGGGCTTTTTTCa                                                                                                            >  1:322687/1‑149 (MQ=255)                                    cccATGGGGATGCTGCAGCAACCTCGCCCATTTTTCATGATCTTTTTTGTCGAGTTATGGGAGCGATTCGGCTACTACGGCGTGCAGGGCGTACTGGC‑GGGTTTCTTCT‑TTTAACATCTTGGGTTCTCCCAATA‑GGA‑GGGTTTTTTTACtt                                                                                                         >  1:539588/1‑151 (MQ=255)                                        tGGGGATGCTGCAGCAACCTCGCCCATTTTTCATGATCTTTTTTGTCGAGTTATGGGAGCGATTCGGCTACTACGGCGTGCAGGGCGTACTGGCGGGTTTTTTTCT‑TTAAA‑ATCTTGGGTTCTCCCCATA‑GAA‑GGGGTTTGTCTCtttttg                                                                                                     >  2:309700/1‑149 (MQ=255)                                        tGGGGATGCTGCAGCAACCTCGCCCATTTTTCATGATCTTTTTTGTCGAGTTATGGGAGCGATTCGGCTACTACGGCGTGCAGGGCGTACTGGC‑GGGTTTTCTCG‑TTTAACAGCCTTGAGTTTCGCAAGA‑GAA‑GGGTTTTTTTACtttttg                                                                                                     >  1:598466/1‑149 (MQ=255)                                        tGGGGATGCTGCAGCAACCTCGCCCATTTTTCATGATCTTTTTTGTCGAGTTATGGGAGCGATTCGGCTACTACGGCGTGCAGGGCGTACTGGC‑GGCTTTTTTCG‑TTTAACAGCTTCGGTTCCCGCAAAA‑CCC‑GGCCTTTGTTCCtttttg                                                                                                     >  2:243652/1‑149 (MQ=255)                                              tgctgcAGCAACCTCGCCCATTTTTCATGATCTTTTTTGTCGAGTTATGGGAGCGATTCGGCTACTACGGCGTGCAGGGCGTACTGGC‑GGTTTTTTTTT‑TTTAAAAACTTGGGTTCTTGCAAGA‑GCA‑GGGGTTTTTTAATTTTTGTGGttt                                                                                               >  1:304669/1‑151 (MQ=255)                                                ctgcAGCAACCTCGCCCATTTTTCATGATCTTTTTTGTCGAGTTATGGGAGCGATTCGGCTACTACGGCGTGCAGGGCGTACTGGC‑GGTTTTTTTTC‑TTTAACAGCTTTGGTTTTCGCAAGA‑ACA‑GGCTTTTTTTCCTTTTTGGGGttttt                                                                                             >  1:890639/1‑150 (MQ=255)                                                           tCGCCCATTTTTCATGATCTTTTTTGTCGAGTTATGGGAGCGATTCGGCTACTACGGCGTGCAGGGCGTACTGGC‑GGTTTTCTTCG‑TTAAACAGCTTGGATTCTCGCAAGA‑GCA‑GGCTTTTGTCACTTTTGGTGCTTTTGCTGCGCTGGTc                                                                                  <  2:422200/151‑1 (MQ=255)                                                           tCGCCCATTTTTCATGATCTTTTTTGTCGAGTTAGGGGAGCGATTCGGCTACTACGGCGTGCAGGGCGTACTGGC‑GGTTTTCTTCG‑TTAAACAGCTTGGATTCTCGCAAGA‑GCA‑GGCTTTTGTCACTTTTGGTGCTTTTGCTGCGCTGGTc                                                                                  <  1:672811/151‑1 (MQ=255)                                                               ccATTTTTCATGATCTTTTTTGTCGAGTTATGGGAGCGATTCGGCTACTACGGCGTGCAGGGCGTACTGGC‑GGTTTTCTTCG‑TTAAACAGCTTGGATTCTCGCAAGA‑GCA‑GGCTTTTGTCACTTTTGGTGCTTTTGCTGCGCTGGTCTA‑Tg                                                                             <  1:328003/151‑1 (MQ=255)                                                                cATTTTTCATGATCTTTTTTGTCGAGTTATGGGAGCGATTCGGCTACTACGGCGTGCAGGGCGTACTGGC‑GGTTTTCTTCG‑TTAAACAGCTTGGATTCTCGCAAGA‑GCA‑GGCTTTTGTCACTTTTGGTGCTTTTGCTGCGCTGGTCTA‑Tgg                                                                            <  1:418676/151‑1 (MQ=255)                                                                cATTTTTCATGATCTTTTTTGTCGAGTTATGGGAGCGATTCGGCTACTACGGCGTGCAGGGCGTACTGGC‑GGTTTTCTTCG‑TTAAACAGCTTGGATTCTCGCAAGA‑GCA‑GGCTTTTGTCACTTTTGGTGCTTTTGCTGCGCTGGTCTA‑Tgg                                                                            <  2:28937/151‑1 (MQ=255)                                                                 aTTTTTCATGATCTTTTTTGTCGAGTTATGGGAGCGATTCGGCTACTACGGCGTGCAGGGCGTACTGGC‑GGTTTTCTTCG‑TTAAACAGCTTGGATTCTCGCAAGA‑GCA‑GGCTTTTGTCACTTTTGGTGCTTTTGCTGCGCTGGTCTA‑TGGc                                                                           <  1:270082/151‑1 (MQ=255)                                                                 aTTTTTCATGATCTTTTTTGTCGAGTTATGGGAGCGATTCGGCTACTACGGCGTGCAGGGCGTACTGGC‑GGTTTTCTTCG‑TTAAACAGCTTGGATTCTCGCAAGA‑GCA‑GGCTTTTGTCACTTTTGGTGCTTTTGCTGCGCTGGTCTA‑TGGc                                                                           <  2:680740/151‑1 (MQ=255)                                                                 aTTTTTCATGATCTTTTTTGTCGAGTTATGGGAGCGATTCGGCTACTACGGCGTGCAGGGCGTACTGGC‑GGTTTTCTTCG‑TTAAACAGCTTGGATTCTCGCAAGA‑GCA‑GGCTTTTGTCACTTTTGGTGCTTTTGCTGCGCTGGTCTA‑TGGc                                                                           <  1:268692/151‑1 (MQ=255)                                                                 aTTTTTCATGATCTTTTTTGTCGAGTTATGGGAGCGATTCGGCTACTACGGCGTGCAGGGCGTACTGGC‑GGTTTTCTTCG‑CTAAACAGCTTGGATTCTCGCAAGA‑GCA‑GGCTTTTGTCACTTTTGGTGCTTTTGCTGCGCTGGTCTA‑TGGc                                                                           <  1:762174/151‑1 (MQ=255)                                                                    tttCATGATCTTTTTTGTCGAGTTATGGGAGCGATTCGGCTACTACGGCGTGCAGGGCGTACTGGC‑GGTTTTTTTTT‑TTTAAAAGCTTGGGTTTTTGCGAGA‑GCA‑GGGGTTTTTTTCTTTTTGGGGTTTTTCCGGGCTGGTCttttggcc                                                                          >  2:388751/1‑143 (MQ=255)                                                                         tGATCTTTTTTGTCGAGTTATGGGAGCGATTCGGCTACTACGGCGTGCAGGGCGTACTGGC‑GGTTTTCTTCG‑TTAAACAGCTTGGATTCTCGCAAGA‑GCA‑GGCTTTTGTCACTTTTGGTGCTTTTGCTGCGCTGGTCTA‑TGGCCTCATTTc                                                                   <  2:180002/151‑1 (MQ=255)                                                                           aTCTTTTTTGTCGAGTTATGGGAGCGATTCGGCTACTACGGCGTGCAGGGCGTACTGGC‑GGTTTTCTTCG‑TTAAACAGCTTGGATTCTCGCAAGA‑GCA‑GGCTTTTGTCACTTTTGGTGCTTTTGCTGCGCTGGTCTA‑TGGCCTCATTTcc                                                                  <  2:496257/150‑1 (MQ=255)                                                                                        aGTTATGGGAGCGATTCGGCTACTACGGCGTGCAGGGCGTACTGGC‑GGGTTTTTTTG‑TTAAAAAGCTTGGGTTTTCGCAAGA‑GCA‑GGCTTTTTTCACTTTTTGTGCTTTTTTTGCGCTGGGCTA‑TGGCCCCCTTTTCATTtgcggggattt                                                    >  2:33054/1‑145 (MQ=255)                                                                                        aGTTATGGGAGCGATTCGGCTACTACGGCGTGCAGGGCGTACTGGC‑GGGTTTTTTTG‑CTTAAAAACCTGGGTTCTCCCAAGA‑GCAGGGCTTTTTTTTTTTTTGG‑GGTTTTTTTTTGCTGGTCTATTGGCCCCTTTTTATTTggcgggtttt                                                     >  2:700513/1‑146 (MQ=255)                                                                                                        cGGCTACTACGGCGTGCAGGGCGTACTGGC‑GGTTTTCTTCG‑TTAAACAGCTTGGATTCTCGCAAGA‑GCA‑GGCTTTTGTCACTTTTGGTGCTTTTGCTGCGCTGGTCTA‑TGGCCTCATTTCCATTGGCGGCTATGTCGGCGACCACctgctg                                    <  1:556732/151‑1 (MQ=255)                                                                                                         ggCTACTACGGCGTGCAGGGCGTACTGGC‑GGTTTTTTTCG‑TTAAACAGCTTGGGTTTTCGCAAGA‑GCA‑GGGGTTTTTCACTTTTTGTGCTTTTTCTGCGCTGGTCTA‑TGGCCTCATTTTCCTTGGCGGCTGTGTTGGCGGCCACCCGCTgg                                   >  1:749113/1‑151 (MQ=255)                                                                                                             actacGGCGTGCAGGGCGTACTGGC‑GGTTTTCTTCG‑TTAAACAGCTTGGATTCTCGCAAGA‑GCA‑GGCTTTTGTCACTTTTGGTGCTTTTGCTGCGCTGGTCTA‑TGGCCTCATTTCCATTGGCGGCTATGTCGGCGACCACCTGCTGGGGAc                               <  2:596871/151‑1 (MQ=255)                                                                                                                  ggCGTGCAGGGCGTACTGGC‑GGTTTTCTTTT‑TTTAACAACTTTGAGTCTTGCAAGA‑GCA‑GGGTTTTTTCTCTTTTTGTGTTTTTGTGGCGTTGGGCTT‑TGGCCTCATTTCCTTTTGTGGGTATGTTGGGGGCCCCCCGGTggggggccaaa                          >  2:45583/1‑144 (MQ=255)                                                                                                                           ggCGTACTGGC‑GGTTTTCTTCG‑TTAAACAGCTTGGATTCTCGCAAGA‑GCA‑GGCTTTTGTCACTTTTGGTGCTTTTGCTGCGCTGGTCTA‑TGGCCTCATTTCCATTGGCGGCTATGTCGGCGACCACCTGCTGGGGACCAAACGCACCAt                   <  2:846653/149‑1 (MQ=255)                                                                                                                             cGTACTGAC‑GGTTTTCTTCG‑TTAAACAGCTTGGATTCTCGCAAGA‑GCA‑GGCTTTTGTCACTTTTGGTGCTTTTGCTGCGCTGGTCTA‑TGGCCTCATTTCCATTGGCGGCTATGTCGGCGACCACCCGCTGGGGGCCAAACCCACCCTTGtt               >  1:509219/1‑151 (MQ=255)                                                                                                                               tACTGGC‑GGTTTTCTTCG‑TTAAACAGCTTGGATTCTCGCAAGA‑GCA‑GGCTTTTGTCACTTTTGGTGCTTTTGCTGCGCTGGTCTA‑TGGCCTCATTTCCATTGGCGGCTATGTCGGCGACCACCTGCTGGGGACCAAACGCACCATTGTTc              <  2:514958/150‑1 (MQ=255)                                                                                                                                 cTGGC‑GGTTTTCTTCG‑TTAAACAGCTTGGATTCTCGCAAGA‑GCA‑GGCTTTTGTCACTTTTGGTGCTTTTGCTGCGCTGGTCTA‑Tgg                                                                            <  2:107986/86‑1 (MQ=255)                                                                                                                                 cTGGC‑GGTTTTCTTCG‑TTAAACAGCTTGGATTCTCGCAAGA‑GCA‑GGCTTTTGTCACTTTTGGTGCTTTTGCTGCGCTGGTCTA‑Tgg                                                                            >  1:107970/1‑86 (MQ=255)                                                                                                                                   ggC‑GGTTTTCTTCG‑TTAAACAGCTTGGATTCTCGCAAGA‑GCA‑GGCTTTTGTCACTTTTGGTGCTTTTGCTGCGCTGGTCTA‑TGGCCTCATTTCCATTGGCGGCTATGTCGGCGACCACCTGCTGGGGACCAAACGCACCATTGTTCTTGGa         >  1:180147/1‑151 (MQ=255)                                                                                                                                     c‑GGTTTTCTTCG‑TTAAACAGCTTGGATTCTCGCAAGA‑GCA‑GGCTTTTGTCACTTTTGGTGCTTTTGCTGCGCTGGTCTA‑TGGCCTCATTTCCATTGGCGGCTATGTCGGCGACCACCTGCTGGGGACCAAACGCACCATTGTTCTTGGAg        <  1:482608/150‑1 (MQ=255)                                                                                                                                     c‑GGTTTTCTTCG‑TTAAACAGCTTGGATTCTCGCAAGA‑GCA‑GGCTTTTGTCACTTTTGGTGCTTTTGCTGCGCTGGTCTA‑TGGCCTCATTTCCATTGGCGGCTATGTCGGCGACCACCTGCTGGGGACCAAACGCACCATTGTTCTTGGAg        <  1:292468/150‑1 (MQ=255)                                                                                                                                       ggTTTTCTTCG‑TTAAACAGCTTGGATTCTCGCAAGA‑GCA‑GGCTTTTGTCACTTTTGGTGCTTTTGCTGCGCTGGTCTA‑TGGCCTCATTTCCATTGGCGGCTTTGTCCGTGCCCCACTTCTGGGGGCCAAACCCACCGTTGGTCTTGGTGCa      >  2:555638/1‑151 (MQ=255)                                                                                                                                        gTTTTCTTCG‑TTAAACAGCTTGGATTCTCGCAAGA‑GCA‑GGCTTTTGTCACTTTTGGTGCTTTTGCTGCGCTGGTCTA‑TGGCCTCATTTCCATTGGCGGCTATGTTGGCGACCACCTGCTGGGGACCAAACGCACCATTTTTCTTGGGGCa      >  2:274009/1‑150 (MQ=255)                                                                                                                                           ttcttcG‑TTAAACAGCTTGGATTCTCGCAAGA‑GCA‑GGCTTTTGTCACTTTTGGTGCTTTTGCTGCGCTGGTCTA‑TGGCCTCATTTCCATTGGCGGCTATGTCGGCGACCACCTGCTGGGGACCAAACGCACCATTGTTCTTGGAGCACttn  <  2:46899/151‑2 (MQ=255)                                                                                                                                           ttcttcG‑TTAAACAGCTTGGATTCTCGCAAGA‑GCA‑GGCTTTTGTCACTTTTGGTGCTTTTGCTGCGCTGGTCTA‑TGGCCTCATTTCCATTGGCGGCTATGTCGGCGACCACCTGCTGGGGACCAAACGCACCATTGTTCTTGGAGCACTtg  <  2:324969/151‑1 (MQ=255)                                                                                                                                             |                                                                                                                                                         GTATTGGAAGGTTAATAAATATGAATACAACAACACCCATGGGGATGCTGCAGCAACCTCGCCCATTTTTCATGATCTTTTTTGTCGAGTTATGGGAGCGATTCGGCTACTACGGCGTGCAGGGCGTACTGGC‑GGTTTTCTTCG‑TTAAACAGCTTGGATTCTCGCAAGA‑GCA‑GGCTTTTGTCACTTTTGGTGCTTTTGCTGCGCTGGTCTA‑TGGCCTCATTTCCATTGGCGGCTATGTCGGCGACCACCTGCTGGGGACCAAACGCACCATTGTTCTTGGAGCACTTG  >  NC_000913/3640842‑3641129

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