Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,097,581	0	A	G	46.2%	144.2 / 27.7	91	noncoding (72/88 nt)	serX	tRNA‑Ser
Reads supporting (aligned to +/- strand):  ref base A (8/41);  new base G (42/0);  total (50/41)
Fisher's exact test for biased strand distribution p-value = 3.50e-18
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CACTGCGTTCACGGTCTGTCCAACTGGCTGTCGCCAGTTGTCGAACCCCGGTCGGGGCTTCTCATCCCCCCGGTGTGTGCAATATACG‑AAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGC‑GGTGAGGGGGGGATTCGAACCCCCGATACGT‑TGCCGTATACACA‑C‑TTTCCAGGCGTGCTCCTTCAGCCACTCGGACACCTCA‑CCAAATTGTCGTTCCTGTC‑TTGC‑TGGAACGGGCGCTAATTTA‑GGGAAATCATGACCTGAGGTCAACAAACTTTTTGAAAAA  >  NC_000913/1097436‑1097731                                                                                                                                                    |                                                                                                                                                              cACTGCGTTCACGGTCTGTCCAACTGGCTGTCGCCAGTTGTCGAACCCCGGTCGGGGCTTCTCATCCCCCCGGTGTGTGCAATATACG‑AAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGC‑GGTGAGGGGGGGGTTCg                                                                                                                                                           >  2:358461/1‑150 (MQ=14)  aCTGCGTTCACGGTCTGTCCAACTGGCTGTCGCCAGTTGTCAAACCCCGGTCGGGGCTTCTCATCCCCCCGGTGTGTGCAATATACG‑AAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGC‑GGTGAGGGGGGGATTCGaa                                                                                                                                                         <  2:576369/151‑1 (MQ=17)      cGTTCACGGTCTGTCCAACTGGCTGTCGCCAGTTGTCGAACCCCGGTCGGGGCTTCTCACCCCCCCGGTGTGTGCAATATACG‑AAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGC‑GGTGAGGGGGGGATTCGAAcccc                                                                                                                                                     <  1:1248983/151‑1 (MQ=25)       gTTCCCGGTCTGTCCAACTGGCTGTCGCCAGTTGTCGAACCCCGGTCGGGGCTTCTCATCCCCCCGGTGTGTGCAATATACG‑AAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGC‑GGTGAGGGGGGGATTCGAAccccc                                                                                                                                                    <  1:568993/151‑1 (MQ=25)       gTTCACGGTCTGTCCAACTGGCTGTCGCTAGTTGTCGAACCCCGGTCGGGGCTTCTCATCCCCCCGGTGTGTGCAATATACG‑‑AAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGC‑GGTGGGGGGGGGGTTTGAAcccccc                                                                                                                                                   >  1:369072/1‑150 (MQ=11)           aCGGTCTGTCCAACTGGCTGTCGCCAGTTGTCGAACCCCGGTCGGGGCTTCTCATCCCCCCGGTGTGTGCAATATACG‑AAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGC‑GGTGAGGGGGGGATTCGAACCCCCGATa                                                                                                                                                <  1:1322732/151‑1 (MQ=39)            cGGTCTGTCCAACTGGCTGTCGCCAGTTGTCGAACCCCGGTCGGGGCTTCTCATCCCCCCGGTGTGTGCAATATACG‑AAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGGGGGGGGGGGGGTTGGAAcccccatac                                                                                                                                                >  2:494013/1‑147 (MQ=11)                  gTCCAACTGGCTGTCGCCAGTTGTCGAACCCCGGTCGGGGCTTCTCATCCCCCCGGTGTGTGCAATATACG‑AAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGC‑GGTGAGGGGGGGATTCGAACCCCCGATACGT‑Tg                                                                                                                                          <  2:925110/149‑1 (MQ=255)                        cTGGCTGTCGCCAGTTGTCGAACCCCGGTCGGGGCTTCTCATCCCCCCGGTGTGTGCAATATACG‑AAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGC‑GGTGGGGGGGGGGGTTGGACCCCCCATAAGT‑TTCCGtata                                                                                                                                   >  2:1230285/1‑150 (MQ=11)                        cTGGCTGTCGCCAGTTGTCGAACCCCGGTCGGGGCTTCTCATCCCCCCGGTGTGTGCAATATACG‑AAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGC‑GGTGAGGGGGGGGTTCGAAACCCCCCTACGT‑TGCCGtata                                                                                                                                   >  2:1365693/1‑150 (MQ=255)                                cGCCAGTTGTCGAACCCCGGTCGGGGCTTCTCATCCCCCCGGTGTGTGCAATATACG‑AAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGC‑GGTGAGGGGGGGGTTTGAAACCCCCATAACG‑TGCCGGATacacaacttt                                                                                                                          >  2:1385638/1‑146 (MQ=11)                                cGCCAGTTGTCGAACCCCGGTCGGGGCTTCTCATCCCCCCGGTGTGTGCAATATACG‑AAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGC‑GGTGAGGGGGGGATTCGAACCCCCGATACGT‑TGCCGTATACACA‑C‑ttt                                                                                                                         <  2:1498813/150‑1 (MQ=255)                                 gCCAGTTGTCGAACCCCGGTCGGGGCTTCTCATCCCCCCGGTGTGTGCAATATACG‑AAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGC‑GGTGAGGGGGGGATTCGAACCCCCGATACGT‑TGCCGTATACACA‑C‑TTTcc                                                                                                                       <  1:1324435/151‑1 (MQ=255)                                  ccAGTTGTCGAACCCCGGTCGGGGCTTCTCATCCCCCCGGTGTGTGCAATATACG‑AAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGC‑GGTGAGGGGGGGGTTCGAACCCCCCATACGT‑TTCCGTATACACA‑C‑ttttcc                                                                                                                      >  2:3431/1‑148 (MQ=255)                                   cAGTTGTCGAACCCCGGTCGGGGCTTCTCATCCCCCCGGTGTGTGCAATATACG‑AAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGC‑GGTGAGGGGGGGGTTCGAACACCCCACACGT‑TGCCGTATACAAA‑C‑ttttccc                                                                                                                     >  1:1563321/1‑147 (MQ=14)                                    aGTTGTCGAACCCCGGTCGTGGCTTCTCATCCCCCCGGTGTGTGCAATATACG‑AAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGC‑GGTGAGGGGGGGATTCGAACCCCCGATACGT‑TGCCGTATACACA‑C‑TTTCCAgg                                                                                                                    <  1:674703/151‑1 (MQ=255)                                      ttGTCGAACCCCGGTCGGGGCTTCTCATCCCCCCGGTGTGTGCAATATACG‑AAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGC‑GGTGAGGGGGGGGTTCGAACCCCCCCAACGT‑TTCCCTATACACC‑C‑CTTcccagcg                                                                                                                  >  1:791334/1‑146 (MQ=14)                                       tGTCGAACCCCGGTCGGGGCTTCTCATCCCCCCGGTGTGTGCAATATACG‑AAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGC‑GGTGAGGGGGGGGTTCGAAACCCCCATACGT‑TGCCGTATACACA‑C‑TTTTCAGGCgg                                                                                                                 >  2:1253306/1‑150 (MQ=255)                                        gTCGAACCCCGGTCGGGGCTTCTCATCCCCCCGGTGTGTGCAATATACG‑AAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGC‑GGTGAGGGGGGGATTCGAACCCCCGATACGT‑TGCCGTATACACA‑C‑TTTCCAGGCGTg                                                                                                                <  1:1607427/151‑1 (MQ=255)                                        gTCGAACCCCGGTCGGGGCTTCTCATCCCCCCGGTGTGTGCAATATACG‑AAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGC‑GGTGAGGGGGGGATTCGAACCCCCGATACGT‑TGCCGTATACACA‑C‑TTTCCAGGCGTg                                                                                                                <  1:275724/151‑1 (MQ=255)                                              ccccGGTCGGGGCTTCTCATCCCCCCGGTGTGTGCAATATACG‑AAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGC‑GGTGAGGGGGGGGGTTGAAACCCCCCTACGT‑TGGCGTATACACA‑C‑TTTccccgggggctccct                                                                                                          >  2:1153957/1‑138 (MQ=14)                                              ccccGGTCGGGGCTTCTCATCCCCCCGGTGTGTGCAATATACG‑AAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGC‑GGTGAGGGGGGGGGTCGAACCCCCCCAACCT‑TTCCGTATACACA‑C‑TTTTCCGGCGTGCTccct                                                                                                          >  2:1259900/1‑149 (MQ=14)                                               cccGGTCGGGGCTTCTCATCCCCCCGGTGTGTGCAATATACG‑AAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGC‑GGTGAGGGGGGGGGTTGAAACCCCCATAAGG‑TGCCGGATAAACA‑A‑TTTTCAGGCGGGGTCCttt                                                                                                         >  2:1127957/1‑150 (MQ=14)                                                 cGGTCGGGGCTTCTCATCCCCCCGGTGTGTGCAATATACG‑AAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGC‑GGTGAGGGGGGGATTCGAACCCCCGATACGT‑TGCCGTATACACA‑C‑TTTCCAGGCGTGCTCCTTCAg                                                                                                       <  2:822430/151‑1 (MQ=255)                                                 cGGTCGGGGCTTCTCATCCCCCCGGTGTGTGCAATATACG‑AAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGC‑GGGGGGGGGGGGGGTTGAACCCCCCCCAAGG‑TTCCCTATAAACACC‑TTTCCAGGctggtcccttct                                                                                                        >  1:813072/1‑140 (MQ=11)                                                     cGGGGCTTCTCATCCCCCCGGTGTGTGCAATATACG‑AAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGC‑GGTGAGGGGGGGGGTTGAACCCCCCCTCCGT‑TGCCGTATACACA‑C‑TTTCCATGCGGGCTCCTTCCGCCAc                                                                                                   >  1:790382/1‑151 (MQ=17)                                                       gggCTTCTCATCCCCCCGGTGTGTGCAATATACG‑AAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGC‑GGTGAGGGGGGGGTTTGAAACCCCCATACGT‑TGCCCGTTACACA‑C‑TTTCCAAGCGTGCTCCTTCCGCCACTc                                                                                                 >  2:349646/1‑151 (MQ=17)                                                        ggCTTCTCATCCCCCCGGTGTGTGCAATATACG‑AAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGC‑GGTGAGGGGGGGATTCGAACCCCCGATACGT‑TGCCGTATACACA‑C‑TTTCCAGGCGTGCTCCTTCAGCCACTCg                                                                                                <  1:1277637/151‑1 (MQ=255)                                                        ggCTTCTCATCCCCCCGGTGTGTGCAATATACG‑AAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGC‑GGTGAGGGGGGGATTCGAACCCCCGATACGT‑TGCCGTATACACA‑C‑TTTCCAGGCGTGCTCCTTCAGCCACTCg                                                                                                <  2:368179/151‑1 (MQ=255)                                                         gCTTCTCATCCCCCCGGTGTGTGCAATATACG‑AAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGC‑GGTGAGGGGGGGATTCGAACCCCCGATACGT‑TGCCGTATACACA‑C‑TTTCCAGGCGTGCTCCTTCAGCCACTCgg                                                                                               <  1:613466/151‑1 (MQ=255)                                                          cTTCTCATCCCCCCGGTGTGTGCAATATACG‑AAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGC‑GGTGAGGGGGGGGGTTGAACCCCCCCTAAGT‑TTCCGGCTAAACA‑C‑CTTCCAGGCGGGCTCCCTCAGCCAccccgc                                                                                              >  2:1519449/1‑146 (MQ=14)                                                            tctcATCCCCCCGGTGTGTGCAATATACG‑AAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGC‑GGTGAGGGGGGGATTCGAACCCCCGATACGT‑TGCCGTATACACA‑C‑TTTCCAGGCGTGCTCCTTCAGCCACTCGGac                                                                                             <  2:1086986/150‑1 (MQ=255)                                                                       cGGTGTGTGCAATATACG‑AAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGC‑GGTGAGGGGGGGATTCGAACCCCCGATACGT‑TGCCGTATACACA‑C‑TTTCCAGGCGTGCTCCTTCAGCCACTCGGACACCTCA‑CCAAAt                                                                                <  1:1575635/151‑1 (MQ=255)                                                                        ggTGTGTGCAATATACG‑AAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGC‑GGTGAGGGGGGGATTCGAACCCCCGATACGT‑TGCCGTATACACA‑C‑TTTCCAGGCGTGCTCCTTCAGCCACTCGGACACCTCA‑CCAAAtt                                                                               <  2:659093/151‑1 (MQ=255)                                                                          tgtgtgCAATATACGAAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGC‑GGTGAGGGGGGGGGTTGAACCCCCCCTACGT‑TTCCGTATACACA‑C‑TTTTCCGGGGTGGTCCTTCAGCCACCCCGACACCCCA‑ccccactt                                                                              >  1:210772/1‑145 (MQ=14)                                                                          tgtgtgCAATATACG‑AAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGC‑GGTGAGGGGGGGGTTTGAACCCCCCATAAGT‑TTCCGTATACACA‑C‑TTTTCAGGCGTGCTCCTTCTGCCACCCCGGCACCCCC‑CCAAAtttt                                                                             >  2:1180311/1‑149 (MQ=14)                                                                          tgtgtgCAATATACG‑AAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGC‑GGTGAGGGGGGGGTTCGGACCCCCGATAAGT‑TTTCGTATACACA‑C‑TTTTCAGGCGGGGTCCTTTCGCCACCCGGACACCCCACCCAAAttt                                                                              >  2:506969/1‑150 (MQ=14)                                                                                cAATATACG‑AAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGC‑GGTGAGGGGGGGGGTTGAACCCCCCATACGT‑TGTCGTATACACA‑C‑CTTCCCAGCGTGCTCCTTCCGCCACCCGGACACCTCA‑CCAAATTTTTGTTcc                                                                       >  2:585631/1‑151 (MQ=14)                                                                                cAATATACG‑AAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGC‑GGTGAGGGGGGGATTCGAACCCCCGATACGT‑TGCCGTATACTCA‑C‑TTTCCAGGCGTGCTCCTTCAGCCACTCGGACACCTCA‑CCAAATTGTCGTTcc                                                                       <  1:376608/151‑1 (MQ=255)                                                                                  atatACG‑AAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGC‑GGTGAGGGGGGGGGTTGAACCCCCCCTACGT‑TTCCGTATACACA‑C‑CTTTCAGGCGTGCTCCCTCCGCCCCCCGGACACCCCC‑CCCAAATTTCGTTCCtt                                                                     >  1:437780/1‑150 (MQ=14)                                                                                  atatACG‑AAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGC‑GGTGAGGGGGGGATTCGAACCCCCGATACGT‑TGCCGTATACACA‑C‑TTTCCAGGCGTGCTCCTTCAGCCACTCGGACACCTCA‑CCAAATTGTCGTTCCTg                                                                     <  2:1557422/151‑1 (MQ=255)                                                                                    ataCG‑AAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGC‑GGTGAGGGGGGGATTCGAACCCCCGATACGT‑TGCCGTATACACA‑C‑TTTCCAGGCGTGCTCCTTCAGCCACTCGGACACCTCA‑CCAAATTGTCGTTCCTGt                                                                    <  1:1615688/150‑1 (MQ=255)                                                                                    ataCG‑AAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGC‑GGTGAGGGGGGGATTCGAACCCCCGATACGT‑TGCCGTATACACA‑C‑TTTCCAGGCGTGCTCCTTCAGCCACTCGGACACCTCA‑CCAAATTGTCGTTCCTGt                                                                    <  1:397054/150‑1 (MQ=255)                                                                                      aCG‑AAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGC‑GGTGAGGGGGGGATTCGAACCCCCGATACGT‑TGCCGTATACACA‑C‑TTTCCAGGCGTGCTCCTTCAGCCACTCGGACACCTCA‑CCAAATTGTCGTTCCTGTC‑tt                                                                <  1:12460/151‑1 (MQ=255)                                                                                      aCG‑AAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGC‑GGTGAGGGGGGGATTCGAACCCCCGATACGT‑TGCCGTATACACA‑C‑TTTCCAGGCGTGCTCCTTCAGCCACTCGGACACCTCA‑CCAAATTGTCGTTCCTGTC‑tt                                                                <  2:134760/151‑1 (MQ=255)                                                                                      aCG‑AAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGC‑GGTGAGGGGGGGATTCGAACCCCCGATACGT‑TGCCGTATACACA‑C‑TTTCCAGGCGTGCTCCTTCAGCCACTCGGACACCTCA‑CCAAATTGTCGTTCCTGTC‑t                                                                 <  2:1580542/150‑1 (MQ=255)                                                                                      aCG‑AAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGC‑GGTGAGGGGGGGATTCGAACCCCCGATACGT‑TGCCGTATACACA‑C‑TTTCCAGGCGTGCTCCTTCAGCCACTCGGACACCTCA‑CCAAATTGTCGTTCCTGTC‑t                                                                 <  2:638293/150‑1 (MQ=255)                                                                                        g‑AAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGC‑GGTGAGGGGGGGGGTCGAACCCCCCCCACGT‑TTCCGTATACAAA‑C‑TTTCCAGGCGGGCTCCTTCAGCCACCCCGACACCTCC‑CCAAAATTTCGTTCCTGTC‑tttt                                                              >  1:694627/1‑149 (MQ=14)                                                                                          aaaaaaaaGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGC‑GGTGAGGGGGGGATTCGAACCCCCGATACGT‑TGCCGTATACACA‑C‑TTTCCAGGCGTGCTCCTTCAGCCACTCGGACACCTCA‑CCAAATTGTCGTTCCTGTC‑TTGC‑t                                                            <  2:367726/151‑1 (MQ=255)                                                                                           aaaaaaaGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGC‑GGAGAAGGGGGGGTTTGAACCCCCCCTACGT‑TTTCGTATACACA‑C‑TTTTCAGGGGGGGTCCTTCCGCCCCTCGGGCACCCCC‑CCCAAATTTCGTTCCTTTC‑tttttt                                                            >  1:277018/1‑147 (MQ=9)                                                                                            aaaaaaGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGC‑GGTGAGGGGGGGGGTTGAACCCCCCATAAGT‑TGTCGTATAAACA‑C‑TTTTCCGGCGTGCTCCTTCAGCCCCTCCGGCACCCCA‑CCCAAATTTCGTTCCTGTC‑TTTC‑Tgg                                                          >  2:1396241/1‑151 (MQ=14)                                                                                             aaaaaGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGC‑GGCGAGGGGGGGGGTCGAAACCCCCCTACGT‑TTTCCTATAAAAA‑CTTTTCCAGG‑GTGGTCCCTTAGCCACCCGGACACCCCA‑CCCAACTTTCGTTCCTGTC‑TTTC‑ttgg                                                         >  1:1288859/1‑148 (MQ=9)                                                                                              aaaaGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGC‑GGTGAGGGGGGGGTTCGAACCCCCCATAAGT‑TGCCGTATTCACA‑C‑TTTCCCGGCGGGCCCCCTCACCCCCCCCGGCACCCCC‑CCCAAATTTCGTTCCTGTC‑TTTC‑ttgg                                                         >  2:655770/1‑147 (MQ=9)                                                                                              aaaaGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGC‑GGTGAGGGGGGGGGTTGAAACCCCCATAAGT‑TTTCGTATAAACA‑C‑TTTTCCGGCGGGCTCCTTCAGCCACCCCGACACCCCC‑CCAAAATTTCGTTTCTGTC‑TTGT‑Tggg                                                         >  1:501770/1‑149 (MQ=9)                                                                                                 aGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGC‑GGTGAGGGGGGGATTCGAACCCCCGATACGT‑TGCCGTATACACA‑C‑TTTCCAGGCGTGCTCCTTCAGCCACTCGGACACCTCA‑CCAAATTGTCGTTCCTGTC‑TTGC‑TGGAACgg                                                     <  2:225067/151‑1 (MQ=255)                                                                                                  gCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGC‑GGTGAGGGGGGGGTTTTAACCCCCCCCAAGT‑TTTCGTTTAAACA‑C‑TTTTCAGGCGGGGTCCCTCAGCCCCCCCGGCACCCCACCCAAATTTTGTTTCCT‑TC‑TTTC‑TGGAAAggg                                                    >  1:1259351/1‑151 (MQ=2)                                                                                                    ccGTACTTTCGTACGAGCTCTTCTTTAAATATGGC‑GGTGAGGGGGGGATTCGAACCCCCGATACGT‑TGCCGTATACACA‑C‑TTTCCAGGCGTGCTCCTTCAGCCACTCGGACACCTCA‑CCAAATTGTCGTTCCTGTC‑TTGC‑TGGAACGGgcg                                                  <  2:1383311/151‑1 (MQ=255)                                                                                                     cGTACTTTCGTACGAGCTCTTCTTTAAATATGGC‑GGTGAGGGGGGGGGTCGAACCCCCCCTACGG‑TTCCGGATACACA‑C‑TTTTCAGGCGTGCTCCTTCCGCCCCTCCGGCACCCCA‑CCCAAATTTCGTTTCTGTC‑TTGC‑TGGGAAGGgcgc                                                 >  2:1544704/1‑151 (MQ=14)                                                                                                        aCTTTCGTACGAGCTCTTCTTTAAATATGGC‑GGTGGGGGGGGGGGTTGGACCCCCGATACGC‑TGCCGTGTACACA‑A‑ATTCCCCGGGTGCGCCCTTCGCCCATCGGAAACCTCC‑CCAAATTTTTGTTCCTGTC‑CTTC‑TTGGACGGgggcca                                               >  2:671974/1‑148 (MQ=2)                                                                                                         cTTTCGTACGAGCTCTTCTTTAAATATGGC‑GGTGGGGGGGGGGTTTGAACCCCCCCTACGT‑TTCCGTATACACA‑C‑TTTTCCGGCGGGCTCCTTCAGCCACTCCGGCACCCCC‑CCCAAATTTCGTTCCTGTC‑TTTC‑TGGAACGGgcgcccat                                             >  1:1387936/1‑147 (MQ=14)                                                                                                         cTTTCGTACGAGCTCTTCTTTAAATATGGC‑GGTGAGGGGGGGGTTCGAACCCCCCCTAAGG‑TGCCGTATACACA‑C‑TTTTCAGGCGGGCTCCTTTAGCCCCTCGGACAACCCA‑CCAAATTTTTGTTCCTGTC‑TTTC‑TTGGAAGGGGGCCAAt                                             >  2:1214384/1‑151 (MQ=11)                                                                                                             cGTACGAGCTCTTCTTTAAATATGGC‑GGTGAGGGGGGGGTTCGAACCCCCGCTACGT‑TGCCGTATACACA‑C‑TTTCCAGGCGTGGTCCTTTAGCCCCTCCGACACCCCC‑CCAAATTTTCGTTCCTGTC‑TTGC‑TGGGACGGGCGCTAAttttg                                         >  2:399375/1‑149 (MQ=17)                                                                                                             cGTACGAGCTCTTCTTTAAATATGGC‑GGTGAGGGGGGGGTTCGAACCCCCCATACGT‑TGCCGTATACACA‑C‑TTTTCCCGCGTGCTCCTTCCGCCACTCGGACACCTCA‑CCCAAATTTTGTTTCTGTC‑TTGC‑TGGAACGGGCGCTAAttttt                                         >  2:184585/1‑149 (MQ=17)                                                                                                             cGTACGAGCTCTTCTTTAAATATGGC‑GGTGAGGGGGGGATTCGAACCCCCGATACGT‑TGCCGTATACACA‑C‑TTTCCAGGCGTGCTCCTTCAGCCACTCGGACACCTCA‑CCAAATTGTCGTTCCTGTC‑TTGC‑TGGAACGGGCGCTAATTTa                                          <  1:480763/150‑1 (MQ=255)                                                                                                              gTACGAGCTCTTCTTTAAATATGGC‑GGTGAGGGGGGGATTCGAACCCCCGATACGT‑TGCCGTATACACA‑C‑TTTCCAGGCGTGCTCCTTCAGCCACTCGGACACCTCA‑CCAAATTGTCGTTCCTGTC‑TTGC‑TGGAACGGGCGCTAATTTA‑gg                                       <  1:17903/151‑1 (MQ=255)                                                                                                               tACGAGCTCTTCTTTAAATATGGC‑GGTGGGGGGGGGGTTCGAACCCCCCATAAGT‑TGTCGTATAAAAA‑C‑CTTTCCGGGGGGGGCCTTCCGCCCCCCCGACACCCCC‑CCAAAATTTCGTTCCTGTC‑TTTCTTGGAACGGG‑GCTAAATTT‑ggg                                      >  1:767162/1‑151 (MQ=9)                                                                                                                aCGAGCTCTTCTTTAAATATGGC‑GGTGAGGGGGGGATTCGAACCCCCGATACGT‑TGCCGTATACACA‑C‑TTTCCAGGCGTGCTCCTTCAGCCACTCGGACACCTCA‑CCAAATTGTCGTTCCTGTC‑TTGC‑TGGAACGGGCGCTAATTTA‑GGGa                                     <  2:879878/151‑1 (MQ=255)                                                                                                                     ctctTCTTTAAATATGGC‑GGTGAGGGGGGGGTTCGAAACCCCCACACCTCTGCC‑TGTATACA‑C‑TTTTCAGGGGGGGTCCTTCTGCCCCCCCGGCACCCCC‑CCAAAATTTCATTCCCTTCTTTGT‑TGGAAGGGGCCCGATTTTA‑GGGAATTc                                 >  1:292343/1‑151 (MQ=2)                                                                                                                     ctctTCTTTAAATATGGC‑GGTGAGGGGGGGGGTCGAAACCCCCCTAAGT‑TTCCGGATAAACA‑C‑CTTCCCAGCGGGCTCCTTCCGCCACTCGGGCACCTCC‑CCCAATTTTCGGTTCTTTC‑TTGC‑TGGGAAGGGCGCTAATTTT‑GGGaaaac                                 >  2:1605680/1‑148 (MQ=9)                                                                                                                      tcttctTTAAATATGGC‑GGTGAGGGGGGGGTTCGAACCCCCCATACGT‑TGCCGTATACACA‑A‑CTTTCAGGCGGGCTCCTTCCGCCCCCCCGACACCTCA‑CCCAAATTTCGTTCCTGTC‑TTGC‑TGGAACGGGCGCCCATTTT‑GGGAAATCAt                               >  2:509876/1‑151 (MQ=17)                                                                                                                        ttcttTAAATATGGC‑GGTGAGGGGGGGGGTTGAAACCCCCATACGT‑TTCCGGATACACA‑C‑TTTTCAGGCGTGCTCCTTCAGCCACCCGGACACCCCA‑CCAAATTTTCGTTCCTTTC‑TTTT‑TGGGACGGGCGCCAATTTT‑GGGGAATCAtt                              >  2:363364/1‑149 (MQ=14)                                                                                                                         tcttTAAATATGGC‑GGTGAGGGGGGGATTCGAACCCCCGATACGT‑TGCCGTATACACA‑C‑TTTCCAGGCGTGCTCCTTCAGCCACTCGGACACCTCA‑CCAAATTGTCGTTCCTGTC‑TTGC‑TGGAACGGGCGCTAATTTA‑GGGAAATCATg                              <  1:1431971/149‑1 (MQ=255)                                                                                                                             tAAATATGGC‑GGTGAGGGGGGGATTCGAACCCCCGATACGT‑TGCCGTATACACA‑C‑TTTCCAGGCGTGCTCCTTCAGCCACTCGGACACCTCA‑CCAAATTGTCGTTCCTGTC‑TTGC‑TGGAACGGGCGCTAATTTA‑GGGAAATCATGACCTg                         <  1:988350/150‑1 (MQ=255)                                                                                                                             tAAATATGGC‑GGTGAGGGGGGGATTCGAACCCCCGATACGT‑TGCCGTATACACA‑C‑TTTCCAGGCGTGCTCCTTCAGCCACTCGGACACCTCA‑CCAAATTGTCGTTCCTGTC‑TTGC‑TGGAACGGGCGCTAATTTA‑GGGAAATCATGACCTg                         <  1:650248/150‑1 (MQ=255)                                                                                                                             tAAATATGGC‑GGTGAGGGGGGGATTCGAACCCCCGATACGT‑TGCCGTATACACA‑C‑TTTCCAGGCGTGCTCCTTCAGCCACTCGGACACCTCA‑CCAAATTGTCGTTCCTGTC‑TTGC‑TGGAACGGGCGCTAATTTA‑GGGAAATCATGACCTg                         <  1:1253104/150‑1 (MQ=255)                                                                                                                              aaaTATGGC‑GGTGAGGGGGGGGTTTGAACCCCCCATACGG‑TTCCGTATACAAA‑C‑TTTTCCGGCGGGGTCCTTCAGCCCCCCCGGCACCCCA‑CCCAAATTTTGTTTCTTTC‑TTTC‑TTGGAAGGGCGCTAATTTATGGGaaaaatgaaccg                          >  1:803346/1‑140 (MQ=2)                                                                                                                              aaaTATGGC‑GGTGAGGGGGGGGGTTGAAACCCCCATACGT‑TGCCGTATAACCA‑C‑CTTTCCGGGGGGGTCCTTCAGCCCCCCGGACACCCCC‑CCCAATTTTCGTTTCTGTC‑TTGC‑TGGAACGGGCGCTCATTTT‑TGGGAAACATGACCTgcg                       >  1:434107/1‑149 (MQ=12)                                                                                                                              aaaTATGGC‑GGTGAGGGGGGGATTCGAACCCCCGATACGT‑TGCCGTATACACA‑C‑TTTCCAGGCGTGCTCCTTCAGCCACTCGGACACCTCA‑CCAAATTGTCGTTCCTGTC‑TTGC‑TGGAACGGGCGCTAATTTA‑GGGAAATCATGACCTGa                        <  2:612828/150‑1 (MQ=255)                                                                                                                              aaaTATGGC‑GGTGAGGGGGGGACTCGAACCCCCGATACGT‑TGCCGTATACACA‑C‑TTTCCAGGCGTGCTCCTTCAGCCACTCGGACACCTCA‑CCAAATTGTCGTTCCTGTC‑TTGC‑TGGAACGGGCGCTAATTTA‑GGGAAATCATGACCTGa                        <  2:499905/150‑1 (MQ=255)                                                                                                                                atatGGC‑GGTGAGGGGGGGATTCGAACCCCCGATACGT‑TGCCGTATACACA‑C‑TTTCCAGGCGTGCTCCTTCAGCCACTCGGACACCTCA‑CCAAATTGTCGTTCCTGTC‑TTGC‑TGGAACGGGCGCTAATTTA‑GGGAAATCATGACCTGAGGt                     >  1:799039/1‑151 (MQ=255)                                                                                                                                    ggC‑GGTGAGGGGGGGATTCGAACCCCCGATACGT‑TGCCGTATACACA‑C‑TTTCCAGGCGTGCTCCTTCAGCCACTCGGACACCTCA‑CCAAATTGTCGTTCCTGTC‑TTGC‑TGGAACGGGCGCTAATTTA‑GGGAAATCATGACCTGAGGTcaac                 <  2:1115633/151‑1 (MQ=255)                                                                                                                                      c‑GGTGAGGGGGGGATTCGAACCCCCGATACGT‑TGCCGTATACACA‑C‑TTTCCAGGCGTGCTCCTTCAGCCACTCGGACACCTCA‑CCAAATTGTCGTTCCTGTC‑TTGC‑TGGAACGGGCGCTAATTTA‑GGGAAATCATGACCTGAGGTcaaca                >  2:561247/1‑150 (MQ=255)                                                                                                                                      c‑GGTGAGGGGGGGATTCGAACCCCCGATACGT‑TGCCGTATACACA‑C‑TTTCCAGGCGTGCTCCTTCAGCCACTCGGACACCTCA‑CCAAATTGTCGTTCCTGTC‑TTGC‑TGGAACGGGCGCTAATTTA‑GGGAAATCATGACCTGAGGTcaaca                >  2:900217/1‑150 (MQ=255)                                                                                                                                      c‑GGTGA‑GGGGGGATTCGAACCCCCGATACGT‑TGCCGTATACACA‑C‑TTTCCAGGCGTGCTCCTTCAGCCACTCGGACACCTCA‑CCAAATTGTCGTTCCTGTC‑TTGC‑TGGAACGGGCGCTAATTTA‑GGGAAATCATGACCTGAGGTCAACaaa              <  1:75623/151‑1 (MQ=255)                                                                                                                                         gTGAGGGGGGGATTCGAACCCCCGATACGT‑TGCCGTATACACA‑C‑TTTCCAGGCGTGCTCCTTCAGCCACTCGGACACCTCA‑CCAAATTGTCGTTCCTGTC‑TTGC‑TGGAACGGGCGCTAATTTA‑GGGAAATCATGACCTGAGGTCAACAAAc             >  2:1015904/1‑151 (MQ=255)                                                                                                                                         gTGAGGGGGGGATTCGAACCCCCGATACGT‑TGCCGTATACACA‑C‑TTTCCAGGCGTGCTCCTTCAGCCACTCGGACACCTCA‑CCAAATTGTCGTTCCTGTC‑TTGC‑TGGAACGGGCGCTAATTTA‑GGGAAATCATGACCTGAGGTCAACAAAc             <  1:1055466/151‑1 (MQ=255)                                                                                                                                           gAGGGGGGGATTCGAACCCCCGATACGT‑TGCCGTATACACA‑C‑TTTCCAGGCGTGCTCCTTCAGCCACTCGGACACCTCA‑CCAAATTGTCGTTCCTGTC‑TTGC‑TGGAACGGGCGCTAATTTA‑GGGAAATCATGACCTGAGGTCAACAAACtt           <  2:1312738/151‑1 (MQ=255)                                                                                                                                             gggggggATTCGAACCCCCGATACGT‑TGCCGTATACACA‑C‑TTTCCAGGCGTGCTCCTTCAGCCACTCGGACACCTCA‑CCAAATTGTCGTTCCTGTC‑TTGC‑TGGAACGGGCGCTAATTTA‑GGGAAATCATGACCTGAGGTCAACAAACtttt         >  2:407678/1‑151 (MQ=255)                                                                                                                                             gggggggATTCGAACCCCCGATACGT‑TGCCGTATACACA‑C‑TTTCCAGGCGTGCTCCTTCAGCCACTCGGACACCTCA‑CCAAATTGTCGTTCCTGTC‑TTGC‑TGGAACGGGCGCTAATTTA‑GGGAAATCATGACCTGAGGTCAACAAACtttt         >  1:194883/1‑151 (MQ=255)                                                                                                                                              ggggggATTCGAACCCCCGATACGT‑TGCCGTATACACA‑C‑TTTCCAGGCGTGCTCCTTCAGCCACTCGGACACCTCA‑CCAAATTGTCGTTCCTGTC‑TTGC‑TGGAACGGGCGCTAATTTA‑GGGAAATCATGACCTGAGGTCAACAAACttttt        >  1:284192/1‑151 (MQ=255)                                                                                                                                                   gATTCGAACCCCCGATACGT‑TGCCGTATACACA‑C‑TTTCCAGGCGTGCTCCTTCAGCCACTCGGACACCTCA‑CCAAATTGTCGTTCCTGTC‑TTGC‑TGGAACGGGCGCTAATTTA‑GGGAAATCATGACCTGAGGTCAACAAACTTTTTGaaa    >  1:929890/1‑150 (MQ=255)                                                                                                                                                    aTTCGAACCCCCGATACGT‑TGCCGTATACACA‑C‑TTTCCAGGCGTGCTCCTTCAGCCACTCGGACACCTCA‑CCAAATTGTCGTTCCTGTC‑TTGC‑TGGAACGGGCGCTAATTTA‑GGGAAATCATGACCTGAGGTCAACAAACTTTTTGaaaaa  <  1:1473710/151‑1 (MQ=255)                                                                                                                                                    |                                                                                                                                                              CACTGCGTTCACGGTCTGTCCAACTGGCTGTCGCCAGTTGTCGAACCCCGGTCGGGGCTTCTCATCCCCCCGGTGTGTGCAATATACG‑AAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGC‑GGTGAGGGGGGGATTCGAACCCCCGATACGT‑TGCCGTATACACA‑C‑TTTCCAGGCGTGCTCCTTCAGCCACTCGGACACCTCA‑CCAAATTGTCGTTCCTGTC‑TTGC‑TGGAACGGGCGCTAATTTA‑GGGAAATCATGACCTGAGGTCAACAAACTTTTTGAAAAA  >  NC_000913/1097436‑1097731

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