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

CCCCCGCAACCAATTAAAATTTGATGAAGTAAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAA‑TC‑AAATTTGATGAAGTAAAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTATTGAACACCCTAACGGGTGTTTT‑TTTGTTTCTGGTCTCCCATAAAAAAGCGCC  >  NC_000913/2947452‑2947740                                                                                                                                              |                                                                                                                                                       ccgggtccaCCAAATAAATTTTGAAGAATTATAGTAATCCAGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAA‑TC‑AAATTTGATGAAGTAAAAGCAGTACg                                                                                                                                             <  3:633573/143‑1 (MQ=14)        aaCCAATTAAAATTTGATGAAGTAAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAA‑TA‑ATATTACATGAAGAAACCGCAGTatgtttaagg                                                                                                                                      >  3:1347815/1‑142 (MQ=14)                                           ggACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAA‑TC‑AAATTTGATGAAGTAAAAGCAGTACGGTGACGAgg                                                                                                                                    >  2:1245497/2‑118 (MQ=255)                                                                  tAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAA‑TG‑AATTTTGAACATGTAAGAGCGGTAGGGTGACGCGGGGTGGAGAAGCCTGGCGGCTCGTCGGGCGCTTCACGCGTAGgtggggcgggcgcct                                                                            >  5:460613/1‑138 (MQ=2)                                                                                             tcgtcgGTTCAAATCCGGCCCCCGCAACCAA‑TG‑AAATTTCAACCAGTAACAGCTGTAGGTTGACGCGGGCTGGCGAAGACTGGTAGCTCGTCGGGCTGATAACCCGAAGGGAGGCGGGGCCACTCCCGCCCCCGCAACAAATTATTGAac                                                  >  7:1238893/1‑150 (MQ=9)                                                                                                            cGGCCCCCGCAACCAA‑TC‑AAATTTGAAGAAGTAAAAGCAGTAAGGTGACGCGGGGTGCAGCACCCTCGTAGCAAGTCGGGCGCATAACCCGAACGGCGAGCGGTCGCATCCGGCCACCGCCACCAATGATTGAACACCCTGACCGGTGatt                                  >  2:675438/1‑151 (MQ=255)                                                                                                            cGGCCCCCGCAACCAA‑TA‑AAATTTGAAGAAGTAAAAGCAGTACGGTGACGAGGGGTGGCGAAGCCTGGTAGCTCGTCGCGCTCTTAACCCGAAGGTAGTACGTTCACATCCCGCCCACGCAACCAATTAATGAACACCCAAACGGCTGttt                                  >  6:824157/1‑151 (MQ=255)                                                                                                             ggCCCCCGCAACCAA‑TC‑AAATTTGATGAAGTAAAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTATTGAACACCCTAACGGGTGtgtt                                 >  6:723041/1‑151 (MQ=255)                                                                                                             ggCCCCCGCAACCAA‑TC‑AAATTTCATGAAGTAAAGGCAGTAGGGTGAACCGGGGTGGCGGAGACTGGTCGCTCGTCGGGCTGATAACCCGTACGGGGTGCGGTCAAGTCCGGCCCGCGCCAACAATGATTGAACACCCAACCGGctgatat                                 >  6:1274835/1‑147 (MQ=255)                                                                                                              gCCCCCGCAACCAA‑TC‑AAATTTGATGAAGTAAAAGCAGTACAGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTATTGAACACCCTAACGGGTGtttt‑t                               >  6:368852/1‑151 (MQ=255)                                                                                                                ccccGCAACCAA‑TC‑AAATTTGATGAAGTAAAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCACCCAATTATTGAACACCCTAACGGGTGtttt‑ttt                             >  5:997264/1‑151 (MQ=255)                                                                                                                ccccGCAACCAA‑TC‑AAATTTGATGAAGTAAAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTATTGAACACCCTAACGGGTGtttt‑ttt                             >  3:1664440/1‑151 (MQ=255)                                                                                                                ccccGCAACCAA‑TC‑AAATTTGATGAAGTAAAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTATTGAACACCCTAACGGGTGtttt‑ttt                             >  5:987699/1‑151 (MQ=255)                                                                                                                ccccGCAACCAA‑TC‑AAATTTGATGAAGTAAAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTATTGAACACCCTAACGGGTGtttt‑ttt                             >  5:135485/1‑151 (MQ=255)                                                                                                                 cccGCAACCAA‑TC‑AAATTTGATGAAGTAAAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTATTGAACACCCTAACGGGTGTTTT‑TTTg                            >  8:1763120/1‑151 (MQ=255)                                                                                                                  ccGCAACCAATTA‑AAATTTGATGAAGT‑AAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTGCAAATCCGGCCCCCGCAACCAATTAATTTACAACCACACGGGTGTGTT‑TTGGt                           >  7:1790816/1‑151 (MQ=11)                                                                                                                  ccGCAACCAA‑TC‑AAATTTGATGAAGTAAAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTATTGAACACCCTAACGGGTGTTTT‑TTTGt                           >  7:859768/1‑151 (MQ=255)                                                                                                                   cGCAACCAA‑TC‑AAATTTGATGAAGTAAAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTATTGAACACCCTAACGGGTGTTTT‑TTTGtt                          >  4:1514898/1‑151 (MQ=255)                                                                                                                    gCAACCAATTA‑AAATTTGATGAAGT‑AAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAAATATTTAACAACCAAACAGGGGTTTtgtttgtt                          >  8:673110/1‑151 (MQ=11)                                                                                                                     cAACCAATTA‑AAATTTGATGAAGT‑AAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAATTAACACCCTATCGGGGGtgtt‑ttgggtgg                        >  1:933200/1‑145 (MQ=11)                                                                                                                     cAACCAATTA‑AAATTTGAAGAAGT‑AAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTATTGAACACCCTGACGGGTGTGTT‑TTTGtttg                        >  1:1596341/1‑150 (MQ=14)                                                                                                                     aaaCCAA‑TAAAAATTTCATCATGT‑AAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTATTGAACACCCTAACGGGTGTTTT‑TTTGTTTc                        <  3:1474307/150‑1 (MQ=14)                                                                                                                      aaCCAA‑TC‑AAATTTGATGAAGTAAAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTATTGAACACCCTAACGGGTGTTTT‑TTTGTTTCt                       >  5:1530755/1‑151 (MQ=255)                                                                                                                      aaCCAA‑TC‑AAATTTGATGAAGTAAAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTATTGAACACCCTAACGGGTGTTTT‑TTTGTTTCt                       >  6:1618523/1‑151 (MQ=255)                                                                                                                      aaCCAA‑TC‑AAATTTGATGAAGTAAAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTATTGAACACCCTAACGGGTGTTTT‑TTTGTTTCt                       >  5:770538/1‑151 (MQ=255)                                                                                                                       aCCAA‑TC‑AAATTTGATGAAGTAAAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTATTGAACACCCTAACGGGTGTTTT‑TTTGTTTCTg                      >  5:1652213/1‑151 (MQ=255)                                                                                                                        ccAA‑TC‑AAATTTGATGAAGTAAAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTATTGAACACCCTAACGGGTGTTTT‑TTTGTTTCTgg                     >  1:1324156/1‑151 (MQ=255)                                                                                                                        ccAA‑TC‑AAATTTGATGAAGTAAAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATACGGCCCCCGCAACCAATTATTGAACACCCTAACGGGTGTTTT‑TTTGTTTCTgg                     >  2:1337478/1‑151 (MQ=255)                                                                                                                        accaattaaaaTTTGATGAAGT‑AAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATCATTTAAAACCCTAAAGGCTGTTGT‑GTGGTTACTg                      >  1:1611912/9‑151 (MQ=11)                                                                                                                        accaattaaaaTTTGATGAAGT‑AAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCAGGCCCCCGCAACCAATTATATAACAACCAAAAGACTCtttt‑tttttaacgg                      >  4:196218/9‑144 (MQ=11)                                                                                                                         cAA‑TC‑AAATTTGATGAAGTAAAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTATTGAACACCCTAACGGGTGTTTT‑TTTGTTTCTGGt                    >  5:184674/1‑151 (MQ=255)                                                                                                                         cAA‑TC‑AAATTTGATGAAGTAAAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTCGTAGCTCGTCGGGCTCATAACCCGAAGGTCCTCGGTTCAAATCCGGCCCCCGCAACCAATTATTGAACACCAAAACGCCTCTTTT‑TTTCATTCTgg                     >  2:221105/1‑150 (MQ=17)                                                                                                                          caattaaaaTTTGATGAAGT‑AAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTATTGAACACCATAAAAGGTGATTT‑GTTGGTTCgggg                    >  2:1592305/7‑150 (MQ=11)                                                                                                                          aa‑TC‑AAATTTGATGAAGTAAAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTATTGAACACCCTAACGGGTGTTTT‑TTTGTTTCTGGtc                   >  1:588691/1‑151 (MQ=255)                                                                                                                           gcttaaaaTTTGATGAAGTAAAAGCAGTCCGGTGCCGCGGGG‑GGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTATTGAACACCCTAACGGGTGTTTT‑TTTGTTTCTGGtc                   <  3:88739/146‑1 (MQ=255)                                                                                                                           aattaaaaTTTGATGAAGT‑AAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTATTTAACACCATCTAGGGTGTTTT‑tgtggtggggggg                   >  3:1332027/6‑142 (MQ=11)                                                                                                                           a‑TC‑AAATTTGATGAAGTAAAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTATTGAACACCCTAACGGGTGTTTT‑TTTGTTTCTGGtct                  >  7:55593/1‑151 (MQ=255)                                                                                                                           a‑TC‑AAATTTGATGAAGTAAAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTATTGAACACCCTAACGGGTGTTTT‑TTTGTTTCTGGtct                  >  5:1669877/1‑151 (MQ=255)                                                                                                                             ttaaaaTTTGATGAAGT‑AAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTATTGAACACCCTAACGGGTGTGTT‑TTTGGTGCTGGtctg                 >  5:1749854/4‑150 (MQ=14)                                                                                                                             ttaaaaTTTGATGAAGT‑AAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATAATTTAACAACCTTACGGGTGTTTT‑TTTGTTTCAGGtctc                 >  7:1399130/4‑151 (MQ=14)                                                                                                                             ttaaaaTTTGATGAAGT‑AAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATAATTGAACAACATAACCGGTGTTTT‑TTTGTTTCGGGtgtc                 >  5:166968/4‑151 (MQ=14)                                                                                                                             tC‑AAATTTGATGAAGTAAAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCACCCAATTATTGAACACCCTAACGGGTGTTTT‑TTTGTTTCTGGtctc                 >  1:1579532/1‑151 (MQ=255)                                                                                                                             tC‑AAATTTGATGAAGTAAAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTATTGAACACCCTAACGGGTGTTTT‑TTTGTTTCTGGtctc                 >  2:1563415/1‑151 (MQ=255)                                                                                                                              taaaaTTTGATGAAGT‑AAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTATTTAACACGATAACGGGTGTGTT‑TGTGTTTGTGGtgtcg                >  1:1679953/3‑150 (MQ=14)                                                                                                                              c‑AAATTTGATGAAGTAAAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTATTGAACACCCTAACGGGTGTTTT‑TTTGTTTCTGGTCTcc                >  4:864830/1‑151 (MQ=255)                                                                                                                              c‑AAATTTGATGAAGTAAAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTATTGAACACCCTAACGGGTGTTTT‑TTTGTTTCTGGTCTcc                >  1:146146/1‑151 (MQ=255)                                                                                                                               aaaaTTTGATGAAGT‑AAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTATTGTTCACCCTCACGGGTGTTTT‑TTTggatgtggggtccc               >  3:1729950/2‑138 (MQ=14)                                                                                                                                aaaTTTGATGAAGTAAAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTATTGAACACCCTAACGGGTGTTTT‑TTTGTTTCTGGTCTccc               >  4:552435/1‑151 (MQ=255)                                                                                                                                aaaTTTGATGAAGTAAAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTATTGAACACCCTAACGGGTGTTTT‑TTTGTTTCTGGTCTccc               >  7:1164599/1‑151 (MQ=255)                                                                                                                                aaaTTTGATGAAGTAAAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTATTGAACACCCTAACGGGTGTTTT‑TTTGTTTCTGGTCTccc               >  4:357983/1‑151 (MQ=255)                                                                                                                                aaaTTTGATGAAGTAAAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTATTGAACACCCTAACGGGTGTTTT‑TTTGTTTCTGGTCTccc               >  4:1447904/1‑151 (MQ=255)                                                                                                                                 aaTTTGATGAAGT‑AAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTATTTAACACCCAAAAAAATCATTA‑TTTGTTTCTGGTGTCCCa              >  4:902840/1‑150 (MQ=11)                                                                                                                                 aaTTTCATGACG‑AAAAGCAGTACGGTGACGCTGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTATTGAACACCCTAACGGGTGTTTT‑TTTGTTTCTGGTCTCCCAt             <  2:936508/151‑1 (MQ=255)                                                                                                                                 aaTTTCATGAAGTAAAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCATCGGTTCAAATCCGGCCCCCGCAACCAATTATTGAACACCCTAACGGGTGTTTT‑TTTGTTTCTGGTCTccc               >  6:611295/1‑150 (MQ=255)                                                                                                                                  aTTTGATGAAGTAAAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTATTGAACACCCTAACGGGTGTTTT‑TTTGTTTCTGGTCTCCCAt             >  1:955778/1‑151 (MQ=255)                                                                                                                                  aTTTGATGAAGTAAAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTATTGAACACCCTAACGGGTGTTTT‑TTTGTTTCTGGTCTCCCAt             >  3:109255/1‑151 (MQ=255)                                                                                                                                  aTTTGATGAAGT‑AAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTATTTAACAACAAAACGCGTCTGTT‑AGTGGTTCTGGTCTCCAAt             >  8:327457/1‑150 (MQ=11)                                                                                                                                   tttGATGAAGTGAAAGCAGTACGGTGACGCGGGGCGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTATTGAACACCCTAACGGGTGTTTT‑TTTGTTTCTGGTCTCCCAt             >  3:446917/1‑150 (MQ=255)                                                                                                                                   tttGATGAAGT‑AAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTATTGAACACGCTATAGGGGGTTTT‑TTTGTTGGTGGTCTCCCATa            >  1:1067713/1‑150 (MQ=14)                                                                                                                                   tttGATGAAGT‑AAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTATTGAACACCCTAACGGGGGTTTT‑GTGGTTGCTGGTGTCCCGTa            >  1:63646/1‑150 (MQ=14)                                                                                                                                    tttcatcacgaaaaGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTATTGAACACCCTAACGGGTGTTTT‑TTTGTTTCTGGTCTCCCATa            <  5:672418/140‑1 (MQ=255)                                                                                                                                    ttGATGAAGTAAAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTATTGAACACCCTAACGGGTGTTTT‑TTTGTTTCTGGTCTCCCATa            >  4:643500/1‑150 (MQ=255)                                                                                                                                    ttGATGAAGTAAAAGCAGTACGGTGACGCGGGGTGGAGCAGCATGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTATTGAACACCCTAACGGGTGTTTT‑TTTGTTTCTGGTCTCCCATa            >  6:520970/1‑150 (MQ=255)                                                                                                                                    ttGATGAAGT‑AAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTATATAACACCCAGACCGGGCTGTT‑ATTGGTGCTGGGCTCCCaaaa           >  6:1854747/1‑150 (MQ=11)                                                                                                                                    ttGATGAAGT‑AAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATAATTGAACACCCTAACGCGTGTTTT‑TTTGTTTCTGGTCTCCCATaa           >  4:881012/1‑150 (MQ=255)                                                                                                                                     tgatgaAGTAAAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTATTGAACACCCTAACGGGTGTTTT‑TTTGTTTCTGGTCTCCCATaa           >  1:1476598/1‑150 (MQ=255)                                                                                                                                     tgatgaAGTAAAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTATTGAACACCCTAACGGGTGTTTT‑TTTGTTTCTGGTCTCCCATaa           >  4:729909/1‑150 (MQ=255)                                                                                                                                     tgatgaAGT‑AAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATAAATGAACAACCTAAAACGTCTTTA‑TTTGTTTCTGGTCTCCCATaaa          >  6:439766/1‑150 (MQ=14)                                                                                                                                      gatgaAGTAAAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTATTGAACACCCTAACGGGTGTTTT‑TTTGTTTCTGGTCTCCCATaaa          >  1:75653/1‑150 (MQ=255)                                                                                                                                      gatgaAGT‑AAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTATTTAACAACAAAACGGCGGTTTT‑TTTGGTGCTGGTGTGCCATatca         >  7:532652/1‑147 (MQ=11)                                                                                                                                       atgaAGTAAAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTATTGAACACCCTAACGGGTGTTTT‑TTTGTTTCTGGTGTCCCATaaaa         >  8:1881347/1‑150 (MQ=255)                                                                                                                                       atgaAGTAAAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTATTGAACACCCTAACGGGTGTTTT‑TTTGTTTCTGGTCTCCCATaaaa         >  6:444063/1‑150 (MQ=255)                                                                                                                                       atgaAGTAAAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTATTGAACACCCTAACGGGTGTATT‑TATGTTTCTGGTCTCCCATaaaa         >  4:763328/1‑150 (MQ=255)                                                                                                                                        tgaAGTAAAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTATTGAACACCCTAACGGGTGTTTT‑TTTGTTTCTGGTCTCCCATaaaaa        >  7:900190/1‑150 (MQ=255)                                                                                                                                         gaAGTAAAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTATTGAACACCCTAACGGGTGTTTT‑TTTGTTTCTGGTCTCCCATaaaaaac      >  1:1338068/1‑150 (MQ=255)                                                                                                                                         gaAGTAAAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTATTGAACACCCTAACGGGTCTTTT‑TTTGTTTCTGGTCTCCCATaaaaa        >  5:670092/1‑149 (MQ=255)                                                                                                                                         gaAGT‑AAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTATTGAACACCCTAACTAGTCATTT‑TTTGTTTCTGGTCTCCAATAAAAAAgc     >  4:472963/1‑151 (MQ=17)                                                                                                                                         gaAGT‑AAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATAATTGAACACCCTAAAAGGTGTTTA‑TTTGTATCAGGTCTCCCATaaaaa        >  4:1208631/1‑148 (MQ=17)                                                                                                                                          aaGTAAAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTATTGAACACCCTAACGGGTGTTTT‑TTTGTTTCTGGTCTCCCATAAAAAAgc     >  4:1544555/1‑151 (MQ=255)                                                                                                                                          aaGTAAAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTATTGAACACCCTAACGGGTGTTTT‑TTTGTTTCTGGTCTCCCATAAAAAAgc     >  4:847303/1‑151 (MQ=255)                                                                                                                                           aGTAAAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTATTGAACACCCTAACGGGTGTTTT‑TTTGTTTCTGGTCTCCCATAAAAAAgcg    >  6:1351916/1‑151 (MQ=255)                                                                                                                                            gTAAAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTATTGAACACCCTAACGGGTGTTTT‑TTTGTTTCTGGTCTCCCATAAAAAAgcgc   >  8:462025/1‑151 (MQ=255)                                                                                                                                             tAAAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTATTGAACACCCTAACGGGTGTTTT‑TTTGTTTCTGGTCTCCCATAAAAAAGCGcc  >  8:1592798/1‑151 (MQ=255)                                                                                                                                             tAAAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCA‑AACCCGCAGGTGGGACGTTCAAATCCGGCCCCCGCAACCAATTATTGAACACCCTAAAGGGGGTGTT‑TTTTTTTCTGGTATACCATAAAAAAGCGcc  >  8:1223276/1‑150 (MQ=18)                                                                                                                                              aaaaGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTATTGAACACCCTAACGGGTGTTTT‑TTTg                            >  5:259559/1‑124 (MQ=255)                                                                                                                                              aaaaGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTATTGAACACCCTAACGGGTGTTTT‑TTTg                            <  8:259199/124‑1 (MQ=255)                                                                                                                                              aaaaGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTATTGAACACCCTAACGGGTGTTTT‑TTTGTTTCTGGTCTCCCATAAAAAAGCGcc  >  3:527098/1‑150 (MQ=255)                                                                                                                                              |                                                                                                                                                       CCCCCGCAACCAATTAAAATTTGATGAAGTAAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAA‑TC‑AAATTTGATGAAGTAAAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTATTGAACACCCTAACGGGTGTTTT‑TTTGTTTCTGGTCTCCCATAAAAAAGCGCC  >  NC_000913/2947452‑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: ‑