Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,947,591	0	A	.	31.5%	214.0 / 108.7	92	intergenic (+18/‑16)	metW/metV	tRNA‑Met/tRNA‑Met
Reads supporting (aligned to +/- strand):  ref base A (58/5);  new base . (28/1);  total (86/6)
Fisher's exact test for biased strand distribution p-value = 6.61e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.47e-01

CGCAACCAATTAAAATTTGATGAAGTAAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATCAAATTT‑GATGAAGTAAAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTATTGAACACCC‑TA‑ACGGGTGTTT‑TTTTGTTTCTGGTCTCCCATAAAAAAGCGCC  >  NC_000913/2947456‑2947740                                                                                                                                         |                                                                                                                                                         ggccagcAATTATAATTTGATGAAGTTTCGGCATACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATCAAATTT‑GATGAAGTAAAAGCAGTACGGTGa                                                                                                                                           <  2:66793/145‑1 (MQ=17)    aaCCAATTAAAATTTGATGAAGTAAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAATTTA‑GATGAAGGAAAAGccgaggggtgcgg                                                                                                                                         >  2:95197/1‑138 (MQ=14)        aaTTAAAATTTGATGAAGTAAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAATGAT‑GATGATCTAAAAGCAGTACGGTGACGCgggg                                                                                                                                    >  3:80754/1‑151 (MQ=17)        aaTTAAAATTTGATGAAGTAAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATGAATTTT‑CAGGAAGTAAAGGCTGTACTGTGACGCgggg                                                                                                                                    >  3:89573/1‑151 (MQ=17)                                         aCGCGGGTTGGAGCAGCCTGGTAGCTCTGCTGGATCTTAACCCGAAGGGGGTCGGTGCCTATCAGGCCCGGGCAACGAATCAACTTTCGTTGTA‑TACTCGCCCTCCCCGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAAc                                                                                                   <  7:548489/151‑1 (MQ=2)                                                       agccttttatgtcggcgGGATCATAACCCGAAGGTCGGCTGTTGTAATCCGGCCCGCGCAAGCAATAAAATTT‑GATTAATTAAACGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCgg                                                                                     <  3:665568/140‑1 (MQ=11)                                                                     tcgGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATGAAATTT‑GATGATGTAAAAGCAGTACGGTGACGCGGGGTGGAGAAGCCTGGTAGCTGGTCGGGCTCATAACGCGAAGGGCGTGGGTTCACCTCCGGccc                                                                       >  7:1106667/1‑151 (MQ=11)                                                                                           gtcgGTTCAAATCCGGCCCCCGCAACCAATAAAATTT‑GATGAAGTAAAAGCAGTACGGTGAcgcg                                                                                                                                       <  2:1184912/65‑1 (MQ=255)                                                                                                 tCAAATCCGGCCCCCGCAACCAATGATATAA‑GATGATGTAGCAGCAGTAAGGTGACGAGGGGTCGCCAAGCAACGAGGCAGGGCGGGCGGATAACCCGAAGGGCGGCGGTTCCAGTCCGGCCCCCGCAACCAATTATTGAACACCC‑Tacac                                          >  4:1633669/1‑148 (MQ=9)                                                                                                  cAAATCCGGCCCCCGCAACCAATCAATTTT‑CAAGAAGTAAAAGCTGTACGTTTAAGGGGGGTGCAGAAGAATGGAGGCTCTTCGCGCGAATAACCCGAAGGTGGTCGGTTAAAATCCGGCCCCCGCAAAAAATTATTGAACACCC‑TA‑AAgg                                        >  4:200632/1‑151 (MQ=12)                                                                                                       ccGGCCCCCGCAACCAATGAAATTT‑GATGAAGTAAAAGCAGTACGGTGACGCGGGGTGGCGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGACGCGCGGCGGTGCAACTCCGGCCCCCGCCACCAATTATTGAACACCC‑TC‑ACGGCTGtt                                   >  1:1406655/1‑151 (MQ=17)                                                                                                        cGGCCCCCGCAACCAATCAAATTT‑GATGAAGTAAAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTATTTAACACCC‑TA‑ACGGGAGttt                                  >  3:551490/1‑151 (MQ=255)                                                                                                          gCCCCCGCAACCAATCAAATTT‑GATGAAGTAAACGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCGCATAACCCGAAGGTCGTCGGGTCACATCCGGCCCCCGCAACCAATTATTGAACACCC‑TA‑CCGGGGGtttgt                                >  2:1633180/1‑149 (MQ=17)                                                                                                          gCCCCCGCAACCAATCAAATTT‑GATGAAGTAAAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTATTGAACACCC‑TA‑ACGGGTGttt‑tt                               >  7:268777/1‑151 (MQ=255)                                                                                                          gCCCCCGCAACCAATCAAATTT‑GATGAAGTAAAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTATTGAACACCC‑TA‑ACGGGTGttt‑tt                               >  5:472476/1‑151 (MQ=255)                                                                                                          gCCCCCGCAACCAATCAAATTT‑GATGAAGTAAAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTATTGAACACCC‑TA‑ACGGGTGttt‑tt                               >  6:1199662/1‑151 (MQ=255)                                                                                                          gCCCCCGCAACCAATCAAATTT‑GATGAAGTAAAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTATTGAACACCC‑TA‑ACGGGTGttt‑tt                               >  4:799864/1‑151 (MQ=255)                                                                                                          gCCCCCGCAACCAATCAAATTT‑GATGAAGTAAAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTATTGAACACCC‑TA‑ACGGGTGttt‑tt                               >  4:910152/1‑151 (MQ=255)                                                                                                           cccccGCAACCAATCAAATTT‑GATGAAGTAAAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCCACCAATTATTGAACACCC‑TA‑ACGGGTGttt‑ttt                              >  2:1223946/1‑151 (MQ=255)                                                                                                           cccccGCAACCAATCAAATTT‑GATGAAGTAAAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTATTGAACACCC‑TA‑ACGGGTGttt‑ttt                              >  7:476663/1‑151 (MQ=255)                                                                                                            ccccGCAACCAATCAAATTT‑GATGAAGTAAAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTATTGAACACCC‑TA‑ACGGGTGttt‑tttt                             >  1:1037984/1‑151 (MQ=255)                                                                                                             cccGCAACCAATCAAATTT‑GATGAAGTAAAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCAGGCCCCCGCAACCAATTATTGAACACCC‑TA‑ACGGGTGTTT‑TTTTg                            >  6:1677559/1‑151 (MQ=255)                                                                                                              ccGCAACCAATCAAATTT‑GATGAAGTAAAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTATTGAACACCC‑TA‑ACGGGTGTTT‑TTTTGt                           >  3:703108/1‑151 (MQ=255)                                                                                                              ccGCAACCAATCAAATTT‑GATGAAGTAAAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTATTGAACACCC‑TA‑ACGGGTGTTT‑TTTTGt                           >  1:332411/1‑151 (MQ=255)                                                                                                               cGCAACCAATCAAATTT‑GATGAAGTAAAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATAATTGAACACCC‑TA‑ACGGGTGTTT‑TTTTGtt                          >  5:205887/1‑151 (MQ=255)                                                                                                                  accaattaaaaTTT‑GATGAAGT‑AAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTATTGAACACCC‑TA‑TCCGGGGTTT‑TGTTGTTTCTg                      >  7:772533/9‑151 (MQ=11)                                                                                                                  accaattaaaaTTT‑GATGAAGT‑AAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTATTGAACACCC‑TA‑ACGCGGGTTT‑TTTTGGTtgtg                      >  5:1315372/9‑151 (MQ=11)                                                                                                                  aaCCAATCAAATTT‑GATGAAGTAAAAGCAGTACGGTGCCGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATACCCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAAGTATTGAACACCC‑TA‑ACGGGGGTTT‑Ttgtggtgcg                       >  6:834205/1‑146 (MQ=18)                                                                                                                  aaCCAATCAAATTT‑GATGAAGTAAAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTATTGAACACCC‑TA‑ACGGGTGTTT‑TTTTGTTTCt                       >  4:137377/1‑151 (MQ=255)                                                                                                                  aaCCAATCAAATTT‑GATGAAGTAAAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGACGTCCGTTGAAATCCGGCCCCCGCAACCAATAATTGAACACCC‑TA‑ACGGCTGTTT‑TATTGTTACt                       >  5:577157/1‑151 (MQ=18)                                                                                                                   ccaattaaaaTTT‑GATGAAGT‑AAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCACCCAATGATTGAACACCC‑TC‑ACCGCTGTGT‑TGTGGTTTGTgg                     >  4:1256200/8‑151 (MQ=1)                                                                                                                   ccaattaaaaTTT‑GATGAAGT‑AAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAATTAACACCA‑TA‑ACCGGTGTTT‑TTTTGTATCTgg                     >  3:1583398/8‑151 (MQ=11)                                                                                                                   ccaattaaaaTTT‑GATGAAGT‑AAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAATGAACAACC‑AA‑ACAGGTGTTT‑TTTTGTTTCTgg                     >  8:988065/8‑151 (MQ=14)                                                                                                                   ccaattaaaaTTT‑GATGAAGT‑AAAGCAGTACGGTGAAGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAATTAACACCA‑TA‑AAAAGTGTTTATTTTGtacagg                      >  2:365149/8‑146 (MQ=11)                                                                                                                   aCCAATCAAATTT‑GATGAAGTAAAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTATTGAACACCC‑TA‑ACGGGTGTTT‑TTTTGTTTCTg                      >  8:267191/1‑151 (MQ=255)                                                                                                                   aCCAATCAAATTT‑GATGAAGTAAAAGCAGTACGGTGACGAGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCAAAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTATTGAACACCC‑TA‑CCGGGTGTTT‑TTTTGTTTCTg                      >  4:949328/1‑151 (MQ=255)                                                                                                                    ccAATCAAATTT‑GATGAAGTAAAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTATTGAACACCC‑TA‑ACGGGTGTTT‑TTTTGTTTCTgg                     >  2:1482745/1‑151 (MQ=255)                                                                                                                    ccAATCAAATTT‑GATGAAGTAAAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTATTGAACACCC‑TA‑ACGGGTGTTT‑TTTTGTTTCTgg                     >  1:878657/1‑151 (MQ=255)                                                                                                                    ccAATCAAATTT‑GATGAAGTAAAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTATTGAACACCC‑TA‑ACGGGTGTTT‑TTTTCTTTCTgg                     >  2:111597/1‑151 (MQ=255)                                                                                                                    caattaaaaTTT‑GATGAAGT‑AAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATGATTGATCACCA‑TA‑TCTCGTGTTT‑TTTTGGTTCTggg                    >  4:1035197/7‑150 (MQ=11)                                                                                                                     aattaaaaTTT‑GATGAAGT‑AAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTATTTTACACCA‑TG‑ACGGGGGTGT‑TTGTGttggtgggg                   >  3:628107/6‑144 (MQ=1)                                                                                                                     aattaaaaTTT‑GATGAAGT‑AAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTATTGAACACCC‑TG‑AAGGGGGTTT‑TGTTGGTTGTgggg                   >  3:327038/6‑149 (MQ=11)                                                                                                                     aattaaaaTTT‑GATGAAGT‑AAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAATGAACAAGC‑TA‑AAGGGTGAGT‑TTTTGTTTCGGGtg                   >  8:1262260/6‑150 (MQ=11)                                                                                                                     aattaaaaTTT‑GATGAAGT‑AAAGCAGTACGGTGACGCGGGGTGGAGAAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGCTTCAAATCCGGCCCCCGCAACCAATAAATGAACAACC‑TA‑ACCTGTGATT‑TTATGTTTcagttc                   >  1:1050329/6‑146 (MQ=11)                                                                                                                      aaTCAAATTT‑GATGAAGTAAAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTATTGAACACCC‑TA‑ACGGGTGTTT‑TTTTGTTTCTGGtc                   >  1:1596350/1‑151 (MQ=255)                                                                                                                      aaTCAAATTT‑GATGAAGTAAAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTATTGAACACCC‑TA‑ACGGGTGTTT‑TTTTGTTTCTGGtc                   >  8:828558/1‑151 (MQ=255)                                                                                                                      aaTCAAATTT‑GATGAAGTAAAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTATTGAACACCC‑TA‑ACGCGTGTTT‑TTTTGTTTCTGGtc                   >  3:110057/1‑151 (MQ=255)                                                                                                                       ttaaaaTTT‑GATGAAGT‑AAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATGATTTAACAACA‑TA‑TATGGGGTTT‑TTTTTTATGTGGtctg                 >  1:154794/4‑150 (MQ=11)                                                                                                                       aTCAAATTT‑GATGAAGTAAAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTATTGAACACCC‑TA‑ACGGGTGTTT‑TTTTGTTTCTGGtct                  >  4:1196677/1‑151 (MQ=255)                                                                                                                       aTCAAATTT‑GATGAAGTAAAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTATTGAACACCC‑TA‑ACGGGTGTTT‑TTTTGTTTCTGGtct                  >  3:759738/1‑151 (MQ=255)                                                                                                                       aTCAAATTT‑GATGAAGTAAAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTATTGAACACCC‑TA‑ACGGGTGTTT‑TTTTGTTTCTGGtct                  >  3:1555674/1‑151 (MQ=255)                                                                                                                        taaaaTTT‑GATGAAGT‑AAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTATTTAACACCC‑TA‑AAGGGTGATT‑TTTTGTTTGAGGtctgc                >  2:179692/3‑149 (MQ=14)                                                                                                                        tCAAATTT‑GATGAAGTAAAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTATTGAACACCC‑TA‑ACGGGTGTTT‑TTTTGTTTCTGGtctc                 >  1:1285356/1‑151 (MQ=255)                                                                                                                         cAAATTT‑GATGAAGTAAAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTATTGAACACCC‑TA‑ACGGGTGTTT‑GTTTGTTTCTGGTCTcc                >  8:1597678/1‑151 (MQ=255)                                                                                                                         aaaaTTT‑GATGAAGT‑AAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTATTGAACACCC‑TC‑AAGGGTGTGT‑TTTGGGTTCTGGTCTccc               >  3:1125795/2‑151 (MQ=14)                                                                                                                         aaaaTTT‑GATGAAGT‑AAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATAATTGAACACCC‑TA‑ACGCGTGTTT‑TATTTTTTCAGGTCTccc               >  1:117795/2‑151 (MQ=14)                                                                                                                          aatttt‑GATGAAGT‑AAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTATTGAACACCC‑TA‑ACGGGTGTTT‑TTTTGTTTCTGGTCTCCCa              <  4:1608498/148‑1 (MQ=255)                                                                                                                          aaaTTT‑GATGAAGTAAAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTATTGAACACCC‑TA‑ACGGGTGTTT‑TTTTGTTTCTGGTCTccc               >  1:664611/1‑151 (MQ=255)                                                                                                                          aaaTTT‑GATGAAGTAAAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTATTGAACACCC‑TA‑ACGGGTGTTT‑TTTTGTTTCTGGTCTccc               >  3:658253/1‑151 (MQ=255)                                                                                                                          aaaTTT‑GATGAAGTAAAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTATTGAACACCC‑TA‑ACGGGTGTTT‑TTTTGTTTCTGGTCTccc               >  3:524058/1‑151 (MQ=255)                                                                                                                          aaaTTT‑GATGAAGTAAAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTATTGAACACCC‑TA‑ACGGGTGTTT‑TTTTGTTTCTGGTCTccc               >  3:1431586/1‑151 (MQ=255)                                                                                                                          aaaTTT‑GATGAAGTAAAAGCAGTACGGTGACGCGGGGTAGAGCAGCCTGGTAGCTCGTCGGGCTCATAGCCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTATTGAACACCC‑TA‑ACGGGTGTTT‑TTTTGTTTCTGGTCTccc               >  1:709889/1‑151 (MQ=255)                                                                                                                          aaaTTT‑GATGAAGT‑AAAGCAGTACGGTGACGCGTGGAGGAGCAGCCTGCTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCACCCGCAACCAATTAATGAACAACA‑AA‑AAAAGAGATT‑TATTTTTTCTGGTGTccc               >  1:1494553/1‑150 (MQ=11)                                                                                                                          aaaTTT‑GATGAAGT‑AAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCCACCAATGATTGATCACCC‑TC‑ACCGCGGTGT‑TTTTGTTTCgggggtgcac              >  4:647207/1‑144 (MQ=11)                                                                                                                          aaaTTT‑GATGAAGT‑AAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCATCGGTTCAAATCCGGCCCCCGCAACCAATAAATGAACAACC‑AC‑AAGAGGGATT‑TTTTGTATCTGGTCTcca               >  1:502698/1‑149 (MQ=11)                                                                                                                           aaTTT‑GATGAAGT‑AAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTATTGAACAGCA‑TG‑AAGGGGGTTT‑TTGTGTTTGTGGTGTCGCa              >  4:1155024/1‑150 (MQ=11)                                                                                                                           aaTTT‑GATGAAGT‑AAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATGATTTAACACCCATT‑ACGGG‑GTTT‑TGTTTGTTCTGGTCTcccg              >  3:225124/1‑149 (MQ=11)                                                                                                                           aaTTT‑GATGAAGT‑AAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATGATTTAACAACC‑TC‑AAGGCTGATT‑ATTTGttcgtgggctgcat              >  5:345442/1‑138 (MQ=11)                                                                                                                            aTTT‑GATGAAGTAAAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAGATCCGGCCCCCGCAACCAATTATTGAACACCC‑TA‑ACGGGTGTTT‑TTTTGTTTCTGGTCTCCCAt             >  8:1339343/1‑151 (MQ=255)                                                                                                                            aTTT‑GATGAAGTAAAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTATTGAACACCC‑TA‑ACGGGTGTTT‑TTTTGTTTCTGGTCTCCCAt             >  5:883730/1‑151 (MQ=255)                                                                                                                            aTTT‑GATGAAGT‑AAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTATTGAACACCC‑TA‑AACGGTGATT‑TTTTGTTTCTGGTCTCCCAt             >  2:1250056/1‑150 (MQ=255)                                                                                                                             ttt‑GATGAACTAAAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGCTAGCTCGTCGGGCTCATAACCCGAAGCTCGTCCGTTCAAATCCGCCCCCCGCAACCAATTATTCAACACCC‑TA‑ACCGGTGTTA‑ATTTGTTTCTGCTCTCCCaa             >  2:637901/1‑149 (MQ=17)                                                                                                                              tt‑GATGAAGTAAAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACTCGAAGGACGTCGGTTCAAATCCGGCCCCCGCAACCAATTATTGAACACCC‑TA‑ACGGGTGTTT‑TTTTGTTTCTGGTCTCCCATa            >  2:211184/1‑150 (MQ=255)                                                                                                                              tt‑GATGAAGT‑AAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCGACCAATGATTTAACACCC‑TATACGGG‑GTGT‑TTTTTTTTCTGGTCTCCCATaa           >  5:547800/1‑150 (MQ=14)                                                                                                                              tt‑GATGAAGT‑AAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTATTGAACACCA‑TA‑AAGCCTGTTT‑TTGTGTTTCTGGTCTCCCATaa           >  6:327060/1‑150 (MQ=17)                                                                                                                               t‑gatgaAGTAAAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTATTGAACACCC‑TA‑ACGGGTGTTT‑TTTTGTTTCTGGTCTCCCATaa           >  3:226179/1‑150 (MQ=255)                                                                                                                                 gatgaAGTAAAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTATTGAACACCC‑TA‑ACGGGTGTTT‑TTTTGTTTCTGGTCTCCCATaaa          >  7:1338941/1‑150 (MQ=255)                                                                                                                                 gatgaAGT‑AAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTATAGAACACCC‑TC‑AAGGGTGTTT‑TTTTTTTTCTGGTCTCCCATataa         >  7:232746/1‑150 (MQ=17)                                                                                                                                  atgaAGTAAAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTATTGAACACCC‑TA‑ACGGGTGTTT‑TTTTGTTTCTg                      <  1:76815/137‑1 (MQ=255)                                                                                                                                  atgaAGTAAAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTATTGAACACCC‑TA‑ACGGGTGTTT‑TTTTGTTTCTg                      >  7:76609/1‑137 (MQ=255)                                                                                                                                  atgaAGTAAAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTATTGAACACCC‑TA‑ACGGGTGTTT‑TTTTGTTTCTGGTATCACATaaaa         >  1:1000947/1‑150 (MQ=255)                                                                                                                                  atgaAGT‑AAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTATTTAACACCC‑TA‑AAAGGGGTTT‑TTTTGTTGCTGGTCTCCCATaaaaac       >  1:1302204/1‑150 (MQ=17)                                                                                                                                  atgaAGT‑AAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTATTGAACAACC‑TA‑AAAGGTGTTT‑TATTGTTACTGGtct                  >  6:376984/1‑140 (MQ=14)                                                                                                                                   tgaAGTAAAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAACTCCGGCCCCCGCACCCAATTATTGAACACCC‑TA‑ACGGGTGTTT‑TTTTGTTTCTGGTCTCCCATaaaaa        >  3:1035375/1‑150 (MQ=255)                                                                                                                                   tgaAGTAAAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTATTGAACACCC‑TA‑ACGGGTGTTT‑TTTTGTTTCTGGTCTCCCATaaaaa        >  4:972499/1‑150 (MQ=255)                                                                                                                                   tgaAGTAAAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTATTGAACACCC‑TA‑ACGGGTGTTT‑TTTTGTTTCTGGTCTCCCATaaaaa        >  1:1404666/1‑150 (MQ=255)                                                                                                                                    gaAGTAAAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAACTCCGGCCCCCGCACCCAATTATTGAACACCC‑TA‑ACGGGTGTTT‑TTTTGTTTCTGGTCTCCCATaaaaa        >  3:251121/1‑149 (MQ=255)                                                                                                                                    gaAGTAAAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTATTGAACACCC‑TA‑ACGGGTGTTT‑TTTTGTTTCTGGTCTCCCATaaaaa        >  4:1651779/1‑149 (MQ=255)                                                                                                                                    gaAGTAAAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTATTGAACACCC‑TA‑ACGGGTGTTT‑TTTTGTTTCTGGTCTCCCATaaaaa        >  3:63294/1‑149 (MQ=255)                                                                                                                                    gaAGTAAAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTATTGAACACCC‑TA‑ACGGGTGTTT‑TTTTGTTTCTGGTCTCCCATaaaaa        >  3:873618/1‑149 (MQ=255)                                                                                                                                    gaAGT‑AAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATGATTGATCACGA‑TG‑AAGCATGTTT‑TGTGGGTAGGGGTGTCCCATATAAAAgc     >  3:887430/1‑151 (MQ=11)                                                                                                                                        tAAAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTATTGAACACCC‑TA‑ACGGGTGTTT‑TTTTGTTTCTGGTCTCCCATAAAAAAGCGcc  >  5:492242/1‑151 (MQ=255)                                                                                                                                         |                                                                                                                                                         CGCAACCAATTAAAATTTGATGAAGTAAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATCAAATTT‑GATGAAGTAAAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTATTGAACACCC‑TA‑ACGGGTGTTT‑TTTTGTTTCTGGTCTCCCATAAAAAAGCGCC  >  NC_000913/2947456‑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: ‑