| Predicted mutation | |||||||
|---|---|---|---|---|---|---|---|
| evidence | seq id | position | mutation | freq | annotation | gene | description |
| RA | minE | 1,530,171 | G→C | 100% | V214V (GTC→GTG) | yfcA ← | conserved inner membrane protein |
| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | minE | 1,530,171 | 0 | G | C | 83.8% | 88.2 / 8.9 | 43 | V214V (GTC→GTG) | yfcA | conserved inner membrane protein |
| Reads supporting (aligned to +/- strand): ref base G (0/7); new base C (36/0); total (36/7) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 3.10e-08 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 | |||||||||||
| Rejected as polymorphism: Variant not supported by required number of reads on each strand. | |||||||||||
TGCTCAACACCAGTCGTGAACCCATGCGCGCGCCGAGGAACTGCCCCACCAGCATGACAAAGCCCGTCGCCCAAATCACTTTGCCGCCGAGAATAAACAGCAGTA > minE/1530116‑1530220 | tGCTCAACACCAGTCTTGACCCCCTGCGCGCGCCGAGGAACTGTCCTACCAACATCACAAATcctg > 1:1695944/1‑64 (MQ=255)tGCTCAACACCAGTCTTGACCCCATGCGCGCGCCGAGGAACTGTCCTACCAACATCACAAATcctg > 1:110387/1‑64 (MQ=255)tGCTCAACACCAGTCTTGACCCCATGCGCGCGCCGAGGAACTGTCCTACCAACATCACAAATcctg > 1:2714545/1‑64 (MQ=255)tGCTCAACACCAGTCTTGACCCCATGCGCGCGCCGAGGAACTGTCCTACCAACATCACAAATcctg > 1:2611841/1‑64 (MQ=255)tGCTCAACACCAGTCTTGACCCCATGCGCGCGCCGAGGAACTGTCCTACCAACATCACAAATcctg > 1:2566363/1‑64 (MQ=255)tGCTCAACACCAGTCTTGACCCCATGCGCGCGCCGAGGAACTGTCCTACCAACATCACAAATcctg > 1:444492/1‑64 (MQ=255)tGCTCAACACCAGTCTTGACCCCATGCGCGCGCCGAGGAACTGTCCTACCAACATCACAAATcctg > 1:2424449/1‑64 (MQ=255)tGCTCAACACCAGTCTTGACCCCATGCGCGCGCCGAGGAACTGTCCTACCAACATCACAAATcctg > 1:2278113/1‑64 (MQ=255)tGCTCAACACCAGTCTTGACCCCATGCGCGCGCCGAGGAACTGTCCTACCAACATCACAAATcctg > 1:2230380/1‑64 (MQ=255)tGCTCAACACCAGTCTTGACCCCATGCGCGCGCCGAGGAACTGTCCTACCAACATCACAAATcctg > 1:2193400/1‑64 (MQ=255)tGCTCAACACCAGTCTTGACCCCATGCGCGCGCCGAGGAACTGTCCTACCAACATCACAAATcctg > 1:1935997/1‑64 (MQ=255)tGCTCAACACCAGTCTTGACCCCATGCGCGCGCCGAGGAACTGTCCTACCAACATCACAAATcctg > 1:1764979/1‑64 (MQ=255)tGCTCAACACCAGTCTTGACCCCATGCGCGCGCCGAGGAACTGTCCTACCAACATCACAAATcctg > 1:1384973/1‑64 (MQ=255)tGCTCAACACCAGTCTTGACCCCATGCGCGCGCCGAGGAACTGTCCTACCAACATCACAAATcctg > 1:132774/1‑64 (MQ=255)tGCTCAACACCAGTCTTGACCCCATGCGCGCGCCGAGGAACTGTCCTACCAACATCACAAATcctg > 1:12825/1‑64 (MQ=255)tGCTCAACACCAGTCTTGACCCCATGCGCGCGCCGAGGAACTGTCCTACCAACATCACAAATcctg > 1:1164705/1‑64 (MQ=255)tGCTCAACACCAGTCTTG‑ACCCATGCGCGCGCCGAGGAACTGTCCGACCAACATCACAAATcctg > 1:3256294/1‑63 (MQ=255) ccAGTCGTGAACCCATGCGCGCGCCGAGGAACTGCCCCACCAGCATGACAAATCCCGTCGCCCAAATCa < 1:1342324/69‑1 (MQ=255) ccAGTCGTGAACCCATGCGCGCGCCGAGGAACTGCCCCACCAGCATGACAAAGCCCGTCGCCCAAATCa < 1:702636/69‑1 (MQ=255) ccAGTCGTGAACCCATGCGCGCGCCGAGGAACTGCCCCACCAGCATGACAAAGCCCGTCGCCCAAATCa < 1:552087/69‑1 (MQ=255) ccAGTCGTGAACCCATGCGCGCGCCGAGGAACTGCCCCACCAGCATGACAAAGCCCGTCGCCCAAATCa < 1:1931418/69‑1 (MQ=255) ccAGTCGTGAACCCATGCGCGCGCCGAGGAACTGCCCCACCAGCATGACAAAGCCCGTCGCCCAAATCa < 1:2464927/69‑1 (MQ=255) ccAGTCGTGAACCCATGCGCGCGCCGAGGAACTGCCCCACCAGCATGACAAAGCCCGTCGCCCAAATCa < 1:3414275/69‑1 (MQ=255) ccAGTCGTGAACCCATGCGCGCGCCGAGGAACTGCCCCACCAGCATGACAAAGCCCGTCGCCCAAATCa < 1:361465/69‑1 (MQ=255) ggAACTGTCCTACCAACATCACAAATCCTGTCGCCCAAATCACTTTGCCCCCgaga > 1:2883605/1‑56 (MQ=255) ggAACTGTCCTACCAACATCACAAATCCTGTCGCCCAAATCACTTTGCCACCGAGAATAAAcagcagca > 1:3342856/1‑67 (MQ=255) ggAACTGTCCTACCAACATCACAAATCCTGTCGCCCAAATCACTTTGCCACCGAGAATAAAcagcagca > 1:348363/1‑67 (MQ=255) ggAACTGTCCTACCAACATCACAAATCCTGTCGCCCAAATCACTTTGCCACCGAGAATAAAcagcagca > 1:995878/1‑67 (MQ=255) ggAACTGTCCTACCAACATCACAAATCCTGTCGCCCAAATCACTTTGCCACCGAGAATAAAcagcagca > 1:39789/1‑67 (MQ=255) ggAACTGTCCTACCAACATCACAAATCCTGTCGCCCAAATCACTTTGCCACCGAGAATAAAcagcagca > 1:1045303/1‑67 (MQ=255) ggAACTGTCCTACCAACATCACAAATCCTGTCGCCCAAATCACTTTGCCACCGAGAATAAAcagcagca > 1:883012/1‑67 (MQ=255) ggAACTGTCCTACCAACATCACAAATCCTGTCGCCCAAATCACTTTGCCACCGAGAATAAAcagcagca > 1:697596/1‑67 (MQ=255) ggAACTGTCCTACCAACATCACAAATCCTGTCGCCCAAATCACTTTGCCACCGAGAATAAAcagcagca > 1:3289835/1‑67 (MQ=255) ggAACTGTCCTACCAACATCACAAATCCTGTCGCCCAAATCACTTTGCCACCGAGAATAAAcagcagca > 1:3226122/1‑67 (MQ=255) ggAACTGTCCTACCAACATCACAAATCCTGTCGCCCAAATCACTTTGCCACCGAGAATAAAcagcagca > 1:3055395/1‑67 (MQ=255) ggAACTGTCCTACCAACATCACAAATCCTGTCGCCCAAATCACTTTGCCACCGAGAATAAAcagcagca > 1:2884694/1‑67 (MQ=255) ggAACTGTCCTACCAACATCACAAATCCTGTCGCCCAAATCACTTTGCCACCGAGAATAAAcagcagca > 1:2638091/1‑67 (MQ=255) ggAACTGTCCTACCAACATCACAAATCCTGTCGCCCAAATCACTTTGCCACCGAGAATAAAcagcagca > 1:2500953/1‑67 (MQ=255) ggAACTGTCCTACCAACATCACAAATCCTGTCGCCCAAATCACTTTGCCACCGAGAATAAAcagcagca > 1:1926696/1‑67 (MQ=255) ggAACTGTCCTACCAACATCACAAATCCTGTCGCCCAAATCACTTTGCCACCGAGAATAAAcagcagca > 1:1905094/1‑67 (MQ=255) ggAACTGTCCTACCAACATCACAAATCCTGTCGCCCAAATCACTTTGCCACCGAGAATAAAcagcagca > 1:1896842/1‑67 (MQ=255) ggAACTGTCCTACCAACATCACAAATCCTGTCGCCCAAATCACTTTGCCACCGAGAATAAAcagcagca > 1:1732022/1‑67 (MQ=255) ggAACTGTCCTACCAACATCACAAATCCTGTCGCCCAAATCACTTTGCCACCGAGAATAAAcagcagca > 1:1162714/1‑67 (MQ=255) | TGCTCAACACCAGTCGTGAACCCATGCGCGCGCCGAGGAACTGCCCCACCAGCATGACAAAGCCCGTCGCCCAAATCACTTTGCCGCCGAGAATAAACAGCAGTA > minE/1530116‑1530220 |
| 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: ‑ |