| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | NC_000913 | 3,434,974 | 0 | T | A | 22.7% | 50.5 / 12.7 | 22 | I254N (ATT→AAT) | fmt | 10‑formyltetrahydrofolate:L‑methionyl‑tRNA(fMet) N‑formyltransferase |
| Reads supporting (aligned to +/- strand): ref base T (6/11); new base A (5/0); total (11/11) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 3.51e-02 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.52e-01 | |||||||||||
| Rejected as consensus: Frequency below/above cutoff threshold. | |||||||||||
GCCAATGAGCTGGCTGGAAATTGAAGGACAGCCGGTTAAAGTCTGGAAAGCATCGGTCATTGATACGGCAACCAACGCTG‑CACCAGGAACGATCCTTGAAGCCAACAAACAAGGC > NC_000913/3434915‑3435029 | gCCAATGAGCTGGCTGGAAATTGAAGGACAGCCGGTTAAAGTCTGGAAAGCATCGGTCAGTGATa > 1:1052153/1‑65 (MQ=255) cAATGAGCTGGCTGGAAATTGAAGGACAGCCGGTTAAAGTCTGGAAAGCATCGGTCATTGATACg < 1:1392940/65‑1 (MQ=255) aTGAGCTGGCTGGAAATTGAAGGACAGCCGGTTAAAGTCTGGAAAGCATCGGTCATTGATACGGc < 1:2391149/65‑1 (MQ=255) aTGAGCTGGCTGGAAATTGAAGGACAGCCGGTTAAAGTCTGGAAAGCATCGGTCATTGATACGGc < 1:2959624/65‑1 (MQ=255) aTGAGCTGGCTGGAAATTGAAGGACAGCCGGTTAAAGTCTGGAAAGCATCGGTCATTGATACGGc < 1:1404/65‑1 (MQ=255) tGAGCTGGCTGGAAATTGAAGGACAGCCGGTTAAAGTCTGGAAAGCATCGGTTAATGATACGGcg > 1:615520/1‑64 (MQ=255) gAGCTGGCTGGAAATTGAAGGACAGCCGGTTAAAGTCTGGAAAGCATCGGTTAATGATACGGcga > 1:256554/1‑63 (MQ=255) gAGCTGGCTGGAAATTGAAGGACAGCCGGTTAAAGTCTGGAAAGCATCGGTCATTGATACGGcaa < 1:1127375/65‑1 (MQ=255) tggctggAAATTGAAGGACAGCCGGTTAAAGTCTGGAAAGCATCGGTCATTGATACGGcaaccaa > 1:813137/1‑65 (MQ=255) ggctggAAATTGAAGGACAGCCGGTTAAAGTCTGGAAAGCATCGGTCATTGATACGGcaaccaac < 1:639165/65‑1 (MQ=255) ctggAAATTGAAGGACAGCCGGTTAAAGTCTGGAAAGCATCGGTCATTGATACGGCAACCAACGc > 1:2333686/1‑65 (MQ=255) aaTTGAAGGACAGCCGGTTAAAGTCTGGAAAGCATCGGTCAATGATACGGCGACCAACGCTG‑CAc > 1:1871867/1‑65 (MQ=255) aTTGAAGGACAGCCGGTTAAAGTCTGGAAAGCATCGGTCATTGATACGGCAACCAACGCTG‑CAcc < 1:2913920/65‑1 (MQ=255) gACAGCCGGTTAAAGTCTGGAAAGCATCGGTCATTGATACGGCAACCAACGCTG‑CACCAGGAACg < 1:378589/65‑1 (MQ=255) aGCCGGTTAAAGTCTGGAAAGCATCGGTCATTGATACGGCAACCAACGCTG‑CACCAGGAACGATc < 1:2139380/65‑1 (MQ=255) aaaGTCTGGAAAGCATCGGTCATTGATACGGCAACCAACGCTG‑CACCAGGAACGATCCTTGAAGc > 1:1392200/1‑65 (MQ=255) aGTCTGGAAAGCATCGGTTAATGATACGGCGACCACCGCGGTCCACACG‑ACGATCCTTGCAGCca > 1:1915750/1‑65 (MQ=255) aGTCTGGAAAGCATCGGTCATTGATACGGCAACCAACGCTG‑CACCAGGAACGATCCTTGAAGCca > 1:2665037/1‑65 (MQ=255) aGTCTGGAAAGCATCGGTCATTGATACGGCAACCAACGCTG‑CACCAGGAACGATCCTTGAAGCca > 1:2254786/1‑65 (MQ=255) aaaGCATCGGTCATTGATACGGCAACCAACGCTG‑CACCAGGAACGATCCTTGAAGCCaacaaaca > 1:1355866/1‑65 (MQ=255) aaGCATCGGTTAATGATACGGCGACCAACGAGG‑CACAAAGGCCGACCCTTGAAGCCAacaaacaa > 1:1015563/1‑65 (MQ=255) cATCGGTCATTGATACGGCAACCAACGCTG‑CACCAGGAACGATCCTTGAAGCCAACAAACAAGGc < 1:2924541/65‑1 (MQ=255) cATCGGTCATTGATACGGCAACCAACGCTG‑CACCAGGAACGATCCTTGAAGCCAACAAACAAGGc < 1:2795433/65‑1 (MQ=255) cATCGGTCATTGATACGGCAACCAACGCTG‑CACCAGGAACGATCCTTGAAGCCAACAAACAAGGc < 1:576565/65‑1 (MQ=255) | GCCAATGAGCTGGCTGGAAATTGAAGGACAGCCGGTTAAAGTCTGGAAAGCATCGGTCATTGATACGGCAACCAACGCTG‑CACCAGGAACGATCCTTGAAGCCAACAAACAAGGC > NC_000913/3434915‑3435029 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 21 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |