Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | REL606 | 4,393,061 | 0 | T | G | 51.4% | 24.1 / 52.2 | 35 | V211G (GTG→GGG) | aidB | isovaleryl CoA dehydrogenase |
Reads supporting (aligned to +/- strand): ref base T (2/15); new base G (18/0); total (20/15) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 4.19e-08 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.55e-08 | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
CAGGGCGGTTCCGATGTTATGAGCAACACCACCCGTGCAGAGCGTCTGGAAGATGGCTCTTATCGGCTGGTGGGGCAT‑AAATGGTTTTTCTCGGTTCCGCAAAGCGATGCGCATCTGGTGCTGGCGCAGACCGCGGGTGGTCTGTC > REL606/4392991‑4393136 | cAGGGCGGTTCCGATGTTATGAGCAACACCACCCGTGCAGAGCGTCTGGAAGATGGCTCTTATCGGCTGggggggg > 2:1099796/1‑75 (MQ=255) gggCGGTTCCGATGTTATGAGCAACACCACCCGTGCAGAGCGTCTGGAAGATGGCTCTTATCGGCTGggggggggt > 1:818028/1‑73 (MQ=255) gggCGGTTCCGATGTTATGAGCAACACCACCCGTGCAGAGCGTCTGGAAGATGGCTCTTATCGGCTGggggggggg > 1:1483188/1‑73 (MQ=255) ggTTCCGATGTTATGAGCAACACCACCCGTGCAGAGCGTCTGGAAGATGGCTCTTATCGGCTGGGGGGGGGT‑aaaa > 2:903842/1‑75 (MQ=255) tCCGATGTTATAAGCAACACCACCCGTGCAGAGCGTCTGGAAGATGGCTCTTATCGGCTGGGGGGGGGT‑AAATggg > 2:415979/1‑75 (MQ=255) gATGTTATGAGCAACACCACCCGTGCAGAGCGTCTGGAAGATGGCTCTTATTGGCTGGTGGGGCAT‑AAATGGtttt < 1:956693/76‑1 (MQ=255) gATGTTATGAGCAACACCACCCGTGCAGAGCGTCTGGAAGATGGCTCTTATCGGCTGGTGGGGCAT‑AAATGGtttt < 1:866162/76‑1 (MQ=255) ttATGAGCAACACCACCCGTGCAGAGCGTCTGGAAGATGGCTCTTATCGGCTGGTGGGGCAT‑AAATGGTTTTtctc < 1:936640/76‑1 (MQ=255) cAACACCACCCGTGCAGAGCGTCTGGAAGATGGCTCTTATCGGCTGGTGGGGCAT‑AAATGGTTTTTCTCGGTTCCg < 2:1433255/76‑1 (MQ=255) aaCACCACCCGTGCAGAGCGTCTGGAAGATGGCTCTTATCGGCTGGGGGGGGGTAAAATGGTTTTTTTGGTTTccc > 2:253116/1‑75 (MQ=255) accgccCGTGCAGAGCGTCTGGAAGATGGCTCTTATCGGCTGGTGGGGCAT‑AAATGGTTTTTCTCGGTTCCGCaaa < 2:448108/76‑1 (MQ=255) ccaccCGTGCAGAGCGTCTGGAAGATGGCTTTTATCGGCTGGGGGGGGAG‑AAATGGGTTTTTTCCGGTCCGCagag > 2:1159505/1‑76 (MQ=255) ccaccCGTGCAGAGCGTCTGGAAGATGGCTCTTATCGGCTGGTGGGGCAT‑AAATGGTTTTTCTCGGTTCCGCAAAg < 2:973821/76‑1 (MQ=255) ccaccCGTGCAGAGCGTCTGGAAGATGGCTCTTATCGGCTGGGGGGGCAT‑AAATGGTTTTTCTCGGTTCCGCAAAg > 2:457135/1‑76 (MQ=255) cGTGCAGAGCGTCTGGAAGATGGCTCTTATCGGCTGGTGGGGCAT‑AAATGGTTTTTCTCGGTTCCGCAAAGCGATg < 2:910011/76‑1 (MQ=255) cGTGCAGAGCGTCTGGAAGATGGCTCTTATCGGCTGGGGGGGGAT‑AAATGGTTTTTTTCCGTTCCGCAAAGCGCTg > 2:963655/1‑76 (MQ=255) gTGCAGAGCGTCTGGAAGATGGCTCTTATCGGCTGGGGGGGGAT‑AAATGGGTTTTCTCGGGTCCGCAAAACGCTgc > 1:404686/1‑76 (MQ=255) tGCAGAGCGTCTGGAAGATGGCTCTTATCGGCTGGTGGGGCAT‑AAATGGTTTTTCTCGGTTCCGCAAAGCGATgcg < 1:497984/76‑1 (MQ=255) gagCGTCTGGAAGATGGCTCTTATCGGCTGGGGGGGGAT‑AAATGGGTTTTTTCGGTTCCGCCAAGCGCTGCGCact > 2:911560/1‑74 (MQ=255) tGGAAGATGGCTCTTATCGGCTGGTGGGGCAT‑AAATGGTTTTTTTCGGTTCCGCAAAACGATGCGCCTCTGGTggc > 2:1036246/1‑74 (MQ=255) tGGAAGATGGCTCTTATCGGCTGGTGGGGCAT‑AAATGGTTTTTCTCGGTTCCGCAAAGCGATGCGCATCTGGTGCt < 2:795426/76‑1 (MQ=255) tGGAAGATGGCTCTTATCGGCTGGGGGGGCAT‑AAATGGTTTTTCTCGGTTCCGCAAAGCGATGCGCCTCTGGTGCt > 2:1471243/1‑76 (MQ=255) ggAAGATGGCTCTTATCGGCTGGGGGGGCAT‑AAATGGTTTTTCTCGGTTCCGCCAAGCGATGCGCATCTGGTGCTg > 2:1175994/1‑76 (MQ=255) aaGATGGCTCTTATCGGCTGGGGGGGCAT‑AAATGGGTTTTCTCGGTTCCGCAAAGCGATGCGCATCTGGGGCTGgc > 2:1373712/1‑76 (MQ=255) aGATGGCTCTTATCGGCTGGTGGGGCAT‑AAAGGGTTTTTCTCGGTTCCGCAAAGCGATGCGCATCTGGTGCTGgcg < 2:836793/76‑1 (MQ=255) tGGCTCTTATCGGCTGGTGGGGCAT‑AAATGGTTTTTCTCGGTTCCGCAAAGCGATGCGCATCTGGTGCTGGCGCAg < 2:288791/76‑1 (MQ=255) tGGCTCTTATCGGCTGGTGGGGCAT‑AAATGGTTTTTCTCGGTTCCGCAAAGCGATGCGCATCTGGTGCTGGCGCAg < 1:632286/76‑1 (MQ=255) tGGCTCTTATCGGCTGGGGGGGGAT‑AAATGGTTTTTCTCGGTTCCGCAAAGCGATGCGCATCTGGTGCTGgcgccg > 2:1130343/1‑74 (MQ=255) tGGCTCTTATCGGCTGGGGGGGGAT‑AAATGGGTTTTTTCGGTTCCGCAAAGCGGTGCGCATCTGGTGCTGGCGCAg > 2:1023611/1‑76 (MQ=255) tGGCTCTTATCGGCTGGGGGGGGAT‑AAATGGGTTTTTTCGGTTCCGCAAAGCGATGCGCATCTGGGGCTGGCGCAg > 2:733421/1‑76 (MQ=255) ggCTCTTATCGGCTGGTGGGGCAT‑AAATGGTTTTTCTCGGTTCCGCAAAGCGATGCGCATCTGGTGCTGGCGCAGa < 2:1544315/76‑1 (MQ=255) ggCTCTTATCGGCTGGTGGGGCAT‑AAATGGTTTTTCTCGGTTCCGCAAAGCGATGCGCATCTGGTGCTGGCGCAGa < 2:625297/76‑1 (MQ=255) tctTATCGGCTGGGGGGGCAT‑AAATGGTTTTTTTCGGTTCCGCAAAGCGATGCGCATCTGGTGCTGGCGCCGACcg > 2:601094/1‑76 (MQ=255) ctTATCGGCTGGGGGGGGAT‑AAATGGGTTTTCTCGGTTCCGCCAAGCGATGCGCGTCTGGGTCTGGCGCAGAcccc > 2:1606143/1‑74 (MQ=255) ttATCGGCTGGTGGGGCAT‑AAATGGTTTTTCTCGGTTCCGCAAAGCGATGCGCATCTGGTGCTGGCGCAGACcgcg < 1:427097/76‑1 (MQ=255) ttATCGGCTGGGGGGGCAT‑AAATGGGTTTTCTCGGTTCCGCAAAGCGATGCGCATCTGGTGCTGGCGCAGACcgcg > 2:1205002/1‑76 (MQ=255) ttATCGGCTGGGGGGGCAT‑AAATGGGTTTTCTCGGGTCCGCAAAGCGATGCGCATCTGGGGCTGGCGCCGGCcgcg > 2:1488254/1‑76 (MQ=255) gCTGGTGGGGCAT‑AAATGGTTTTTCTCGGTTCCGCAAAGCGATGCGCATCTGGTGCTGGCGCAGACCGCGggtggt > 1:203529/1‑76 (MQ=255) tggtggGGCAT‑AAATGGTTTTTCTCGGTTCCGCAAAGCGATGCGCATCTGGTGCTGGCGCAGACCGCGGGTGGTCt < 2:504461/76‑1 (MQ=255) tggtggGGCAT‑AAATGGTTTTTCTCGGTTCCGCAAAGCGATGCGCATCTGGTGCTGGCGCAGACCGCGGGTGGTCt < 2:415630/76‑1 (MQ=255) tggtggGGCAT‑AAATGGTTTTTCTCGGTTCCGCAAAGCGATGCGCATCTGGTGCTGGCGCAGACCGCGGGTGGTCt < 2:1547303/76‑1 (MQ=255) tggtggGGCAT‑AAATGGTTTTTCTCGGTTCCGCAAAGCGATGCGCATCTGGTGCTGGCGCAGACCGCGGGTGGTCt < 2:987222/76‑1 (MQ=255) ggtggGGCAT‑AAATGGTTTTTCTCGGTTCCGCAAAGCGATGCGCATCTGGTGCTGGCGCAGACCGCGGGTGGTCTg < 1:39833/76‑1 (MQ=255) tggGGCAT‑AAATGGTTTTTCTCGGTTCCGCAAAGCGATGCGCATCTGGTGCTGGCGCAGACCGCGGGTGGTCTGTc > 1:1091437/1‑76 (MQ=255) | CAGGGCGGTTCCGATGTTATGAGCAACACCACCCGTGCAGAGCGTCTGGAAGATGGCTCTTATCGGCTGGTGGGGCAT‑AAATGGTTTTTCTCGGTTCCGCAAAGCGATGCGCATCTGGTGCTGGCGCAGACCGCGGGTGGTCTGTC > REL606/4392991‑4393136 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 19 ≤ ATCG/ATCG < 29 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |