Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,947,591	0	A	.	44.4%	176.3 / 161.3	90	intergenic (+18/‑16)	metW/metV	tRNA‑Met/tRNA‑Met
Reads supporting (aligned to +/- strand):  ref base A (45/5);  new base . (35/5);  total (80/10)
Fisher's exact test for biased strand distribution p-value = 7.46e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.31e-01

GCAACCAATTAAAATT‑TGATGAAGTAAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAA‑TC‑AAATTTGATGAAGTAAAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTATTGAACACCCT‑AACGGGTGTTTTTTTGTTTCTGGTCTCCCATAAAAAAGCGCC  >  NC_000913/2947457‑2947740                                                                                                                                          |                                                                                                                                                       gggaagcaaaaaaaTTGTGATGAATTAAAGCACTCCGCTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAA‑TC‑AAATTTGATGAAGTAAAAGCAGTACGGTGa                                                                                                                                         <  3:302628/141‑1 (MQ=17)                           aacggattACGGGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGTCTCATAACCCGAAGGTGGTCGGTTGTAATCCGGCCCCGGCAACCAA‑TA‑AAATTTGATGAAGTAAAAGGAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCtc                                                                                                               <  2:107874/144‑1 (MQ=11)                            acccactACGCTGACGAGGTGGGGAGCAGCCTGGTCGCTCGGCGGGCTCAAAACCCGAAGGTCGTCGGTGCAAAACCGCCCCCCGCAACCAA‑TC‑AAATTTGATGAACTAAAACCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCtcg                                                                                                              <  2:717270/145‑1 (MQ=14)                                                               aGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAA‑TC‑AAATTTCAGGAAGTAAAAGCAGTAAGGTGAAGCGGGGTGGAGCACCATGGTAGCTCGTCGGGCTCATAAACCGAAGGTCGTCGGTTCAAATc                                                                           >  2:403193/1‑151 (MQ=14)                                                                       gggCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAA‑TC‑AAATTTGATGAAGTAAAAGCAGTACGGTGAcg                                                                                                                                       <  2:1471653/83‑1 (MQ=255)                                                                           tCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAA‑TC‑AAATTTGAGGAAGTAAACGCAGTAGTGTGACGCGGGGTGGAGCAGCATGGTAGCTCGTCGGGCTGTTAACGCGAAGGTCGGCGGGGCCCCTCCGGCCCgcgcc                                                                >  2:707118/1‑149 (MQ=11)                                                                                   cGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAA‑TT‑AAATTTTATGTATTAAAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGgtc                                                                                        <  1:1125415/118‑1 (MQ=255)                                                                                    gAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAA‑TC‑AAATTTCAAGAAGTAAAAGCAGTAAGTTGAAGCGGTGTCCAGAAGACTCGAACCTATTAGGGCTCATAAAACGAAGGTCTAACGGACCAATCCGCCCCCCGCAACAAattatt                                                      >  3:300944/1‑151 (MQ=9)                                                                                                   gtgctcttccgaggtaaCCAA‑TTAAAATTTGATCAAGT‑AAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTATTGAACACCCT‑AACggg                                      <  8:1274070/136‑1 (MQ=11)                                                                                                     aTCCGGCCCCCGCAACCAA‑TC‑AAATTTCATGAAGTAAAAGCAGTACGGTTACGAGGGGTGGATCAGAATCGTAGCTCGTCGGGCTCATAACCCGAAGGAAGTCGGTTAAAATCCGGCCAACGCAAACAATAATTGAACACCCT‑aaacgatg                                    >  6:588473/1‑145 (MQ=17)                                                                                                         ggCCCCCGCAACCAA‑TC‑AAATTTGATGAAGTAAAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCCACTCCGGCCCCCGCAACCAATTATTGAACACCCT‑AACGGGTGgttt                                >  3:1486209/1‑151 (MQ=255)                                                                                                         ggCCCCCGCAACCAA‑TC‑AAATTTCATCAAGTAAAAGCAGTACGGTGACCAGGGGTGGAGCAGCATCGAAGCTCGACGGGCTCATAACCCGAACGACGGACGTACAAATCCGGCCCACGCAACCAATAATTGAACACCCT‑AAAGGgtgata                                 >  4:361368/1‑147 (MQ=255)                                                                                                           cccccGCAACCAA‑TC‑AAATTTGATGAAGTAAAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTATTGAACACCCT‑AACGGGTGtttttt                              >  5:284377/1‑151 (MQ=255)                                                                                                            ccccGCAACCAA‑TC‑AAATTTGATGAAGTAAAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTATTGAACACCCT‑AACGGGTGttttttt                             >  6:1797995/1‑151 (MQ=255)                                                                                                            ccccGCAACCAA‑TC‑AAATTTGATGAAGTAAAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTATTGAACACCCT‑AACGGGTGttttttt                             >  8:1762891/1‑151 (MQ=255)                                                                                                            caccGCAACCAA‑TC‑AAATTTGATGAAGTAAAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCCAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTATTGAACACCCT‑AACGGGTGttttttt                             >  7:1605058/3‑151 (MQ=255)                                                                                                             cccGCAACCAA‑TC‑AAATTTGATGAAGTAAAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTATTGAACACCCT‑AACGGGTGTTTTTTTg                            >  3:382787/1‑151 (MQ=255)                                                                                                             cccGCAACCAA‑TC‑AAATTTGATGAAGTAAAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTATTGAACACCCT‑AACGGGTGTTTTTTTg                            >  2:721646/1‑151 (MQ=255)                                                                                                               cGCAACCAA‑TC‑AAATTTGATGAAGTAAAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTATTGAACACCCT‑AACGGGTGTTTTTTTGtt                          >  1:1188125/1‑151 (MQ=255)                                                                                                                gCAACCAATTA‑AAATTTGATGAAGT‑AAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATAAATTAACACCCA‑ATACGGTGtttttttttt                          >  6:49111/1‑150 (MQ=11)                                                                                                                gCAACCAA‑TC‑AAATTTGATGAAGTAAAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTAGGTTCAAATCCTGCACCCGCAACCAATTATTGAACACCCT‑AACGGGTGTTTTTTTGttt                         >  2:184769/1‑151 (MQ=255)                                                                                                                 cAACCAATTA‑AAATTTGATGAAGT‑AAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATAAATTAACACGAT‑CACGGCTGTTTAtttttta                         >  2:580809/1‑149 (MQ=11)                                                                                                                 cAACCAA‑TC‑AAATTTGATGAAGTAAAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATACGGCCCCCGCAACCAATTATTGAACACCCT‑AACGGGTGTTTTTTTGTTTc                        >  7:1479640/1‑151 (MQ=255)                                                                                                                  aaCCAA‑TTAAAATTTGATGAAGT‑AAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTATTGTACACGCT‑GACGGGGGTGTTTTTGgttgg                       >  2:1808967/1‑149 (MQ=1)                                                                                                                  aaCCAA‑TTAAAATTTGATGAAGT‑AAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAATTAACAACAA‑AAACGGTCTGTTTTTGTTTcg                       >  1:110692/1‑150 (MQ=11)                                                                                                                  aaCCAA‑TTAAAATTTGATGAAGT‑AAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGGCGGTTCAAATCCGGCCCCCGCAACCAATAATTGAACACCCT‑AAACCGTGTTTTTTGGtttgt                       >  8:15591/1‑149 (MQ=11)                                                                                                                  aaCCAA‑TC‑AAATTTGATGAAGTAAAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTATTGAACACCCT‑AACGGGTGTTTTTTTGTTTCt                       >  4:382540/1‑151 (MQ=255)                                                                                                                  aaCCAA‑TC‑AAATTTGATGAAGTAAAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTATTGAACACCCT‑AACGGGTGTTTTTTTGTTTCt                       >  4:7673/1‑151 (MQ=255)                                                                                                                  aaCCAA‑TC‑AAATTTGATGAAGTAAAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTATTGAACACCCT‑AACGGGTGTTTTTTTGTTTCt                       >  4:950449/1‑151 (MQ=255)                                                                                                                  aaCCAA‑TC‑AAATTTGATGAAGTAAAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTATTGAACACCCT‑AACGGGTGTTTTTTTGTTTCt                       >  5:168296/1‑151 (MQ=255)                                                                                                                  aaCCAA‑TC‑AAATTTGATGAAGTAAAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTATTGAACACCCT‑AACGGGTGTTTTTTTGTTTCt                       >  1:1584635/1‑151 (MQ=255)                                                                                                                  aaCCAA‑TC‑AAATTTGATGAAGTAAAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTATTGAACACCCT‑AACGGGTGTTTTTTTGTTTCt                       >  1:1032608/1‑151 (MQ=255)                                                                                                                    ccAA‑TC‑AAATTTGATGAAGTAAAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTATTGAACACCCT‑AACGGGTGTTTTTTTGTTTCTgg                     >  4:430986/1‑151 (MQ=255)                                                                                                                    ccAA‑TC‑AAATTTGATGAAGTAAAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTATTGAACACCCT‑AACCGGTGTTTTTTTGTTTCTgg                     >  8:1190769/1‑151 (MQ=255)                                                                                                                     gcaattatatttttATTAAGT‑AAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTATTGAACACCCT‑AACGGGTGTTTTTTTGtt                          <  1:960727/136‑1 (MQ=17)                                                                                                                     ccaattaaaaTTTGATGAAGT‑AAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTATTGAACACCCA‑AACGGGTGTTTTTTTGGTTCTgg                     >  3:359439/8‑151 (MQ=14)                                                                                                                     ccaattaaaaTTTGATGAAGT‑AAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAATTAACACCCA‑AACCGGGGTTTTTTTGtttgagg                     >  1:708046/8‑147 (MQ=11)                                                                                                                     ccaattaaaaTTTGATGAAGT‑AAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAATTAACAACAT‑AAAAGGTGTTTTATTGTTTCTgg                     >  1:1584395/8‑151 (MQ=11)                                                                                                                     ccaattaaaaTTTGATGAAGT‑AAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATAATAGAACACCCT‑AACACGTGATTATTAGtt                          >  7:962270/8‑146 (MQ=11)                                                                                                                     cAA‑TC‑AAATTTGATGAAGTAAAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTATTGAACACCCT‑AACGGGTGTTTTTTTGTTTCTGGt                    >  1:1388823/1‑151 (MQ=255)                                                                                                                     cAA‑TC‑AAATTTGATGAAGTAAAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTATTGAACACCCT‑AACGGGTGTTTTTTTGTTTCTGGt                    >  7:1683264/1‑151 (MQ=255)                                                                                                                     cAA‑TC‑AAATTTGATGAAGTAAAAGCAGCACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTATTGAACACCCT‑AACGGGTGTTTTTTTGTTTCTGGt                    >  8:858421/1‑151 (MQ=255)                                                                                                                      caattaaaaTTTGATGAAGT‑AAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTATTTAACAACCT‑AACGGGTGTTTATTTGtttgaggt                    >  1:978236/7‑146 (MQ=14)                                                                                                                      caattaaaaTTTGATGAAGT‑AAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTATTGAACAGCCT‑GACGGGTGTGTTTTGGTTTCgggg                    >  2:1240329/7‑150 (MQ=11)                                                                                                                      caattaaaaTTTGATGAAGT‑AAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAATGAACACCCT‑TACGCGGGTGTTTgtgtggggggg                    >  7:901023/7‑144 (MQ=11)                                                                                                                      aa‑TC‑AAATTTGATGAAGTAAAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCTACTCCGGCCCCCGCAACCTATTATTGAACACCCT‑AACGGGGGTTGTTTTGTTTCTGGtg                   >  2:192064/1‑150 (MQ=255)                                                                                                                      aa‑TC‑AAATTTGATGAAGTAAAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTATTGAACACCCT‑AACGGGTGTTTTTTTGTTTCTGGtc                   >  3:753611/1‑151 (MQ=255)                                                                                                                       aattaaaaTTTGATGAAGT‑AAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTATTGATCACCCT‑ATCGGGTGTTTTGTTGTTTGTgggg                   >  3:1410089/6‑149 (MQ=14)                                                                                                                        attaaaaTTTTATGAAGT‑AAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTATTGAACACCCT‑AACGGGTGTTTTTTTGttt                         <  6:1546221/140‑1 (MQ=255)                                                                                                                        attaaaaTTTGATGAAGT‑AAAGCAGTACGGTGACGCGGGGTGGCGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCACATCCGGCCCCCGCAACCAATGAATGTACAGCCT‑GTAGCGTGTTTTGTTGGATCTGGtgt                  >  8:884549/5‑149 (MQ=11)                                                                                                                        attaaaaTTTGATGAAGT‑AAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTATTTAACACCAT‑AACCCGGGTTTATGTGGTTCAGGtct                  >  7:1326154/5‑151 (MQ=11)                                                                                                                        attaaaaTTTGATGAAGT‑AAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATGATTGATCACCCT‑AACGGGTGTTTTGTTGTTTCTGGtct                  >  8:400887/5‑151 (MQ=14)                                                                                                                        attaaaaTTTGATGAAGT‑AAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATAAAAGAACACCCA‑AACGCGTCTTTTATTGTTTCTGGtct                  >  3:1031282/5‑151 (MQ=11)                                                                                                                        attaaaaTTTGATGAAGT‑AAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATACGGCCCCCGCAACCAATTAATGAACACCCA‑AAAggg                                      >  5:1731627/5‑131 (MQ=11)                                                                                                                        attaaaaTTTGATGAAGT‑AAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAAACCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTATTGAACACCCT‑AACGGGTGTGTTTTTGttt                         >  8:1549180/5‑144 (MQ=255)                                                                                                                         ttaaaaTTTGATGAAGT‑AAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAATTAACAACCT‑AAAGGAGCTTTTGTTGTTTGAGGtctc                 >  6:682397/4‑151 (MQ=11)                                                                                                                         ttaaaaTTTGATGAAGT‑AAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAATGAACACCCT‑GAAGCGGGTTTAATTGTTACTGGtctg                 >  2:668193/4‑150 (MQ=11)                                                                                                                         tC‑AAATTTGATGAAGTAAAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTATTGAACACCCT‑AACGGGTGTTTTTTTGTTTCTGGtctc                 >  3:950643/1‑151 (MQ=255)                                                                                                                         tC‑AAATTTGATGAAGTAAAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTATTGAACACCCT‑AACGGGTCTTTTTTTGTTTATGGtctc                 >  5:435219/1‑151 (MQ=255)                                                                                                                          c‑AAATTTGATGAAGTAAAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTATTGAACACCCT‑AACGGGTGTTTTTTTGTTTCTGGTCTcc                >  6:512115/1‑151 (MQ=255)                                                                                                                           aaaaTTTGATGAGGT‑AAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTATTGAACACCCT‑AACGGGTGTTTTTTTGTTTCTGGt                    <  7:38331/145‑1 (MQ=255)                                                                                                                           aaaaTTTGATGAAGT‑AAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATGATTGAACACGCTGAACGCG‑CATTTTTTGTTTCAggggtcgc               >  7:655356/2‑145 (MQ=11)                                                                                                                           aaaaTTTGATGAAGT‑AAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGCTCGTCGGTTCAAATCCGGCCCCCGCAACCAATAATATAACACCCT‑AAAGGGGCTTTTATTGTATCAGGt                    >  1:38257/2‑146 (MQ=11)                                                                                                                            aaaTTTGATGAAGTAAAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTATTGAACACCCT‑AACGGGTGTTTTTTTGTTTCTGGTCTccc               >  4:1098679/1‑151 (MQ=255)                                                                                                                            aaaTTTGATGAAGTAAAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTATTGAACACCCT‑AACGGGTGTTTTTTTGTTTCTGGTCTccc               >  6:757207/1‑151 (MQ=255)                                                                                                                            aaaTTTGATGAAGTAAAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTATTGAACACCCT‑AACGGGTGTTTTTTTGTTTCTGGTCTccc               >  6:1413123/1‑151 (MQ=255)                                                                                                                             aaTTTGATGAAGTAAAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTATTGAACACCCT‑AACGGGTGTTTTTTTGTTTCTGGTCTccc               >  3:1236002/1‑150 (MQ=255)                                                                                                                             aaTTTGATGAAGTAAAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTATTGAACACCCT‑AACGGGTGTTTTTTTGTTTCTGGTCTccc               >  8:550834/1‑150 (MQ=255)                                                                                                                             aaTTTGATGAAGT‑AAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTATTGAACACCCA‑AAAGCGTGTTTTTTTTTTTCTGGTCTCCCa              >  4:1695226/1‑150 (MQ=17)                                                                                                                              aTTTGATGAAGTAAAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTATTGAACACCCT‑AACGGGTGTTTTTTTGTTTCTGGTCTCCCAt             >  6:1396977/1‑151 (MQ=255)                                                                                                                              aTTTGATGAAGTAAAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCCCAACCAATTATTGAACACCCT‑AAAGGGTGTTTTTTTGTTTCTGGTCTCCCAt             >  2:491333/1‑151 (MQ=255)                                                                                                                              aTTTGATGAAGT‑AAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTATTGAACACGCT‑AACGGGTGTGTTTTTGTTTCTGGTCTCCCAt             >  8:1510234/1‑150 (MQ=255)                                                                                                                              aTTTGATGAAGT‑AAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATAATTGAACACCCT‑AACTGGTCTGTTTTTGTTACGGGTCTCCCAt             >  7:1026598/1‑150 (MQ=14)                                                                                                                               tttGATGAAGTAAAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTATTGAACACCCT‑AACGGGTGTTTTTTTGTTTCTGGTCTCCCAt             >  1:378677/1‑150 (MQ=255)                                                                                                                               tttGATGAAGT‑AAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTATTGAACACCCT‑AACGGGTGTTTTTTTGttt                         <  5:119318/137‑1 (MQ=255)                                                                                                                               tttGATGAAGT‑AAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTATATAACAACCT‑AACGCCTGATTATTTGtttagatcggagagcac           >  3:119571/1‑137 (MQ=11)                                                                                                                               tttGATGAAGT‑AAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATAATTTAACAACCT‑AAAGAGTGTTTTATTGTTTCAGGTCTCCAATa            >  6:1060611/1‑150 (MQ=14)                                                                                                                               tttGATGAAGT‑AAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATAATTGAACACCCT‑AACACGGCAGTATTTGTTACAGGTCTCCCATa            >  5:626687/1‑150 (MQ=11)                                                                                                                                ttGATGAAGTAAAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTATTGAACACCCT‑AACGGGTGTTTTTTTGTTTCTGGTCTCCCATa            >  6:10528/1‑150 (MQ=255)                                                                                                                                ttGATGAAGT‑AAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTATTGAACACCCT‑AACGGGTGTTTTTTTTTTTCTGGTCTCCCATaa           >  1:1173009/1‑150 (MQ=255)                                                                                                                                ttGATGAAGT‑AAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATAATTGAACACCAT‑AAAGGGTCTTTTTTTGTTTCTGGTCTCCCATaa           >  6:1150898/1‑150 (MQ=17)                                                                                                                                 tgatgaAGTAAAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTATTGAACACCCT‑AACGGGTGTTTTTTTGTTTCTGGTCTCCCATaa           >  4:737448/1‑150 (MQ=255)                                                                                                                                  gatgaAGT‑AAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATGAATGAACAACAA‑CACCGGGGTGTTGTGGTATCTGGTCTCCCaaaaaa         >  7:303323/1‑150 (MQ=11)                                                                                                                                  gaggaAGT‑AAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATAATTGAACAACAT‑AAAAAATCATTTTTTGTATCAGGTCTCCCATaaaa         >  5:1518387/4‑150 (MQ=11)                                                                                                                                    tgaAGTAAAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTATTGAACACCCT‑AACGGGTGTTTTTTTGTTTCTGGTCTCCCATaaaaa        >  3:40533/1‑150 (MQ=255)                                                                                                                                     gaAGTAAAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGACCCCGCAACCAATTATTGAACACCCT‑AACGGGTGTTTTTTTGTTTCTGGTCTCCCATaaaaaac      >  1:1582905/1‑150 (MQ=255)                                                                                                                                     gaAGT‑AAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTATTGAACACCCT‑AACGGGTGTTTTTTTGTTTCTGGTCTCCCATAAAAACgc     >  5:1577434/1‑151 (MQ=255)                                                                                                                                     gaAGT‑AAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATGATTGAACACCCT‑AAACCCTGATTTTTTGTAACAGGTCTCCCATAAAAAAgc     >  6:520427/1‑151 (MQ=14)                                                                                                                                      aaGTAAAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTATTGAACACCCT‑AACGGGTGTTTTTTTGTTTCTGCTCTCCCATaaaaaacc     >  6:463529/1‑149 (MQ=255)                                                                                                                                      aaGTAAAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTATTGAACAACAT‑AACGGGTGTTTTTTTGTTTCTGGTCTCCCATAAAAAAgc     >  4:1471076/1‑151 (MQ=255)                                                                                                                                          aaaaGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTATTGAACACCCT‑AACGGGTGTTTTTTTGTTTCTGGTCTCCCATAAAAAAGCGcc  >  1:360852/1‑150 (MQ=255)                                                                                                                                          |                                                                                                                                                       GCAACCAATTAAAATT‑TGATGAAGTAAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAA‑TC‑AAATTTGATGAAGTAAAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTATTGAACACCCT‑AACGGGTGTTTTTTTGTTTCTGGTCTCCCATAAAAAAGCGCC  >  NC_000913/2947457‑2947740

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