Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | Exported | 2,045,818 | C→T | 6.5% | intergenic (+580/+215) | yqiK → / ← rfaE | conserved hypothetical protein/fused heptose 7‑phosphate kinase and heptose 1‑phosphate adenyltransferase |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | Exported | 2,045,818 | 0 | C | T | 6.5% | 210.4 / 3.9 | 76 | intergenic (+580/+215) | yqiK/rfaE | conserved hypothetical protein/fused heptose 7‑phosphate kinase and heptose 1‑phosphate adenyltransferase |
Reads supporting (aligned to +/- strand): ref base C (31/40); new base T (2/3); total (33/43) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 1.00e+00 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.59e-01 |
AAAGCATTGCGCGGGGTCCGGCGCTCATCAACAATCGGGGGGCAGCAAGGGGCTGAAACGGGAAAGCCCCTCCCGAAGAAG‑GGGCCTTGTATAAGGAAAGGGTTATGATGAAGCTCGTCATCATACTGGTTGTGTTGT > Exported/2045750‑2045887 | aaaGCATTGCGCGGGGTCCGGCGCTCATCAACAATCGGGGGGCAGCAAGGGGCTGAAACGGGAAAGCCCCt < 1:1302570/71‑1 (MQ=255) ttGCGCGGGGTCCGGCGCTCATCAACAATCGGGGGGCAGCAAGGGGCTGAAACGGGAAAGCCCCTCCCgaa < 1:4422296/71‑1 (MQ=255) tGCGCGGGGTCCGGCGCTCATCAACAATCGGGGGGCAGCAAGGGGCTGAAACGGGAAAGCCCCTCCCg > 1:565428/1‑68 (MQ=255) cgcgGGGTCCGGCGCTCATCAACAATCGGGGGGCAGCAAGGGGCTGAAACGGGAAAGCCCCTCCCgaagaa < 1:3636028/71‑1 (MQ=255) cgcgGGGTCCGGCGCTCATCAACAATCGGGGGGCAGCAAGGGGCTGAAACGGGAAAGCCCCTCCCgaag > 1:1446895/1‑69 (MQ=255) gTCCGGCGCTCATCAACAATCGGGGGGCAGCAAGGGGCTGAAACGGGAAAGCCCCTCCCGAAGAAG‑GGGcc > 1:4650558/1‑71 (MQ=255) cGGCGCTCATCAACAATCGGGGGGCAGCAAGGGGCTGAAACGGGAAAGCCCCTCCCgaaga > 1:4429500/1‑61 (MQ=255) gcTCATCAACAATCGGGGGGCAGCAAGGGGCTGAAACGGGAAAGCCCCTCCCGAAGAAG‑GGGCCTTGta > 1:5442407/1‑69 (MQ=255) cTCATCAACAATCGGGGGGCAGCAAGGGGCTGAAACGGGAAAGCCCCTCCCGAAGAAG‑GGGCCTTGtata < 1:1887124/70‑1 (MQ=255) atcaACAATCGGGGGGCAGCAAGGGGCTGAAACGGGAAAGCCCCTCCCGAAGAAG‑GGGCCTTGtat < 1:2333700/66‑1 (MQ=255) atcaACAATCGGGGGGCAGCAAGGGGCTGAAACGGGAAAGCCCCTCCCGAAGAAG‑GGGCCTTGTATAAgg > 1:3343442/1‑70 (MQ=255) atcaACAATCGGGGGGCAGCAAGGGGCTGAAACGGGAAAGCCCCTCCCGAAGAAG‑GGGCCTTGTATAAg < 1:5466060/69‑1 (MQ=255) tcaACAATCGGGGGGCAGCAAGGGGCTGAAACGGGAAAGCCCCTCCCGAAGAAG‑GGGCCTTg < 1:4827154/62‑1 (MQ=255) tcaACAATCGGGGGGCAGCAAGGGGCTGAAACGGGAAAGCCCCTCCCGAAGAAG‑GGGCCTTg < 1:1552405/62‑1 (MQ=255) tcaACAATCGGGGGGCAGCAAGGGGCTGAAACGGGAAAGCCCCTCCCGAAGAAG‑GGGCCTTGTATAAGGaa < 1:2717963/71‑1 (MQ=255) aacaaTCGGGGGGCAGCAAGGGGCTGAAACTGGAAAGCCCCTCCCGAAGAAG‑GGGCCTTGTATAAGGa < 1:892892/68‑1 (MQ=255) aacaaTCGGGGGGCAGCAAGGGGCTGAAACGGGAAAGCCCCTCCCgaag > 1:822613/1‑49 (MQ=255) aacaaTCGGGGGGCAGCAAGGGGCTGAAACGGGAAAGCCCCTCCCGAAGAAG‑GGGCCTTGTATAAGGAAAg > 1:2466718/1‑71 (MQ=255) aacaaTCGGGGGGCAGCAAGGGGCTGAAACGGGAAAGCCCCTCCCGAAGAAG‑GGGCCTTGTATAAGGAAAg > 1:5640262/1‑71 (MQ=255) ctaaTCGGGGGGCAGCAAGGGGCTGAAACGGGAAAGCCCCTCCCGAAGAAG‑GGGCCTTGTATAAGGaaa < 1:4132819/67‑1 (MQ=255) acaaTCGGGGGGCAGCAAGGGGCTGAAACGGGAAAGCCCCTCCCGAAGAAG‑GGGCCTTg < 1:5090951/59‑1 (MQ=255) acaaTCGGGGGGCAGCAAGGGGCTGAAACGGGAAAGCCCCTCCCGAAGAAG‑GGGCCTTg < 1:711585/59‑1 (MQ=255) acaaTCGGGGGGCAGCAAGGGGCTGAAACGGGAAAGCCCCTCCCGAAGAAG‑GGGCCTTg < 1:1271467/59‑1 (MQ=255) gggggCAGCAAGGGGCTGAAACGGGAAAGCCCCTCCCgaag < 1:914440/41‑1 (MQ=255) gggggCAGCAAGGGGCTGAAACGGGAAAGCCCCTCCCgaag < 1:3182772/41‑1 (MQ=255) ggCAGCAAGGGGCTGAAACGGGAAAGCCCCTCCCGAAGAAG‑GGGCCTTGTATAAgg < 1:4051813/56‑1 (MQ=255) ggCAGCAAGGGGCTGAAACGGGAAAGCCCCTCCCGAAGAAG‑GGGCCTTGTATAAgg < 1:276982/56‑1 (MQ=255) ggCAGCAAGGGGCTGAAACGGGAAAGCCCCTCCCGAAGAAG‑GGGCCTTGTATAAgg < 1:346786/56‑1 (MQ=255) ggCAGCAAGGGGCTGAAACGGGAAAGCCCCTCCCGAAGAAG‑GGGCCTTGTATAAgg < 1:2562837/56‑1 (MQ=255) gcagcaAGGGGCTGAAACGGGAAAGCCCCTCCCGAAGAAG‑GGGCCTTGTATAAgg < 1:165441/55‑1 (MQ=255) gcagcaAGGGGCTGAAACGGGAAAGCCCCTCCCGAAGAAG‑GGGCCTTGTATAAGGa < 1:733399/56‑1 (MQ=255) cagcaAGGGGCTGAAACGGGAAAGCCCCTCCCGAAGAAg‑g > 1:2605999/1‑40 (MQ=40) cagcaAGGGGCTGAAACGGGAAAGCCCCTCCCGAAGAAG‑GGGCCTTGTATAAGGAAAGGGTTATGATGAAg > 1:830053/1‑71 (MQ=255) agcaAGGGGCTGAAACGGGAAAGCCCCTCCCGAAGAAG‑GGGCCTTGTATAAGGAAAGGGTTa < 1:4329956/62‑1 (MQ=255) agcaAGGGGCTGAAACGGGAAAGCCCCTCCCGAAGAAG‑GGGCCTTGTATAAGGAAAGGGTTa < 1:2531398/62‑1 (MQ=255) gcaAGGGGCTGAAACGGGAAAGCCCCTCCCgaaga < 1:4260276/35‑1 (MQ=40) gcaAGGGGCTGAAACGGGAAAGCCCCTCCCgaaga < 1:3957118/35‑1 (MQ=40) gcaAGGGGCTGAAACGGGAAAGCCCCTCCCgaaga < 1:4963063/35‑1 (MQ=40) gcaAGGGGCTGAAACGGGAAAGCCCCTCCCGAAGAAG‑GGGCCTTGTATAAgg < 1:4661619/52‑1 (MQ=255) gcaAGGGGCTGAAACGGGAAAGCCCCTCCCGAAGAAG‑GGGCCTTGTATAAGGAAAgg < 1:3000704/57‑1 (MQ=255) gcaAGGGGCTGAAACGGGAAAGCCCCTCCCGAAGAAG‑GGGCCTTGTATAAGGAAAGGGTTatg > 1:4328361/1‑63 (MQ=255) aaGGGGCTGAAACGGGAAAGCCCCTCCCGAAGAAG‑GGGCCTTGTATAAGGAAAGGGt > 1:781119/1‑57 (MQ=255) aGGGGCTGAAACGGGAAAGCCCCTCCCGAAGAAG‑GGGCCTTGTATAAGGAAAGGGTTa > 1:3736880/1‑58 (MQ=255) aGGGGCTGAAACGGGAAAGCCCCTCCCGAAGAAG‑GGGCCTTGTATAAGGAAAGGGTTATGATGAAGCTc < 1:657273/69‑1 (MQ=255) gggCTGAAACGGGAAAGCCCCTCCCGAAGAAG‑GGTCCTTGTATAAGGAAAGGGTTa < 1:1977594/56‑1 (MQ=255) acgCTGAAACGGGAAAGCCTCTCCCGGAGAAGAGGGCTTTTAATAAGGAAAGGGTTATGATGAAGCACGt > 1:303292/3‑70 (MQ=0) ggCTGAAACGGGAAAGCCCCTCCCGAAGAAG‑GGGCCTTGtat < 1:57129/42‑1 (MQ=255) ggCTGAAACGGGAAAGCCCCTCCCGAAGAAG‑GGGCCTTGTATAAGGAAAGGGTTATGATGAAGCTCGtcat > 1:552372/1‑71 (MQ=255) ggCTGAAACGGGAAAGCCCCTCCCGAAGAAG‑GGGCCTTGTATAAGGAAAGGGTTATGATGAAGCTCGtcat > 1:5627124/1‑71 (MQ=255) ggCTGAAACGGGAAAGCCCCTCCCGAAGAAG‑GGGCCTTGTATAAGGAAAGGGTTATGATGAAGCTCGtcat > 1:3187602/1‑71 (MQ=255) cgCTGAAACGGGAAAGCCTCTCCCGGAGAAGAGGGCTTTTAATAAGGAAAGGGTTATGATGAAGCACGtca < 1:1786500/70‑1 (MQ=1) cTGAAACGGGAAAGCCCCTCCCGAGGAAG‑GGGCCATAAATAAGGaaa > 1:3653992/1‑47 (MQ=14) cTGAAACGGGAAAGCCCCTCCCGAGGAAG‑GGGCCATAAATAAGGAAAGGGTCatg < 1:144998/55‑1 (MQ=14) cTGAAACGGGAAAGCCCCTCCCGAAGAAG‑GGGCCTTGTATaa > 1:442617/1‑42 (MQ=255) gAAACGGGAAAGCCCCTCCCGAGGAAG‑GGGCCATAAATAAGGaa > 1:1492151/1‑44 (MQ=14) gAAACGGGAAAGCCCCTCCCGAGGAAG‑GGGCCATAAATAAGGAAAGGGTCATGATGAAGCTACtcat > 1:1889915/1‑67 (MQ=14) gAAACGGGAAAGCCCCTCCCGAGGAAG‑GGGCCATAAATAAGGAAAGGGTCATGATGAAGCTACtcat > 1:1896669/1‑67 (MQ=14) gAAACGGGAAAGCCCCTCCCGAAGAAG‑GGGCCTTGTATAAGGAAAGGGtt < 1:1190765/50‑1 (MQ=255) aaaCGGGAAAGCCCCTCCCGAGGAAG‑GGGCCATAAATAAGGAAAGGGTCatgat < 1:4403377/54‑1 (MQ=14) aaaCGGGAAAGCCCCTCCCGAGGAAG‑GGGCCATAAATAAGGAAAGGGTCatgat < 1:3787364/54‑1 (MQ=14) aaaCGGGAAAGCCCCTCCCGAGGAAG‑GGGCCATAAATAAGGAAAGGGTCATGATGAAGCTACTcatcatcg < 1:565443/71‑3 (MQ=11) aaaCGGGAAAGCCCCTCCCGAAGAAG‑GGGCCTTGTATAAGGAAAGGGTTATGATGAAGCTCGTCATCATAc > 1:4899675/1‑71 (MQ=255) aaCGGGAAAGCCCCTCCCGAAGAAG‑GGGCCTTGTATAAGGAAAGGGTTATGATGAAg > 1:3535738/1‑57 (MQ=255) aCGGGAAAGCCCCTCCCGAGGAAG‑GGGCCATAAATAAGGAAAGGGTCAtgatga > 1:3633700/1‑54 (MQ=14) aCGGGAAAGCCCCTCCCGAGGAAG‑GGGCCATAAATAAGGAAAGGGTCAtgatga > 1:3590871/1‑54 (MQ=14) aCGGGAAAGCCCCTCCCGAAGAAG‑GGGCCTTGTATAATGAAAGGGTTATGATGAAGCTCGt < 1:5193989/61‑1 (MQ=255) ggAAAGCCCCTCCCGAGGAAG‑GGGCCATAAATAAGGAAAGGGTCatgat < 1:2622584/49‑1 (MQ=14) ggAAAGCCCCTCCCGAGGAAG‑GGGCCATAAATAAGGAAAGGGTCATGATGAAGCTACTCATCATCGTGgtg > 1:1155474/1‑70 (MQ=2) ggAAAGCCCCTCCCGAAGAAG‑GGGCCTTGTATAAGGAAAGGGTTATGATGAAGCTCGTCATCATACTGgtt > 1:3793091/1‑71 (MQ=255) gAAAGCCTCTCCCGGAGAAGAGGGCTTTTAATAAGGAAAGGGTTATGATGAAGCACGTCATCATACTGgt > 1:3680797/1‑70 (MQ=11) gAAAGCCCCTCCCGAGGAAG‑GGGCCATAAATAAGGAAAGGGTCATGATGAAGCTACTCATCATCGTGgtgc > 1:3332284/1‑69 (MQ=2) gAAAGCCCCTCCCGAAGAAG‑GGGCCTTGTATAAGGAAAGGGTTATGATGAAGCTCGTCATCATACTGgttg < 1:2081641/71‑1 (MQ=255) aaaGCCCCTCCCGAGGAAG‑GGGCCATAAATAAGGAAAGGGTCATGATGAAGCTACTCATCATCGTGgtgc > 1:4164927/1‑68 (MQ=2) aaaGCCCCTCCCGAGGAAG‑GGGCCATAAATAAGGAAAGGGTCATGATGAAGCTACTCATCATCGTGgtgc > 1:3132334/1‑68 (MQ=2) aGCCTCTCCCGGAGAAGAGGGCTTTTAATAAGGAAAGGGt < 1:3505689/40‑1 (MQ=0) aGCCTCTCCCGGAGAAGAGGGCTTTTAATAAGGAAAGGGt < 1:22867/40‑1 (MQ=0) cccTCCCGAAGAAG‑GGGCCTTGTATAAGGAAAGGGTTATGATGAAGCTCGTCATCATACt < 1:1253870/60‑1 (MQ=255) ccTCCCGAAGAAG‑GGGCCTTGTATAAGGAAAGGGTTATGATGAAGCTCGTCATCATACTGgttgtgttgt > 1:3907027/1‑70 (MQ=255) | AAAGCATTGCGCGGGGTCCGGCGCTCATCAACAATCGGGGGGCAGCAAGGGGCTGAAACGGGAAAGCCCCTCCCGAAGAAG‑GGGCCTTGTATAAGGAAAGGGTTATGATGAAGCTCGTCATCATACTGGTTGTGTTGT > Exported/2045750‑2045887 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |