Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,571,036	0	T	G	49.2%	75.0 / 22.1	67	N239T (AAC→ACC)	eutE	aldehyde oxidoreductase, ethanolamine utilization protein
Reads supporting (aligned to +/- strand):  ref base T (11/23);  new base G (33/0);  total (44/23)
Fisher's exact test for biased strand distribution p-value = 7.13e-10
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CATCAGTTCATCGGCTACGCTATCAACAACAATCAGTACCTTTTCGTCGGCACAAATGATGTTGTTATCGAAAGAAGCGCCTTTGACGATGGACTGAGCGGCACGGGCGAGGTCGGCGGTTTCATCCACCACTACCGGCGGGTTGCCAGCGCCTGCGGCAATCAGACGTTTA‑TTGGTGT‑GTTTACGCGCCGCTTCTACTACCGCTTCGCCGCCGGTTACCACCAGCA‑GGCCGAT‑A‑CCCGGAAACTTGAACAAGCGTTGCGCGGTTTCGATATCCGGATTTGCCACAGTAACC  >  NC_000913/2570894‑2571185                                                                                                                                               |                                                                                                                                                           cATCAGTTCATCGGCTACGCTATCAACAACAATCAGTACCTTTTCGTCGGCACAAATGATGTTGTTATCGAAAGAAGCGCCTTTGACGATGGACTGAGCGGCACGGGCGAGGTCGGCGGTTTCATCCACCACTACCGGCGGGGTGCCcgcg                                                                                                                                                    >  2:1535279/1‑151 (MQ=255) cATCAGTTCATCGGCTACGCTATCAACAACAATCAGTACCTTTTCGTCGGCACAAATGATGTTGTTATCGAAAGAAGCGCCTTTGACGATGGACTGAGCGGCACGGGCGAGGTCGGCGGTTTCATCCACCACTACCGGCGGGGTGCCAgcg                                                                                                                                                    >  2:1676165/1‑151 (MQ=255)   tCAGTTCATCGGCTACGCTATCAACAACAATCAGTACCTTTTCGTCGGCACAAATGATGTTGTTATCGAAAGAAGCGCCTTTGACGATGGACTGAGCGGCACGGGCGAGGTCGGCGGTTTCATCCACCACTACCGGCGGGTTGCCAGCGcc                                                                                                                                                  <  2:232970/151‑1 (MQ=255)   ncAGTTCATCGGCTACGCTATCAACAACAATCAGTACCTTTTCGTCGGCACAAATGATGTTGTTATCGAAAGAAGCGCCTTTGACGATGGACTGAGCGGCACGGGCGAGGTCGGCGGTTTCATCCACCACTACCGGCGGGGTGCCCgcgc                                                                                                                                                   >  2:1264828/2‑150 (MQ=255)    cAGTTCATCGGCTACGCTATCAACAACAATCAGTACCTTTTCGTCGGCACAAATGATGTTGTTATCGAAAGAAGCGCCTTTGACGATGGACTGAGCGGCACGGGCGAGGTCGGCGGTTTCATCCACCACTACCGGCGGGGTTCCcgcggc                                                                                                                                                  >  1:63966/1‑148 (MQ=255)      gTTCATCGGCTACGCTATCAACAACAATCAGTACCTTTTCGTCGGCACAAATGATGTTGTTATCGAAAGAAGCGCCTTTGACGATGGACTGATCGGCACGGGCGAGGTCGGCGGTTTCATCCACCACTACCGGCGGGGTggccgccccccg                                                                                                                                               >  1:1473035/1‑140 (MQ=255)                  cGCTATCAACAACAATCAGTACCTTTTCGTCGGCACAAATGATGTTGTTATCGAAAGAAGCGCCTTTGACGATGGACTGAGCGGCACGGGCGAGGTCGGCGGTTTCCTCCACCACTACCGGCGGGGTGGCAGCGcccccgggcattagagg                                                                                                                                   >  1:1279935/1‑136 (MQ=255)                   gCTATCAACAACAATCAGTACCTTTTCGTCGGCACAAATGATGTTGTTATCGAAAGAAGCGCCTTTGACGATGGACTGAGCGGCACGGGCGAGGTCGGCGGTTTCATCCACCACTACCGGCGGGGTGCCCGCGCCCCCGGCAATCatccgt                                                                                                                                  >  2:944464/1‑146 (MQ=255)                      aTCAACAACAATCAGTACCTTTTCGTCGGCACAAATGATGTTGTTATCGAAAGAAGCGCCTTTGACGATGGACTGAGCGGCACGGGCGAGGTCGGCGGTTTCATCCACCACTACCGGCGGGTTGCCAGGGCCTGCGGCAATCAGACGTTTa                                                                                                                               >  2:1104261/1‑151 (MQ=255)                      aTCAACAACAATCAGTACCTTTTCGTCGGCACAAATGATGTTGTTATCGAAAGAAGCGCCTTTGACGATGGACTGAGCGGCACGGGCGAGGTCGGCGGTTTCATCCACCACTACCGGCGGGGTTGCAGGGCCCGCGGCCAACAGACGgttt                                                                                                                               >  1:1193151/1‑150 (MQ=255)                      aTCAACAACAATCAGTACCTTTTCGTCGGCACAAATGATGTTGTTATCGAAAGAAGCGCCTTTGACGATGGACTGAGCGGCACGGGCGAGGTCGGCGGTTTCATCCACCACTACCGGCGGGGTTCCAGCGCCCGCGGCCATCAGACGtttt                                                                                                                               >  1:1101029/1‑150 (MQ=255)                       tCAACAACAATCAGTACCTTTTCGTCGGCACAAATGATGTTGTTATCGAAAGAAGCGCCTTTGACGATGGACTGAGCGGCACGGGCGAGGTCGGCGGTTTCATCCACCACTACCGGCGGGTTGCCAGCGCCTGCGGCAATCAGACGTTTA‑t                                                                                                                             <  1:956608/151‑1 (MQ=255)                       tCAACAACAATCAGTACCTTTTCGTCGGCACAAATGATGTTGTTATCGAAAGAAGCGCCTTTGACGATGGACTGAGCGGCACGGGCGAGGTCGGCGGTTTCATCCACCACTACCGGCGGGTTGCCAGCGCCTGCGGCAATCAGACGTTTA‑t                                                                                                                             <  1:950020/151‑1 (MQ=255)                        caacaacaaTCAGTACCTTTTCGTCGGCACAAATGATGTTGTTATCGAAAGAAGCGCCTTTGACGATGGACTGAGCGGCACGGGCGAGGTCGGCGGTTTCATCCACCACTACCGGCGGGTTGCCCGCGCCCGCGGCCATCAGAAGtttt‑tt                                                                                                                            >  2:964338/1‑151 (MQ=255)                        caacaacaaTCAGTACCTTTTCGTCGGCACAAATGATGTTGTTATCGAAAGAAGCGCCTTTGACGATGGACTGAGCGGCACGGGCGAGGTCGGCGGTTTCATCCACCACTACCGGCGGGGTGCCAGCGCCCGCGGGGATCAGACGGTtt‑tt                                                                                                                            >  2:384744/1‑151 (MQ=255)                          acaacaaTCAGTACCTTTTCGTCGGCACAAATGATGTTGTTATCGAAAGAAGCGCCTTTGACGATGGACTGAGCGGCACGGGCGAGGTCGGCGGTTTCATCCACCACTACCGGCGGGTTGCCCGGGCCCGCGGGCATCAGACGCTtt‑tttg                                                                                                                          >  1:1234112/1‑149 (MQ=255)                          acaacaaTCAGTACCTTTTCGTCGGCACAAATGATGTTGTTATCGAAAGAAGCGCCTTTGACGATGGACTGAGCGGCACGGGCGAGGTCGGCGGTTTCATCCACCACTACCGGCGGGTTGCCAGCGCCTGCGGCAATCAGACGTTTA‑TTgg                                                                                                                          <  1:1560543/151‑1 (MQ=255)                          acaacaaTCAGTACCTTTTCGTCGGCACAAATGATGTTGTTATCGAAAGAAGCGCCTTTGACGATGGACTGAGCGGCACGGGCGAGGTCGGCGGTTTCATCCACCACTACCGGCGGGTTGCCAGCGCCTGCGGCAATCAGACGTTTA‑TTgg                                                                                                                          <  1:872589/151‑1 (MQ=255)                          acaacaaTCAGTACCTTTTCGTCGGCACAAATGATGTTGTTATCGAAAGAAGCGCCTTTGACGATGGACTGAGCGGCACGGGCGAGGTCGGCGGTTTCATCCACCACTACCGGCGGGGTGCCAGCGCCTGCGGCAATCAGACGTTTA‑TTgg                                                                                                                          >  1:17941/1‑151 (MQ=255)                          acaacaaTCAGTACCTTTTCGTCGGCACAAATGATGTTGTTATCGAAAGAAGCGCCTTTGACGATGGACTGAGCGGCACGGGCGAGGTCGGCGGTTTCATCCACCACTACCGGCGGGGTGCCAGCGCCCGCGGCAATCAGACGGTtt‑tttg                                                                                                                          >  2:490466/1‑149 (MQ=255)                               aaTCAGTACCTTTTCGTCGGCACAAATGATGTTGTTATCGAAAGAAGCGCCTTTGACGATGGACTGAGCGGCACGGGCGAGGTCGGCGGTTTCATCCACCACTACCGGCGGGGTGCCAGCGCCCGCGGCCAACACACGCTTT‑TTGgttt‑gt                                                                                                                    >  1:845474/1‑151 (MQ=255)                               aaTCAGTACCTTTTCGTCGGCACAAATGATGTTGTTATCGAAAGAAGCGCCTTTGACGATGGACTGAGCGGCACGGGCGAGGTCGGCGGTTTCATCCACCACTACCGGCGGGGTGCCACCGCCCGCGGGAATCAGGCGGTTT‑TTGgggt‑gt                                                                                                                    >  2:793818/1‑151 (MQ=255)                               aaTCAGTACCTTTTCGTCGGCACAAATGATGTTGTTATCGAAAGAAGCGCCTTTGACGATGGACTGAGCGGCACGGGCGAGGTCGGCGGCTTCATCCACCACTAACGGCGGGTTTCCCGGGCCTGCGCCGAACAGACGGTTA‑TTtgtgt‑gg                                                                                                                    >  1:1464261/1‑150 (MQ=255)                                aTCAGTACCTTTTCGTCGGCACAAATGATGTTGTTATCGAAAGAAGCGCCTTTGACGATGGACTGAGCGGCACGGGCGAGGTCGGCGGTTTCATCCACCACTACCGGCGGGTGGCCAGCGCCGGCGGCAATCAGACGGTTA‑TTGgtggggt                                                                                                                    >  1:1061101/1‑147 (MQ=255)                                 tCAGTACCTTTTCGTCGGCACAAATGATGTTGTTATCGAAAGAAGCGCCTTTGACGATGGACTGAGCGGCACGGGCGAGGTCGGCGGTTTCATCCACCACTACCGGCGGGTTGCCAGCGCCTGCGGCAATCAGACGTTTA‑TTGGTGT‑Gttt                                                                                                                  <  1:1435123/151‑1 (MQ=255)                                    gTACCTTTTCGTCGGCACAAATGATGTTGTTATCGAAAGAAGCGCCTTTGACGATGGACTGAGCGGCACGGGCGAGGTCGGCGGTTTCATCCACCACTACCGGCGGGGTGCCAGCGCCCGCGGGCATCAGACGGTTT‑TTGGGGT‑GTTTTcg                                                                                                               >  2:1335793/1‑151 (MQ=255)                                    gTACCTTTTCGTCGGCACAAATGATGTTGTTATCGAAAGAAGCGCCTTTGACGATGGACTGAGCGGCACGGGCGAGGTCGGCGGTTTCATCCACCACTACCGGCGGGGTGCCAGCGCCCGCGGGAATTAGACGGTTT‑TTTGGGT‑GGTTTcg                                                                                                               >  2:1639075/1‑151 (MQ=255)                                     tACCTTTTCGTCGGCACAAATGATGTTGTTATCGAAAGAAGCGCCTTTGACGATGGACTGAGCGGCACGGGCGAGGTCGGCGGTTTCATCCACCACTACCGGCGGGTTGCCAGCGCCTGCGGCAATCAGACGTTTA‑TTGGTGT‑GTTTAcgc                                                                                                              <  1:1731716/151‑1 (MQ=255)                                     tACCTTTTCGTCGGCACAAATGATGTTGTTATCGAAAGAAGCGCCTTTGACGATGGACTGAGCGGCACGGGCGAGGTCGGCGGTTTCATCCACCACTACCGGCGGGTTGCCAGCGCCTGCGGCAATCAGACGTTTA‑TTGGTGT‑GTTTAcg                                                                                                               <  1:1335195/150‑1 (MQ=255)                                                  gCACAAATGATTTTGTTATCGAAAGAAGCGCCTTTGACGATGGACTGAGCGGCACGGGCGAGGTCGGCGGTTTCATCCACCACTACCGGCGGGGTGGCCGCGCCCGCGGGAAACCGACGGTTT‑TTGGGGT‑GTTTTCGCGCCGCCTcttcta                                                                                                 >  1:1728986/1‑151 (MQ=255)                                                  gCACAAATGATGTTGTTATCGAAAGAAGCGCCTTTGACGATGGACTGAGCGGCACGGGCGAGGTCGGCGGTTTCATCCACCACTACCGGCGGGTTGCCAGCGCCTGCGGCAATCAGACGTTTA‑TTGGTGT‑GTTTACGCGCCGCTTctact                                                                                                  <  2:1103457/150‑1 (MQ=255)                                                  gCACAAATGATGTTGTTATCGAAAGAAGCGCCTTTGACGATGGACTGAGCGGCACGGGCGAGGTCGGCGGTTTCATCCACCACTACCGGCGGGTTGCCAGCGCCTGCGGCAATCAGACGTTTA‑TTGGTGT‑GTTTACGCGCCGCTTctact                                                                                                  <  2:86829/150‑1 (MQ=255)                                                  gCACAAATGATGTTGTTATCGAAAGAAGCGCCTTTGACGATGGACTGAGCGGCACGGGCGAGGTCGGCGGTTTCATCCACCACTACCGGCGGGGTGCCACCGCCCGCGGCCATCAGACGTTTT‑TTGGTGT‑GTTTTAGCGCCGCCTctacta                                                                                                 >  1:1349172/1‑151 (MQ=255)                                                  gCACAAATGATGTTGTTATCGAAAGAAGCGCCTTTGACGATGGACTGAGCGGCACGGGCGAGGGCGGCGGTTTCAACCACCACCACCGGGGGGTTGCCAAGCCCAGCCGCGATCTAACCTCTATTTGGTG‑‑GTTTGCCCGCCGCCTctactt                                                                                                 >  2:373432/1‑150 (MQ=255)                                                   cacaAATGATGTTGTTATCGAAAGAAGCGCCTTTGACGATGGACTGAGCGGCACGGGCGAGGTCGGCGGTTTCATCCACCACTACCGGCGGGTTGCCAGCGCCTGCGGCAATCAGACGTTTA‑TTGGTGT‑GTTTACGCGCCGCTTCtactac                                                                                                <  1:903074/151‑1 (MQ=255)                                                     caAATGATGTTGTTATCGAAAGAAGCGCCTTTGACGATGGACTGAGCGGCACGGGCGAGGTCGGCGGTTTCATCCACCACTACCGGCGGGTTGCCAGCGCCTGCGGCAATCAGACGTTTA‑TTGGTGT‑GTTTACGCGCCGCTTCTACTACCg                                                                                              <  2:1280524/151‑1 (MQ=255)                                                      aaaTGATGTTGTTATCGAAAGAAGCGCCTTTGACGATGGACTGAGCGGCACGGGCGAGGTCGGCGGTTTCATCCACCACTACCGGCGGGTTGGCAGCCCCCGCGGGCATCAAACCGTTT‑TTGGTGT‑GTTTTCGCGCCGCCTCTaccaccc                                                                                              >  2:1282418/1‑149 (MQ=255)                                                        atgatgTTGTTATCGAAAGAAGCGCCTTTGACGATGGACTGAGCGGCACGGGCGAGGTCGGCGGTTTCATCCACCACTACCGGCGGGGTGCCCGCGCCCGCGGCCATCAGACGGTTT‑TTTGTGT‑GTTTTCGCGCCCCTTCTACTAcccc                                                                                             >  2:37037/1‑147 (MQ=255)                                                             gttgttATCGAAAGAAGCGCCTTTGACGATGGACTGAGCGGCACGGGCGAGGTCGGCGGTTTCATCCACCACTACCGGCGGGGTGCCAGCGCCCGCGGGAAACAGACGGTTT‑TTGGGGG‑GTTTTCGCGCCCCCTTTACTCCCGCCTcccc                                                                                       >  2:1419726/1‑150 (MQ=255)                                                                   aTCGAAAGAAGCGCCTTTGACGATGGACTGAGCGGCACGGGCGAGGTCGGCGGTTTCATCCACCACTACCGGCGGGGTGCCCGCGCCTTCGGCCATCAGACGTTTT‑TTTGTGG‑GTTTTCGCGCCGCCTCTACTACCGCCTCGCCGCCgggt                                                                                >  1:1491556/1‑149 (MQ=255)                                                                          gaagCGCCTTTGACGATGGACTGAGCGGCACGGGCGAGGTCGGCGGTTTCATCCACCACTACCGGCGGGTTGCCCGCGCCTGCGGCAATCAGACGTTTT‑TTGGTGT‑GTTTTCGCGCCGCCTCTACTACCCCCTCCCCGCCgggttccccc                                                                          >  2:877302/1‑142 (MQ=255)                                                                           aagCGCCTTTGACGATGGACTGAGCGGCACGGGCGAGGTCGGCGGTTTCATCCACCACTACCGGCGGGTTGCCAGCGCCTGCGGCAATCAGACGTTTA‑TTGGTGT‑GTTTACGCGCCGCTTCTACTACCGCTTCGCCGCCGGTTACCACcag                                                                        <  2:557224/151‑1 (MQ=255)                                                                               gcCTTTGACGATGGACTGAGCGGCACGGGCGAGGTCGGCGGCTTCATCCACCACTACCGGCGGGTTGCCAGCGCCTGCGGCAATCAGACGTTTA‑TTGGGGT‑GTTTACGCGCCGCTTCTACTACCGCTTCGCCGCCGGTTACCACCAGCA‑gg                                                                   >  2:33640/1‑151 (MQ=255)                                                                                     gACGATGGACTGAGCGGCACGGGCGAGGTCGGCGGTTTCATCCACCACTACCGGCGGGGTGCCAGCGCCTGCGGGAATCCGACGTTTA‑TTGGTGT‑GTTTTCGCGCCGCTTCTACTACCGCCTCGCCGCCGGGTTCCCCCAGCA‑GGCCGCT‑a                                                            >  2:127842/1‑151 (MQ=255)                                                                                      aCGATGGACTGAGCGGCACGGGCGAGGTCGGCGGTTTCATCCACCACTACCGGCGGGTTGCCAGCGCCTGCGGCAATCAGACGTTTA‑TTGGTGT‑GTTTACGCGCCGCTTCTACTACCGCTTCGCCGCCGGTTACCACCAGCA‑GGCCGAT‑A‑c                                                          <  2:1492246/151‑1 (MQ=255)                                                                                      aCGATGGACTGAGCGGCACGGGCGAGGTCGGCGGTTTCATCCACCACTACCGGCGGGGTGCCCGCGCCCGCGGCAATCAGACGTTTT‑TTGGTGT‑GTTTTCGCGCCGCCTCTACTACCCCTTCCCCGCCGGTTACCACCAGCA‑GGCCGAT‑a                                                            >  1:159947/1‑150 (MQ=255)                                                                                      aCGATGGACTGAGCGGCACGGGCGAGGTCGGCGGTTTCATCCACCACTACCGGCGGGGTGCCCGCGCCCGCGGCAATCAGACGTTTA‑TTTGGGT‑GTTTTCGCGCCGCTTCTCCTACCGCCTCGCCGCCGGGTTCCCCCCGCA‑AGcccgta                                                             >  1:1474252/1‑146 (MQ=255)                                                                                           ggACTGAGCGGCACGGGCGAGGTCGGCGGTTTCATCCCCCACTACCGGCGGGTTGCCCGCGCCCGCGGCGCACAGAGGTTTA‑TTGGGGC‑GGTTTCGCGCCCCCTCTCCCCCCCCCTCGCCCCCCGTTTCCCCCCGCC‑GGCCTAT‑T‑CCCCGa                                                     >  2:1341980/1‑151 (MQ=255)                                                                                           ggACTGAGCGGCACGGGCGAGGTCGGCGGTTTCATCCCCCACCACCGGCGGGGTGCCCGCGGCCGCGGCCATCAGGCGCTTT‑TTTGTGG‑GGTTGCGCGCCGCTTCTACTACCCCCTCGCCGCCGGCTACCACCCGCC‑GGCCGCT‑A‑ccccgg                                                     >  1:106752/1‑148 (MQ=255)                                                                                           ggACTGAGCGGCACGGGCGAGGTCGGCGGTTTCATCCACCACTACCGGCGGGGTTCCCGCGGCCGCGGCCATCAAGCGGTTT‑TTGGTGG‑GTTTATGCGCCGCTTCTTCTATCGCCTCGCCGCCGGGTACCACCAgcc‑ggccgggaacccgg                                                       >  1:497531/1‑142 (MQ=255)                                                                                           ggACTGAGCGGCACGGGCGAGGTCGGCGGTTTCATCCACCACTACCGGCGGGGTGCCAGCGCCCTCGGGAATCAGAAGGTTT‑TTTGGGT‑GTTTTCGCGCCGCTTCTACTACCGCCTCCCCGCCGGGTAACACCAACACGGCCG‑T‑A‑CCCggg                                                     >  1:1414734/1‑150 (MQ=255)                                                                                           ggACTGAGCGGCACGGGCGAGGTCGGCGGTTTCATCCACCACTACCGGCGGGGTGCCAGCGCCCGCGGCAATCAGACGGTTT‑TTGGGGG‑GTTTACGCGCCGCTTCTACTACCCCTTCGCCCCCGGGTTACACCAGCA‑GGCCGAG‑A‑CCCGGa                                                     >  2:833936/1‑151 (MQ=255)                                                                                            gACTGAGCGGCACGGGCGAGGTCGGCGGTTTCATCCACCACTACCGGCGGGGTGCCCGCGCCCGCGGGAATTCGACGCTTT‑TTGGTGTGGTTTAC‑CGGCGCCTCTACTACCGCCTTCCCCCCGGGTTCCACCAACA‑GGcccctaacccgga                                                      >  1:247583/1‑140 (MQ=255)                                                                                              cTGAGCGGCACGGGCGAGGTCGGCGGTTTCATCCACCACTACCGGCGGGTTGCCAGCGCCTGCGGCAATCAGACGTTTA‑TTGGTGT‑GTTTACGCGCCGCTTCTACTACCGCTTCGCCGCCGGTTACCACCAGCA‑GGCCGAT‑A‑CCCGGAAAc                                                  <  1:1540295/151‑1 (MQ=255)                                                                                                gAGCGGCACGGGCGAGGTCGGCGGTTTCATCCACCACTACCGGCGGGGTGCCCGCGCCTGCGGCCATCCGACGGTTT‑TTTGTGG‑GTTTTCTCGCCCCTTCTTCCACCGCTTCGCCGCCCGTTTCCACCAGCA‑GGCCGGT‑A‑CCCCGGAACtt                                                >  1:720936/1‑151 (MQ=255)                                                                                                     gCACGGGCGAGGTCGGCGGTTTCATCCACCACTACCGGCGGGGTGCCCGCGCCTGCGGGCATCAGACGTTTT‑TTTGTGT‑GTTTACGCCCCCCCTCTACCACCCCCTCCCCCCCCGGTTCCCCCAGCA‑GGCCCAT‑A‑CCCCGGAACTTTAACa                                           >  2:826952/1‑151 (MQ=255)                                                                                                       aCGGGCGAGGTCGGCGGTTTCATCCACCACTACCGGCGGGTTGCCAGCGCCTGCGGCAATCAGACGTTTA‑TTGGTGT‑GTTTACGCGCCGCTTCTACTACCGCTTCGCCGCCGGTTACCACCAGCA‑GGCCGAT‑A‑CCCGGAAACTTGAACAAg                                         <  2:1124508/151‑1 (MQ=255)                                                                                                          ggCGAGGTCGGCGGTTTCATCCACCACTACCGGCGGGGTGCCCGCGCCCGCGGCAATTCGACGTTTA‑TTTGGGG‑GGTTTCGCGCCGCTTCTACTTCCCCTTCGCCCCCCGGTTCCCCCAGCA‑GGGCCATAA‑CCCGGAAACTTAACCAACCg                                       >  1:946445/1‑151 (MQ=255)                                                                                                                gTCGGCGGTTTCATCCACCACTACCGGCGGGTTGCCAGCGCCTGCGGCAATCAGACGTTTA‑TTGGTGT‑GTTTACGCGCCGCTTCTACTACCGCTTCGCCGCCGGTTACCACCAGCA‑GGCCGAT‑A‑CCCGGAAACTTGAACAAGCGTTgcgc                                 <  2:1414246/150‑1 (MQ=255)                                                                                                                  cggcggTTTCATCCACCACTACCGGCGGGGTGCCCGCGCCCCCGGCCATCACACGGTTT‑TTTGTGT‑GTTTTTGCGCCCCCTCTTCTTCCCCCTCCCCCCCCGGTTCCCCCCCCA‑GGCCCGT‑ACCCCGGAACCTTG‑ACAAACCCTGCGCggg                              >  1:38845/1‑150 (MQ=255)                                                                                                                     cggTTTCATCCACCACTACCGGCGGGTTGCCAGCGCCTGCGGCAATCAGACGTTTA‑TTGGTGT‑GTTTACGCGCCGCTTCTACTACCGCTTCGCCGCCGGTTACCACCAGCA‑GGCCGAT‑A‑CCCGGAAACTTGAACAAGCGTTGCGCGGttt                            <  2:1532379/150‑1 (MQ=255)                                                                                                                       gTTTCATCCACCACTACCGGCGGGTTGCCAGCGCCTGCGGCAATCAGACGTTTA‑TTGGTGG‑GGTTACGCGCCGCTTCTACTACCGCTTCGCCGCCGGTTACCACCAGCA‑GGCCGAT‑A‑CCCGGGAACTTGAACAAGCGTTGCGCGGGTTCGa                         >  2:170103/1‑151 (MQ=255)                                                                                                                        tttCATCCACCACTACCGGCGGGTTGCCAGCGCCTGCGGCAATCAGACGTTTA‑TTGGTGT‑GTTTACGCGCCGCTTCTACTACCGCTTCGCCGCCGGTTACCACCAGCA‑GGCCGAT‑A‑CCCGGAAACTTGAACAAGCGTTGCGCGGTTTCGa                         <  2:1387172/150‑1 (MQ=255)                                                                                                                         ttCATCCACCACTACCGGCGGGTTGCCAGCGCCTGCGGCAATCAGACGTTTA‑TTGGTGT‑GTTTACGCGCCGCTTCTACTACCGCTTCGCCGCCGGTTACCACCAGCA‑GGCCGAT‑A‑CCCGGAAACTTGAACAAGCGTTGCGCGGTTTCGata                       <  2:288792/151‑1 (MQ=255)                                                                                                                         ttCATCCACCACTACCGGCGGGTTGCCAGCGCCTGCGGCAATCAGACGTTTA‑TTGGTGT‑GTTTACGCGCCGCTTCTACTACCGCTTCGCCGCCGGTTACCACCAGCA‑GGCCGAT‑A‑CCCGGAAACTTGAACAAGCGTTGCGCGGTTTCGata                       <  1:180947/151‑1 (MQ=255)                                                                                                                                cccccTACCGGCGGGTTGCCAGCGCCTGCGGCAATCAGACGTTTA‑TTGGTGT‑GTTTACGCGCCGCTTCTACTACCGCTTCGCCGCCGGTTACCACCAGCA‑GGCCGAT‑A‑CCCGGAAACTTGAACAAGCGTTGCGCGGTTTCGATATCCGGAt                <  2:1313641/147‑1 (MQ=255)                                                                                                                                     tACCGGCGGGTTGCCAGCGCCTGCGGCAATCAGACGTTTA‑TTGGTGT‑GTTTACGCGCCGCTTCTACTACCGCTTCGCCGCCGGTTACCACCAGCA‑GGCCGAT‑A‑CCCGGAAACTTGAACAAGCGTTGCGCGGTTTCGATATCCGGATTTGcc           <  2:788415/151‑1 (MQ=255)                                                                                                                                               ttGCCAGCGCCTGCGGCAATCAGACGTTTA‑TTGGTGT‑GTTTACGCGCCGCTTCTACTACCGCTTCGCCGCCGGTTACCACCAGCA‑GGCCGAT‑A‑CCCGGAAGCTTGAACAAGCGGTGCGCGGTTTCGATATCCGGATTTGCCACAGTAAcc  <  1:1184998/150‑1 (MQ=255)                                                                                                                                               ttGCCAGCGCCTGCGGCAATCAGACGTTTA‑TTGGTGT‑GTTTACGCGCCGCTTCTACTACCGCTTCGCCGCCGGTTACCACCAGCA‑GGCCGAT‑A‑CCCGGAAACTTGAACAAGCGTTGCGCGGTTTCGATATCCGGATTTGCCACAGTAAcc  <  2:63988/150‑1 (MQ=255)                                                                                                                                               |                                                                                                                                                           CATCAGTTCATCGGCTACGCTATCAACAACAATCAGTACCTTTTCGTCGGCACAAATGATGTTGTTATCGAAAGAAGCGCCTTTGACGATGGACTGAGCGGCACGGGCGAGGTCGGCGGTTTCATCCACCACTACCGGCGGGTTGCCAGCGCCTGCGGCAATCAGACGTTTA‑TTGGTGT‑GTTTACGCGCCGCTTCTACTACCGCTTCGCCGCCGGTTACCACCAGCA‑GGCCGAT‑A‑CCCGGAAACTTGAACAAGCGTTGCGCGGTTTCGATATCCGGATTTGCCACAGTAACC  >  NC_000913/2570894‑2571185

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