| Predicted mutation | |||||||
|---|---|---|---|---|---|---|---|
| evidence | seq id | position | mutation | freq | annotation | gene | description |
| RA | minE | 749,764 | G→A | 57.1% | intergenic (‑135/+9) | ycfS ← / ← mfd | conserved hypothetical protein/transcription‑repair coupling factor |
| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | minE | 749,764 | 0 | G | A | 57.1% | 10.5 / 42.2 | 35 | intergenic (‑135/+9) | ycfS/mfd | conserved hypothetical protein/transcription‑repair coupling factor |
| Reads supporting (aligned to +/- strand): ref base G (9/6); new base A (8/12); total (17/18) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 3.15e-01 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 | |||||||||||
AGTGTCGGATAGTGCAGGGAAGTGCAAAGAATTTGTAAATGTTGCAGATGGGGGCGCAGAAACGCCCCCGATTTACCATTAAGCGATCGCGTTCTCTTCCAGTTCACGCATAAACTGGCGTACCCATTCGA > minE/749695‑749825 | aGTGTCGGATAGTTCAGGGAAGTGCAAAGAATTTGTAAATGTTGCAGATGGGGGCGCAGAAACGCCCCCGa < 1:2889961/71‑1 (MQ=255)aGTGTCGGATAGTGCAGGGAAGTGCAAAGAATTTGTAAATGTTGCAGATGGGGGCGCAGAAACGCCCCCg < 1:1565440/70‑1 (MQ=255) tgtCGGCTAGTGCAGGGAAGTGCAAAGAATTTGTAAATGTTGCAGATGGGGGCGCAGAAACGCCCCCAAtt < 1:1324106/71‑1 (MQ=255) tCGGATAGTGCAGGGAAGTGCAAAGAATTTGTAAATGTTGCAGATGGGGGCGCAGAAACGCCCCCAAttt < 1:222691/70‑1 (MQ=255) ggATAGTGCAGGGAAGTGCAAAGAATTTGTAAATGTTGCAGATGGGGGCGCAGAAACGCCCCCAATTTAcc < 1:2244896/71‑1 (MQ=255) aTAGTGCAGGGAAGTGCAAAGAATTTGTAAATGTTGCAGATGGGGGCGCAGAAACGCCCCCGATTTACCa < 1:2838895/70‑1 (MQ=255) aTAGTGCAGGGAAGTGCAAAGAATTTGTAAATGTTGCAGATGGGGGCGCAGAAACGCCCCCAATTTACCAt < 1:3798528/71‑1 (MQ=255) tAGTGCAGGGAAGTGCAAAGAATTTGTAAATGTTGCAGATGGGGGCGCAGAAACGCCCCCAATTTACCa > 1:3124085/1‑69 (MQ=255) cAGGGAAGTGCAAAGAATTTGTAAATGTTGCAGATGGGGGCGCAGAAACGCCCCCAATTTACCATTAAGCg < 1:2086022/71‑1 (MQ=255) aGGGAAGTGCAAAGAATTTGTAAATGTTGCAGATGGGGGCGCAGAAACGCCCCCGAtt < 1:1685333/58‑1 (MQ=255) aaGTGCAAAGAATTTGTAAATGTTGCAGATGGGGGCGCAGAAACGCCCCCAATTTACCATTAAGCGATcgc > 1:889284/1‑71 (MQ=255) aaGTGCAAAGAATTTGTAAATGTTGCAGATGGGGGCGCAGAAACGCCCCCAATTTACCATTAAGCGATcgc > 1:901457/1‑71 (MQ=255) aaGTGCAAAGAATTTGTAAATGTTGCAGATGGGGGCGCAGAAACGCCCCCAATTTACCATTAAGCGATc < 1:1828656/69‑1 (MQ=255) tGCAAAGAATTTGTAAATGTTGCAGATGGGGGCGCAGAAACGCCCCCGATTTACCATTAAGCGATCGCGtt > 1:718054/1‑71 (MQ=255) aaaGAATTTGTAAATGTTGCAGATGGGGGCGCAGAAACGCCCCCGATTTACCATTa < 1:3507109/56‑1 (MQ=255) aaGAATTTGTAAATGTTGCAGATGGGGGCGCAGAACCGCCCCCAATTTACCATTAAGCGATCGCGTTctct > 1:3521812/1‑71 (MQ=255) aaGAATTTGTAAATGTTGCAGATGGGGGCGCAGAAACGCCCCCGATTTACCATTAAGCGATCGCGTtctc > 1:52514/1‑70 (MQ=255) gAATTTGTAAATGTTGCAGATGGGGGCGCAGAAACGCCCCCGAttt < 1:360426/46‑1 (MQ=255) aaTTTGTAAATGTTGCAGATGGGGGCGCAGAAACGCCCCCAATTTACCATTAAGCGATCGCGTTCTCTTcc < 1:3467912/71‑1 (MQ=255) gTAAATGTTGCAGATGGGGGCGCAGAAACGCCCCCAATTTACCATTAAGCGATCGCGTTCTCTTCCAGtt < 1:2135004/70‑1 (MQ=255) aaTGTTGCAGATGGGGGCGCAGAAACGCCCCCAATTTACCATTAAGCGATCGCGTTCTCTTCCAGTTCAc < 1:3455584/70‑1 (MQ=255) tGTTGCAGATGGGGGCGCAGAAACGCCCCCGATTTACCATTAAGCGATCGCGTTCTCTTCCAGTTCACGCa < 1:1077153/71‑1 (MQ=255) ttGCAGATGGGGGCGCAGAAACGCCCCCAATTTACCATTAAGCGATCGCGTTCTCTTCCAGTTCACGCATa < 1:428805/71‑1 (MQ=255) gCAGATGGGGGCGCAGAAACGCCCCCGATTTACCATTAAGCGATCGCGTTCTCTTCCAGTTCACGCATaaa > 1:84/1‑71 (MQ=255) gggggCGCAGAAACGCCCCCAATTTACCATTAAGCGATCGCGTTCTCTTCCAGTTCACGCATAAACTGGCg > 1:2268152/1‑71 (MQ=255) gAAACGCCCCCGATTTACCATTAAGCGATCGCGTTCTCTTCCAg > 1:1216617/1‑44 (MQ=255) gAAACGCCCCCGATTTACCATTAAGCGATCGCGTTCTCTTCCAGTTCACGCATAAACTGGCGTa > 1:3741677/1‑64 (MQ=255) gAAACGCCCCCGATTTACCATTAAGCGATCGCGTTCTCTTCCAGTTCACGCATAAACTGGCGTACCCAtt > 1:645491/1‑70 (MQ=255) aCGCCCCCGATTTACCATTAAGCGATCGCGTTCTCTTCCAGTTCACGCATAAACTGGCGTACCCATTCGa > 1:3750610/1‑70 (MQ=255) aCGCCCCCAATTTACCATTAAGCGATCGCGTTCTCTTCCAGTTCACGCATAAACTgg > 1:3009881/1‑57 (MQ=255) aCGCCCCCAATTTACCATTAAGCGATCGCGTTCTCTTCCAGTTCACGCATAAACTGGCGTACCCATTCGa > 1:1573760/1‑70 (MQ=255) aCGCCCCCAATTTACCATTAAGCGATCGCGTTCTCTTCCAGTTCACGCATAAACTGGCGTACCCATTCGa > 1:149728/1‑70 (MQ=255) cccccGATTTACCATTAAGCGATCGCGTTCTCTTCCAGTTCACGCATaaa > 1:1357541/1‑50 (MQ=255) ccccAATTTACCATTAAGCGATCGCGTTCTCTTCCAGTTCACGGATa < 1:670159/47‑1 (MQ=39) ccccAATTTACCATTAAGCGATCGCGTTCTCTTCCAGTTCACGCATa < 1:2276993/47‑1 (MQ=255) ccGATTTACCATTAAGCGATCGCGTTCTCTTCCAGTTCACGCATAAACTGGCGTa > 1:2990847/1‑55 (MQ=255) | AGTGTCGGATAGTGCAGGGAAGTGCAAAGAATTTGTAAATGTTGCAGATGGGGGCGCAGAAACGCCCCCGATTTACCATTAAGCGATCGCGTTCTCTTCCAGTTCACGCATAAACTGGCGTACCCATTCGA > minE/749695‑749825 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |