Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,911,754	0	T	G	42.1%	125.9 / 13.1	76	T275P (ACC→CCC)	htpX	putative endopeptidase
Reads supporting (aligned to +/- strand):  ref base T (7/37);  new base G (32/0);  total (39/37)
Fisher's exact test for biased strand distribution p-value = 5.87e-15
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GTGTCAGTGGTTTACGTATCACTCAGCCACGATCCAGGGCATAAAAAAAGCGCGTCGATCAGGACGCGCTTTTTAGTATTTACTTCATATTACTTCAGGTATTCACCCGTACGCAGAGCTTCAATTCGTTTATCCAGCGGCGGGTGGGTCATGAACAACTCACTGAGCGATTTCGACTTACCGTTAATGCAGAGAGCCATCATGCTGGTTGCTTCTTGCGGTTCATAG‑CTGGTTTTCAGGCGCTGCAGCGCGGCAATCATTTTCTCGCGACCAACCAGTTTTGCCGAACCGGCATC  >  NC_000913/1911606‑1911901                                                                                                                                                     |                                                                                                                                                     gtgtCAGTGGTTTACGTATCACTCAGCCACGATCCAGGGCATAAAAAAAGCGCGTCGATCAGGACGCGCTTTTTAGTATTTACTTCATATTACTTCAGGTATTCACCCGTACGCAGAGCTTCAATTCGTTTATCCAGCGGCGGGTGGGTc                                                                                                                                                     <  2:897714/150‑1 (MQ=255)   gtCAGTGGTTTACGTATCACTCAGCCACGATCCAGGGCATAAAAAAAGCGCGTCGATCAGGACGCGCTTTTTAGTATTTACTTCATATTACTTCAGGTATTCACCCGTACGCAGAGCTTCAATTCGTTTATCCAGCGGCGGGTGGGTCAt                                                                                                                                                   <  2:859665/150‑1 (MQ=255)     cAGTGGTTTACGTATCACTCAGCCACGATCCAGGGCATAAAAAAAGCGCGTCGATCAGGACGCGCTTTTTAGTATTTACTTCATATTACTTCAGGTATTCACCCGTACGCAGAGCTTCAATTCGTTTATCCAGCGGCGGGTGGGTCATGaa                                                                                                                                                <  1:85572/151‑1 (MQ=255)     cAGTGGTTTACGTATCACTCAGCCACGATCCAGGGCATAAAAAAAGCGCGTCGATCAGGACGCGCTTTTTAGTATTTACTTCATATTACTTCAGGTATTCACCCGTACGCAGAGCTTCAATTCGTTTATCCAGCGGCGGGTGGGTCATGa                                                                                                                                                 <  1:785532/150‑1 (MQ=255)     cAGTGGTTTACGTATCACTCAGCCACGATCCAGGGCATAAAAAAAGCGCGTCGATCAGGACGCGCTTTTTAGTATTTACTTCATATTACTTCAGGTATTCACCCGTACGCAGAGCTTCAATTCGTTTATCCAGCGGCGGGTGGGTCATGa                                                                                                                                                 <  1:1221811/150‑1 (MQ=255)     cAGTGGTTTACGTATCACTCAGCCACGATCCAGGGCATAAAAAAAGCGCGTCGATCAGGACGCGCTTTTTAGTATTTACTTCATATTACTTCAGGTATTCACCCGTACGCAGAGCTTCAATTCGTTTATCCAGCGGCGGGTGGGTCATGa                                                                                                                                                 <  1:957840/150‑1 (MQ=255)           tttACGTATCACTCAGCCACGATCCAGGGCATAAAAAAAGCGCGTCGATCAGGACGCGCTTTTTAGTATTTACTTCATATTACTTCAGGTATTCACCCGTACGCAGAGCTTCAATTCGTTTATCCAGCGGCGGGTGGGGCATGaacaa                                                                                                                                             >  2:720296/1‑148 (MQ=255)                  aTCACTCAGCCACGATCCAGGGCATAAAAAAAGCGCGTCGATCAGGACGCGCTTTTTAGTATTTACTTCATATTACTTCAGGTATTCACCCGTACGCAGAGCTTCAATTCGTTTATCCAGCGGCGGGTGGGTCATGAACAACTCACTGAg                                                                                                                                    <  1:244472/150‑1 (MQ=255)                     aCTCAGCCACGATCCAGGGCATAAAAAAAGCGCGTCGATCAGGACGCGCTTTTTAGTATTTACTTCATATTACTTCAGGTATTCACCCGTACGCAGAGCTTCAATTCGTTTATCCAGCGGCGGGGGGGTCATGAACAACTCACTTAGCgg                                                                                                                                 >  1:1130192/1‑149 (MQ=255)                       tCAGCCACGATCCAGGGCATAAAAAAAGCGCGTCGATCAGGACGCGCTTTTTAGTATTTACTTCATATTACTTCAGGTATTCACCCGTACGCAGAGCTTCAATTCGTTTATCCAGCGGCGGGGGGGGCATGAACAACTCACTGAGCGAtt                                                                                                                               >  1:1146364/1‑150 (MQ=255)                        cAGCCACGATCCAGGGCATAAAAAAAGCGCGTCGATCAGGACGCGCTTTTTAGTATTTACTTCATATTACTTCAGGTATTCACCCGTACGCAGAGCTTCAATTCGTTTATCCAGCGGCGGGTGGGTCATGAACAACTCACTGAGCGATTTc                                                                                                                             <  1:1062341/151‑1 (MQ=255)                         aGCCACGATCCAGGGCATAAAAAAAGCGCGTCGATCAGGACGCGCTTTTTAGTATTTACTTCATATTACTTCAGGTATTCACCCGTACGCAGAGCTTCAATTCGTTTATCCAGCGGCGGGTGGGTCATGAACAACTCACTGAGCGATTTCg                                                                                                                            <  1:1210309/151‑1 (MQ=255)                            cACGATCCAGGGCATAAAAAAATCGCGTCGATCAGGACGCGCTTTTTAGTATTTACTTCATATTACTTCAGGTATTCACCCGTACGCAGAGCTTCAATTCGTTTATCCAGCGGCGGGTGGGGCATGAACAACCCACTGAGGGGttttcact                                                                                                                         >  2:835419/1‑146 (MQ=255)                              cGATCCAGGGCATAAAAAAAGCGCGTCGATCAGGACGCGCTTTTTAGTATTTACTTCATATTACTTCAGGTATTCACCCGTACGCAGAGCTTCAATTCGTTTATCCAGCGGCGGGGGGGGCATGAACAAACCACTGAGCGGTTTTGACtt                                                                                                                        >  1:465868/1‑150 (MQ=255)                               gATCCAGGGCATAAAAAAAGCGCGTCGATCAGGACGCGCTTTTTAGTATTTACTTCATATTACTTCAGGTATTCACCCGTACGCAGAGCTTCAATTCGTTTATCCAGCGGCGGGTGGGGCATGAACAACTCACTGAGCGATTTCGACtt                                                                                                                        >  2:468614/1‑149 (MQ=255)                               gATCCAGGGCATAAAAAAAGCGCGTCGATCAGGACGCGCTTTTTAGTATTTACTTCATATTACTTCAGGTATTCACCCGTACGCAGAGCTTCAATTCGTTTATCCAGCGGCGGGGGGGGCATGAACAACTCACTGAGCGATTTTGACtt                                                                                                                        >  2:431365/1‑149 (MQ=255)                                 tCCAGGGCATAAAAAAAGCGCGTCGATCAGGACGCGCTTTTTAGTATTTACTTCATATTACTTCAGGTATTCACCCGTACGCAGAGCTTCAATTCGTTTATCCAGCGGCGGGTGGGTCATGAACAACTCACTGAGCGATTTCTACTTAcc                                                                                                                     >  1:126483/1‑150 (MQ=255)                                    aGGGCATAAAAAAAGCGCGTCGATCAGGACGCGCTTTTTAGTATTTACTTCATATTACTTCAGGTATTCACCCGTACGCAGAGCTTCAATTCGTTTATCCAGCGGCGGGGGGGGCCAGAAAAAATCACTGAGGGGGTTTGACTTACCGt                                                                                                                   >  2:1297927/1‑149 (MQ=255)                                      ggCATAAAAAAAGCGCGTCGATCAGGACGCGCTTTTTAGTATTTACTTCATATTACTTCAGGTATTCACCCGTACGCAGAGCTTCAATTCGTTTATCCAGCGGCGGGTGGGTCATGAACAACTCACTGAGCGATTTCGACTTACCGTTAAt                                                                                                               <  1:1023693/151‑1 (MQ=255)                                                aaGCGCGTCGATCAGGACGCGCTTTTTAGTATTTACTTCATATTACTTCAGGTATTCACCCGTACGCAGAGCTTCAATTCGTTTATCCAGCGGCGGGTGGGTCATTAACAACTCACTGAGCGATTTCGACTTACCGTTAATGCAGAGAGcc                                                                                                     >  2:736524/1‑151 (MQ=255)                                                 aGCGCGTCGATCAGGACGCGCTTTTTAGTATTTACTTCATATTACTTCAGGTATTCACCCGTACGCAGAGCTTCAATTCGTTTATCCAGCGGCGGGTGGGTCATGAACAACTCACTGAGCGATTTCGACTTACCGTTAATGCAGAGAGCca                                                                                                    <  1:1534/151‑1 (MQ=255)                                                    gcgTCGATCAGGACGCGCTTTTTAGTATTTACTTCATATTACTTCAGGTATTCACCCGTACGCAGAGCTTCAATTCGTTTATCCAGCGGCGGGTGGGTCATGAACAACTCACTGAGCGATTTCGACTTACCGTTAATGCAGAGAGCcatca                                                                                                 <  2:414917/151‑1 (MQ=255)                                                     cgTCGATCAGGACGCGCTTTTTAGTATTTACTTCATATTACTTCAGGTATTCACCCGTACGCAGAGCTTCAATTCGTTTATCCAGCGGCGGGTGGGTCATGAACAACTCGCTGAGCGATTTCGACTTACCGTTAATGCAGAGAGCcatcat                                                                                                <  2:122919/151‑1 (MQ=255)                                                      gTCGATCAGGACGCGCTTTTTAGTATTTACTTCATATTACTTCAGGTATTCACCCGTACGCAGAGCTTCAATTCGTTTATCCAGCGGCGGGTGGGTCATGAACAACTCACTGAGCGATTTCGACTTACCGTTAATGCAGAGAGCCATCATg                                                                                               <  2:1130264/151‑1 (MQ=255)                                                      gTCGATCAGGACGCGCTTTTTAGTATTTACTTCATATTACTTCAGGTATTCACCCGTACGCAGAGCTTCAATTCGTTTATCCAGCGGCGGGTGGGGCATGAACAACTCACTGAGCGATTTCGACTTACCGTTAATGCAGAGAGCcatc                                                                                                  >  2:737269/1‑148 (MQ=255)                                                         gATCAGGACGCGCTTTTTAGTATTTACTTCATATTACTTCAGGTATTCACCCGTACGCAGAGCTTCAATTCGTTTATCCAGCGGCGGGTGGGTTATTAACAAAACACACACAGATTTCTAGTTACCCTTAAAGCAGAGAGGcaaca                                                                                                 >  1:983055/1‑146 (MQ=255)                                                          aTCAGGACGCGCTTTTTAGTATTTACTTCATATTACTTCAGGTATTCACCCGTACGCAGAGCTTCAATTCGTTTATCCAGCGGCGGGTGGGTCATGAACAACTCACTGAGCGATTTCGACTTACCGTTAATGCAGAGAGCCATCATGCTgg                                                                                           <  2:833289/151‑1 (MQ=255)                                                              ggACGCGCTTTTTAGTATTTACTTCATATTACTTCAGGTATTCACCCGTACGCAGAGCTTCAATTCGTTTATCCAGCGGCGGGTGGGTCATTAAAAACTCACTGAGCGATTTCGACTTACCGTTTATGCAGAGAGCCATCATGCTgg                                                                                           >  1:819938/1‑147 (MQ=255)                                                                 cgcgCTTTTTAGTATTTACTTCATATTACTTCAGGTATTCACCCGTACGCAGAGCTTCAATTCGTTTATCCAGCGGCGGGTGGGGCATGAACAACCCACTGAGCGATTTCGACTTACCGTTAATGCAGAGAGCCATCATGCTTGGTGctt                                                                                     >  1:1271412/1‑150 (MQ=255)                                                                     cTTTTTAGTATTTACTTCATATTACTTCAGGTATTCACCCGTACGCAGAGCTTCAATTCGTTTATCCAGCGGCGGGTGGGTCATGAACAACTCACTGAGCGATTTCGACTTACCGTTAATGCAGAGAGCCATCATGCTGGTTGcttctt                                                                                  >  2:26922/1‑149 (MQ=255)                                                                      tttttAGTATTTACTTCATATTACTTCAGGTATTCACCCGTACGCAGAGCTTCAATTCGTTTATCCAGCGGCGGGGGGGGCATGAACAACTCACTGAGCGATTTTGACTTACCGTTAATGCAGAGAGCCATCATGCTGGTTTCTTCTTTCg                                                                               >  2:923672/1‑151 (MQ=255)                                                                         ttAGTATTTACTTCATATTACTTCAGGTATTCACCCGTACGCAGAGCTTCAATTCGTTTATCCAGCGGCGGGTGGGGCCTGAACAACTCACTTAGCGATTTTGACTTACCGTTAATGCAGAGAGCCATCATGCTGGTTGCTTCTTTCGGtt                                                                            >  1:480590/1‑151 (MQ=255)                                                                          tAGTATTTACTTCATATTACTTCAGGTATTCACCCGTACGCAGAGCTTCAATTCGTTTATCCAGCGGCGGGGGGGGCATGAACAACTCACTGAGCGATTTCGACTTACCGTTAATGCAGAGAGCCATCATGCTGGTTGCTTCTTTCGGtt                                                                            >  2:728647/1‑150 (MQ=255)                                                                           aGTATTTACTTCATATTACTTCAGGTATTCACCCGTACGCAGAGCTTCAATTCGTTTATCCAGCGGCGGGTGGGTCATGAACAACTCACTGAGCGATTTCGACTTACCGTTAATGCAGAGAGCCATCATGCTGGTTGCTTCTTGCGGTTCa                                                                          <  1:1132179/151‑1 (MQ=255)                                                                           aGTATTTACTTCATATTACTTCAGGTATTCACCCGTACGCAGAGCTTCAATTCGTTTATCCAGCGGCGGGTGGGTCATGAACAACTCACTGAGCGATTTCGACTTACCGTTAATGCAGAGAGCCATCATGCTGGTTGCTTCTTGCGGTTCa                                                                          <  2:685950/151‑1 (MQ=255)                                                                            gTATTTACTTCATATTACTTCAGGTATTCACCCGTACGCAGAGCTTCAATTCGTTTATCCAGCGGCGGGTGGGGCATGAACAAATCACCTAGCGATTTTTACTTACCGGTTATGGAGAGAGCCATCATGCTGGGTGCTTCTTTCGGTTCAt                                                                         >  1:335336/1‑151 (MQ=255)                                                                             tATTTACTTCATATTACTTCAGGTATTCACCCGTACGCAGAGCTTCAATTCGTTTATCCAGCGGCGGGTGGGGCATGAACAAAACACTGAGCGATTTTTACCTTTCGTTAATGGAGAGAGCCATCATGCTGGTTGCTTTTTTCGGGTCAt                                                                         >  1:637460/1‑150 (MQ=255)                                                                               tttACTTCATATTACTTCAGGTATTCACCCGTACGCAGAGCTTCAATTCGTTTATCCAGCGGCGGGTGGGGCATGAACAAATCACCGAGCGATTTTGACTTACCGCTAATGCAGAGAGCCATCATGCTGGTTGCTTCTTTCGGGTCATa                                                                        >  1:1051151/1‑149 (MQ=255)                                                                                ttACTTCATATTACTTCAGGTATTCACCCGTACGCAGAGCTTCAATTCGTTTATCCAGCGGCGGGTGGGGCATGAACAACTCACTGAGCGATTTTGTCTTTCCGTTAATGCAGAGAGCCATCATGCTGGTTGCTTCTTTCGGGTTATAG‑c                                                                     >  1:947826/1‑150 (MQ=255)                                                                                     tCATATTACTTCAGGTATTCACCCGTACGCAGAGCTTCAATTCGTTTATCCAGCGGCGGGTGGGTCATGAACAACTCACTGAGCGATTTCGACTTACCGTTAATGCAGAGAGCCATCATGCTGGTTGCTTCTTGCGGTTCATAG‑CTGGtt                                                                <  1:720241/150‑1 (MQ=255)                                                                                     tCATATTACTTCAGGTATTCACCCGTACGCAGAGCTTCAATTCGTTTATCCAGCGGCGGGTGGGTCATGAACAACTCACTGAGCGATTTCGACTTACCGTTAATGCAGAGAGCCATCATGCTGGTTGCTTCTTGCGGTTCATAG‑CTGGt                                                                 <  2:814599/149‑1 (MQ=255)                                                                                     tCATATTACTTCAGGTATTCACCCGTACGCAGAGCTTCAATTCGTTTATCCAGCGGCGGGGGGGGCATGAACAACTCACTGAGCGATTTCGACTTACCGTTAATTCAGAGAGCCCTCATGCTGGTTTGTTCTTGCGGTTCATAG‑CTGGgtt                                                               >  1:1219822/1‑151 (MQ=255)                                                                                      cATATTACTTCAGGTATTCACCCGTACGCAGAGCTTCAATTCGTTTATCCAGCGGCGGGGGGGGCCTGAACAAATCACCGAGCGGTTTTGACTTACCGTTAATGCAGAGAGCCATCATGCTGGGTGCTTTTTTCGGTTCATAG‑CTGGgttt                                                              >  1:903291/1‑151 (MQ=255)                                                                                         atTACTTCAGGTATTCACCCGTACGCAGAGCTTCAATTCGTTTATCCAGCGGCGGGTGGGGCATGAACAACTCACTGAGCGATTTCGACTTACCGTTAATGCAGAGAGCCATCATGCTGGGTGCTTCTTTCGGTTCATAG‑CTGGtttt                                                              >  2:1249331/1‑148 (MQ=255)                                                                                           tACTTCAGGTATTCACCCGTACGCAGAGCTTCAATTCGTTTATCCAGCGGCGGGGGGGGGATTAACACATAACTCAGCGGTTTTTTCTTACCGGTGATGCAGGGAGGCATCATGCTGGGTGCTTTTTTCGGGTCATAGACGGGtttta                                                             >  1:382347/1‑147 (MQ=255)                                                                                            aCTTCAGGTATTCACCCGTACGCAGAGCTTCAATTCGTTTATCCAGCGGCGGGTGGGGCATGCACAAATCCCTGAGCGGTTTTGACTTTACGCTAATGCAGAGGGCCATCATGCTGGTTGGTTCTTGCGGGTCCTAA‑GTGGTTTTTAGgc                                                         >  1:917612/1‑150 (MQ=255)                                                                                            aCTTCAGGTATTCACCCGTACGCAGAGCTTCAATTCGTTTATCCAGCGGCGGGGGGGGCAAGAACAACTCACTGAGGGATTTCGACTTACCGTTAATGCAGAGAGCCATCATGCTGGTTGCTTTTTGCGGTTTATTG‑CTGGGTTTTAGgcg                                                        >  2:840831/1‑151 (MQ=255)                                                                                                  ggTATTCACCCGTACGCAGAGCTTCAATTCGTTTATCCAGCGGCGGGTGGGGCATGAAAAACTCACTGAGCGATTTTGACTTACCGTTAATGCAGAGAGCCATCATGCTGGTTGCTTCTTTCGGTTTATAG‑CTGGTTTTCAGGCGCTGc                                                    >  2:1196232/1‑149 (MQ=255)                                                                                                   gTATTCACCCGTACGCAGAGCTTCAATTCGTTTATCCAGCGGCGGGTGGGTCATGAACAAATCACTGAGCGATTTTGACTTACCGTTAATTCAGAGAGCCATCAAGCTGGTTGCTTCTTTCGGTTCATAG‑ATGGGTTTCAGGCGCTGCAgc                                                 >  1:1119272/1‑151 (MQ=255)                                                                                                     aTTCACCCGTACGCAGAGCTTCAATTCGTTTATCCAGCGGCGGGTGGGTCATGAACAACTCACTGAGCGATTTCGACTTACCGTTAATGCAGAGAGCCATCATGCTGGTTGCTTCTTGCGGTTCATAG‑CTGGTTTTCAGGCGCTGCAgcgc                                               <  2:652147/151‑1 (MQ=255)                                                                                                     aTTCACCCGTACGCAGAGCTTCAATTCGTTTATCCAGCGGCGGGTGGGTCATGAACAACTCACTGAGCGATTTCGACTTACCGTTAATGCAGAGAGCCATCATGCTGGTTGCTTCTTGCGGTTCATAG‑CTGGTTTTCAGGCGCTGCAgcg                                                <  1:201072/150‑1 (MQ=255)                                                                                                     aTTCACCCGTACGCAGAGCTTCAATTCGTTTATCCAGCGGCGGGTGGGTCATGAACAACTCACTGAGCGATTTCGACTTACCGTTAATGCAGAGAGCCATCATGCTGGTTGCTTCTTGCGGTTCATAG‑CTGGTTTTCAGGCGCTGCAgcg                                                <  2:330442/150‑1 (MQ=255)                                                                                                      ttCACCCGTACGCAGAGCTTCAATTCGTTTATCCAGCGGCGGCGGTGGCCGGAACTACACACTGTTAGAGTTTGAATTACCGTGCAGGGGGGGAGACAACAACATTGTTTCTTCTTTCGGTTCTTAG‑GTGGTTTTAAAgcgccac                                                    >  1:875978/1‑142 (MQ=255)                                                                                                        cACCCGTACGCAGAGCTTCAATTCGTTTATCCAGCGGCGGGTGGGTCATGAACAACTCACTGAGCGATTTCGACTTACCGTTAATGCAGAGAGCCATCATGCTGGTTGCTTCTTGCGGTTCATAG‑CTGGTTTTCAGGCGCTGCAGCGCgg                                             <  2:590575/150‑1 (MQ=255)                                                                                                        cACCCGTACGCAGAGCTTCAATTCGTTTATCCAGCGGCGGGGGGGGCATGAACAAAACACTGAGCGATTTCGACTTAACGGTAAAGCAGAGAGGCAACAAGGTGGGTGCTTCTTGCGGGTCATAA‑CTGGGTTTTAGGGGCTGCAGCGCgg                                             >  1:555226/1‑150 (MQ=255)                                                                                                          cccGTACGCAGAGCTTCAATTCGTTTATCCAGCGGCGGGTGGGTCATGAACAACTCACTGAGCGATTTCGACTTACCGTTAATGCAGAGAGCCATCATGCTGGTTGCTTCTTGCGGTTCATAG‑CTGGTTTTCAGGCGCTGCAGCGCGGCa                                           <  1:931123/150‑1 (MQ=255)                                                                                                           ccGTACGCAGAGCTTCAATTCGTTTATCCAGCGGCGGGTGGGTCATGAACAACTCACTGAGCGATTTCGACTTACCGTTAATGCAGAGAGCCATCATGCTGGTTGCTTCTTGCGGTTCATAG‑CTGGTTTTCAGGCGCTGCAGCGCGGCAAt                                         <  1:44654/151‑1 (MQ=255)                                                                                                           ccGTACGCAGAGCTTCAATTCGTTTATCCAGCGGCGGGTGGGGCATGAACAACTCCCTGAGCGAATTCGGCCTTCCCATAAAGCAAAGAGACATCCTGCTGGTTGGTTCTTGCGGTTCCTAG‑ATGGGTTTCACGCGCTGCAGCGCGGGaa                                          >  1:957449/1‑150 (MQ=255)                                                                                                           ccGTACGCAGAGCTTCAATTCGTTTATCCAGCGGCGGGGGGGGCATGAACAAATCACTGAGCGATTTTGACTTACCGTTAATGCAGAGAGGCCTCATGCTGGGTGCTTCTTTCGGGTCATAG‑CTGGGTTTTAGGCGCTGCAGCGCgg                                             >  1:368358/1‑147 (MQ=255)                                                                                                               aCGCAGAGCTTCAATTCGTTTATCCAGCGGCGGGTGGGTCATGAACAACTCACTGAGCGATTTCGACTTACCGTTAATGCAGAGAGCCATCATGCTGGTTGCTTCTTGCGGTTCATAG‑CTGGTTTTCAGGCGCTGCAGCGCGGCAATCAtt                                     <  2:65081/151‑1 (MQ=255)                                                                                                                 gCAGAGCTTCAATTCGTTTATCCAGCGGCGGGGGGGGCATGAAAAACTCACTGAGGGATTTCGACTTACCGTTAATGCAGAGAGCCATCATGCTGGTTGCTTCTTTCGGGTCATAG‑CTGGTTTTCAGGCGCTGCAGCGCGGCCATCAtttt                                   >  2:496137/1‑151 (MQ=255)                                                                                                                  cAGAGCTTCAATTCGTTTATCCAGCGGCGGGTGGGTCATGAACAACTCACTGAGCGATTTCGACTTACCGTTAATGCAGAGAGCCATCATGCTGGTTGCTTCTTGCGGTTCATAG‑CTGGTTTTCAGGCGCTGCAGCGCGGCAATCATTTtc                                  <  1:873567/151‑1 (MQ=255)                                                                                                                    gagCTTCAATTCGTTTATCCAGCGGCGGGTGGGTCATGAACAACTCACTGAGCGATTTCGACTTACCGTTAATGCAGAGAGCCATCATGCTGGTTGCTTCTTGCGGTTCATAG‑CTGGTTTTCAGGCGCTGCAGCGCGGCAATCATTTtctc                                <  1:222087/151‑1 (MQ=255)                                                                                                                    gagCTTCAATTCGTTTATCCAGCGGCGGGTGGGTCATGAACAACTCACTGAGCGATTTCGACTTACCGTTAATGCAGAGAGCCATCATGCTGGTTGCTTCTTGCGGTTCATAG‑CTGGTTTTCAGGCGCTGCAGCGCGGCAATCATTTtctc                                <  2:355955/151‑1 (MQ=255)                                                                                                                     agCTTCAATTCGTTTATCCAGCGGCGGGTGGGTCATGAACAACTCACTGAGCGATTTCGACTTACCGTTAATGCAGAGAGCCATCATGCTGGTTGCTTCTTGCGGTTCATAG‑CTGGTTTTCAGGCGCTGCAGCGCGGCAATCATTTTCTcg                               <  2:781368/151‑1 (MQ=255)                                                                                                                        ttCAATTCGTTTATCCAGCGGCGGGTGGGGCATGAACAACTCACTGAGCGATTTCGCCTTACCGTTAATGCAGAGAGCCCACATGCTGGTTGCTTTTTGCGGGTCATAG‑CTGGTTTTCAGGCGCTGCAGCGCGGCAAACAtttttttg                               >  2:424378/1‑144 (MQ=255)                                                                                                                        ttCAATTCGTTTATCCAGCGGCGGGTGGGGCATGAACAACTCACTGAGCGATTTCGACTTACCGTTAATGCAGAGAGCCATCATGCTGGTTGCTTCTTTCGGTTCATAG‑CTGGGTTTCAGGCGCTGCAGCGCGGCAATCATTTtct                                 >  2:310718/1‑146 (MQ=255)                                                                                                                          cAATTCGTTTATCCAGCGGCGGGTGGGTCATGAACAACTCACTGAGCGATTTCGACTTACCGTTAATGCAGAGAGCCATCATGCTGGTTGCTTCTTGCGGTTCATAG‑CTGGTTTTCAGGCGCTGCAGCGCGGCAATCATTTTCTCGCGacc                          <  1:704246/151‑1 (MQ=255)                                                                                                                            aTTCGTTTATCCAGCGGCGGGTGGGTCATGAACAACTCACTGAGCGATTTCGACTTACCGTTAATGCAGAGAGCCATCATGCTGGTTGCTTCTTGCGGTTCATAG‑CTGGTTTTCAGGCGCTGCAGCGCGGCAATCATTTTCTCGCGaccaa                        <  2:250376/151‑1 (MQ=255)                                                                                                                               cGTTTATCCAGCGGCGGGTGGGTCATGAACAACTCACTGAGCGATTTCGACTTACCGTTAATGCAGAGAGCCATCATGCTGGTTGCTTCTTGCGGTTCATAG‑CTGGTTTTCAGGCGCTGCAGCGCGGCAATCATTTTCTCGCGaccaacca                     <  1:1054717/151‑1 (MQ=255)                                                                                                                               cGTTTATCCAGCGGCGGGTGGGTCATGAACAACTCACTGAGCGATTTCGACTTACCGTTAATGCAGAGAGCCATCATGCTGGTTGCTTCTTGCGGTTCATAG‑CTGGTTTTCAGGCGCTGCAGCGCGGCAATCATTTTCTCGCGaccaacc                      <  2:1271090/150‑1 (MQ=255)                                                                                                                               cGTTTATCCAGCGGCGGGTGGGTCATGAACAACTCACTGAGCGATTTCGACTTACCGTTAATGCAGAGAGCCATCATGCTGGTTGCTTCTTGCGGTTCATAG‑CTGGTTTTCAGGCGCTGCAGCGCGGCAATCATTTTCTCGCGaccaacc                      <  1:908618/150‑1 (MQ=255)                                                                                                                                gTTTATCCAGCGGCGGGTGGGGCATGAACAACTCACTGAGCGATTTCGACTTACCGTTAATGCAGAGAGCCATCATGCTGGTTGCTTCTTTCGGGTCATAG‑CTGGTTTTCAGGCGCTGCAGCGCGGCAATCATTTTTTCGCGaccaacc                      >  2:19869/1‑149 (MQ=255)                                                                                                                                    aTCCAGCGGCGGGTGGGGCATGAACAACTCACTGAGCGATTTCGACTTACCGTTAATGCAGAGAGGCATCATGCTGGTTGCTTCTTTCGGGTCATAG‑CTGGTTTTCAGGCGCTGCAGCGCGGCAATCATTTTTTCGCGACCAACCCGtttt                >  1:145164/1‑151 (MQ=255)                                                                                                                                        aGCGGCGGGTGGGTCATGAACAACTCACTGAGCGATTTCGACTTACCGTTAATGCAGAGAGCCATCATGCTGGTTGCTTCTTGCGGTTCATAG‑CTGGTTTTCAGGCGCTGCAGCGCGGCAATCATTTTCTCGCGACCAACCAGTTTTGCCg            <  1:627879/151‑1 (MQ=255)                                                                                                                                             cggGTGGGTCATGAACAACTCACTGAGCGATTTCGACTTACCGTTAATGCAGAGAGCCATCATGCTGGTTGCTTCTTGCGGTTCATAG‑CTGGTTTTCAGGCGCTGCAGCGCGGCAATCATTTTCTCGCGACCAACCAGTTTTGCCGAACCg       <  2:651056/151‑1 (MQ=255)                                                                                                                                                gtgggtCATGAACAACTCACTGAGCGATTTCGACTTACCGTTAATGCAGAGAGCCATCATGCTGGTTGCTTCTTGCGGTTCATAG‑CTGGTTTTCAGGCGCTGCAGCGCGGCAATCATTTTCTCGCGACCAACCAGTTTTGCCGAACCGGCa    <  2:626702/151‑1 (MQ=255)                                                                                                                                                    gtCATGAACAACTCACTGAGCGATTTCGACTTACCGTTAATGCAGAGAGCCATCATGCTGGTTGCTTCTTGCGGTTCATAG‑CTGGTTTTCAGGCGCTGCAGCGCGGCAATCATTTTCTCGCGACCAACCAGTTTTGCCGAACCGGCATc  >  1:1001430/1‑149 (MQ=255)                                                                                                                                                     |                                                                                                                                                     GTGTCAGTGGTTTACGTATCACTCAGCCACGATCCAGGGCATAAAAAAAGCGCGTCGATCAGGACGCGCTTTTTAGTATTTACTTCATATTACTTCAGGTATTCACCCGTACGCAGAGCTTCAATTCGTTTATCCAGCGGCGGGTGGGTCATGAACAACTCACTGAGCGATTTCGACTTACCGTTAATGCAGAGAGCCATCATGCTGGTTGCTTCTTGCGGTTCATAG‑CTGGTTTTCAGGCGCTGCAGCGCGGCAATCATTTTCTCGCGACCAACCAGTTTTGCCGAACCGGCATC  >  NC_000913/1911606‑1911901

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: ‑