Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | minE | 1,516,244 | C→T | 77.3% | intergenic (‑229/+30) | cvpA ← / ← dedD | membrane protein required for colicin V production/conserved hypothetical protein |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | minE | 1,516,244 | 0 | C | T | 77.3% | 36.1 / 7.3 | 35 | intergenic (‑229/+30) | cvpA/dedD | membrane protein required for colicin V production/conserved hypothetical protein |
Reads supporting (aligned to +/- strand): ref base C (3/5); new base T (13/14); total (16/19) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 7.00e-01 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.35e-01 |
CGTGTCGGATGCGGCGTAAACGCCTTATCCGACCTACGGTTCTACCCCTGCGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCAAGACCGTATTAATTCGGCGTATAGCCCATTACCACGCCACTTAA > minE/1516178‑1516309 | cGTGTCGGATGCGGCGTAAACGCCTTATCCGACCTACGGTTCTACCCCTGCGTAGGCCTGATAAGAcgc > 1:389857/1‑69 (MQ=255) cGTGTCGGATGCGGCGTAAACGCCTTATCCGACCTACGGTTCTACCCCTGCGTAGGCCTGATAAGAcgc > 1:434621/1‑69 (MQ=255) tGCGGCGTAAACGCCTTATCCGACCTACGGTTCTACCCCTGCGTAGGCCTGATAAGACGCGCCAGCGt > 1:2731005/1‑68 (MQ=255) tGCGGCGTAAACGCCTTATCCGACCTACGGTTCTACCCCTGCGTAGGCCTGATAAGACGCGCCAGCGt > 1:2585682/1‑68 (MQ=255) gcgTAAACGCCTTATCCGACCTACGGTTCTACCCCTGCGTAGGCCTGATAAGACGCGCCAGCGTCGCa < 1:1205899/68‑1 (MQ=255) gcgTAAACGCCTTATCCGACCTACGGTTCTACCCCTGCGTAGGCCTGATAAGACGCGCCAGCGTCGCa < 1:1731852/68‑1 (MQ=255) gCCTTATCCGACCTACGGCTCTACCCCTGCACAGGCCTGATAAGATGCGCCAGCGTCGCATCAGGCAAg > 1:3168003/1‑69 (MQ=21) gCCTTATCCGACCTACGGCTCTACCCCTGCACAGGCCTGATAAGATGCGCCAGCGTCGCATCAGGCAAg > 1:2873095/1‑69 (MQ=21) gCCTTATCCGACCTACGGCTCTACCCCTGCACAGGCCTGATAAGATGCGCCAGCGTCGCATCAGGCAAg > 1:3326974/1‑69 (MQ=21) gCCTTATCCGACCTACGGCTCTACCCCTGCACAGGCCTGATAAGATGCGCCAGCGTCGCATCAGGCAAg > 1:374672/1‑69 (MQ=21) gCCTTATCCGACCTACGGCTCTACCCCTGCACAGGCCTGATAAGATGCGCCAGCGTCGCATCAGGCAAg > 1:1885099/1‑69 (MQ=21) gCCTTATCCGACCTACGGCTCTACCCCTGCACAGGCCTGATAAGATGCGCCAGCGTCGCATCAGGCAAg > 1:1830008/1‑69 (MQ=21) gCCTTATCCGACCTACGGCTCTACCCCTGCACAGGCCTGATAAGATGCGCCAGCGTCGCATCAGGCAAg > 1:561681/1‑69 (MQ=21) gCCTTATCCGACCTACGGCTCTACCCCTGCACAGGCCTGATAAGATGCGCCAGCGTCGCATCAGGCAAg > 1:1651708/1‑69 (MQ=21) gCCTTATCCGACCTACGGCTCTACCCCTGCACAGGCCTGATAAGATGCGCCAGCGTCGCATCAGGCAAg > 1:612660/1‑69 (MQ=21) gCCTTATCCGACCTACGGCTCTACCCCTGCACAGGCCTGATAAGATGCGCCAGCGTCGCATCAGGCAAg > 1:147323/1‑69 (MQ=21) gCCTTATCCGACCTACGGCTCTACCCCTGCACAGGCCTGATAAGATGCGCCAGCGTCGCATCAGGCAAg > 1:1457175/1‑69 (MQ=21) gCCTTATCCGACCTACGGCTCTACCCCTGCACAGGCCTGATAAGATGCGCCAGCGTCGCATCAGGCAAg > 1:1443978/1‑69 (MQ=21) gCCTTATCCGACCTACGGCTCTACCCCTGCACAGGCCTGATAAGATGCGCCAGCGTCGCATCAGGCAAg > 1:1279956/1‑69 (MQ=21) cTACCCCTGCACAGGCCTGATAAGATGCGCCAGCGTCGCATCAGGCAAGACCGTATTAATTCGGTGtat < 1:51339/69‑1 (MQ=21) cTACCCCTGCACAGGCCTGATAAGATGCGCCAGCGTCGCATCAGGCAAGACCGTATTAATTCGGTGtat < 1:3170951/69‑1 (MQ=21) cTACCCCTGCACAGGCCTGATAAGATGCGCCAGCGTCGCATCAGGCAAGACCGTATTAATTCGGTGtat < 1:663647/69‑1 (MQ=21) cTACCCCTGCACAGGCCTGATAAGATGCGCCAGCGTCGCATCAGGCAAGACCGTATTAATTCGGTGtat < 1:961474/69‑1 (MQ=21) cTACCCCTGCACAGGCCTGATAAGATGCGCCAGCGTCGCATCAGGCAAGACCGTATTAATTCGGTGtat < 1:1254729/69‑1 (MQ=21) cTACCCCTGCACAGGCCTGATAAGATGCGCCAGCGTCGCATCAGGCAAGACCGTATTAATTCGGTGtat < 1:3050464/69‑1 (MQ=21) cTACCCCTGCACAGGCCTGATAAGATGCGCCAGCGTCGCATCAGGCAAGACCGTATTAATTCGGTGtat < 1:3017230/69‑1 (MQ=21) cTACCCCTGCACAGGCCTGATAAGATGCGCCAGCGTCGCATCAGGCAAGACCGTATTAATTCGGTGtat < 1:257830/69‑1 (MQ=21) cTACCCCTGCACAGGCCTGATAAGATGCGCCAGCGTCGCATCAGGCAAGACCGTATTAATTCGGTGtat < 1:2397845/69‑1 (MQ=21) cTACCCCTGCACAGGCCTGATAAGATGCGCCAGCGTCGCATCAGGCAAGACCGTATTAATTCGGTGtat < 1:1872313/69‑1 (MQ=21) cTACCCCTGCACAGGCCTGATAAGATGCGCCAGCGTCGCATCAGGCAAGACCGTATTAATTCGGTGtat < 1:184224/69‑1 (MQ=21) cTACCCCTGCACAGGCCTGATAAGATGCGCCAGCGTCGCATCAGGCAAGACCGTATTAATTCGGTGtat < 1:1713975/69‑1 (MQ=21) cTACACCTGCACAGGCCTGATAAGATGCGCCAGCGTCGCATCAGGCAAGACCGTATTAATTCGGTGtat < 1:1341058/69‑1 (MQ=21) ccTGCGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCAAGACCGTATTAATTCGGCg < 1:886418/61‑1 (MQ=255) ccTGCGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCAAGACCGTATTAATTCGGCGTATAGc < 1:2503991/67‑1 (MQ=255) ccTGCGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCAAGACCGTATTAATTCGGCGTATAGc < 1:1612176/67‑1 (MQ=255) aGGCCTGATAAGATGCGCCAGCGTCGCATCAGGCAAGACCGTATTAATTCGGTGtat < 1:1814456/57‑1 (MQ=21) aGACGCGCCAGCGTCGCATCAGGCAAGACCGTATTAATTCGGCGTATAGCCCATTACCACGCCACTTaa > 1:1279962/1‑69 (MQ=255) | CGTGTCGGATGCGGCGTAAACGCCTTATCCGACCTACGGTTCTACCCCTGCGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCAAGACCGTATTAATTCGGCGTATAGCCCATTACCACGCCACTTAA > minE/1516178‑1516309 |
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: ‑ |