Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | REL606 | 4,393,061 | 0 | T | G | 52.0% | 40.2 / 81.1 | 50 | V211G (GTG→GGG) | aidB | isovaleryl CoA dehydrogenase |
Reads supporting (aligned to +/- strand): ref base T (4/20); new base G (26/0); total (30/20) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 2.25e-10 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.45e-11 | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
AAAGCAGGGCGGTTCCGATGTTATGAGCAACACCACCCGTGCAGAGCGTCTGGAAGATGGCTCTTATCGGCTGGTGGGGCATAAATGGTTTTTCTCGGTTCCGCAAAGCGATGCGCATCTGGTGCTGGCGCAGACCGCGGGTGGTCT > REL606/4392987‑4393133 | aaaGCAGGGCGGTTCCGATGTTATGAGCAACACCCCCCGTGCAGAGCGTCTGGAAGATGGCTCTTATCGGCtggtg < 1:521273/76‑1 (MQ=255) aaaGCAGGGCGGTTCCGATGTTATGAGCAACACCACCCGTGCAGAGCGTCTGGAAGATGGCTCTTATCGGCtggtg < 1:249139/76‑1 (MQ=255) aaaGCAGGGCGGTTCCGATGTTATGAGCAACACCACCCGTGCAGAGCGTCTGGAAGATGGCTCTTATCGGCtggtg < 2:398276/76‑1 (MQ=255) aaGCAGGGCGGTTCCGATGTTATGAGCAACACCACCCGTGCAGAGCGTCTGGAAGATGGCTCTTATCGGCtggtgg > 1:2017068/1‑76 (MQ=255) gCAGGGCGGTTCCGATGTTATGAGCAACACCACCCGTGCAGAGCGTCTGGAAGATGGCTCTTATCGGCTGGTgggg < 2:1181930/76‑1 (MQ=255) cAGGGCGGTTCCGATGTTATGAGCAACACCACCCGTGCAGAGCGTCTGGAAGATGGCTCTTATCGGCTGggggggg > 2:67864/1‑75 (MQ=255) cAGGGCGGTTCCGATGTTATGAGCAACACCACCCGTGCAGAGCGTCTGGAAGATGGCTCTTATCGGCTGGTGGGGc < 2:150619/76‑1 (MQ=255) aGGGCGGTTCCGATGTTATGAGCAACACCACCCGTGCAGAGCGTCTGGAAGATGGCTCTTATCGGCTGgggggggg > 1:543672/1‑74 (MQ=255) cGGTTTCGATGTTATGAGCAACACCACCCGTGCAGAGCGTCTGGAAGATGGCTCTTATCGGCTGGGGGGGGGTaaa > 2:1977901/1‑76 (MQ=255) cGGTTCCGATGTTATGAGCAACACCACCCGTGCAGAGCGTCTGGAAGATGGCTCTTATCGGCTGGTGGGGCATaaa < 2:1546866/76‑1 (MQ=255) cGGTTCCGATGTTATGAGCAACACCACCCGTGCAGAGCGTCTGGAAGATGGCTCTTATCGGCTGGTGGGGCATaaa < 2:552950/76‑1 (MQ=255) tCCGATGTTATGAGCAACACCACCCGTGCAGAGCGTCTGGAAGATGGCTCTTATCGGCTGGGGGGGGATAAATGGt > 1:277495/1‑76 (MQ=255) ccGATGTTATGAGCAACACCACCCGTGCAGAGCGTCTGGAAGATGGCTCTTATCGGCTGGTGGGGCATAAATGGtt < 1:466394/76‑1 (MQ=255) gATGTTATGAGCAACACCACCCGTGCAGAGCGTCTGGAAGATGGCTCTTATCGGCTGGTGGGGCATAAATGGtttt < 1:1683432/76‑1 (MQ=255) gATGTTATGAGCAACACCACCCGTGCAGAGCGTCTGGAAGATGGCTCTTATCGGCTGGTGGGGCATAAATGGtttt < 1:337003/76‑1 (MQ=255) gATGTTATGAGCAACACCACCCGTGCAGAGCGTCTGGAAGATGGCTCTTATCGGCTGGGGGGGCCTAAAAGGgttt > 1:1921391/1‑76 (MQ=255) aTGTTATGAGCAACACCACCCGTGCAGAGCGTCTGGAAGATGGCTCTTATCGGCTGGTGGGGCATAAATGGttttt < 2:2019771/76‑1 (MQ=255) aTGTTATGAGCAACACCACCCGTGCAGAGCGTCTGGAAGATGGCTCTTATCGGCTGGGGGGGGATAAAAGGttttt > 2:440047/1‑76 (MQ=255) gTTATGAGCAACACCACCCGTGCAGAGCGTCTGGAAGATGGCTCTTATCGGCTGGGGGGGGGTAAAAGGTTTTtct > 2:1664160/1‑76 (MQ=255) gTTATGAGCAACACCACCCGTGCAGAGCGTCTGGAAGATGGCTCTTATCGGCTGGGGGGGGATAAATGGttttttt > 1:1735236/1‑74 (MQ=255) ttATGAGCAACACCACCCGTGCAGAGCGTCTGGAAGATGGCTCTTATCGGCTGGTGGGGCATAAATGGTTTTtctc < 1:1286104/76‑1 (MQ=255) tATGAGCAACACCACCCGTGCAGAGCGTCTGGAAGATGGCTCTTATCGGCTGGGGGGGGATAAATGGTTTTTTTCg > 1:1892868/1‑76 (MQ=255) aTGAGCAACACCACCCGTGCAGAGCGTCTGGAAGATGGCTCTTATCGGCTGGGGGGGGATAAATGGGTTTTTTCgg > 1:1584793/1‑76 (MQ=255) acacCACCCGTGCAGAGCGTCTGGAAGATGGCTCTTATCGGCTGGTGGGGCATAAATGGTTTTTCTCGGTTCCGCa < 2:1548704/76‑1 (MQ=255) acacCACCCGTGCAGAGCGTCTGGAAGATGGCTCTTATCGGCTGGTGGGGCATAAATGGTTTTTCTCGGTTCCGCa < 1:413576/76‑1 (MQ=255) accaccCGTGCAGAGCGTCTGGAAGATGGCTCTTATCGGCTGGTGGGGCATAAATGGTTTTTCTCGGTTCCGCaaa < 1:32604/76‑1 (MQ=255) caccCGTGCAGAGCGTCTGGAAGATGGCTCTTATCGGCTGGTGGGGCATAAATGGTTTTTCTCGGTTCCGCAAAGc < 1:1390813/76‑1 (MQ=255) caccCGTGCAGAGCGTCTGGAAGATGGCTCTTATCGGCTGGGGGGGGATAAATGGGTTTTTTCGGTTCCGCaaaac > 1:1581858/1‑74 (MQ=255) accCGTGCAGAGCGTCTGGAAGATGGCTCTTATCGGCTGGTGGGGCATAAATGGTTTTTCTCGGTTCCGCAAAGCg < 2:767999/76‑1 (MQ=255) accCGTGCAGAGCGTCTGGAAGATGGCTCTTATCGGCTGGGGGGGGATAAATGGGTTTTTTCGGTTCCGCAAAGCg > 1:1040693/1‑76 (MQ=255) cGTGCAGAGCGTCTGGAAGATGGCTCTTATCGGCTGGGGGGGGGTAAATGGTTTTTTTCGGTTCCGCAAAGCGATg > 2:1673406/1‑76 (MQ=255) gTGCAGAGCGTCTGGAAGATGGCTCTTATCGGCTGGGGGGGGATAAATGGTTTTTCTTGGTTCCGCAAAGCGATgc > 2:194792/1‑76 (MQ=255) tGCAGAGCGTCTGGAAGATGGCTCTTATCGGCTGGGGGGGGATAAATGGGTTTTTTCGGTTTCGCCAAGCGATgcg > 1:1622280/1‑76 (MQ=255) gCAGAGCGTCTGGAAGATGGCTCTTATCGGCTGGGGGGGCGTAAATGGTTTTTTTCTGTTCCGCACAGCGATgcgc > 1:1219473/1‑76 (MQ=255) cAGAGCGTCTGGAAGATGGCTCTTATCGGCTGGGGGGGGCTAAATGGGTTTTTTCGGTTCCCCAAAGCGATgggcc > 2:539590/1‑75 (MQ=255) agagCGTCTGGAAGATGGCTCTTATCGGCTGGTGGGGCATAAATGGTTTTTCTCGGTTCCGCAAAGCGATGCGCAt < 1:1782336/76‑1 (MQ=255) agagCGTCTGGAAGATGGCTCTTATCGGCTGGGGGGGGATAAATGGTTTTTTTCGGTTCCGCAAAGCGATGCGCAt > 2:654037/1‑76 (MQ=255) agagCGTCTGGAAGATGGCTCTTATCGGCTGGGGGGGGATAAATGGGTTTTTTCGGTTTCGCCAAACGATGCGCAt > 1:1404771/1‑76 (MQ=255) agagCGTCTGGAAGATGGCTCTTATCGGCTGGGGGGGCCTAAATGGTTTTTCTCGGTTCCGCAAAGCGATGCGCAt > 1:1877322/1‑76 (MQ=255) gagCGTCTGGAAGATGGCTCTTATCGGCTGGGGGGGGATAAATGGGTTTTTTCGGTTTCGCCAAGCGATGCGCATc > 1:1143218/1‑76 (MQ=255) gTCTGGAAGATGGCTCTTATCGGCTGGGGGGGCCTAAATGGTTTTTCTCGGTTCCGCCAAGCGATGCGCATCTggg > 2:633057/1‑75 (MQ=255) tCTGGAAGATGGCTCTTATCGGCTGGTGGGGCATAAATGGGTTTTTTTGGTTCCGCAAAGCGATGCGCATCTgggg > 2:1355274/1‑74 (MQ=255) tCTGGAAGATGGCTCTTATCGGCTGGGGGGGGGTAAATGGGTTTTTTCGGTTCCGCCAAGCGAGGCGCACCTgggg > 2:1328996/1‑74 (MQ=255) cTGGAAGATGGCTCTTATCGGCTGGGGGGGGATAAATGGTTTTTTTCGGGTCCGCAAAGCGATGCGCATCTGGGGc > 2:1687199/1‑76 (MQ=255) cTGGAAGATGGCTCTTATCGGCTGGGGGGGGATAAATGGTTTTTCTCGGTTCCGCCAAGCGATGCGCATCTGGTGc > 1:1790721/1‑76 (MQ=255) ggAAGATGGCTCTTATCGGCTGGGGGGGCATAAATGGTTTTTCTCGGTTCCGCCAAGCGCTGCGCATCTGGTGCTg > 1:1405049/1‑76 (MQ=255) ggAAGATGGCTCTTATCGGCTGGGGGGGCATAAATGGTTTTTCTCGGTTCCGCCAAGCGATGCGCATCTGGGGCTg > 2:148940/1‑76 (MQ=255) aaGATGGCTCTTATCGGCTGGGGGGGCATAAATGGGTTTTTTCGGGTCCGCAAAGCGATGCGCATCTGGTGCTggg > 1:674777/1‑75 (MQ=255) gATGGCTCTTATCGGCTGGTGGGGCATAAATGGTTTTTCTCGGTTCCGCAAAGCGATGCGCATCTGGTGCTGgcgc < 1:573359/76‑1 (MQ=255) ggCTCTTATCGGCTGGTGGGGCATAAATGGTTTTTCTCGGTTCCGCAAAGCGATGCGCATCTGGTGCTGGCGCAGa < 2:481128/76‑1 (MQ=255) ggCTCTTATCGGCTGGTGGGGCATAAATGGTTTTTCTCGGTTCCGCAAAGCGATGCGCATCTGGTGCTGGCGCAGa < 1:910725/76‑1 (MQ=255) ggCTCTTATCGGCTGGTGGGGCATAAATGGTTTTTCTCGGTTCCGCAAAGCGATGCGCATCTGGTGCTGGCGCAGa < 1:211411/76‑1 (MQ=255) ggCTCTTATCGGCTGGGGGGGCATAAATGGTTTTTCTCGGTTCCGCAAAGCGATGCGCATCTGGTGCTGGCGCAgc > 1:1990442/1‑75 (MQ=255) tgtCGGCTGGTGGGGCATAAATGGTTTTTCTCGGTTCCGCAAAGCGATGCGCATCTGGTGCTGGCGCAGACCGCgg < 1:980930/74‑1 (MQ=255) tATCGGCTGGTGGGGCATAAATGGTTTTTCTCGGTTCCGCAAAGCGATGCGCATCTGGTGCTGGCGCAGACcgcgn < 1:666789/76‑2 (MQ=255) cGGCTGGTGGGGCATAAATGGTTTTTCTCGGTTCCGCAAAGCGATGCGCATCTGGTGCTGGCGCAGACCGCGggtg > 2:152497/1‑76 (MQ=255) cGGCTGGTGGGGCATAAATGGTTTTTCTCGGTTCCGCAAAGCGATGCGCATCTGGTGCTGGCGCAGACCGCGggtg > 2:966901/1‑76 (MQ=255) ggCTGGTGGGGCATAAATGGTTTTTCTCGGTTCCGCAAAGCGATGCGCATCTGGTGCTGGCGCAGACCGCgggggg > 1:478133/1‑76 (MQ=255) tggtggGGCATAAATGGTTTTTCTCGGTTCCGCAAAGCGATGCGCATCTGGTGCTGGCGCAGACCGCGGGTGGTCt > 1:670160/1‑76 (MQ=255) | AAAGCAGGGCGGTTCCGATGTTATGAGCAACACCACCCGTGCAGAGCGTCTGGAAGATGGCTCTTATCGGCTGGTGGGGCATAAATGGTTTTTCTCGGTTCCGCAAAGCGATGCGCATCTGGTGCTGGCGCAGACCGCGGGTGGTCT > REL606/4392987‑4393133 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 19 ≤ ATCG/ATCG < 28 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |