Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,947,591	0	A	.	38.6%	175.1 / 122.7	83	intergenic (+18/‑16)	metW/metV	tRNA‑Met/tRNA‑Met
Reads supporting (aligned to +/- strand):  ref base A (49/2);  new base . (30/2);  total (79/4)
Fisher's exact test for biased strand distribution p-value = 6.37e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.33e-01

TCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAA‑TC‑AAATTTGATGAAGTAAAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTATTGAACACCCTAACG‑GGTGTTTTTTTGTTTCTGGTCTCCCATAAAAAAGCGCC  >  NC_000913/2947530‑2947740                                                                |                                                                                                                                                       tCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAA‑TC‑AAATTTGAACAAGTAAAAGCAGTAAGGTGACGCGGTGTCGAGCAGCCTCGAAGCTCGTCGGGCTCATAACCCGAAGCACGACGGTTCAAATCCGGCCACCCCa                                                                >  5:664046/1‑150 (MQ=14)                       gcccctCCCGCCCCCGCAACCAA‑TC‑AAAATTGCAGCACTAAAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTATTGAACACCCTAACg                                         <  5:1182822/146‑1 (MQ=17)                                 cccccGCAACCAA‑TC‑AAATTTGATGAAGTAAAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTATTGAACACCCTAACG‑GGTGtttttt                              >  7:927739/1‑151 (MQ=255)                                  ccccGCAACCAA‑TC‑AAATTTGATGAAGTAAAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTATTGAACACCCTAACG‑GGTGttttttt                             >  4:1076998/1‑151 (MQ=255)                                   cccGCAACCAA‑TC‑AAATTTGATGAAGTAAAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTATTGAACACCCTAACG‑GGTGTTTTTTTg                            >  7:1563353/1‑151 (MQ=255)                                    ccGCAACCAA‑TC‑AAATTTGATGAAGTAAAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTATTGAACACCCTAACG‑GGTGTTTTTTTGt                           >  3:1656119/1‑151 (MQ=255)                                     cGCAACCAA‑TC‑AAATTTGATGAAGTAAAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTGGTCGGTTCAAATCCGGCCCCCGCAACCAATTATTGAACACCCTAACG‑GGTGTTTTTTTGtt                          >  5:453347/1‑151 (MQ=255)                                      gCAACCAA‑TC‑AAATTTGATGAAGTAAAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTATTGAACACCCTAACG‑GGTGTTTTTTTGttt                         >  1:1378596/1‑151 (MQ=255)                                       cAACCAATTA‑AAATTTGATGAAGT‑AAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATGATAGAACAGCCTAACT‑CGTGTGTTGTGGttgg                        >  8:238935/1‑149 (MQ=11)                                       aaaCCAA‑TTAAAATTTCATCAAG‑AAAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTATTGAACACCCTAACG‑GGTGTTTTTTTGTTTc                        <  3:1296026/150‑1 (MQ=14)                                        aaCCAA‑TTAAAATTTGATGAAGT‑AAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCACCCAATTATTGAACACCCAGACG‑GGTGTTTTGTttttacg                       >  7:726706/1‑148 (MQ=11)                                        aaCCAA‑TC‑AAATTTGATGAAGTAAAAGCAGTGCGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTATTGAACACCCTAACG‑GGTGTTTTTTTGTTTCt                       >  4:871700/1‑151 (MQ=255)                                        aaCCAA‑TC‑AAATTTGATGAAGTAAAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTATTGAACACCCTAACG‑GGTGTTTTTTTGTTTCt                       >  5:355872/1‑151 (MQ=255)                                         aCCAA‑TC‑AAATTTGATGAAGTAAAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTATTGAACACCCTAACG‑GGTGTTTTTTTGTTTCTg                      >  5:601366/1‑151 (MQ=255)                                          ccAA‑TC‑AAATTTGATGAAGTAAAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTATTGAACACCCTAACG‑GGTGTTTTTTTGTTTCTgt                     >  6:819721/1‑150 (MQ=255)                                          ccAA‑TC‑AAATTTGATGAAGTAAAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTATTGAACACCCTAACG‑GGTGTTTTTTTGTTTCTgg                     >  6:442726/1‑151 (MQ=255)                                          ccAA‑TC‑AAATTTGATGAAGTAAAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTATTGAACACCCTAACG‑GGTGTTTTTTTGTTTCTgg                     >  7:100697/1‑151 (MQ=255)                                           ccaattaaaaTTTGATGAAGT‑AAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTATTTAACAACCTAACG‑GGTGTTTTTTTGTTTc                        >  2:1303878/8‑148 (MQ=14)                                           ccaattaaaaTTTGATGAAGT‑AAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTATTGAACACCCTAACG‑GGTGTTTTTTTGTTTCTgg                     >  8:1583148/8‑151 (MQ=255)                                           cAA‑TC‑AAATTTGATGAAGTAAAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTATTGAACACCCTAACG‑GGTGTTTTTTTGTTTCTggg                    >  1:165831/1‑150 (MQ=255)                                           cAA‑TC‑AAATTTGATGAAGTAAAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTATTGAACACCCTAACG‑GGTGTTTTTTTGTTTCTGGt                    >  2:717964/1‑151 (MQ=255)                                            catcctataTTTTATTAAGT‑AAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTATTGAACACCCTAACG‑GGTGTTTTTTTGTTTCTGGt                    <  1:1220229/143‑1 (MQ=17)                                            aa‑TC‑AAATTTGATGAAGTAAAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATACCCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCACCCAATTATTGAACACCCTAACG‑GGTGTTTTTTTGTTTCTGGtc                   >  1:1138470/1‑151 (MQ=255)                                            aa‑TC‑AAATTTGATGAAGTAAAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTATTGAACACCCTAACG‑GGTGTTTTTTTGTTTCTGGtc                   >  5:1410243/1‑151 (MQ=255)                                             aattaaaaTTTGATGAAGT‑AAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAATGAACAACCAAAAG‑AGTCAGTTTTTGTTTCTgggg                   >  4:445921/6‑149 (MQ=11)                                             a‑TC‑AAATTTGATGAAGTAAAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACACGAAGGTCGTCGGTTCAAAACCGGCCCCCGCAACCAATTATTGAACACCCTAACC‑GGAGTTATTTTGTTTCTGGtct                  >  6:299638/1‑151 (MQ=255)                                             a‑TC‑AAATTTCATGAAGTAAAAGCAGTACGCTGACGCCGCGTGCAGAAGCCAGGTAGCTCGACGGCCTCATAACCCGAACGTCGACGCATCAAATCCGGCCCCCGCAACCAAATATTGAACACCCTAACG‑GGAGTATTTTTGTAACTGGta                   >  5:1155112/1‑149 (MQ=17)                                              attaaaaTTTGATGAAGT‑AAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTATTTAACAACCTAAAA‑GGTGTTTTTTTGTTGCAGGtct                  >  4:137046/5‑151 (MQ=14)                                              attaaaaTTTGATGAAGT‑AAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATAATTTAACACCCTAACG‑GGGGTTTTATTTTTTCAGGtct                  >  2:1250624/5‑151 (MQ=14)                                               ttaaaaTTTGATGAAGT‑AAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTATTTAACACCCTCACC‑GGGGTGTTTGTGTTGCTGGtctg                 >  8:105145/4‑150 (MQ=11)                                               ttaaaaTTTGATGAAGT‑AAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTATTGAACAACCTAAAG‑GGTCTTTTTTTGTTTCTGGtctc                 >  4:1326068/4‑151 (MQ=14)                                               tC‑AAATTTGATGAAGTAAAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTATTGAACACCCTAACG‑GGTGTTTTTTTGTTTCTGGtctc                 >  1:1302829/1‑151 (MQ=255)                                                taaaaTTTGATGAAGT‑AAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATAATTGAACACCCTAACA‑GCTGTTTTTTTTTTACGGGtctgc                >  2:1073675/3‑149 (MQ=14)                                                c‑AAATTTGATGAAGTAAAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTATTGAACACCCTAACG‑GGTGTTTTTTTGTTTCTGGTCTcc                >  7:1635972/1‑151 (MQ=255)                                                c‑AAATTTGATGAAGTAAAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTATTGAACACCCTAACG‑GGTGTTTTTTTGTTTCTGGTCTcc                >  4:341529/1‑151 (MQ=255)                                                c‑AAATTTGATGAAGTAAAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTATTGAACACCCTAAAG‑GGTGTTTTTTTGTTTCTGGTCTcc                >  1:944840/1‑151 (MQ=255)                                                 aaaaTTTGATGAAGT‑AAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCACCCAATGATTGAACACCCAAACG‑CGGCTTGTAGTGTTCCCGGTCTccc               >  8:260174/2‑151 (MQ=11)                                                 aaaaTTTGATGAAGT‑AAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTATTGAACACCCTAACG‑GGTGTTTTTTTGTTTGTGGTCTccc               >  8:1404399/2‑151 (MQ=255)                                                 aaaaTTTGATGAAGT‑AAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTATTGAACACCCTAACG‑GGTGTTTTTTTGTTTCTGGTCTccc               >  1:813407/2‑151 (MQ=255)                                                  aaaTTTGATGAAGTAAAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTATTGAACACCCTAACG‑GGTGTTTTTTTGTTTCTGGTCTccc               >  5:576296/1‑151 (MQ=255)                                                  aaaTTTGATGAAGTAAAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTATTGAACACCCTAACG‑GGTGTTTTTTTGTTTCTGGTCTccc               >  8:968852/1‑151 (MQ=255)                                                  aaaTTTGATGAAGTAAAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTATTGAACACCCTAACG‑GGTGTTTTTTTGTTTCTGGTCTccc               >  5:691027/1‑151 (MQ=255)                                                  aaaTTTGATGAAGTAAAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTATTGAACACCCTAACG‑GGTGTTTTTTTGTTTCTGGTCTccc               >  8:1065313/1‑151 (MQ=255)                                                  aaaTTTGATGAAGTAAAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTATTGAACACCCTAACG‑GGTGTTTTTTTGTTTCTGGTCTccc               >  6:892391/1‑151 (MQ=255)                                                  aaaTTTGATGAAGTAAAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTATTGAACACCCTAACG‑GGTGTTTTTTTGTTTCTGGTCTccc               >  3:323556/1‑151 (MQ=255)                                                  aaaTTTGATGAAGTAAAAGAAGTACGCTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCCAAGGTCGTCGTTTCAACTCCCCCCCACGCAACCAATTATTGAACACACAAACC‑GGTGTTTTTTTGTATCTGGTCTccc               >  1:1669893/1‑151 (MQ=18)                                                  aaaTTTGATGAAGT‑AAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATGATTGAACACCCTAAAC‑GCTGATTTTTTGGTTCTGGTCTccc               >  8:229034/1‑150 (MQ=14)                                                  aaaTTTGATGAAGT‑AAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATGATATAACACCCAGACG‑GGGGTTTTTGTGTTTGGGGGCTccc               >  5:63341/1‑150 (MQ=11)                                                  aaaTTCGATGAAGT‑AAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTATTGAACAACATAACG‑GGTGATTTGTTGTTGCTGTTGTccc               >  2:1402532/1‑150 (MQ=14)                                                   acgttcttgaAGT‑AAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTATTGAACACCCTAACG‑GGTGTTTTTTTGTTTCTGGTCTCCCAt             <  1:1079227/144‑1 (MQ=255)                                                   aaTTTGATGAAGTAAAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTATTGAACACCCTAACG‑GGTGTTTTTTTGTTTCTGGTCTccc               >  6:1577355/1‑150 (MQ=255)                                                   aaTTTGATGAAGTAAAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTATTGAACACCCTAACG‑GGTGTTTTTTTGTTTCTGGTCTccc               >  1:220476/1‑150 (MQ=255)                                                    aTTTGATGAAGT‑AAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAATGAACACCCTAACA‑CGGGTGTTTTTGTTTCTGGTCTCCCAt             >  7:1274954/1‑150 (MQ=17)                                                    aTTTGATGAAGT‑AAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATAATTGAACACCCTAACG‑GGTGTTTTTTTGTTTCTGGTCTCCCAt             >  4:1541817/1‑150 (MQ=255)                                                    aTTTGATGAAGT‑AAAGCAGTACGGAGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTATTGAACACCCTAACG‑CGTGTTTTTTTGTTTCTGGTCTCCCAt             >  8:1036397/1‑150 (MQ=255)                                                     tttGATGAAGTAAAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTATTGAACACCCTAACG‑GGTGTTTTTTTGTTTCTGGTCTCCCAt             >  5:1577472/1‑150 (MQ=255)                                                     tttGATGAAGTAAAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTATTGAACACCCTAACG‑GGTGTTTTTTTGTTTCTGGTCTCCCAt             >  2:569760/1‑150 (MQ=255)                                                     tttGATGAAGTAAAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTATTGAACACCCTAACG‑GGTGTTTTTTTGTTTCTGGTCTCCCAt             >  3:1228284/1‑150 (MQ=255)                                                     tttGATGAAGT‑AAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTATTGATCAACCTCACG‑GGTGTTTTTTTGTTTCTGGTCTCCCATa            >  5:1295750/1‑150 (MQ=255)                                                     tttGATGAAGT‑AAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAATTAACACCATAAAC‑GGTGTTTTTTTGTTTCTGGTCTCCCATa            >  2:226285/1‑150 (MQ=17)                                                      ttGATGAAGTAAAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTATTGAACACCCTAACG‑GGGGTTTTTTTGTTTCTGGTCTCCCATa            >  3:499086/1‑150 (MQ=255)                                                      ttGATGAAGT‑AAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTATAGAACACCCTAACC‑GGGGATTTGGTGTTCCTGGTCTCCCTTaa           >  7:245133/1‑150 (MQ=14)                                                      ttGATGAAGT‑AAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTATAGAACAACCTAACGCGGCATTAGTTTG‑AACTGGTCTCCCATaa           >  2:849147/1‑150 (MQ=11)                                                       tgatgaAGTAAAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTATTGAACACCCTAACG‑GGTGTTTTTTTGTTTCTGGTCTCCCATaa           >  6:757902/1‑150 (MQ=255)                                                       tgatgaAGTAAAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTATTGAACACCCTAACG‑GGTGTTTTTGTGTTTCTGGTCTCCCATaa           >  4:192286/1‑150 (MQ=255)                                                       tgatgaAGT‑AAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTATTGAACAACATGACT‑GGTGTTTTTTTGTTTCTGGTCTCCAATaaa          >  3:1458623/1‑150 (MQ=17)                                                         atgaAGT‑AAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTATTGAACACCCTCACG‑GGTGTGTTTTTGTTTGTGGGCTCCCATaaaaac       >  7:499068/1‑150 (MQ=17)                                                         atgaAGT‑AAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTATTGAACAACATAAAG‑GGTGTTTTTTTGTTTCTGGTCTCCCAt             >  5:1272464/1‑145 (MQ=255)                                                          tgaAGTAAAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTATTGAACACCCTAATG‑GGTGTTTTTTTGTTTCTGGTCTCCCATaaaaa        >  1:1261748/1‑150 (MQ=255)                                                          tgaAGTAAAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTATTGAACACCCTAACG‑GGTGTTTTTTTGTTTCTGGTCTCCCATaaaaa        >  7:1377683/1‑150 (MQ=255)                                                          tgaAGTAAAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTATTGAACACCCTAACG‑GGTGTTTTTTTGTTTCTGGTCTCCCATaaaaa        >  3:52718/1‑150 (MQ=255)                                                          tgaAGTAAAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTATTGAACACCCTAACG‑GGTGTTTTTTTGTTTCTGGTCTCCCATaaaaa        >  8:207572/1‑150 (MQ=255)                                                          tgaAGT‑AAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTATTGAACACCCTCACG‑GCTGTTTTTTGGTTTGTGGGCTCCCCTaacaacg      >  8:982138/1‑149 (MQ=14)                                                          tgaAGT‑AAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTATTGAACACCCTAACG‑GGTGTGTTTTTGTTTCTGGTCTCCCagagagcgg      >  5:1076688/1‑143 (MQ=255)                                                          tgaAGT‑AAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTATTGAACAACCTAACC‑GGTGTTTTTTTGTTTCTGGTCTCCCATaaaaacg      >  7:1518658/1‑149 (MQ=255)                                                          tgaAGT‑AAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATAATTGAACACCCAAACG‑GGGGATTTGTTGTTGCGGGTCTccaaaaacaa        >  8:656803/1‑141 (MQ=14)                                                           gaAGTAAAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTATTGAACACCCTAACG‑GGTGTTTTTTTGTTTCTGGTCTCCCATaaaaa        >  7:191277/1‑149 (MQ=255)                                                            aaGTAAAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTATTGAACACCCTAACG‑GGTGTTTTTTTGTTTCTGGTCTCCCATAAAAAAgc     >  8:1161337/1‑151 (MQ=255)                                                              gTAAAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTATTGAACACCCTAACG‑GGTGTTTTTTTGTTTCTGGTCTCCCATAAAAAAgcgc   >  6:634940/1‑151 (MQ=255)                                                              gTAAAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTATTGAACACCCTAACG‑GGTGTTTTTTTGTTTCTGGTCTCCCATAAAAAAgcgc   >  5:795929/1‑151 (MQ=255)                                                               tAAAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTATTGAACACCCTAACG‑GGTGTTTTTTTGTTTCTGGTCTCCCATAAAAAAGCGcc  >  7:498442/1‑151 (MQ=255)                                                               tAAAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTATTGAACACCCTAACG‑GGTGTTTTTTTGTTTCTGGTCTCCCATAAAAAAGCGcc  >  5:496380/1‑151 (MQ=255)                                                               tAAAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCCCAACCAATTATTGAACACCCTAACG‑GGTGTTTTTTTGTTTCTTGTCTCCCATAAAAAAgcgc   >  5:1567681/1‑150 (MQ=255)                                                                aaaaGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCCACCAATTATTGAACACCCTAACG‑GGTGTTTTTTTGTTTCTGGTCTCCCATAAAAAAGCGcc  >  4:1256221/1‑150 (MQ=255)                                                                aaaaGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTATTGAACACCCTAACG‑GGTGTTTTTTTGTTTCTGGTCTCCCATAAAAAAGCGcc  >  7:867258/1‑150 (MQ=255)                                                                aaaaGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTATTGAACACCCTAACG‑GGTGTTTTTTTGTTTCTGGTCTCCCATAAAAAAGCGcc  >  1:903075/1‑150 (MQ=255)                                                                aaaaGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTATTGAACACCCTAACG‑GGTGTTTTTTTGTTTCTGGTCTCCCATAAAAAAGAGcc  >  2:882314/1‑150 (MQ=255)                                                                aaaaGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTATTGAACACCATAACG‑GGTGTTTTTTTGTTTCTGGTCTCCCATAAAAAAGCGcc  >  1:275658/1‑150 (MQ=255)                                                                |                                                                                                                                                       TCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAA‑TC‑AAATTTGATGAAGTAAAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTATTGAACACCCTAACG‑GGTGTTTTTTTGTTTCTGGTCTCCCATAAAAAAGCGCC  >  NC_000913/2947530‑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: ‑