| Predicted mutation | |||||||
|---|---|---|---|---|---|---|---|
| evidence | seq id | position | mutation | freq | annotation | gene | description |
| RA | minE | 1,707,359 | A→G | 100% | L6L (TTA→CTA) | yphB ← | conserved hypothetical protein |
| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | minE | 1,707,359 | 0 | A | G | 82.5% | 51.7 / 3.1 | 34 | L6L (TTA→CTA) | yphB | conserved hypothetical protein |
| Reads supporting (aligned to +/- strand): ref base A (2/4); new base G (28/0); total (30/4) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 3.23e-04 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.93e-01 | |||||||||||
| Rejected as polymorphism: Variant not supported by required number of reads on each strand. | |||||||||||
GCGCCAGAAGCCTTCGATTACGCCGCCCTGGTCAGAAACGTCCAGCTTTAGCGACCCGTGGGATAAGGTATAGATGGTCATCCGCTAACTCCTTAATCCGGGAAGTTAATCACAACTTTCCCGCATTTGCCGCT > minE/1707294‑1707427 | gcgcCAGAAGCCTTCGATTACGCCGCCCTGGTCAGAAACGTCCAGCTTTAGCGACCCGTGGGATAAgg < 1:930669/68‑1 (MQ=255) aGAAGCCTTCGATTACGCCGCCCTGGTCAGAAACGTCCAGCTTTAGCGACCCGTGGGATAAGGTATAGa < 1:3055133/69‑1 (MQ=255) gATTACGCCGCCCTGGTCAGAAACGTCCAGCTTTAGCGCCCCGTGGGATAAGGTATAGATGGTCATcc > 1:2442061/1‑68 (MQ=255) gATTACGCCGCCCTGGTCAGAAACGTCCAGCTTTAGCGACCCGTGGGATAAGGTATAGATGGTCATcc > 1:1288022/1‑68 (MQ=255) aGAAACGTCCAGCTTTAGCGACCCGTGGGATAAGGTATAGa < 1:1474089/41‑1 (MQ=255) cTTTAGCGACCCGTGGGATAAGGTATAGATGGTCATCCGCTAACTCCTTAATCCGGGAAGTTAATcaca < 1:1327975/69‑1 (MQ=255) ttAGCAGCCCGTGAGATAGGGTATAGATGGTCATCCGCAGACTCCTTAATCCGGGAAATTAATCACaa > 1:735501/1‑68 (MQ=255) ttAGCAGCCCGTGAGATAGGGTATAGATGGTCATCCGCAGACTCCTTAATCCGGGAAATTAATCACAAc > 1:99819/1‑69 (MQ=255) ttAGCAGCCCGTGAGATAGGGTATAGATGGTCATCCGCAGACTCCTTAATCCGGGAAATTAATCACAAc > 1:697598/1‑69 (MQ=255) ttAGCAGCCCGTGAGATAGGGTATAGATGGTCATCCGCAGACTCCTTAATCCGGGAAATTAATCACAAc > 1:593360/1‑69 (MQ=255) ttAGCAGCCCGTGAGATAGGGTATAGATGGTCATCCGCAGACTCCTTAATCCGGGAAATTAATCACAAc > 1:3270142/1‑69 (MQ=255) ttAGCAGCCCGTGAGATAGGGTATAGATGGTCATCCGCAGACTCCTTAATCCGGGAAATTAATCACAAc > 1:2991002/1‑69 (MQ=255) ttAGCAGCCCGTGAGATAGGGTATAGATGGTCATCCGCAGACTCCTTAATCCGGGAAATTAATCACAAc > 1:2870759/1‑69 (MQ=255) ttAGCAGCCCGTGAGATAGGGTATAGATGGTCATCCGCAGACTCCTTAATCCGGGAAATTAATCACAAc > 1:2462973/1‑69 (MQ=255) ttAGCAGCCCGTGAGATAGGGTATAGATGGTCATCCGCAGACTCCTTAATCCGGGAAATTAATCACAAc > 1:2387455/1‑69 (MQ=255) ttAGCAGCCCGTGAGATAGGGTATAGATGGTCATCCGCAGACTCCTTAATCCGGGAAATTAATCACAAc > 1:1641441/1‑69 (MQ=255) ttAGCAGCCCGTGAGATAGGGTATAGATGGTCATCCGCAGACTCCTTAATCCGGGAAATTAATCACAAc > 1:2049920/1‑69 (MQ=255) ttAGCAGCCCGTGAGATAGGGTATAGATGGTCATCCGCAGACTCCTTAATCCGGGAAATTAATCACAAc > 1:2091430/1‑69 (MQ=255) ttAGCAGCCCGTGAGATAGGGTATAGATGGTCATCCGCAGACTCCTTAATCCGGGAAATTAATCACAAc > 1:2150600/1‑69 (MQ=255) ttAGCAGCCCGTGAGATAGGGTATAGATGGTCATCCGCAGACTCCTTAATCCGGGAAATTAATCACAAc > 1:2380195/1‑69 (MQ=255) ttAGCAGCCCGTGAGAAAGGGTATAGATGGTCATCCGCAGACTCCTTAATCCGGGAAATTAATCACAAc > 1:2193808/1‑69 (MQ=255) gagATAGGGTATAGATGGTCATCCGCAGACTCCTTAATCCGGGAAATTAAt > 1:2848917/3‑51 (MQ=255) gagATAGGGTATAGATGGTCATCCGCAGACTCCTTAATCCGGGAAATTAATCACAACTTTCCCGCAttt > 1:2868707/3‑69 (MQ=255) gagATAGGGTATAGATGGTCATCCGCAGACTCCTTAATCCGGGAAATTAATCACAACTTTCCCGCAttt > 1:3174924/3‑69 (MQ=255) gagATAGGGTATAGATGGTCATCCGCAGACTCCTTAATCCGGGAAATTAATCACAACTTTCCCGCAttt > 1:35569/3‑69 (MQ=255) gagATAGGGTATAGATGGTCATCCGCAGACTCCTTAATCCGGGAAATTAATCACAACTTTCCCGCAttt > 1:501808/3‑69 (MQ=255) gagATAGGGTATAGATGGTCATCCGCAGACTCCTTAATCCGGGAAATTAATCACAACTTTCCCGCAttt > 1:582570/3‑69 (MQ=255) gagATAGGGTATAGATGGTCATCCGCAGACTCCTTAATCCGGGAAATTAATCACAACTTTCCCGCAttt > 1:2618723/3‑69 (MQ=255) gagATAGGGTATAGATGGTCATCCGCAGACTCCTTAATCCGGGAAATTAATCACAACTTTCCCGCAttt > 1:2029221/3‑69 (MQ=255) gagATAGGGTATAGATGGTCATCCGCAGACTCCTTAATCCGGGAAATTAATCACAACTTTCCCGCAttt > 1:712269/3‑69 (MQ=255) gagATAGGGTATAGATGGTCATCCGCAGACTCCTTAATCCGGGAAATTAATCACAACTTTCCCGCAttt > 1:718526/3‑69 (MQ=255) gagATAGGGTATAGATGGTCATCCGCAGACTCCTTAATCCGGGAAATTAATCACAACTTTCCCGCAttt > 1:1429728/3‑69 (MQ=255) gagATAGGGTATAGATGGTCATCCGCAGACTCCTTAATCCGGGAAATTAATAACAACTTTCCCGCAttt > 1:30294/3‑69 (MQ=255) gagATAGGGTATAGATAGTCATCCGCAGACTCCTTAATCCGGGAAATTAATCACAACTTTCCCGCAt > 1:2664595/3‑67 (MQ=255) aGGTATAGATGGTCATCCGCTAACTCCTTAATCCGGGAAGTTAATCACAACTTTCCCGCATTTGCCGCt > 1:2451650/1‑69 (MQ=255) aGGTATAGATGGTCATCCGCTAACTCCTTAATCCGGGAAGTTAATCACAACTTTCCCGCATTTGCCGCt < 1:1487824/69‑1 (MQ=255) | GCGCCAGAAGCCTTCGATTACGCCGCCCTGGTCAGAAACGTCCAGCTTTAGCGACCCGTGGGATAAGGTATAGATGGTCATCCGCTAACTCCTTAATCCGGGAAGTTAATCACAACTTTCCCGCATTTGCCGCT > minE/1707294‑1707427 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 21 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |