Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | minE | 1,530,162 | C→T | 100% | V217V (GTG→GTA) | yfcA ← | conserved inner membrane protein |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | minE | 1,530,162 | 0 | C | T | 76.1% | 53.6 / 12.8 | 47 | V217V (GTG→GTA) | yfcA | conserved inner membrane protein |
Reads supporting (aligned to +/- strand): ref base C (2/9); new base T (35/0); total (38/9) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 4.99e-08 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.47e-01 | |||||||||||
Rejected as polymorphism: Variant not supported by required number of reads on each strand. |
ATCAGCTTTTGTCCTTTGCTCAACACCAGTCGTGAACCCATGCGCGCGCCGAGGAACTGCCCCACCAGCATGACAAAGCCCGTCGCCCAAATCACTTTGCCGCCGAGAATAAACAGCAGTA > minE/1530100‑1530220 | aTCAGCTTTTGTCCTTTGCTCAACACCAGTCGTGAACCCATGCGCGCGCCGAGGAACTGCCCCACCAGc < 1:3179328/69‑1 (MQ=255) tCAGCTTTTGTCCTTTGCTCAACACCAGTCGTGAACCCATGCGCGCGCCGAGGAACTGCCCCACCAGc > 1:1837521/1‑68 (MQ=255) tCAGCTTTTGTCCTTTGCTCAACACCAGTCGTGAACCCATGCGCGCGCCGAGGAACTGCCCCACCAGc > 1:2041217/1‑68 (MQ=255) tGCTCAACACCAGTCTTGACCCCCTGCGCGCGCCGAGGAACTGTCCTACCAACATCACAAATcctg > 1:1695944/1‑64 (MQ=255) tGCTCAACACCAGTCTTGACCCCATGCGCGCGCCGAGGAACTGTCctacca > 1:3209488/1‑51 (MQ=255) tGCTCAACACCAGTCTTGACCCCATGCGCGCGCCGAGGAACTGTCCTACCAACATCACAAATcctg > 1:110387/1‑64 (MQ=255) tGCTCAACACCAGTCTTGACCCCATGCGCGCGCCGAGGAACTGTCCTACCAACATCACAAATcctg > 1:444492/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: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:1164705/1‑64 (MQ=255) tGCTCAACACCAGTCTTGACCCCATGCGCGCGCCGAGGAACTGTCCTACCAACATCACAAATcctg > 1:1764979/1‑64 (MQ=255) tGCTCAACACCAGTCTTGACCCCATGCGCGCGCCGAGGAACTGTCCTACCAACATCACAAATcctg > 1:1384973/1‑64 (MQ=255) tGCTCAACACCAGTCTTGACCCCATGCGCGCGCCGAGGAACTGTCCTACCAACATCACAAATcctg > 1:12825/1‑64 (MQ=255) tGCTCAACACCAGTCTTGACCCCATGCGCGCGCCGAGGAACTGTCCTACCAACATCACAAATcctg > 1:132774/1‑64 (MQ=255) tGCTCAACACCAGTCTTG‑ACCCATGCGCGCGCCGAGGAACTGTCCGACCAACATCACAAATcctg > 1:3256294/1‑63 (MQ=255) ccAGTCGTGAACCCATGCGCGCGCCGAGGAACTGCCCCACCAGCATGACAAATCCCGTCGCCCAAATCa < 1:1342324/69‑1 (MQ=255) ccAGTCGTGAACCCATGCGCGCGCCGAGGAACTGCCCCACCAGCATGACAAAGCCCGTCGCCCAAATCa < 1:552087/69‑1 (MQ=255) ccAGTCGTGAACCCATGCGCGCGCCGAGGAACTGCCCCACCAGCATGACAAAGCCCGTCGCCCAAATCa < 1:361465/69‑1 (MQ=255) ccAGTCGTGAACCCATGCGCGCGCCGAGGAACTGCCCCACCAGCATGACAAAGCCCGTCGCCCAAATCa < 1:3414275/69‑1 (MQ=255) ccAGTCGTGAACCCATGCGCGCGCCGAGGAACTGCCCCACCAGCATGACAAAGCCCGTCGCCCAAATCa < 1:2464927/69‑1 (MQ=255) ccAGTCGTGAACCCATGCGCGCGCCGAGGAACTGCCCCACCAGCATGACAAAGCCCGTCGCCCAAATCa < 1:702636/69‑1 (MQ=255) ccAGTCGTGAACCCATGCGCGCGCCGAGGAACTGCCCCACCAGCATGACAAAGCCCGTCGCCCAAATCa < 1:1931418/69‑1 (MQ=255) tCGTGAACCCATGCGCGCGCCGAGGAACTGCCCCACCAGc < 1:3230814/40‑1 (MQ=255) ggAACTGTCCTACCAACATCACAAATCCTGTCGCCCAAATCACTTTGCCCCCgaga > 1:2883605/1‑56 (MQ=255) ggAACTGTCCTACCAACATCACAAATCCTGTCGCCCAAATCACTTTGCCACCGAGAATAAAcagcagca > 1:3342856/1‑67 (MQ=255) ggAACTGTCCTACCAACATCACAAATCCTGTCGCCCAAATCACTTTGCCACCGAGAATAAAcagcagca > 1:995878/1‑67 (MQ=255) ggAACTGTCCTACCAACATCACAAATCCTGTCGCCCAAATCACTTTGCCACCGAGAATAAAcagcagca > 1:348363/1‑67 (MQ=255) ggAACTGTCCTACCAACATCACAAATCCTGTCGCCCAAATCACTTTGCCACCGAGAATAAAcagcagca > 1:39789/1‑67 (MQ=255) ggAACTGTCCTACCAACATCACAAATCCTGTCGCCCAAATCACTTTGCCACCGAGAATAAAcagcagca > 1:1045303/1‑67 (MQ=255) ggAACTGTCCTACCAACATCACAAATCCTGTCGCCCAAATCACTTTGCCACCGAGAATAAAcagcagca > 1:697596/1‑67 (MQ=255) ggAACTGTCCTACCAACATCACAAATCCTGTCGCCCAAATCACTTTGCCACCGAGAATAAAcagcagca > 1:883012/1‑67 (MQ=255) ggAACTGTCCTACCAACATCACAAATCCTGTCGCCCAAATCACTTTGCCACCGAGAATAAAcagcagca > 1:3289835/1‑67 (MQ=255) ggAACTGTCCTACCAACATCACAAATCCTGTCGCCCAAATCACTTTGCCACCGAGAATAAAcagcagca > 1:3226122/1‑67 (MQ=255) ggAACTGTCCTACCAACATCACAAATCCTGTCGCCCAAATCACTTTGCCACCGAGAATAAAcagcagca > 1:1162714/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:1732022/1‑67 (MQ=255) ggAACTGTCCTACCAACATCACAAATCCTGTCGCCCAAATCACTTTGCCACCGAGAATAAAcagcagca > 1:1896842/1‑67 (MQ=255) ggAACTGTCCTACCAACATCACAAATCCTGTCGCCCAAATCACTTTGCCACCGAGAATAAAcagcagca > 1:1905094/1‑67 (MQ=255) ggAACTGTCCTACCAACATCACAAATCCTGTCGCCCAAATCACTTTGCCACCGAGAATAAAcagcagca > 1:1926696/1‑67 (MQ=255) | ATCAGCTTTTGTCCTTTGCTCAACACCAGTCGTGAACCCATGCGCGCGCCGAGGAACTGCCCCACCAGCATGACAAAGCCCGTCGCCCAAATCACTTTGCCGCCGAGAATAAACAGCAGTA > minE/1530100‑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: ‑ |