Predicted mutation | ||||||
---|---|---|---|---|---|---|
evidence | seq id | position | mutation | annotation | gene | description |
RA | NC_000913 | 4,616,491 | T→C | D192G (GAT→GGT) | yjjI ← | DUF3029 family protein, putative glycine radical enzyme |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_000913 | 4,616,491 | 0 | T | C | 89.1% | 154.1 / 9.0 | 55 | D192G (GAT→GGT) | yjjI | DUF3029 family protein, putative glycine radical enzyme |
Reads supporting (aligned to +/- strand): ref base T (5/1); new base C (28/21); total (33/22) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 3.84e-01 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.52e-01 |
ATCATGCACCGGACCGTTGGCGATGTGCGGTTTGCTACATTCACAGATGTTCTTCGCCACTTCCAGCAGCAGGTCATCAGGGGTGATTTCAGGATCGTAGATAAAGGTCAGGTTCGGTGAAACCTGCTTCAACTCTGCATCTGCACGTAAGATCGCGCGGGTAATGGGCGAATCAGACGGGCCGATATTGG > NC_000913/4616397‑4616587 | aTCATGCACCGGACCGTTGGCGATGTGCGGTTTGCTACATTCACAGATGTTCTTCGCCACTTCCAGCAGCAGGTCATCAGGGGTGATTTCAGGATCGTAg < 2:649102/100‑1 (MQ=255) tCATGCACCGGACCGTTGGCGATGTGCGGTTTGCTACATTCACAGATGTTCTTCGCCACTTCCAGCAGCAGGTCATCAGGGGTGATTTCAGGACCGTAGa < 1:198460/100‑1 (MQ=255) aTGCACCGGACCGTTGGCGATGTGCGGTTTGCTACATTCACAGATGTTCTTCGCCACTTCCAGCAGCAGGTCATCAGGGGTGATTTCAGGACCGTAGATa > 1:670664/1‑100 (MQ=255) gCACCGGACCGTTGGCGATGTGCGGTTTGCTACATTCACAGATGTTCTTCGCCACTTCCAGCAGCAGGTCATCAGGGGTGATTTCAGGACCGTAGATaaa > 2:457783/1‑100 (MQ=255) cACCGGACCGTTGGCGATGTGCGGTTTGCTACATTCACAGATGTTCTTCGCCACTTCCAGCAGCAGGTCATCAGGGGTGATTTCAGGACCGTAGATAAAg < 1:1129251/100‑1 (MQ=255) cACCGGACCGTTGGCGATGTGCGGTTTGCTACATTCACAGATGTTCTTCGCCACTTCCAGCAGCAGGTCATCAGGGGTGATTTCAGGACCGTAGATAAAg < 1:1129315/100‑1 (MQ=255) aCCGGACCGTTGGCGATGTGCGGTTTGCTACATTCACAGATGTTCTTCGCCACTTCCAGCAGCAGGTCATCAGGGGTGATTTCAGGACCGTAGATAAAgg < 2:725141/100‑1 (MQ=255) aCCGGACCGTTGGCGATGTGCGGTTTGCTACATTCACAGATGTTCTTCGCCACTTCCAGCAGCAGGTCATCAGGGGTGATTTCAGGACCGTAGATAAAgg < 2:725134/100‑1 (MQ=255) cGGACCGTTGGCGATGTGCGGTTTGCTACATTCACAGATGTTCTTCGCCACTTCCAGCAGCAGGTCATCAGGGGTGATTTCAGGACCGTAGATAAAGGt > 2:402777/1‑99 (MQ=255) ggACCGTTGGCGATGTGCGGTTTGCTACATTCACAGATGTTCTTCGCCACTTCCAGCAGCAGGTCATCAGGGGTGATTTCAGGACCGTAGATAAAGGTCa > 1:924121/1‑100 (MQ=255) ggACCGTTGGCGATGTGCGGTTTGCTACATTCACAGATGTTCTTCGCCACTTCCAGCAGCAGGTCATCAGGGGTGATTTCAGGACCGTAGATAAAGGTCa < 1:689235/100‑1 (MQ=255) ccGTTGGCGATGTGCGGTTTGCTACATTCACAGATGTTCTTCGCCACTTCCAGCAGCAGGTCATCAGGGGTGATTTCAGGACCGTAGATAAAGGTCAGGt < 1:651776/100‑1 (MQ=255) ggCGATGTGCGGTTTGCTACATTCACAGATGTTCTTCGCCACTTCCAGCAGCAGGTCATCAGGGGTGATTTCAGGACCGTAGATAAAGGTCAGGTTCGGt < 1:1276264/100‑1 (MQ=255) gCGATGTGCGGTTTGCTACATTCACAGATGTTCTTCGCCACTTCCAGCAGCAGGTCATCAGGGGTGATTTCAGGACCGTAGATAAAGGTCAGGTTCGGTg > 2:104781/1‑100 (MQ=255) tgCGGTTTGCTACATTCACAGATGTTCTTCGCCACTTCCAGCAGCAGGTCATCAGGGGTGATTTCAGGATCGTAGATAAAGGTCAGGTTCGGTGAAACCt > 1:1298532/1‑100 (MQ=255) tgCGGTTTGCTACATTCACAGATGTTCTTCGCCACTTCCAGCAGCAGGTCATCAGGGGTGATTTCAGGACCGTAGATAAAGGTCAGGTTCGGTGAAACCt > 2:705324/1‑100 (MQ=255) tgCGGTTTGCTACATTCACAGATGTTCTTCGCCACTTCCAGCAGCAGGTCATCAGGGGTGATTTCAGGACCGTAGATAAAGGTCAGGTTCGGTGAAACCt > 1:1117576/1‑100 (MQ=255) gTTTGCTACATTCACAGATGTTCTTCGCCACTTCCAGCAGCAGGTCATCAGGGGTGATTTCAGGACCGTAGATAAAGGTCAGGTTCGGTGAAACCTGCtt > 1:1290754/1‑100 (MQ=255) gTTTGCTACATTCACAGATGTTCTTCGCCACTTCCAGCAGCAGGTCATCAGGGGTGATTTCAGGACCGTAGATAAAGGTCAGGTTCGGTGAAACCTGCtt > 2:656977/1‑100 (MQ=255) tttGCTACATTCACAGATGTTCTTCGCCACTTCCAGCAGCAGGTCATCAGGGGTGATTTCAGGACCGTAGATAAAGGTCAGGTTCGGTGAAACCTGCTTc < 2:939124/100‑1 (MQ=255) tttGCTACATTCACAGATGTTCTTCGCCACTTCCAGCAGCAGGTCATCAGGGGTGATTTCAGGACCGTAGATAAAGGTCAGGTTCGGTGAAACCTGCTTc > 2:749659/1‑100 (MQ=255) ttGCTACATTCACAGATGTTCTTCGCCACTTCCAGCAGCAGGTCATCAGGGGTGATTTCAGGACCGTAGATAAAGGTCAGGTTCGGTGAAACCTGCTTcn < 1:540407/100‑2 (MQ=255) ttGCTACATTCACAGATGTTCTTCGCCACTTCCAGCAGCAGGTCATCAGGGGTGATTTCAGGACCGTAGATAAAGGTCAGGTTCGGTGAAACCTGCTTCa > 2:1275523/1‑100 (MQ=255) ttGCTACATTCACAGATGTTCTTCGCCACTTCCAGCAGCAGGTCATCAGGGGTGATTTCAGGACCGTAGATAAAGGTCAGGTTCGGTGAAACCTGCTTCa > 1:62131/1‑100 (MQ=255) cTACATTCACAGATGTTCTTCGCCACTTCCAGCAGCAGGTCATCAGGGGTGATTTCAGGACCGTAGATAAAGGTCAGGTTCGGTGAAACCTGCTTCAAct > 2:1152174/1‑100 (MQ=255) tACATTCACAGATGTTCTTCGCCACTTCCAGCAGCAGGTCATCAGGGGTGATTTCAGGACCGTAGATAAAGGTCAGGTTCGGTGAAACCTGCTTCAActc > 1:300265/1‑100 (MQ=255) aCATTCACAGATGTTCTTCGCCACTTCCAGCAGCAGGTCATCAGGGGTGATTTCAGGACCGTAGATAAAGGTCAGGTTCGGTGAAACCTGCTTCAActct > 2:1047222/1‑100 (MQ=255) aTTCACAGATGTTCTTCGCCACTTCCAGCAGCAGGTCATCAGGGGTGATTTCAGGATCGTAGATAAAGGTCAGGTTCGGTGAAACCTGCTTCAACtctgc > 2:831689/1‑100 (MQ=255) aTTCACAGATGTTCTTCGCCACTTCCAGCAGCAGGTCATCAGGGGTGATTTCAGGACCGTAGATAAAGGTCAGGTTCGGTGAAACCTGCTTCAACtctgc > 1:106011/1‑100 (MQ=255) ttCACAGATGTTCTTCGCCACTTCCAGCAGCAGGTCATCAGGGGTGATTTCAGGACCGTAGATAAAGGTCAGGTTCGGTGAAACCTGCTTCAACtctgca > 2:413993/1‑100 (MQ=255) ttCACAGATGTTCTTCGCCACTTCCAGCAGCAGGTCATCAGGGGTGATTTCAGGACCGTAGATAAAGGTCAGGTTCGGTGAAACCTGCTTCAACtctgca < 1:1033771/100‑1 (MQ=255) gTTCTTCGCCACTTCCAGCAGCAGGTCATCAGGGGTGATTTCAGGATCGTAGATAAAGGTCAGGTTCGGTGAAACCTGCTTCAACTCTGCATCTGCACGt > 1:585203/1‑100 (MQ=255) cttcGCCACTTCCAGCAGCAGGTCATCAGGGGTGATTTCAGGACCGTAGATAAAGGTCAGGTTCGGTGAAACCTGCTTCAACTCTGCATCTGCACGTAAg > 2:263805/1‑100 (MQ=255) gCCACTTCCAGCAGCAGGTCATCAGGGGTGATTTCAGGACCGTAGATAAAGGTCAGGTTCGGTGAAACCTGCTTCAACTCTGCATCTGCACGTAAGATcg > 2:3659/1‑100 (MQ=255) cTTCCAGCAGCAGGTCATCAGGGGTGATTTCAGGACCGTAGATAAAGGTCAGGTTCGGTGAAACCTGCTTCAACTCTGCATCTGCACGTAAGATcgcgcg < 2:1003791/100‑1 (MQ=255) tCCAGCAGCAGGTCATCAGGGGTGATTTCAGGACCGTAGATAAAGGTCAGGTTCGGTGAAACCTGCTTCAACTCTGCATCTGCACGTAAGATCGCGCggg > 2:1150936/1‑100 (MQ=255) tCCAGCAGCAGGTCATCAGGGGTGATTTCAGGACCGTAGATAAAGGTCAGGTTCGGTGAAACCTGCTTCAACTCTGCATCTGCACGTAAGATCGCGCggg < 2:260905/100‑1 (MQ=255) ccAGCAGCAGGTCATCAGGGGTGATTTCAGGACCGTAGATAAAGGTCAGGTTCGGTGAAACCTGCTTCAACTCTGCATCTGCACGTAAGATCGCGCGGGt > 1:754401/1‑100 (MQ=255) gcagcagGTCATCAGGGGTGATTTCAGGACCGTAGATAAAGGTCAGGTTCGGTGAAACCTGCTTCAACTCTGCATCTGCACGTAAGATCGCGCGGGTAAt < 1:522472/100‑1 (MQ=255) gcagGTCATCAGGGGTGATTTCAGGATCGTAGATAAAGGTCAGGTTCGGTGAAACCTGCTTCAACTCTGCATCTGCACGTAAGATCGCGCGGGTAATggg > 2:644244/1‑100 (MQ=255) cagGTCATCAGGGGTGATTTCAGGACCGTAGATAAAGGTCAGGTTCGGTGAAACCTGCTTCAACTCTGCATCTGCACGTAAGATCGCGCGGGTAATGGGc < 2:168014/100‑1 (MQ=255) cagGTCATCAGGGGTGATTTCAGGACCGTAGATAAAGGTCAGGTTCGGTGAAACCTGCTTCAACTCTGCATCTGCACGTAAGATCGCGCGGGTAATGGGc < 2:167652/100‑1 (MQ=255) agGTCATCAGGGGTGATTTCAGGACCGTAGATAAAGGTCAGGTTCGGTGAAACCTGCTTCAACTCTGCATCTGCACGTAAGATCGCGCGGGTAATGGGCg > 1:577082/1‑100 (MQ=255) agGTCATCAGGGGTGATTTCAGGACCGTAGATAAAGGTCAGGTTCGGTGAAACCTGCTTCAACTCTGCATCTGCACGTAAGATCGCGCGGGTAATGGGCg < 1:635465/100‑1 (MQ=255) tcatcaGGGGTGATTTCAGGACCGTAGATAAAGGTCAGGTTCGGTGAAACCTGCTTCAACTCTGCATCTGCACGTAAGATCGCGCGGGTAATGGGCGAAt < 1:543821/100‑1 (MQ=255) catcaGGGGTGATTTCAGGACCGTAGATAAAGGTCAGGTTCGGTGAAACCTGCTTCAACTCTGCATCTGCACGTAAGATCGCGCGGGTAATGGGCGAATc > 2:91597/1‑100 (MQ=255) gggTGATTTCAGGATCGTAGATAAAGGTCAGGTTCGGTGAAACCTGCTTCAACTCTGCATCTGCACGTAAGATCGCGCGGGTAATGGGCGAATCAGACgg > 1:238528/1‑100 (MQ=255) tGATTTCAGGACCGTAGATAAAGGTCAGGTTCGGTGAAACCTGCTTCAACTCTGCATCTGCACGTAAGATCGCGCGGGTAATGGGCGAATCAGACGGGcc < 2:55885/100‑1 (MQ=255) tGATTTCAGGACCGTAGATAAAGGTCAGGTTCGGTGAAACCTGCTTCAACTCTGCATCTGCACGTAAGATCGCGCGGGTAATGGGCGAATCAGACGGGcc < 1:1150947/100‑1 (MQ=255) aTTTCAGGACCGTAGATAAAGGTCAGGTTCGGTGAAACCTGCTTCAACTCTGCATCTGCACGTAAGATCGCGCGGGTAATGGGCGAATCAGACGGGCCGa > 1:522895/1‑100 (MQ=255) aTTTCAGGACCGTAGATAAAGGTCAGGTTCGGTGAAACCTGCTTCAACTCTGCATCTGCACGTAAGATCGCGCGGGTAATGGGCGAATCAGACGGGCCGa > 1:1219804/1‑100 (MQ=255) tttCAGGACCGTAGATAAAGGTCAGGTTCGGTGAAACCTGCTTCAACTCTGCATCTGCACGTAAGATCGCGCGGGTAATGGGCGAATCAGACGGGCCGat > 1:126365/1‑100 (MQ=255) tCAGGACCGTAGATAAAGGTCAGGTTCGGTGAAACCTGCTTCAACTCTGCATCTGCACGTAAGATCGCGCGGGTAATGGGCGAATCAGACGGGCCGatat > 2:245722/1‑100 (MQ=255) aGGACCGTAGATAAAGGTCAGGTTCGGTGAAACCTGCTTCAACTCTGCATCTGCACGTAAGATCGCGCGGGTAATGGGCGAATCAGACGGGCCGATATTg < 1:749664/100‑1 (MQ=255) ggACCGTAGATAAAGGTCAGGTTCGGTGAAACCTGCTTCAACTCTGCATCTGCACGTAAGATCGCGCGGGTAATGGGCGAATCAGACGGGCCGATATTgg > 1:1232328/1‑100 (MQ=255) | ATCATGCACCGGACCGTTGGCGATGTGCGGTTTGCTACATTCACAGATGTTCTTCGCCACTTCCAGCAGCAGGTCATCAGGGGTGATTTCAGGATCGTAGATAAAGGTCAGGTTCGGTGAAACCTGCTTCAACTCTGCATCTGCACGTAAGATCGCGCGGGTAATGGGCGAATCAGACGGGCCGATATTGG > NC_000913/4616397‑4616587 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |