| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | NC_000913 | 2,947,591 | 0 | A | . | 21.3% | 132.5 / 32.8 | 47 | intergenic (+18/‑16) | metW/metV | tRNA‑Met/tRNA‑Met |
| Reads supporting (aligned to +/- strand): ref base A (26/11); new base . (7/3); total (33/14) | |||||||||||
| 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.50e-01 | |||||||||||
| Rejected as consensus: Frequency below/above cutoff threshold. | |||||||||||
GACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATCAAATTTG‑ATGAAGTAAAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTATTGAACACCCTAACGGGTGTTTTTTTGTTTCTGGTCTCCCATAAAAA > NC_000913/2947495‑2947734 | gACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATCAAATTTG‑ATGAAGTa < 1:2613802/97‑1 (MQ=39)gACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATCAAATTTG‑ATGAAGTAAAAGCa < 1:2613785/103‑1 (MQ=255)gACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATCAAATTTG‑ATGAAGTAAAAGCa < 1:2613786/103‑1 (MQ=255)gACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATCAAATTTG‑ATGAAGTAAAAGCAGTa < 1:2613779/106‑1 (MQ=255)gACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATCAAATTTG‑ATGAAGTAAAAGCAGTAc < 1:2613812/107‑1 (MQ=255)gACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATCAAATTTG‑ATGAAGTAAAAGCAGTACgg < 1:2613816/109‑1 (MQ=255)gACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATCAAATTTG‑ATGAAGTAAAAGCAGTACGGt < 1:2613792/110‑1 (MQ=255)gACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATCAAATTTG‑ATGAAGTAAAAGCAGTACGGTg < 1:2613797/111‑1 (MQ=255)gACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATCAAATTTG‑ATGAAGTAAAAGCAGTACGGTGa < 1:2613814/112‑1 (MQ=255)gACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATCAAATTTG‑ATGAAGTAAAAGCAGTACGGTGAcgc < 1:2613781/115‑1 (MQ=255) aCGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATCAAATTTG‑ATGAAGTAAAAGCAGTACGGt < 1:2613822/109‑1 (MQ=255) gcCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATCAAATTTG‑ATGAAGTAAAAGCAGTAc < 1:2613825/91‑1 (MQ=255) ttCAAATCCGGCCCCCGCAACCAATAAAATTTG‑ATGAAGTAAAAGCAGTAAGGTGAAGAGGGGTGGCGAAGCCTGGAAGCTCGTCGGGCTCATAACCCGAAGGGCGGCGGGACAAATCCGGCCCCCGCCACCATTTAATGAACAcccaa > 2:2613916/1‑147 (MQ=14) cccccGCAACCAATCAAATTTG‑ATGAAGTAAAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCTTAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTATTGAACACCCTAACGGGTGtttt > 2:2613952/1‑149 (MQ=255) cccccGCAACCAATAAATTTTC‑AAGAAGTAAAAGCAGTAATTTTAAGAGGGGTGGAGAAGAATGGTAGCTCGTCGGGCTCATAAACCGAAGGTCGTCGGTTCAAATACGGCCCCCGCAACCAATTATTGAACACCCTAACGGGTGtttt > 2:2613911/1‑149 (MQ=17) ccccGCAACCAATCAAATTTG‑ATGAAGTAAAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTATTGAACACCCTAACGGGTGttttt > 2:2613908/1‑149 (MQ=255) ccGCAACCAATCAAATTTG‑ATGAAGTAAAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGGTCAAATCCGGCCCCCGCAACCAATTATTGAACACCCTAACGGGTGttttttt > 2:2613866/1‑149 (MQ=255) cGCAACCAATCAAATTTG‑ATGAAGTAAAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTATTGAACACCCTAACGGGTGttttttt > 1:2613956/1‑148 (MQ=255) aCCAATTAAAATTGTATTAAGT‑AAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTATTGAACACCCTAACGGGTGTTTTTTTGTTTc < 1:2613835/149‑1 (MQ=14) aCCAATCAAATTTG‑ATGAAGTAAAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTATTGAACACCCTAACGGGTGTTTTTTTGTTTc > 2:2613982/1‑149 (MQ=255) aattaaaaTTTG‑ATGAAGT‑AAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTATTGAacac > 2:2613831/6‑123 (MQ=0) aattaaaaTTTG‑ATGAAGT‑AAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTATTGAacac < 1:2613831/118‑1 (MQ=0) aaTCAAATTTG‑ATGAAGTAAAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTATTGAACACCCTAACGGGTGTTTTTTTGTTTCTgg > 1:2613955/1‑149 (MQ=255) aaTCAAATTTG‑ATGAAGTAAAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTATTGAACACCCTAACGGGTGTTTTTTTGTTTCTgg > 2:2613978/1‑149 (MQ=255) aaTCAAATTTG‑ATGAAGTAAAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTAGTCGGGCTCATCACCCGAAGTGCGTCGGTTCAAATCCGGCGCCCGCAACCAATTATTGAACACCCTAACGGGAGTTTTTTTGTTTCTgg > 1:2613914/1‑149 (MQ=18) ttaaaaTTTG‑ATGAAGT‑AAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTATTGAACAcc > 2:2613832/4‑122 (MQ=12) ttaaaaTTTG‑ATGAAGT‑AAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTATTGAACAcc < 1:2613832/119‑1 (MQ=12) aTCAAATTTG‑ATTAAGTAAAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTATCTCGTCGGGCTCATAGCCCGAAGGTCGTCTGATCAATTCCGGGCCCCGCAAGCAATTATTGAACACCATAACGGGTGTTTTTTTGTTTCTGGt > 2:2613968/1‑149 (MQ=18) aTCAAATTTG‑ATGAAGTAAAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTATTGAACACCCTAACGGGTGTTTTTTTGTTTCTGGt > 2:2613969/1‑149 (MQ=255) aTCAAATTTG‑ATGAAGTAAAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTATTGAACACCCTAACGGGTGTTTTTTTGTTTCTGGt > 1:2613944/1‑149 (MQ=255) aTCAAATTTG‑ATGAAGTAAAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTATTGAACACCCTAACGGGTGTTTTTTTGTTTCTGGt > 2:2613925/1‑149 (MQ=255) tCAAATTTG‑ATGAAGTAAAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTATTGAACACCCTAACGGGTGTTTTTTTGTTTCTGTtc > 2:2613939/1‑149 (MQ=255) tCAAATTTG‑ATGAAGTAAAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTATTGAACACCCTAACGGGTGTTTTTTTGTTTCTGGtc > 2:2613994/1‑149 (MQ=255) tCAAATTTG‑ATGAAGTAAAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTATTGAACACCCTAACGGGTGTTTTTTTGTTTCTGGtc > 1:2613998/1‑149 (MQ=255) tCAAATTTG‑ATGAAGTAAAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTATTGAACACCCTAACGGGTGTTTTTTTGTTTCTGGtc > 1:2613937/1‑149 (MQ=255) aaaTTTG‑ATGAAGTAAAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTATTGAACACCCTAACGGGTGTTTTTTTGTTTCTGGtctc > 2:2613972/1‑149 (MQ=255) aaaTTTG‑ATGAAGT‑AAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACACGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTATTGAACAACCTTACGGGTGTTTTTTTGTTGCTGTTCTcc > 2:2613920/1‑149 (MQ=17) aaTTTG‑ATGAAGTAAAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTATTGAACACCCTAACGGGTGTTTTTTTGTTTCTGGtctc > 1:2613897/1‑148 (MQ=255) aTTTG‑ATGAAGTAAAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTGCAAATCCGGCCCCCGCAACCAATTATTGAacac > 1:2613993/1‑117 (MQ=255) aTTTG‑ATGAAGTAAAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACACGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTATTGAACACCCTAACGGGTGTTTTTTTGTTTCTGGTCTccc > 2:2613899/1‑149 (MQ=255) tttG‑ATGAAGTAAAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTATTGAACACCCTAACGGGTGTTTTTTTGTTTCTGGTCTcc > 1:2613990/1‑147 (MQ=255) ttG‑ATGAAGT‑AAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTATTGAACACCCTAACGGGTGTTTTTTTGTTTCTGGTCTCCCATa > 2:2613841/1‑149 (MQ=255) g‑atgaAGTAAAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTATTGAACACCCTAACGGGTGTTTTTTTGTTTCTGGTCTCCCATaa > 2:2613913/1‑149 (MQ=255) g‑atgaAGT‑AAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTATTGaa > 2:2613838/1‑109 (MQ=14) g‑atgaAGT‑AAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTATTGAACACCCTAACGGGTGTTTTTTTGTTTCTGGTCTCCCATaaa > 2:2613865/1‑149 (MQ=255) atgaAGTAAAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTATTGAACACCCTAACGGGTGTTTTTTTGTTTCTGGTCTCCCATaaa > 2:2613895/1‑149 (MQ=255) tgaAGT‑AAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTATTGAACACCCTAACGGGTGTTTTTTTGt > 2:2613834/1‑130 (MQ=255) gaAGTAAAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCCACCAATTATTGAACACCCTAACGGGTGTTTTTTTGTTTCTGGTCTCCCATaaaaa > 2:2613981/1‑149 (MQ=255) aaaaGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCgg < 2:2613829/47‑1 (MQ=32) | GACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATCAAATTTG‑ATGAAGTAAAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTATTGAACACCCTAACGGGTGTTTTTTTGTTTCTGGTCTCCCATAAAAA > NC_000913/2947495‑2947734 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 21 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |