| Predicted mutation | |||||||
|---|---|---|---|---|---|---|---|
| evidence | seq id | position | mutation | freq | annotation | gene | description |
| RA | NC_000913 | 3,313,207 | G→A | 40.0% | intergenic (‑29/+135) | rbfA ← / ← infB | 30s ribosome binding factor/translation initiation factor IF‑2 |
| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | NC_000913 | 3,313,207 | 0 | G | A | 40.0% | 0.4 / 7.8 | 10 | intergenic (‑29/+135) | rbfA/infB | 30s ribosome binding factor/translation initiation factor IF‑2 |
| Reads supporting (aligned to +/- strand): ref base G (3/3); new base A (2/2); total (5/5) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 1.00e+00 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.18e-01 | |||||||||||
GCCCAGGCGAGGATCTTTAATTTCACGCTGCAGGATGAGAGCGATCTCTTTTTGCATTTCCTGCGCTACGCGCTGCGGGCGACCAAATTCTTTCGCCATAATAAATTCTCCTGACAAAAAAGGGGCTGTTAGCCCCTTTTTAAAATTAATTTCAGGTGGAAGGGCT > NC_000913/3313080‑3313245 | gCCCAGGCGAGGATCTTTAATTTCACGCTGCAGGATGAGAGCGATCTCTTTTTGCATTTCCTGCGCTACGCGCTGCGGGCGACCAAATTCTTTCGCCATAATAAATTCTCCTGACAAAAAAGGGGCTGTTAGc > 1:1427221/1‑133 (MQ=255)gCCCAGGCGAGGATCTTTAATTTCACGCTGCAGGATGAGAGCGATCTCTTTTTGCATTTCCTGCGCTACGCGCTGCGGGCGACCAAATTCTTTCGCCATAATAAATTCTCCTGACAAAAAAGGGGCTGTTAGc < 2:1427221/133‑1 (MQ=255) ccAGGCGAGGATCTTTAATTTCACGCTGCAGGATGAGAGCGATCTCTTTTTGCATTTCCTGCGCTACGCGCTGCGGGCGACCAAATTCTTTCGCCATAATAAATTCTCCTGACAAAAAAGGGGCTAACAGCCCCttt < 1:1302261/137‑1 (MQ=255) ccAGGCGAGGATCTTTAATTTCACGCTGCAGGATGAGAGCGATCTCTTTTTGCATTTCCTGCGCTACGCGCTGCGGGCGACCAAATTCTTTCGCCATAATAAATTCTCCTGACAAAAAAGGGGCTAACAGCCCCttt > 2:1302261/1‑137 (MQ=255) aCGCTGCAGGATGAGAGCGATCTCTTTTTGCATTTCCTGCGCTACGCGCTGCGGGCGACCAAATTCTTTCGCCATAATAAATTCTCCTGACAAAAAAGGGGCTAACAGCCCCtttt > 1:1264526/1‑116 (MQ=255) aCGCTGCAGGATGAGAGCGATCTCTTTTTGCATTTCCTGCGCTACGCGCTGCGGGCGACCAAATTCTTTCGCCATAATAAATTCTCCTGACAAAAAAGGGGCTAACAGCCCCtttt < 2:1264526/116‑1 (MQ=255) cTACGCGCTGCGGGCGACCAAATTCTTTCGCCATAATAAATTCTCCTGACAAAAAAGGGGCTGtt > 1:1282459/1‑65 (MQ=255) cTACGCGCTGCGGGCGACCAAATTCTTTCGCCATAATAAATTCTCCTGACAAAAAAGGGGCTGtt < 2:1282459/65‑1 (MQ=255) cATAATAAATTCTCCTGACAAAAAAGGGGCTGTTAGCCCCTTTTTAAAATTAATTTCAGGTGGAAGGGCt > 1:702373/1‑70 (MQ=255) cATAATAAATTCTCCTGACAAAAAAGGGGCTGTTAGCCCCTTTTTAAAATTAATTTCAGGTGGAAGGGCt < 2:702373/70‑1 (MQ=255) | GCCCAGGCGAGGATCTTTAATTTCACGCTGCAGGATGAGAGCGATCTCTTTTTGCATTTCCTGCGCTACGCGCTGCGGGCGACCAAATTCTTTCGCCATAATAAATTCTCCTGACAAAAAAGGGGCTGTTAGCCCCTTTTTAAAATTAATTTCAGGTGGAAGGGCT > NC_000913/3313080‑3313245 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |