Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,947,591	0	A	.	36.5%	168.6 / 103.8	74	intergenic (+18/‑16)	metW/metV	tRNA‑Met/tRNA‑Met
Reads supporting (aligned to +/- strand):  ref base A (45/2);  new base . (27/0);  total (72/2)
Fisher's exact test for biased strand distribution p-value = 5.30e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.40e-01

AAATCCGGCCCCCGCAACCAATTAAAATTTGATGAAGTAAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATCAAATTTGATGAAGTAAAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTATTGAACACCCT‑AACGGGTGTTTTTTTGTTTCTGGTCTCCCATAAAAAAGCGCC  >  NC_000913/2947444‑2947740                                                                                                                                                    |                                                                                                                                                       aaaTACTGCCGCCGGAATCAATTAAAATTTGATGAATTAAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATCAAATTTGATGAAGTaaaa                                                                                                                                                     <  7:686573/151‑1 (MQ=17)                aaCGAATTATATTTTGTTGAAGTTTCGCATTACGGGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATCAAATTTGATGAAGTAAAAGCAGTACGGTGAcgc                                                                                                                                      <  1:283879/151‑1 (MQ=14)                   cAATTAAAATTTGATGAAGTAAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGCTCGTCGCTTCAAATCCGGCCCCCGCAACCAATCAAATAAGAAGAACAAAAAGCAGTAATGTGACGAggg                                                                                                                                   >  7:390548/1‑151 (MQ=17)                                                                                      cTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATCAAATTTGATGAAGTAAAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGc                                                                 >  3:1269821/1‑150 (MQ=40)                                                                                                      cgtcgGTTCAAATCCGGCCCCCGCAACCAATCATATTTGATGATGTAAAAGCAGTACTGTGACGGGGGGTGGAGAAGCCTGGTAGCTCTTCGGGCGCATTACCCGAAGGTCTTCGGGTCCCATCGGGCCCGCGCCACCATTTATTGAacac                                                >  3:358061/1‑151 (MQ=14)                                                                                                                   ccGGCCCCCGCAACCAATCAAATAACAAGAAGTAAAATCAGTACGGTGACCCGGGGTGCAGAAGAATAGAAGCTATTCGGGCACATAAAACGAACGGCGGGCGGACAAATCCCCCCCCCGCAACAAAAAATTGAACACCCA‑AACgctagtt                                  >  2:1477611/1‑145 (MQ=255)                                                                                                                      gaccccGCAACCAATCAAATTTGATGAAGTAAAATCAGTACGGTGACGCGGGGTGGAGCAGACTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTATTGAACACCCT‑AACGGGTGttttt                               >  5:1174318/3‑151 (MQ=255)                                                                                                                      gCCCCCGCAACCAATCAAATTTGATGAAGTAAAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTATTGAACACCCT‑AACGGGTGttttt                               >  7:79690/1‑151 (MQ=255)                                                                                                                        ccccGCAACCAATCAAATTTGATGAAGTAAAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCACCCGACAACAATTATTGAACACCCT‑AACCGGTGttttttt                             >  6:854445/1‑151 (MQ=255)                                                                                                                         cccGCAACCAATCAAATTTGATGAAGTAAAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTATTGAACACCCT‑AACGGGTGTTTTTTTg                            >  8:554674/1‑151 (MQ=255)                                                                                                                          ccGCAACCAATCAAATTTGATGAAGTAAAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATTACCCGAAGGTCGTGGGTTCAAATCCGGCGCCCGCAACCAATTTTTGAACACCCG‑AACGGGTGTTTTTTTGt                           >  5:854412/1‑151 (MQ=255)                                                                                                                          ccGCAACCAATCAAATTTGATGAAGTAAAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGACGTCGGTTCAAATCCGGCCCCCGCAACCAATTATTGAACACCCT‑AACGGGTGTTTTTTTGt                           >  4:657545/1‑151 (MQ=255)                                                                                                                            gCAACCAATCAAATTTGATGAAGTAAAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTATTGAACACCCT‑AACGGGTGTTTTTTTGttt                         >  6:1079487/1‑151 (MQ=255)                                                                                                                            gCAACCAATCAAATTTGATGAAGTAAAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTATTGAACACCCT‑AACGGGTGTTTTTTTGttt                         >  1:1348175/1‑151 (MQ=255)                                                                                                                            gCAACCAATCAAATTTGATGAAGTAAAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTATTGAACACCCT‑AACGGGTGTTTTTTTGttt                         >  7:1159534/1‑151 (MQ=255)                                                                                                                             cAACCAATCAAATTTGATGAAGTAAAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTATTGAACACCCT‑AACGGGTGTTTTTTTGttt                         >  1:1291953/1‑150 (MQ=255)                                                                                                                             cAACCAATCAAATTTGATGAAGTAAAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTATTGAACACCCT‑AACGGGTGTTTTTTTGTTTc                        >  2:1404929/1‑151 (MQ=255)                                                                                                                             cAACCAATCAAATTTGATGAAGTAAAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTATTGAACACCCT‑AACGGGTGTTTTTTTGTTTc                        >  6:549839/1‑151 (MQ=255)                                                                                                                              aaCCAATCAAATTTGATGAAGTAAAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTATTGAACACCCT‑AACGGGTGTTTTTTTGTTTCt                       >  5:537306/1‑151 (MQ=255)                                                                                                                              aaCCAATCAAATTTGATGAAGTAAAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTATTGAACACCCT‑AACGGGTGTTTTTTTGTTTCt                       >  6:664045/1‑151 (MQ=255)                                                                                                                              aaCCAATCAAATTTGATGAAGTAAAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTATTGAACACCCT‑AACGGGTGTTTTTTTGTTTCt                       >  7:1110246/1‑151 (MQ=255)                                                                                                                               ccaattaaaaTTTGATGAAGT‑AAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTATTGATCACCCT‑AACAGGTGTTTTATTGTATCTgg                     >  6:27841/8‑151 (MQ=14)                                                                                                                               ccaattaaaaTTTGATGAAGT‑AAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTATTGAACAACCT‑AAAAGGTGTTTTTTTGTTTCAgg                     >  1:97556/8‑151 (MQ=14)                                                                                                                               ccaattaaaaTTTGATGAAGT‑AAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATGAATGAACACCCT‑AACCGGTGTGTTTTTGTTTGTgg                     >  2:1319532/8‑151 (MQ=11)                                                                                                                               aCCAATCAAATTTGATGAAGTAAAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATACCCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCCACCAATTATTGAACACCCT‑AACGGGTGGTTTTTTGTTTCTg                      >  1:1471518/1‑151 (MQ=255)                                                                                                                               aCCAATCAAATTTGATGAAGTAAAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTATTGAACACCCT‑AACGGGTGTTTTTTTGTTTCTg                      >  3:1184921/1‑151 (MQ=255)                                                                                                                               aCCAATCAAATTTGATGAAGTAAAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTATTGAACACCCT‑AACGGGTGTTTTTTTGTTTCTg                      >  7:589799/1‑151 (MQ=255)                                                                                                                               aCCAATCAAATTTGATGAAGTAAAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTATTGAACACCCT‑AACGGGTGTTTTTTTGTTTCTg                      >  3:1742445/1‑151 (MQ=255)                                                                                                                                caattaaaaTTTGATGAAGT‑AAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTATTGAACACCCTGAACGGGGTTttgtttgtgctggg                     >  3:197872/7‑145 (MQ=11)                                                                                                                                 cAATCAAATTTGATGAAGTAAAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAGCCAATTATTGAACACCCT‑AACGGGTGTTTTTTTGTTTCTGGt                    >  2:765792/1‑151 (MQ=255)                                                                                                                                 cAATCAAATTTGATGAAGTAAAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTATTGAACACCCT‑AACGGGTGTTTTTTTGTTTCTGGt                    >  5:1643910/1‑151 (MQ=255)                                                                                                                                 cAATCAAATTTGATGAAGTAAAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTATTGAACACCCT‑AACGGGTGTTTTTTTGTTTCTGGt                    >  4:534202/1‑151 (MQ=255)                                                                                                                                 aattaaaaTTTGATGAAGT‑AAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTATTTAACACCCT‑AACGCGTCTTTTTTTTTTTCTgggg                   >  1:933125/6‑149 (MQ=14)                                                                                                                                 aattaaaaTTTGATGAAGT‑AAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTATTTAACAACAT‑AACGTGGGATTTTTTTTTTCAGGtg                   >  6:1268136/6‑150 (MQ=11)                                                                                                                                 aattaaaaTTTGATGAAGT‑AAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTATTGAACACCCT‑CACGGGGGTTTTgttgttggtgggg                   >  8:583508/6‑144 (MQ=14)                                                                                                                                  aaTCAAATTTGATGAAGTAAAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTATTGAACACCCT‑AACGGGTGTTTTTTTGTTTCTGGtc                   >  4:1324447/1‑151 (MQ=255)                                                                                                                                   aTCAAATTTGATGAAGTAAAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTATTGAACACCCT‑AACGGGTGTTTTTTTGTTTCTGGtct                  >  8:1269167/1‑151 (MQ=255)                                                                                                                                    taaaaTTTGATGAAGT‑AAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTATTGAACACCCT‑AACGCGTGTTTTATTGTTTCTGGTCTcc                >  1:1213490/3‑151 (MQ=255)                                                                                                                                    taaaaTTTGATGAAGT‑AAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATAATAGAACACCCT‑CAAGCGGGTTTTTTTTTATCGGGTGTcc                >  1:1501893/3‑151 (MQ=11)                                                                                                                                    taaaaTTTGATGAAGT‑AAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATAAATTAACACCCA‑AACGCGTCATTATTTGGTTCTGGTGTcc                >  5:141332/3‑151 (MQ=11)                                                                                                                                    taaaaTTTCATCAAG‑AAAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCGCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTATTGAACACCCT‑AACGGGTGTTTTTTTGTTTCTGGtctc                 <  8:296837/148‑1 (MQ=18)                                                                                                                                     aaaaTTTGATGAAGT‑AAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTATTGAACACCCT‑AACGGGTGTTTTTTTGTTTGTGGTCTccc               >  8:1585516/2‑151 (MQ=255)                                                                                                                                      aaaTTTGATGAAGTAAAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTATTGAACACCCT‑AACGGGTGTTTTTTTGTTTCTGGTCTccc               >  4:645475/1‑151 (MQ=255)                                                                                                                                      aaaTTTGATGAAGTAAAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTATTGAACACCCT‑AACGGGAGTTTATTTGTTTCTGGTCTccc               >  6:761204/1‑151 (MQ=255)                                                                                                                                      aaaTTTGATGAAGT‑AAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGTGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTATAGAACAACAT‑ATATTATCTTTTATTGTATATGGTCTccc               >  5:1026238/1‑150 (MQ=11)                                                                                                                                      aaaTTTGATGAAGT‑AAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATGATTGAACACCCT‑AACGGCTGTTTTTTTGTTTCTGGTCTcca               >  1:330621/1‑149 (MQ=255)                                                                                                                                       aaTTTGATGAAGT‑AAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATAAATTATCACCCT‑ATAAGGTGTTTTTTGTTTTCTGGTCTCCCa              >  5:379655/1‑150 (MQ=11)                                                                                                                                        aTTTGATGAAGT‑AAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTATTGAACACCCT‑AACGGGTGTTTTTTTGTTACTGGTCTCCCAt             >  4:1729728/1‑150 (MQ=255)                                                                                                                                        aTTTGATGAAGT‑AAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTATTGAACACCCT‑AACGGGGGTTTTTTTGTTTCTGGTCTcccct             >  8:1093727/1‑148 (MQ=255)                                                                                                                                         tttGATGAAGTAAAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTATTGAACACCCT‑AACGGGTGTTTTTTTGTTTCTGGTCTCCCAt             >  5:1397969/1‑150 (MQ=255)                                                                                                                                         tttGATGAAGT‑AAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTATTTAACACCCT‑AACGGGTGTTTTTTTTTTTCTGGTCTCCCATa            >  5:56413/1‑150 (MQ=255)                                                                                                                                         tttGATGAAGT‑AAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTATTGAACACCCT‑AACGGGTGTTTTTTTGTTTCTGGGCTCCCATa            >  3:1050306/1‑150 (MQ=255)                                                                                                                                         tttGATGAAGT‑AAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATAAATGAACACCCA‑AAAACGTGATTTTTTGTTTCTGGTCTCCCATa            >  7:759352/1‑150 (MQ=14)                                                                                                                                           tgatgaAGTAAAAGCAGTACGGTGCCGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTATTGAACACCCT‑AACGGGTGTTTTTTTGTTTCTGGTCTCCCATaa           >  5:891181/1‑150 (MQ=255)                                                                                                                                           tgatgaAGTAAAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTATTGAACACCCT‑AACGGGTGTTTTTTTGTTTCTGGTCTCCCATaa           >  6:1703354/1‑150 (MQ=255)                                                                                                                                           tgatgaAGTAAAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTATTGAACACCCT‑AACGGGTGTTTTTTTGTTTCTGGTCTCCCATaa           >  7:126605/1‑150 (MQ=255)                                                                                                                                           tgatgaAGTAAAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTATTGAACACCCT‑AACGGGTGTTTTTTTGTTTCTGGTCTCCCATaa           >  5:1445524/1‑150 (MQ=255)                                                                                                                                           tgatgaAGT‑AAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTATTGAACACCCT‑AAAGGGTGTTTTTTTGTTTCTGGTCTCCCATaaa          >  2:163561/1‑150 (MQ=255)                                                                                                                                            gatgaAGT‑AAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATAAATGAACACCCT‑AAAAGCGCTTTTTTTTTTACTGGTCTCCCaaaaaa         >  3:1574944/1‑150 (MQ=14)                                                                                                                                             atgaAGT‑AAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTATTGAACACCCT‑CACGGGTGTGTTGTGGTTGCTGGGCTCCCATaaaaac       >  8:1118232/1‑150 (MQ=17)                                                                                                                                             atgaAGT‑AAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATAATTGAACACCCT‑AAAGCGTGTGTTGTTGGATCTGGTCTCCCATaaaaa        >  8:1316739/1‑150 (MQ=17)                                                                                                                                             atgaAGT‑AAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATAATTGAACAACAT‑AAAGCGTCTTTTATTTTAACAGGTCTCCaaaaaaaa        >  7:479381/1‑150 (MQ=11)                                                                                                                                              tgaAGT‑AAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTATTGAACACCCT‑CACGGGTGTGTTTTTGTTTCTGGTCTCCCATaaaaacg      >  3:1163165/1‑149 (MQ=255)                                                                                                                                               gaAGTAAAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTATTGAACACCCT‑AACGGGTGTTTTTTTGTTTCTGGTCTCCCATaaaaa        >  6:1026051/1‑149 (MQ=255)                                                                                                                                               gaAGT‑AAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTATTGAACACCCT‑AAAGGGTGATTTTTTGTTTCTGGTCTCCCATAAAAAAgc     >  1:1372657/1‑151 (MQ=255)                                                                                                                                               gaAGT‑AAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACACGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTATTGAACACCCT‑AACGGGTCttttttt                             >  5:1598688/1‑127 (MQ=255)                                                                                                                                                aaGTAAAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTATTGAACACCCT‑AACGGGTGTTTTTTTGTTTCTGGTCTCCCATAAAAAAgc     >  4:785591/1‑151 (MQ=255)                                                                                                                                                aaGTAAAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTATTGAACACCCT‑AACGGGTGTTTTTTTGTTTCTGGTCTCCCATAAAAAAgc     >  2:44364/1‑151 (MQ=255)                                                                                                                                                aaGTAAAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTATTGAACACCCT‑AACGCGTGTTTTTTTGTTTCTGGTCTCCCATAAAAAAgc     >  4:665307/1‑151 (MQ=255)                                                                                                                                                 aGTAAAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTATTGAACACCCT‑AACGGGTGTTTTTTTGTTTCTGGTCTCCCATAAAAAAgcg    >  2:301377/1‑151 (MQ=255)                                                                                                                                                 aGTAAAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTATTGAACACCCT‑AACGGGTGTTTTTTTGTTTCTGGTCTCCCATAAAAAAgcg    >  3:1705336/1‑151 (MQ=255)                                                                                                                                                 aGTAAAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTATTGAACACCCT‑AACGGGTGTTTTTTTGTTTCTGGTCTCCCATAAAAAAgcg    >  3:777629/1‑151 (MQ=255)                                                                                                                                                 aGTAAAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTATTGAACACCCT‑AACGGGTGTTTTTTTGTTTCTGGTCTCCCATAAAAAAgcg    >  4:1124730/1‑151 (MQ=255)                                                                                                                                                  gTAAAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTATTGAACACCCT‑AACGGGTGTTTTTTTGTTTCTGGTCTCCCATAAAAAAgcgc   >  4:592695/1‑151 (MQ=255)                                                                                                                                                  gTAAAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTACCTCGTCGGTCTCATAACACTAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTATTGAACACACT‑AACGGGTGTTATTTTGTTTCTGGTCTCCCATAAAAAAgacc   >  3:112080/1‑148 (MQ=21)                                                                                                                                                    aaaaGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTATTGAACACCCT‑AACGGGTGTTTTTTTGTTTCTGGTCTCCCATAAAAAAGCGcc  >  3:1235732/1‑150 (MQ=255)                                                                                                                                                    aaaaGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTATTGAACACCCT‑AACGGGTGTTTTTTTGTTTCTGGTCTCCCATAAAAAAGCGcc  >  3:487726/1‑150 (MQ=255)                                                                                                                                                    |                                                                                                                                                       AAATCCGGCCCCCGCAACCAATTAAAATTTGATGAAGTAAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATCAAATTTGATGAAGTAAAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTATTGAACACCCT‑AACGGGTGTTTTTTTGTTTCTGGTCTCCCATAAAAAAGCGCC  >  NC_000913/2947444‑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: ‑