Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | minE | 2,187,308 | T→C | 34.6% | P183P (CCA→CCG) | nanT ← | sialic acid transporter |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | minE | 2,187,308 | 0 | T | C | 34.6% | 22.3 / 23.1 | 26 | P183P (CCA→CCG) | nanT | sialic acid transporter |
Reads supporting (aligned to +/- strand): ref base T (14/3); new base C (4/5); total (18/8) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 7.81e-02 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-01 |
CGTGTTTCTCTTTCCAGTCTTCCGCTTCCGGGATGTTTTTACGCAGCCAGAGAGCAAAGATGATTGGCAAAATGCCGATAAAGAACAGCGCACGCCAGCCCCAGACCGGAACCACCAGGCTATAGAC > minE/2187244‑2187370 | cGTGTTTCTCTTTCCAGTCTTCCGCTTCCGGGATGTTTTTACGCAGCCAGAGAGCAAAGATGATTGGCaaa < 1:466545/71‑1 (MQ=255) cAGTCTTCCGCTTCCGGGATGTTTTTACGCAGCCAGAGAGCAAAGATGATTGGCAAAATGcc > 1:93162/1‑62 (MQ=255) cttccgcttccgGGATGTTTTTACGCAGCCAGAGAGCAAAGATGATTGGCAAAATGc > 1:781662/1‑57 (MQ=255) ccgcttccgGGATGTTTTTACGCAGCCAGAGAGCAAAGATGATCGGCaaaa > 1:1096012/1‑51 (MQ=255) gggATGTTTTTACGCAGCCAGAGAGCAAAGATGATTGGCAAAATGCCGATAAAGAACAGCGCa > 1:1884972/1‑63 (MQ=255) aTGTTTTTACGCAGCCAGAGAGCAAAGATGATTGGCaa < 1:2112608/38‑1 (MQ=255) aTGTTTTTACGCAGCCAGAGAGCAAAGATGATCGGCAAAATGCCGATAAAGAACAGCGCACGCCAGccc > 1:232810/1‑69 (MQ=255) gTTTTTACGCAGCCAGAGAGCAAAGATGATTGGCAAAATGCCGATAAAGAACAGCGCACGCCAGCCCCAGa > 1:413820/1‑71 (MQ=255) ttttACGCAGCCAGAGAGCAAAGATGATTGGCaaaa > 1:55073/1‑36 (MQ=255) tttACGCAGCCAGAGAGCAAAGATGATCGGCAAAATGCCGATAAAGAACAGCGCACGCCAGCCCCAGACCg > 1:2089105/1‑71 (MQ=255) gCCAGAGAGCAAAGATGATCGGCAAAATGCCGATAAAGAACAGCGCACGc < 1:142262/50‑1 (MQ=255) agagCAAAGATGATCGGCAAAATGCCGATAAAGAACAGCGCACGCCAGCCCCAGACCGGAaccac < 1:2143354/65‑1 (MQ=255) agagCAAAGATGATCGGCAAAATGCCGATAAAGAACAGCGCACGCCAGCCCCAGACCGGAaccac < 1:431798/65‑1 (MQ=255) gagCAAAGATGATTGGCAAAATGCCGATAAAGAACAGCGCACGCCAGCCCCAGACCGGAACCACCAGGCta < 1:666745/71‑1 (MQ=255) gagCAAAGATGATCGGCAAAATGCCGATAAAGAACAGCGCACGCCAGCCCCAGACCGGAACCACCAGGCta < 1:221450/71‑1 (MQ=255) cAAAGATGATCGGCAAAATGCCGATAAAGAACAGCGCACGCCAGCCCCAGACCGGAACCACCAg < 1:590402/64‑1 (MQ=255) aaGATGATTGGCAAAATGCCGATAAAGAACAGCGCACGCCAGCCCCAGACCGGAACCACCAgcc > 1:999982/1‑62 (MQ=255) aaGATGATTGGCAAAATGCCGATAAAGAACAGCGCACGCCAGCCCCAGACCGGAACCACCAGGc > 1:592097/1‑64 (MQ=255) aaGATGATTGGCAAAATGCCGATAAAGAACAGCGCACGCCAGCCCCAGACCGGAACCACCAGGc > 1:652922/1‑64 (MQ=255) aaGATGATTGGCAAAATGCCGATAAAGAACAGCGCACGCCAGCCCCAGACCGGAACCACCAGGc > 1:136974/1‑64 (MQ=255) aaGATGATCGGCAAAATGCCGATAAAGAACAGCGCACGCCAGCCCCAGACCGGAACCACCAGGCTATAGAc > 1:2095945/1‑71 (MQ=255) gatgatTGGCAAAATGCCGATAAAGAACAGCGCACGCCAGCCCCAGACCGGAac > 1:129231/1‑54 (MQ=255) tgatTGGCAAAATGCCGATAAAGAACAGCGCACGc > 1:1408791/1‑35 (MQ=255) gatTGGCAAAATGCCGATAAAGAACAGCGCACGCCAGCCCCa > 1:1616168/1‑42 (MQ=255) gatTGGCAAAATGCCGATAAAGAACAGCGCACGCCAGCCCCa > 1:1291602/1‑42 (MQ=255) atTGGCAAAATGCCGATAAAGAACAGCGCACGCCAGCCCCa > 1:543916/1‑41 (MQ=255) | CGTGTTTCTCTTTCCAGTCTTCCGCTTCCGGGATGTTTTTACGCAGCCAGAGAGCAAAGATGATTGGCAAAATGCCGATAAAGAACAGCGCACGCCAGCCCCAGACCGGAACCACCAGGCTATAGAC > minE/2187244‑2187370 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |