| Predicted mutation | |||||||
|---|---|---|---|---|---|---|---|
| evidence | seq id | position | mutation | freq | annotation | gene | description |
| RA | minE | 1,716,835 | G→A | 7.2% | intergenic (+93/+105) | yphH → / ← glyA | predicted DNA‑binding transcriptional regulator/serine hydroxymethyltransferase |
| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | minE | 1,716,835 | 0 | G | A | 7.2% | 212.9 / 5.5 | 70 | intergenic (+93/+105) | yphH/glyA | predicted DNA‑binding transcriptional regulator/serine hydroxymethyltransferase |
| Reads supporting (aligned to +/- strand): ref base G (31/29); new base A (2/3); total (33/37) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 6.72e-01 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.38e-01 | |||||||||||
CCCGAAGGAAGCCTTGATGATCATAAACGAAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACACGGAG‑ATTGCAATATATTGAATTTGCAAAGTTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAA > minE/1716767‑1716905 | cccGAAGGAAGCCTTGATGATCATAAACGAAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACACGga < 1:3585983/70‑1 (MQ=255)cccGAAGGAAGCCTTGATGATCATAAACGAAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACACGga < 1:1290254/70‑1 (MQ=255)cccGAAGGAAGCCTTGATGATCATAAACGAAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACACGga < 1:1661947/70‑1 (MQ=255)cccGAAGGAAGCCTTGATGATCATAAACGAAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACACGga < 1:3790146/70‑1 (MQ=255) ccGAAGGAAGCCTTGATGATCATAAACGAAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACACGgag < 1:3879586/70‑1 (MQ=255) cGAAGGAAGCCTTGATGATCATAAACGAAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACACGga < 1:3261819/68‑1 (MQ=255) cGAAGGAAGCCTTGATGATCATAAACGAAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACACGGAG‑At < 1:1044909/71‑1 (MQ=255) agCCTTGATGATCATAAACGAAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACACGGAG‑ATTGCAata > 1:2977423/1‑71 (MQ=255) gCCTTGATGATCATAAACGAAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACACGga < 1:4078590/60‑1 (MQ=255) ccTTGATGATCATAAACGAAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACACGGAG‑ATTGCAatata > 1:3051359/1‑71 (MQ=255) ttGATGATCATAAACGAAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACACGga < 1:2090217/57‑1 (MQ=255) tgatgaTCATAAACGAAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACACGGAG‑ATTGCAATATAtt > 1:1846672/1‑70 (MQ=255) gatgatCATAAACGAAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACACGGAG‑ATTGCAATATATTg < 1:3964185/70‑1 (MQ=255) atgatCATAAACGAAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACACGGAG‑ATTGCaa > 1:1685762/1‑62 (MQ=255) atgatCATAAACGAAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACACGGAG‑ATTGCaa > 1:3674266/1‑62 (MQ=255) tgatCATAAACGAAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACACGGAG‑ATTGCAATATATTGAAt > 1:1039740/1‑71 (MQ=255) gatCATAAACGAAAAATTGCCTGATGCGCTGCGCTTATCAGGCCTACACGgag < 1:301944/53‑1 (MQ=255) tCATAAACGAAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACACGGAG‑ATTGCAATATATTGAATTTg < 1:1006996/71‑1 (MQ=255) cATAAACGAAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACACGGAG‑ATTGCAATATATTGAATTTg > 1:2701873/1‑70 (MQ=255) cATAAAAGAAAAATTGCCTGATGCTCTACGCTTATCAGGCCTACACGGAG‑ATTGCAATATATTGAATTTg > 1:1212670/1‑70 (MQ=255) aTAAACGAAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACACTga < 1:3571552/48‑1 (MQ=255) aTAAACGAAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACACGGAG‑ATTGCAata > 1:1481500/1‑58 (MQ=255) tAAACGAAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACACGGAG‑ATTGCAatata < 1:3940519/59‑1 (MQ=255) aaaCGAAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACACGGAG‑ATTGCAATATATTGAATTTGCa > 1:2461888/1‑69 (MQ=255) aaaCGAAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACACGGAG‑ATTGCAATATATTGAATTTGCa > 1:2118127/1‑69 (MQ=255) aaCGAAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACACGgag‑a > 1:194831/1‑47 (MQ=255) aaCGAAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACACGGAG‑ATTGCAATATATTGAATTTGCAAAg > 1:1661596/1‑71 (MQ=255) aaCGAAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACACGGAG‑ATTGCAATATATTGAATTTGCAAAg > 1:3265183/1‑71 (MQ=255) aCGAAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACACGGAG‑ATTGCAATATATTGAATTTGCAAAGt < 1:1174803/71‑1 (MQ=255) aCGAAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACACGGAG‑ATTGCAATATATTGAATTTGCAAAGt < 1:2285852/71‑1 (MQ=255) aCGAAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACACGGAG‑ATTGCAATATATTGAATTTGCAAAGt < 1:1419646/71‑1 (MQ=255) gAAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACACGGAG‑ATTGCAATATATTGAATTTGCAAAGtt < 1:1172683/70‑1 (MQ=255) caaaTTGCCTGATGCGCTACGCTTATCGGGCCTACAAGTGC‑ATTGCAATCTATTGAATTTGCGATGTTTTg < 1:585976/70‑1 (MQ=11) caaaTTGCCTGATGCGCTACGCTTATCAGGCCTACAAGTGC‑ATTGCAATCTATTGAATTTGCGATGTTTTg < 1:2980459/70‑1 (MQ=11) caaaTTGCCTGATGCGCTACGCTTATCAGGCCTACAAGTGC‑ATTGCAATCTATTGAATTTGCGATGTTTTg < 1:1180270/70‑1 (MQ=11) aaaaTTGCCTGATGCGCTACGCTTATCAGGCCTACACGGAG‑ATTGCAATATATTGAATTTGCAAAGTTTTg < 1:2767673/71‑1 (MQ=255) aaaaTTGCCTGATGCGCTACGCTTATCAGGCCTACACGGAG‑ATTGCAATATATTGAATTTGCAAAGTTTTg < 1:3799106/71‑1 (MQ=255) aaaaTTGCCTGATGCGCTACGCTTATCAGGCCTACACGGAG‑ATTGCAATATATTGAATTTGCAAAGTTTTg < 1:2501861/71‑1 (MQ=255) aaaTTGCCTGATGCGCTGCGCTTATCAGGCCTACAAGTTGAATTGCAATTTATTGAATTTGCACATTTTTg < 1:1206041/71‑1 (MQ=2) aaaTTGCCTGATGCGCTACGCTTATCAGGCCTACACGGAG‑ATTTCAATATATTGAATTTGCAAAGTTTTGt < 1:3970829/71‑1 (MQ=255) aaaTTGCCTGATGCGCTACGCTTATCAGGCCTACACGGAG‑ATTGCAATATATTGAATTTGCAAAGTTTTg > 1:1864276/1‑70 (MQ=255) aaaTTGCCTGATGCGCTACGCTTATCAGGCCTACACGGAG‑ATTGCAATATATTGAATTTGCAAAGTTTTGt < 1:811366/71‑1 (MQ=255) gCCTGATGCGCTACGCTTATCAGGCCTACACGGAG‑ATTGCAATATATTGAATTTGCAAAGTTTTGTAGGcc > 1:911080/1‑71 (MQ=255) ccTGATGCGCTACGCTTATCAGGCCTACACGGAG‑ATTGCAATATATTGAATTTGCAAAGTTTTGTAGGCCg < 1:61138/71‑1 (MQ=255) cTGATGCGCTACGCTTATCAGGCCTACACGGAG‑ATTGCAATATATTGAATTTGCAAAGTTTTGTAGGCCgg < 1:1305034/71‑1 (MQ=255) gATGCGCTACGCTTATCAGGCCTACACGGAG‑ATTGCAATATATTGAATTTGCAAAGTTTTGTAGGCCGGa > 1:3293141/1‑70 (MQ=255) gATGCGCTACGCTTATCAGGCCTACACGGAG‑ATTGCAATATATTGAATTTGCAAAGTTTTGTAGGCCGGa > 1:1815686/1‑70 (MQ=255) gATGCGCGATGCTTATCAGGCCTACCAGAAG‑ATTGCAATATATTGAATTTGCACTGTTTTGTAGGCCGGa > 1:834745/1‑70 (MQ=11) aTGCGCTACGCTTATCAGGCCTACACGGAG‑ATTGCAATATATTGAATTTGCAAAGTTTTGTAGGCCGGATa < 1:1412157/71‑1 (MQ=255) aTGCGCGATGCTTATCAGGCCTACCAGAAG‑ATTGCAATATATTGAATTTGCACTGTTTTGTAGGCCGGAt < 1:1877370/70‑1 (MQ=11) cTACGCTTATCAGGCCTACACGGAG‑ATTGCAATATATTGAATTTGCAAAGTTTTGTAGGCCGGATAAGGCg < 1:346244/71‑1 (MQ=255) aCGCTTATCAGGCCTACACGGAG‑ATTGCAATATATTGAATTTGCaaa > 1:1564827/1‑47 (MQ=255) aCGCTTATCAGGCCTACACGGAG‑ATTGCAATATATTGAATTTGCAAAGtttt > 1:1861805/1‑52 (MQ=255) aCGCTTATCAGGCCTACACGGAG‑ATTGCAATATATTGAATTTGCAAAGTTTTGTAGGCCGGATAAGGCGtt > 1:3359635/1‑71 (MQ=255) aCGCTTATCAGGCCTACACGGAG‑ATTGCAATATATTGAATTTGCAAAGTTTTGTAGGCCGGATAAGGCGtt > 1:2876386/1‑71 (MQ=255) tgCTTATCAGGCCTACCAGAAG‑ATTGCAATATATTGAATTTGCACTGTTTTg < 1:3181163/51‑1 (MQ=11) tgCTTATCAGGCCTACCAGAAG‑ATTGCAATATATTGAATTTGCACTGTTTTg < 1:1662940/51‑1 (MQ=11) cGCTTATCAGGCCTACACGGAG‑ATTGCAATATATTGAATTTGCAAAGTTTTGTAGGCCGGATAAGGCGTTc < 1:843926/71‑1 (MQ=255) cTTATCAGGCCTACACGGAG‑ATTGCAATATATTGAATTTGCAAAGttt > 1:2490392/1‑48 (MQ=255) cTTATCAGGCCTACACGGAG‑ATTGCAATATATTGAATTTGCAAAGTTTTGTAGGCCGGATAAGGCGTTCAc > 1:2571736/1‑71 (MQ=255) tATCAGGCCTACCAGAAG‑ATTGCAATATATTGAATTTGCact > 1:3949887/1‑40 (MQ=11) aTCAGGCCTACACGGAG‑ATTGCAATATATTGAATTTGCAAAGTTTTGTAGGCCGGATAAGGCGTTCAcgc < 1:2465711/70‑1 (MQ=255) aTCAGGCCTACAAGTGC‑ATTGCAATCTATTGAATTTGCGATGTTTTGTAGGACGGATAAGGCGCTcgcg < 1:2009239/69‑1 (MQ=2) aGGCCTACACGGAG‑ATTGCAATATATTGAATTTGCAAAGTTTTGTAGGCCGGATAAGGCGTTCAcgccg > 1:181541/1‑69 (MQ=255) gCCTACACGGAG‑ATTGCAATATATTGAATTTGCAAAGTTTTGTAGGCCGGATAAGGCGTTCACGCCGCATc > 1:528223/1‑71 (MQ=255) gCCTACACGGAG‑ATTGCAATATATTGAATTTGCAAAGTTTTGTAGGCCGGATAAGGCGTTCACGCCGCATc < 1:2128009/71‑1 (MQ=255) cTACACGGAG‑ATTGCAATATATTGAATTTGCAAAGttt < 1:1282849/38‑1 (MQ=255) cTACACGGAG‑ATTGCAATATATTGAATTTGCAAAGTTTTGTAGGCCGGATAAGGCGtt > 1:939185/1‑58 (MQ=255) cTACACGGAG‑ATTGCAATATATTGAATTTGCAAAGTTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCg > 1:1275493/1‑71 (MQ=255) tACACGGAG‑ATTGCAATATATTGAATTTGCAAAGttt < 1:3633273/37‑1 (MQ=255) tACACGGAG‑ATTGCAATATATTGAATTTGCAAAGTTTTGTAGGCCgg < 1:1101214/47‑1 (MQ=255) tACACGGAG‑ATTGCAATATATTGAATTTGCAAAGTTTTGTAGGCCg < 1:3525913/46‑1 (MQ=255) tACACGGAG‑ATTGCAATATATTGAATTTGCAAAGTTTTGTAGGCCGGAt < 1:1665877/49‑1 (MQ=255) tACACGGAG‑ATTGCAATATATTGAATTTGCAAAGTTTTGTAGGCCGGATAAg < 1:798626/52‑1 (MQ=255) acacGGAG‑ATTGCAATATATTGAATTTGCAAAGTTTTGTAGGCCGGAt < 1:575918/48‑1 (MQ=255) acacGGAG‑ATTGCAATATATTGAATTTGCAAAGTTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCgg > 1:2866669/1‑70 (MQ=255) cacGGAG‑ATTGCAATATATTGAATTTGCAAAGTTTTGTAGGCCGGATAAGGCGTTCa > 1:2585801/1‑57 (MQ=255) cacGGAG‑ATTGCAATATATTGAATTTGCAAAGTTTTGTAGGCCGGATAAGGCGTTCAcgccgc > 1:3393738/1‑63 (MQ=255) cacGGAG‑ATTGCAATATATTGAATTTGCAAAGTTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCgg < 1:307968/69‑1 (MQ=255) cacGGAG‑ATTGCAATATATTGAATTTGCAAAGTTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCa < 1:413119/71‑1 (MQ=255) cacGGAG‑ATTGCAATATATTGAATTTGCAAAGTTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCa > 1:829181/1‑71 (MQ=255) gag‑aTTGCAATATATTGAATTTGCAAAGTTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGaa > 1:935588/1‑71 (MQ=255) | CCCGAAGGAAGCCTTGATGATCATAAACGAAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACACGGAG‑ATTGCAATATATTGAATTTGCAAAGTTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAA > minE/1716767‑1716905 |
| 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: ‑ |