Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | NC_000913 | 3,210,839 | A→T | 17.2% | K20M (AAG→ATG) | rpsU → | 30S ribosomal subunit protein S21 |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_000913 | 3,210,839 | 0 | A | T | 17.2% | 64.2 / 8.2 | 35 | K20M (AAG→ATG) | rpsU | 30S ribosomal subunit protein S21 |
Reads supporting (aligned to +/- strand): ref base A (8/20); new base T (4/2); total (13/22) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 1.54e-01 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.99e-01 |
AATTGCAGTAAAATTCCGCACCATTTTGAAATAAGCTGGCGTTGATGCCAGCGGCAAACCGAATTAATCAAAGGTGAGAGGCACATGCCGGTAATTAAAGTACGTGAAAACGAGCCGTTCGACGTAGCTCTGCGTCGCTTCAAGCGTTCCTGCGAAAAAGCAGGTGTTCTGGCGGAAGTTCGTCGTCGTGAGTTCTATGAAAAACCGACTACCGAACGTAAGCGCGCTAAAGCTTCTGCAGTGAAACGTCACGCGAAGAAACTGGCTCGCGAAAACGCACGC > NC_000913/3210697‑3210978 | aaTTGCAGGAAAACTCCGCACCAGTTTGAAATAAGCTGGCGTTGATGCCAGCGGCAAACCGAATTAATCCAACATGAGACGCACATGCCGGTAATTAAAGTACGTGAAAACGAGCCGTTCGACGTAGCTCTGCGTCGCTTCaa < 2:290213/143‑1 (MQ=255) ttGCAGTAAAATTCCGCACCATTTTGAAATAAGCTGGCGTTGATGCCAGCGGCAAACCGAATTAATCAAAGGTGAGAGGCACATGCCGGTAATTAAAGTACGTGAAAACGAGCCGTTCGACGTAGCTCTGCGTCGCTTCAAGc > 1:320802/1‑143 (MQ=255) agggaaaTTGCGCACCAGTTTGAAGTAAGCTGGCGGTGCTGCCAGCGGGAAACCGAGTTAATCAAAGGTGAGAGGCACATGCCGGTAATTAAAGTACGTGAAAACGAGCCGTTCGACGTAGCTCTGCGTCGCTTCAAGCGTTc < 2:594978/139‑1 (MQ=255) tAAAATTCCGCACCATTTTGAAATAAGCTGGCGTTGATGCCAGCGGCAAACCGACTTAATCAAAGGTGAGAGGCACATGCCGGCAAGTAACGTACGTGAAAACGAGCCGTTCGACGTAGCTCTGCGTCGCTTCCCGCGTTCCt > 1:774802/1‑143 (MQ=255) gCACCATTTTGAAATAAGCTGGCGTTGATGCCAGCGGCAAACCGAATTAATCAAAGGTGAGAGGCACATGCCGGTAATTAAAGTACGTGAAAACGAGCCGTTCGACGTAGCTCTGCGTCGCTTCAAGCGTTCATGCGAAAAAg > 1:383740/1‑143 (MQ=255) cACCATTTTGAAATAAGCTGGCGTTGATGCCAGCGGCAAACCGAATTAATCAAAGGTGAGAGGCACATGCCGGTAATTAAAGTACGTGAAAACGAGCCGTTCGACGTAGCTCTGCGTCGCTTCATGCGTTCCTGCGAAAAAGc > 1:96922/1‑143 (MQ=255) cACCATTTTGAAATAAGCGGGCGTTGATGCCAGCGGCAAACCGAATTAATCAAAGGTGAGAGGCACATGCCGGTAATTAAAGTACGTGAAAACGAGCCGTTCGACGTAGCTCTGCGTCGCTTCAAGCGTTCCTGCGAAAAAGc < 2:759295/143‑1 (MQ=255) ccATTTTGAAATAAGCTGGCGTTGATGCCAGCGGCAAACCGAATTAATCAAAGGTGAGAGGCACATGCCGGTAATTAAAGTACGTGAAAACGAGCCGTTCGACGTAGCTCTGCGTCGCTTCAAGCGTTCCTGCGAAAAAGCAg < 1:510647/143‑1 (MQ=255) ccATTTTGAAATAAGCTGGCGTTGATGCCAGCGGCAAACCGAATTAACCAAAGGTGAGAGGCACATGCCGGTAATTAAAGTACGTGAAAACGAGCCGTTCGACGTAGCTCTGCGTCGCTTCAAGCGTTCCTGCGAAAAAGCAg < 1:290375/143‑1 (MQ=255) aTAAGCTGGCGTTGATGCCAGCGGCAAACCGAATTAATCAAAGGTGAGAGGCACATGCCGGTAATTAAAGTACGTGAAAACGAGCCGTTCGACGTAGCTCTTCGTCGCTTCAAGCGTTCCTGCGAAAAAGCAGGTGTTCTGGc > 1:672459/1‑143 (MQ=255) ggCGTTGATGCCAGCGGCAAACCGAATTAATCAAAGGTGAGAGGCACATGCCGGTAATTAAAGTACGTGAAAACGAGCCGTTCGACGTAGCTCTGCGTCGCTTCAAGCGTTCCTGCGAAAAAGCAGGTGTTCTGGCGGAAGtt > 2:775412/1‑143 (MQ=255) ggCGTTGATGCCAGCGGCAAACCGAATTAATCAAAGGTGAGAGGCACATGCCGGTAATTAAAGTACGTGAAAACGAGCCGTTCGACGTAGCTCTGCGTCGCTTCAAGCGTTCCTGCGAAAAAGCAGGTGTTCTGGCGGAAGtt < 1:597582/143‑1 (MQ=255) ttGATGCCAGCGGCAAACCGAATTAATCAAAGGTGAGAGGCACATGCCGGTAATTAAAGTACGTGAAAACGAGCCGTTCGACGTAGCTCTGCGTCGCTTCAAGCGTTCCTGCGAAAAAGCAGGTGTTCTGGCGGAAGTtcgtc < 2:444687/143‑1 (MQ=255) tGATGCCAGCGGCAAACCGAATTAATCAAAGGTGAGGGGCTCCGCCCGGTAATTAAAGTACGTGAAAACGAGCCGTTCGACGTAGCTCTGCGTCGCTTCAAGCGTTCCTGCGAAAAAGCAGGTGTTCTGGCGGAAGTtcgtcg < 2:580591/143‑1 (MQ=255) cagcagCAAAGCTAATCAATCAAAGCTGAGAGCCACATGCCGGTAAGTAAACTACGTGCAAACGAGCCGTTCGACGTAGCTCTCCGTCGCTTCAAGCGTTCCTGCGAAAAAGCAGGTGTTCTGGCGGAAGTTCGTCGTCGTGa < 1:56472/143‑1 (MQ=255) cAAACAGAATTAATCAAAGGTGAGAGGCACATGCCGGTAATTAAAGTACGTGAAAACGAGCCGTTCGACGTAGCTCTGCGTCGCTTCAAGCGTTCCTGCGAAAAAGCAGGTGTTCTGGCGGAAGTTCGTCGTCGTGAGTTCTa < 2:619444/143‑1 (MQ=255) aaaCCGAATTCATCAAAGTGGAGGGGCACATGCCGGTAATTAAAGTACGGGAAAACGAGCCGGTCGACGGAGCTCTGCGTCGCTTCAAGCGTTCCTGCGAGAAAGCAGGTGTTCTGGCGGAAGTTCGTCGTCGTGAGTTCTAt < 2:308141/143‑1 (MQ=255) cccGACTTAAACAAAGGCGACAGGCACATGCCGGTAATTAAAGTACGTGAAAACGAGCCGTTCGACGTAGCTCTGCGTCGCTTCAAGCGTTCCTGCGAAAAAGCAGGTGTTCTGGCGGAAGTTCGTCGTCGTGAGTTCTATGa < 1:11039/142‑1 (MQ=255) aCCGAATTAATCAAAGGGGAGAGGCCCATGCCGGTAATTAAAGTACGTGAAAACGAGCCGTTCGACGGAGCTCTGCGTCGCTGCAAGCGTTTCTGCGAAAAAGCAGGTGTTCTGGCGGAAGTTCGTCGTCGTGAGTTCTATGa < 1:39353/143‑1 (MQ=255) aaTGAATCCAAGGTGAGAGGCACATGCCGGTAATTAAAGTACGTGAAAACGAGCCGTTCGACGTAGCTCTGCGTCGCTTCAAGCGTTCCTGCGAAAAAGCAGGTGTTCTGGCGGAAGTGCGTCGTCGTGAGTTGTATGagaaa > 2:504463/1‑143 (MQ=255) aTTAATCAAAGGTGAGAGGCACATGCCGGTAATTAAAGTACGTGAAAACGAGCCGTTCGACGTAGCTCTGCGTCGCTTCAAGCGTTCCTGCGAAAAAGCAGGTGTTCTGGCGGAAGTTCGTCGTCGTGAGTTCTATGAAAAAc < 1:750444/143‑1 (MQ=255) agGCACCTGCCGGTAATTAAAGTACGTGAAAACGAGCCGTTCGACGTAGCTCTGCGTCGCTTCATGCGTTCCTGCGAAAAAGCAGGTGTTCTGGCGGAAGTTCGTCGTCGTGAGTTCTATGAAAAACCGACTACCGAACGTaa < 2:734203/143‑1 (MQ=255) aTGCCGGTAATTAAAGTACGTGAAAACGAGCCGTTCGACGTAGCTATGCGTCGCTTCAAGCGTTCCTGCGAAAAAGCAGGTGTTCTGGCGGAAGTTCGTCGTCGTGAGTTCTATGAAAAACCGACTACCGAACGTAAgcgcgc < 2:692775/143‑1 (MQ=255) aTTAAAGTACGTGAAAACGAGCCGTTCGACGTAGCTCTGCGTCGCTTCATGCGTTCCTGCGAAAATGCAGGAGTTCTGGCGTAAGGTCGTCGTCGTGAGTTCTATGAAAAACCGACTACCGAACGTACGCGCGCTAAAGCtta > 2:445262/1‑142 (MQ=255) aTTAAAGTACGTGAAAACGAGCCGTTCGACGTAGCTCTGCGTCGCTTCAAGCGTTCCTGCGAAAAAGCAGGTGTTCTGGCGGAAGTTCGTCGTCGTGAGTTCTATGAAAAACCGACTACCGAACGTAAGCGCGCTAAAGCTTc < 2:654957/143‑1 (MQ=255) tGAAAACGAGCCGTTCGACGTAGCTCTGCGTCGCTTCATGCGTTCCTGCGAAAAAGCAGGTGGTCTGGAGGAAGTTCGTCGTCGTGAGTTCTATGAAAAACCGACTACCGAACGTAAGCGCGCTAAAGCTTCTGCAGTGaagc > 1:39405/1‑141 (MQ=255) tGAAAACGAGCCGTTCGACGTAGCTCTGCGTCGCTTCAAGCGTTCCTGCGACAAAGCAGGTGTTCTGGCGGAAGTTCGTCGTCGTGAGTTCTATGAAAAACCGACTACCGCACGTAAGCGCGCTAAAGCTTCTGCAGTGaagc > 2:749169/1‑141 (MQ=255) ttCGACGTAGCTCTGCGTCGCTTCAAGCGTTCCTGCGAAAAAGCAGGTGTTCTGGCGGAAGTTCGTCGTCGTGAGTTCTATGAAAAACCGACTACCGAACGTAAGCGCGCGAAAGCTTCTGCAGTGAAACGACACGCgaagaa > 1:141778/1‑143 (MQ=255) gtggcGTAGCTCTGTGTAGCTTCAAGCGGTCCTGCGAAAAAGCAGCTGTTCTGGCGGAAGTTTGTCGTCGGGAGTCCTATGAAAGACCGACTACCGAACGTAAGCGCTCTAAAGCTTCTGCAGTGAAACGTCACGCGAAGaaa < 2:633866/139‑1 (MQ=255) aactcgGCGTCGCTTCAAGCGTTCCTGCGAAAAAGCAGGTGTTCTGGCGGAAGTTCGTCGTCGTGAGTTCTATGAAAAACCGACTACCGAACGTAAGCGCGCTAAAGCTTCTGCAGTGAAACGTCACGCGAAGAAACTGGCTc < 1:725718/141‑1 (MQ=255) aGCGCTGCGGCGCTTCTAGCGTTCCTGCGAAAAATGAGGTGTTCTGGCGGAAGTTCGTCGTCGTGAGTTCTATGAAAAACCGACTACCGATCGTAAGCGCGCTAAAGCTTCTGCAGTGAAACGTCACGCGAAGAAACTGGCTc < 2:570384/143‑1 (MQ=255) ctctGCGTCGCTTCATGCGTTCCTGCGAAAAAGCAGATGTTCTGGCGGAAGTTCGCCGGCGAGAGGTCTATGAAAAACCGACTACCGAACGTACGGGCGCTACAGCTGCTGCAGTGAAACGTGACGCGAAGAAACTggagaac > 2:78608/1‑138 (MQ=255) cGTCGCTTCATGCGTTCCTGCGAAAAAGCAGGTGTTCTGGCGGAAGTTCGTCGTCGTGAGTTCTATGAAAAACCGACTACCGAACGTAAGCGCGCTAAAGCTTCTGCAGTGAAACGTCACGCGAAGAAACTGGCTCGCGaaaa < 1:492633/143‑1 (MQ=255) cGTCGCGTCAAGCGTTCCTGCGAAAAAGCAGGTGTTCTGGCGGAAGTTCGTCGTCGTGAGTTCTATGAAAAACCGACTACCGAACGTAAGCGCGCTAAAGCTTCTGCAGTGAAACGTCACGCGAAGAAACTGGCTCGCGaaaa < 2:740502/143‑1 (MQ=255) tCAAGCGTTCCTGCGAAAAAGCAGGTGTTCTGGCGGAAGTTCGTCGTCGTGAGTTCTGTGAAAAACCGACTACCGAACGTAAGCGCGCTAAAGCTTCTGCAGTGAACCGTCACGCGAAGCAACTGGCTCGCGAAAaagcacgc > 2:603098/1‑143 (MQ=255) tCAAGCGTGCCTGCGAAAAAGCAGGTGGTCTGGCGGCAGTTCGTCGTCGTGAGTTCTATGAAAAACCGACTACCGAACGTAAGCGCGCTAAAGCTTCTGCAGTGAAACGTCACGCGAAGAAACTGGCTCGCGAAAacgcacgc < 2:732289/143‑1 (MQ=255) | AATTGCAGTAAAATTCCGCACCATTTTGAAATAAGCTGGCGTTGATGCCAGCGGCAAACCGAATTAATCAAAGGTGAGAGGCACATGCCGGTAATTAAAGTACGTGAAAACGAGCCGTTCGACGTAGCTCTGCGTCGCTTCAAGCGTTCCTGCGAAAAAGCAGGTGTTCTGGCGGAAGTTCGTCGTCGTGAGTTCTATGAAAAACCGACTACCGAACGTAAGCGCGCTAAAGCTTCTGCAGTGAAACGTCACGCGAAGAAACTGGCTCGCGAAAACGCACGC > NC_000913/3210697‑3210978 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 8 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |