Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,097,581	0	A	G	42.6%	101.9 / 20.1	63	noncoding (72/88 nt)	serX	tRNA‑Ser
Reads supporting (aligned to +/- strand):  ref base A (5/30);  new base G (26/0);  total (33/30)
Fisher's exact test for biased strand distribution p-value = 1.62e-12
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CACTGCGTTCACGGTCTGTCCAACTGGCTGTCGCCAGTTGTCGAACCCCGGTCGGGGCTTCTCATCCCCCCGGTGTGTGCAATATACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGT‑GA‑GGGGGGGATTCG‑AACCCCCGATACGTTGCCGTATACACACTTTCCAGGCGTGCTCCTTCAGCCACTCGGACACCTCACCAAATTGTCGTTCCTGTCTTGCTGGAACGGGCG‑CTAATTTAGGGAAATCATGACCTGAGGTCAACAAACTTTTTGAAAA  >  NC_000913/1097436‑1097730                                                                                                                                                    |                                                                                                                                                        cACTGCGTTCCCGGTCTGTCCAACTGGCTGTCGCCAGTTGTCGAACCCCGGTCGGGGCTTCTCATCCCCCCGGTGTGTGCAATATACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATCGCGGT‑GA‑GGGGGGGATTc                                                                                                                                                      <  2:155768/149‑1 (MQ=11)  aCTGCGTTCACGGTCTGTCCAACTGGCTGTCGCCAGTTGTCGAACCCCGGTCGGGGCTTCTCATCCCCCCGGTGTGTGCAATATACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGT‑GA‑GGGGGGGATTCG‑aa                                                                                                                                                  <  1:1123062/151‑1 (MQ=32)      cGTTCACGGTCTGTCCAACTGGCTGTCGCCAGTTGTCGAACCCCGGTCGGGGCTTCTCATCCCCCCGGTGTGTGCAATATACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGT‑GA‑GGGGGGGGGTCG‑AAcccc                                                                                                                                              >  1:869621/1‑151 (MQ=18)      cGTTCACGGTCTGTCCAACTGGCTGTCGCCAGTTGTCGAACCCCGGTCGGGGCTTCTCATCCCCCCGGTGTGTGCAATATACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGT‑GA‑GGGGGGGATTCG‑AAcccc                                                                                                                                              <  2:1009121/151‑1 (MQ=35)       gTTCACGGTCTGTCCAACTGGCTGTCGCCAGTTGTCGAACCCCGGTCGGGGCTTCTCATCCCCCCGGTGTGTGCAATATACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGT‑GA‑GGGGGGGATTCG‑AAccccc                                                                                                                                             <  1:137918/151‑1 (MQ=35)        ttCACGGTCTGTCCAACTGGCTGTCGCCAGTTGTCGAACCCCGGTCGGGGCTTCTCATCCCCCCGGTGTGTGCAATATACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGT‑GA‑GGGGGGGGGTCG‑AAccccc                                                                                                                                             >  1:643225/1‑150 (MQ=18)         tCACGGTCTGTCCAACTGGCTGTCGCCAGTTGTCGAACCCCGGTCGGGGCTTCTCATCCCCCCGGTGTGTGCAATATACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGT‑GA‑GGGGGGGATTCG‑AACCCCCg                                                                                                                                            <  1:324626/150‑1 (MQ=37)             ggTCTGTCCAACTGGCTGTCGCCAGTTGTCGAACCCCGGTCGGGGCTTCTCATCCCCCCGGTGTGTGCAATATACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGC‑GG‑AGGGGGGGGTCG‑AAccccc                                                                                                                                             >  2:428100/1‑145 (MQ=255)                   tCCAACTGGCTGTCGCCAGTTGTCGAACCCCGGTCGGGGCTTCTCATCCCCCCGGTGTGTGCAATATACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGT‑GA‑GGGGGGGATTCG‑AACCCCCGATACGTTGc                                                                                                                                   <  1:867986/149‑1 (MQ=255)                    ccAACTGGCTGTCGCCAGTTGTCGAACCCCGGTCGGGGCTTCTCATCCCCCCGGTGTGTGCAATATACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGT‑GA‑GGGGGGGGTTTG‑GAccccccctaccttttc                                                                                                                                  >  1:922412/1‑138 (MQ=11)                      aaCTGGCTTTCGCCAGTTGTCGAACCCCGGCCGGGGCTTCTCATCCCCCCGGTGTGTGCAATATACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGT‑GA‑GGGGGGGATTCG‑AACCCCCGATACGTTGCCGtat                                                                                                                              <  1:234668/151‑1 (MQ=255)                        cTGGCTGTCGCCAGTTGTCGAACCCCGGTCGGGGCTTCTCATCCCCCCGGTGTGTGCAATATACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGT‑GA‑GGGGGGGATTCG‑AACCCCCGATACGTTGCCGtata                                                                                                                             <  2:873229/150‑1 (MQ=255)                           gCTGTCGCCAGTTGTCGAACCCCGGTCGGGGCTTCTCATCCCCCCGGTGTGTGCAATATACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGG‑GA‑GGGGGGGGTTTG‑AGCCCCCCCCACGTTTCCGCATtca                                                                                                                           >  2:530503/1‑149 (MQ=11)                             tGTCGCCAGTTGTCGAACCCCGGTCGGGGCTTCTCATCCCCCCGGTGTGTGCAATATACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGT‑GA‑GGGGGGGGGTCG‑AACCCCCCACACGGTTCCGTATaca                                                                                                                           >  1:473417/1‑147 (MQ=14)                              gTCGCCAGTTGTCGAACCCCGGTCGGGGCTTCTCATCCCCCCGGTGTGTGCAATATACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGT‑GG‑GGGGGGGGTTTG‑AAACCCCCCTACGTTTCCGCATAAACACtt                                                                                                                      >  1:1015922/1‑151 (MQ=11)                              gTCGCCAGTTGTCGAACCCCGGTCGGGGCTTCTCATCCCCCCGGTGTGTGCAATATACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGT‑GA‑GGGGGGGATTCG‑AACCCCCGATACGTTGCCGTATacacac                                                                                                                        <  2:218723/149‑1 (MQ=255)                              gTCGCCAGTTGTCGAACCCCGGTCGGGGCTTCTCATCCCCCCGGTGTGTGCAATATACG‑AAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGT‑GA‑GGGGGGGATTCG‑AACCCCCGATACGTTGCCGTATACACACtt                                                                                                                      <  2:243457/150‑1 (MQ=255)                                  ccAGTTGTCGAACCCCGGTCGGGGCTTCTCATCCCCCCGGTGTGTGCAATATACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGT‑GA‑GGGGGGGGGTCG‑AACCCCCCCCACCTTTCCGCCTACACACCTTcc                                                                                                                   >  2:355053/1‑150 (MQ=11)                                     gTTGTCGAACCCCGGTCGGGGCTTCTCATCCCCCCGGTGTGTGCAATATACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGT‑GA‑GGGGGGGATTCG‑AACCCCCGATACGTTGCCGTATACACACTTTCCAGGc                                                                                                               <  1:789581/151‑1 (MQ=255)                                       tGTCGAACCCCGGTCGGGGCTTCTCATCCCCCCGGTGTGTGCAATATACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGT‑GA‑GGGGGGGATTCG‑AACCCCCGATACGTTGCCGTATACACACTTTCCAGGCg                                                                                                              <  1:404372/150‑1 (MQ=255)                                        gTCGAACCCCGGTCGGGGCTTCTCATCCCCCCGGTGTGTGCAATATACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGT‑GA‑GGGGGGGATTCG‑AACCCCCGATACGTTGCCGTATACACACTTTCCAGGCGTg                                                                                                            <  1:882630/151‑1 (MQ=255)                                        gTCGAACCCCGGTCGGGGCTTCTCATCCCCCCGGTGTGTGCAATATACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGT‑GA‑GGGGGGGATTCG‑AACCCCCGATACGTTGCCGTATACACACTTTCCAGGCGTg                                                                                                            <  2:243425/151‑1 (MQ=255)                                        gTCGAACCCCGGTCGGGGCTTCTCATCCCCCCGGTGTGTGCAATATACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGC‑GG‑GGGGGGGGGTCG‑AAACCCCCAATCGCTTCCGCTTACACACttttcccgccgg                                                                                                            >  2:756955/1‑142 (MQ=11)                                                     cGGGGCTTCTCATCCCCCCGGTGTGTGCAATATACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGT‑GA‑GGGGGGGGTTCG‑AACCCCCGCTACCTTTCCGTATACACACTTTCCCCGCGTGCTcccttcgcc                                                                                                 >  2:857413/1‑142 (MQ=17)                                                       gggCTTCTCATCCCCCCGGTGTGTGCAATATACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGT‑GA‑GGGGGGGGTTCG‑AACCCCCCGCACGTTGTCGTCTAAACACTTTCCAGGGGTGCTCCTTCCGcccccc                                                                                             >  2:462062/1‑147 (MQ=17)                                                                aTCCCCCCGGTGTGTGCAATATACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGT‑GA‑GGGGGGGATTCG‑AACCCCCGATACGTTGCCGTATACACACTTTCCAGGCGTGCTCCTTCAGCCACTCGGACACCTc                                                                                    <  1:876124/151‑1 (MQ=255)                                                                   cccccGGTGTGTGCAATATACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGT‑GA‑GGGGGGGGTTCG‑AACCCCCCCGACCTTGCCCCATACACACCTTCCCGGCGTGCTCCCTCCGCCCCCCCGACACCCCAcc                                                                                 >  1:572068/1‑151 (MQ=17)                                                                     cccGGTGTGTGCAATATACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGT‑GA‑GGGGGGGATTCG‑AACCCCCGATACGTTGCCGTATACACACTTTCCAGGCGTGCTCCTTCAGCCACTCGGACACCTCACCaa                                                                               <  2:269660/151‑1 (MQ=255)                                                                        ggTGTGTGCAATATACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGT‑GA‑GGGGGGGATTCG‑AACCCCCGATACGTTGCCGTATACACACTTTCCAGGCGTGCTCCTTCAGCCACTCGGACACCTCACCAAAtt                                                                            <  1:833366/151‑1 (MQ=255)                                                                        ggTGTGTGCAATATACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGA‑GA‑GGGGGGGGTTTT‑AAACCCCCAGAACTTTCCGCATAAACACATTTCCGGGGGGCTCCTTCTCCCACTCGGCGACCCCCCCaa                                                                               >  2:24186/1‑148 (MQ=2)                                                                                cAATATACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGC‑GA‑GGGGGGGGTTTC‑AAACCCCCATCCCTTTTCCTATACACACATTTCAAGGGTGCGCCTTCCCCCCCCCCGACACCCCCCCCAAATATCGtt                                                                      >  2:1110363/1‑149 (MQ=2)                                                                                 aaTATACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGT‑GA‑GGGGGGGGGTTG‑AACCCCCCCTACGTTGCCGGATACAAACTTTTCAGGCGGGCTCCTTCAGCCACCCCGGCACCCCACCCAAATGTCGtt                                                                      >  1:805560/1‑148 (MQ=14)                                                                                  atatACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGT‑GA‑GGGGGGGGTTCG‑AACCCCCCCTACGTTTCCGTATACACACCTTCCAGGCGTGCCCCTTCAGCCACCCGGACACCCCCCCAAATTTTCGTTcc                                                                    >  1:105523/1‑149 (MQ=14)                                                                                   tataCGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGG‑GA‑GGGGGGGGTTCC‑AACCCCCGATACGTTTCCGTATACACCCTTTCCACGGGTGCTCCTTCAACCACCCGGAAACCCCCCCAAATTTTCGTTTCtt                                                                  >  1:449012/1‑149 (MQ=14)                                                                                       cGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGT‑GA‑GGGGGGGGGTCG‑CACCCCCCCTAAGTTGCCGTATACAAACCTTTCAGGCGTGCTCCCTCCGCCCCCCCGGCACCCCACCCAAATGTCGTTTCTGTCtt                                                              >  1:939612/1‑150 (MQ=14)                                                                                       cGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGA‑GA‑GGGGGGGGTTTG‑AACCCCCCCGACCTTTGCGTATACACACTTTTCCGGGGGGGGCCTTCAGCCACCCCGGCACCTCCCCCAAATTTCGTTCCTGTCttt                                                             >  2:472112/1‑150 (MQ=9)                                                                                        gAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGT‑GA‑GGGGGGGGTTCG‑AACCCCCCCTACGTTGCCGGATAAACACTTTCCCGGGGTGCTCCTTCAGCCACCCCGACACCCCCCCAAAATTTCGTTCCCGTCtt                                                              >  1:554032/1‑149 (MQ=14)                                                                                        gAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGA‑GAGGGGGGGGTTTCG‑AACCCCC‑ATACGTTTTCGTTTAAACAAATTCCCGGGGGGGGCCTTCAACCACCCCGGCACCCCCCCACATTTTTCTTCCCGTCtt                                                              >  2:312992/1‑149 (MQ=2)                                                                                              aaaGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGT‑GT‑GGAGGGGGTACG‑TAACCCCGATTCCTTTTCACGTAAAGAATTTTCTGGCGTGCTGCTTCCGCCACTTCGGCACCttccccacatatc                                                                         >  2:524494/1‑120 (MQ=2)                                                                                                   ccGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGT‑GA‑GGGGGGGGGTTT‑AACCCCCCCTCCGTTGCCGTAGACACACCTTTCAGGCGTGCGCCTTCAGCCCCCCGGCCACCCTCCCAAAATATCGTTTCTGTCTTTCTGGGACggg                                                   >  2:649/1‑149 (MQ=9)                                                                                                      tACTTTCGTACGAGCTCTTCTTTAAATATGGCGGT‑GA‑GGGGGGGGTTTG‑AACCCCCGCTAAGTTTCCGTATACACACCTTTCAGGCGGGCTCCTTCACCCACCCGGACACCCCCCCAAAATTTCGTTTCTGTCttttttggacggggg                                                 >  2:738208/1‑135 (MQ=2)                                                                                                          ttCGTACGAGCTCTTCTTTAAATATGGCGGTGGG‑GGGGGGGTTTCG‑AACCCCC‑ATACCTTTCCGTATACACACTTTTCCGGCGGGGTCCTTTCCCCCCCCCGACACCCCACCAAAATTTCGTTTCTTTCTTTCTGGGACGGGCG‑CCAAttt                                         >  2:177457/1‑151 (MQ=9)                                                                                                          ttCGTACGAGCTCTTCTTTAAATATGGCGGT‑GA‑GGGGGGGATTCG‑AACCCCCGATACGTTGCCGTATACACACTTTCCAGGCGTGCTCCTTCAGCCACTCGGACACCTCACCAAATTGTCGTTCCTGTCTTGCTGGAACGGGCG‑CTAAt                                           <  2:354297/149‑1 (MQ=255)                                                                                                             gTACGAGCTCTTCTTTAAATATGGCGGT‑GA‑GGGGGGGGTTCGAAACCCCCCACAGGTTG‑CGTATACACACCTTTCAGGCGGGCCCCTTCAGCCACCCCGACAACCCCCCCAATTTTCGTTTCTGTCTTGGTGGGAAGGGCG‑CCAATTTTgg                                      >  1:728983/1‑151 (MQ=9)                                                                                                               aCGAGCTCTTCTTTAAATATGGCGGT‑GA‑GGGGGGGATTCG‑AACCCCCGATACGTTGCCGTATACACACTTTCCAGGCGTGCTCCTTCAGCCACTCGGACACCTCACCAAATTGTCGTTCCTGTCTTGCTGGAACGGGCG‑CTAATTTAGGGa                                    <  2:457463/151‑1 (MQ=255)                                                                                                                   gCTCTTCTTTAAATATGGCGGT‑GA‑GGGGGGGGTTCG‑AAACCCCCATACGTTGCCGTATACACACTTTTCACGCGGGCTCCTTCAGCCCCCCCGACACCCCCCCCAAATTTCGTTTCCTTCTTTCTGGGAAGGGCGCCAAATTTTGGGaa                                   >  1:891135/1‑149 (MQ=2)                                                                                                                        tcttTAAATATGGCGGT‑GA‑GGGGGGGGTTCG‑AAACCCCGCTACGTTGCCGTATAAACACTTTCCAGGCGTGGTCCTTCCGCCACTCGGGCACCCCCCCCAATTTTCGTTCCTTTCTTTCTGGAACGGGCG‑CCAATTTTGGGaaa                                  >  1:866288/1‑144 (MQ=17)                                                                                                                        tcttTAAATATGGCGGT‑GA‑GGGGGGGATTCG‑AACCCCCGATACGTTGCCGTATACACACTTTCCAGGCGTGCTCCTTCAGCCACTCGGACACCTCACCAAATTGTCGTTCCTGTCTTGCTGGAACGGGCG‑CTAATTTAGGGAAATCAt                              <  2:672700/148‑1 (MQ=255)                                                                                                                        tcttTAAATATGGCGGT‑GA‑GGGGGGGATTCG‑AACCCCCGATACGTTGCCGTATACACACTTTCCAGGCGTGCTCCTTCAGCCACTCGGACACCTCACCAAATTGTCGTTCCTGTCTTGCTGGAACGGGCG‑CTAATTTAGGGAAATCATGa                            <  1:234563/150‑1 (MQ=255)                                                                                                                            tAAATATGGCGGT‑GA‑GGGGGGGATTCG‑AACCCCCGATACGTTGCCGTATACACACTTTCCAGGCGTGCTCCTTCAGCCACTCGGACACCTCACCAAATTGTCGTTCCTGTCTTGCTGGAACGGGCG‑CTAATTTAGGGAAATCATGACCTg                        <  2:1015662/150‑1 (MQ=255)                                                                                                                               atatGGCGGT‑GA‑GGGGGGGATTCG‑AACCCCCGATACGTTGCCGTATACACACTTTCCAGGCGTGCTCCTTCAGCCACTCGGACACCTCACCAAATTGTCGTTCCTGTCTTGCTGGAACGGGCG‑CTAATTTAGGGAAATCATGACCTGAgg                     >  1:682153/1‑150 (MQ=255)                                                                                                                               atatGGCGGT‑GA‑GGGGGGGATTCG‑AACCCCCGATACGTTGCCGTATACACACTTTCCAGGCGTGCTCCTTCAGCCACTCGGACACCTCACCAAATTGTCGTTCCTGTCTTGCTGGAACGGGCG‑CTAATTTAGGGAAATCATGACCTGAGGt                    <  2:836063/151‑1 (MQ=255)                                                                                                                                  tGGCGGT‑GA‑GGGGGGGATTCG‑AACCCCCGATACGTTGCCGTATACACACTTTCCAGGCGTGCTCCTTCAGCCACTCGGACACCTCACCAAATTGTCGTTCCTGTTTTGCTGGAACGGGGG‑CTAATTTAGGGACATCATGACCTGAGGTca                  >  1:1079026/1‑150 (MQ=255)                                                                                                                                      ggT‑GA‑GGGGGGGATTCG‑AACCCCCGATACGTTGCCGTATACACACTTTCCAGGCGTGCTCCTTCAGCCACTCGGACACCTCACCAAATTGTCGTTCCTGTCTTGCTGGAACGGGCG‑CTAATTTAGGGAAATCATGACCTGAGGTCAACaaa             <  1:899356/151‑1 (MQ=255)                                                                                                                                          gA‑GGGGGGGATTCG‑AACCCCCGATACGTTGCCGTATACACACTTTCCAGGCGTGCTCCTTCAGCCACTCGGACACCTCACCAAATTGTCGTTCCTGTCTTGCTGGAACGGGCG‑CTAATTTAGGGAAATCATGACCTGAGGTCAACAAACt           <  1:464022/150‑1 (MQ=255)                                                                                                                                          gA‑GGGGGGGATTCG‑AACCCCCGATACGTTGCCGTATACACACTTTCCAGGCGTGCTCCTTCAGCCACTCGGACACCTCACCAAATTGTCGTTCCTGTCTTGCTGGAACGGGCG‑CTAATTTAGGGAAATCATGACCTGAGGTCAACAAACt           <  2:869405/150‑1 (MQ=255)                                                                                                                                              ggggggATTCG‑AACCCCCGATACGTTGCCGTATACACACTTTCCAGGCGTGCTCCTTCAGCCACTCGGACACCTCACCAAATTGTCGTTCCTGTCTTGCTGGAACGGGCG‑CTAATTTAGGGAAATCATGACCTGAGGTCAACAAACttttt       <  2:243070/151‑1 (MQ=255)                                                                                                                                                ggggATTCG‑AACCCCCGATACGTTGCCGTATACACACTTTCCAGGCGTGCTCCTTCAGCCACTCGGACACCTCACCAAATTGTCGTTCCTGTCTTGCTGGAACGGGCG‑CTAATTTAGGGAAATCATGACCTGAGGTCAACAAACTTTTTGa     <  1:50314/151‑1 (MQ=255)                                                                                                                                                 gggATTCG‑ACCCCCCGATACGTTGCCGTATACACACTTTCCAGGCGTGCTCCTTCAGCCACTCGGACACCTCACCAAATTGTCGTTCCTGTCTTGCTGGAACGGGCG‑CTAATTTAGGGAAATCATGACCTGAGGTCAACAAACTTTTTGaa    <  2:755734/151‑1 (MQ=255)                                                                                                                                                 gggATTCG‑AACCCCCGATACGTTGCCGTATACACACTTTCCAGGCGTGCTCCTTCAGCCACTCGGACACCTCACCAAATTGTCGTTCCTGTCTTGCTGGAACGGGCG‑CTAATTTAGGGAAATCATGACCTGAGGTCAACAAACTTTTTGa     >  2:418727/1‑150 (MQ=255)                                                                                                                                                  ggATTCG‑AACCCCCGATACGTTGCCGTATACACACTTTCCAGGCGTGCTCCTTCAGCCACTCGGACACCTCACCAAATTGTCGTTCCTGTCTTGCTGGAACGGGCG‑CTAATTTAGGGAAATCATGACCTGAGGTCAACAAACTTTTTGaa    >  2:823600/1‑150 (MQ=255)                                                                                                                                                   gATTCG‑AACCCCCGATACGTTGCCGTATACACACTTTCCAGGCGTGCTCCTTCAGCCACTCGGACACCTCACCAAATTGTCGTTCCTGTCTTGCTGGAACGGGCG‑CTAATTTAGGGAAATCATGACCTGAGGTCAACAAACTTTTTGaaaa  <  1:68095/151‑1 (MQ=255)                                                                                                                                                   gATTCG‑AACCCCCGATACGTTGCCGCATACACACTTGCCAGGCGTGCTCCTTCAGCCACTCGGACACCTCACCAAATTGTCGTTCCTGTCTTGCTGGAACGGGCG‑CTAATTTAGGGAACTCATGACCTGAGGTCAACAAACTTTTTGaaa   >  2:607290/1‑150 (MQ=255)                                                                                                                                                    |                                                                                                                                                        CACTGCGTTCACGGTCTGTCCAACTGGCTGTCGCCAGTTGTCGAACCCCGGTCGGGGCTTCTCATCCCCCCGGTGTGTGCAATATACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGT‑GA‑GGGGGGGATTCG‑AACCCCCGATACGTTGCCGTATACACACTTTCCAGGCGTGCTCCTTCAGCCACTCGGACACCTCACCAAATTGTCGTTCCTGTCTTGCTGGAACGGGCG‑CTAATTTAGGGAAATCATGACCTGAGGTCAACAAACTTTTTGAAAA  >  NC_000913/1097436‑1097730

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