Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | REL606 | 4,393,061 | 0 | T | G | 43.8% | 42.9 / 12.8 | 32 | V211G (GTG→GGG) | aidB | isovaleryl CoA dehydrogenase |
Reads supporting (aligned to +/- strand): ref base T (0/18); new base G (14/0); total (14/18) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 2.12e-09 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
AAAGCAGGGCGGTTCCGATGTTATGAGCAACACCACCCGTGCAGAGCGTCTGGAAGATGGCTCTTATCGGCTGGTGGGGCATAAATGGTTTTTCTCGGTTCCGCAAAGCGATGCGCATCTGGTGCTGGCGCAGACCGCGGGTGGTCTGTC > REL606/4392987‑4393136 | aaaGCAGGGCGGTTCCGATGTTATGAGCAACACCACCCGTGCAGAGCGTCTGGAAGATGGCTCTTATCGGCtggtg < 2:287852/76‑1 (MQ=255) aaaGCAGGGCGGTTCCGATCTTATGAGCAACACCGCCCGTGCAGAGCGTCTGGAAGATGGCTCTTATCGGCtggtg < 2:1189398/76‑1 (MQ=255) cAGGGCGGTTCCGATGTTATGAGCAACACCACCCGTGCAGAGCGTCTGGAAGATGGCTCTTATCGGCTGGTGGGGc < 1:770605/76‑1 (MQ=255) aGGGCGGTTCCGATGTTATGAGCAACACCACCCGTGCAGAGCGTCTGGAAGATGGCTCTTATCGGCTGgggggggg > 1:1251849/1‑74 (MQ=255) gggCGGTTCCGATGTTATGAGCAACACCACCCGTGCAGAGCGTCTGGAAGATGGCTCTTATCGGCTGgggggggaa > 2:20627/1‑73 (MQ=255) gggCGGTTCCGATGTTATGAGCAACACCACCCGTGCAGAGCGTCTGGAAGATGGCTCTTATCGGCTGGGGGGGcgg > 2:234608/1‑74 (MQ=255) ggTTCCGATGTTATGAGCAACACCACCCGTGCAGAGCGTCTGGAAGATGGCTCTTATCGGCTGGTGGGGCATAAAt < 1:333899/76‑1 (MQ=255) ttCCGATGTTATGAGCAACACCACCCGTGCAGAGCGTCTGGAAGATGGCTCTTATCGGCTGGTGGGGCAt < 1:543396/70‑1 (MQ=255) ttCCGATGTTATGAGCAACACCACCCGTGCAGAGCGTCTGGAAGATGGCTCTTATCGGCTGGGGGGGGAt > 2:543396/1‑70 (MQ=255) ttCCGATGTTATGAGCAACACCACCCGTGCAGAGCGTCTGGAAGATGGCTCTTATCGGCTGGGGGGGGATAAAAgg > 2:1587621/1‑76 (MQ=255) tCCGATGTTATGAGCAACACCACCCGTGCAGAGCGTCTGGAAGATGGCTCTTATCGGCTGGGGGGGGGTAAATGGt > 1:1876664/1‑76 (MQ=255) cGATGTTATGAGCAACACCACCCGTGCAGAGCGTCTGGAAGATGGCTCTTATCGGCTGGTGGGGCATAAATGGttt < 2:1069033/76‑1 (MQ=255) gATGTTATGAGCAACACCACCCGTGCAGAGCGTCTGGAAGATGGCTCTTATCGGCTGGTGGGGCATAAATGGtttt < 2:1538218/76‑1 (MQ=255) gATGTTATGAGCAACACCACCCGTGCAGAGCGTCTGGAAGATGGCTCTTATCGGCTGGTGGGGCATAAATGGtttt < 1:868921/76‑1 (MQ=255) aTGTTATGAGCAACACCACCCGTGCAGAGCGTCTGGAAGATGGCTCTTATCGGCTGGTGGGGCATAAATGGttttt < 1:2051437/76‑1 (MQ=255) tGTTATGAGCAACACCACCCGTGCAGAGCGTCTGGAAGATGGCTCTTATCGGCTGGGGGGGGATAAATAGGTTTtc > 2:1990679/1‑76 (MQ=255) tATGAGCAACACCACCCGTGCAGAGCGTCTGGAAGATGGCTCTTATCGGCTGGGGGGGGATAAAAGGGTTTTTTCg > 2:198777/1‑76 (MQ=255) aGCAACACCACCCGTGCAGAGCGTCTGGAAGATGGCTCTTATCGGCTGGTGGGGCATAAATGGTTTTTCTCGGTTc < 2:691442/76‑1 (MQ=255) gCAACACCACCCGTGCAGAGCGTCTGGAAGATGGCTCTTATCGGCTGGTGGGGCATAAATGGTTTTTCTCGGTTcc < 2:282084/76‑1 (MQ=255) acacCACCCGTGCAGAGCGTCTGGAAGATGGCTCTTATCGGCTGGTGGGGCATAAATGGTTTTTCTCGGTTCCGCa < 1:1986616/76‑1 (MQ=255) caccacCCGTGCAGAGCGTCTGGAAGATGGCTCTTATCGGCTGGGGGGGGGTAAATGGTTTTTTTTTGTTTCGCaa > 1:703488/1‑76 (MQ=255) ccaccCGTGCAGAGCGTCTGGAAGATGGCTCTTATCGGCTGGTGGGGCATAAATGGTTTTTCTCGGTTCCGCAAAg < 2:1909817/76‑1 (MQ=255) ccaccCGTGCAGAGCGTCTGGAAGATGGCTCTTATCGGCTGGTGGGGCATAAATGGTTTTTCTCGGTTCCGCAAAg < 2:1292234/76‑1 (MQ=255) caccCGTGCAGAGCGTCTGGAAGATGGCTCTTATCGGCTGGGGGGGGGTAAATGGGGTTTTTCGGGTTCGCAAAGc > 1:1283236/1‑76 (MQ=255) gCAGAGCGTCTGGAAGATGGCTCTTATCGGCTGGGGGGGGATAAATGGGTTTTTTCGGGTCCCCCAAGCGATgcgc > 2:1760330/1‑76 (MQ=255) cAGAGCGTCTGGAAGATGGCTCTTATCGGCTGGTGGGGCATAAATGGTTTTTCTCGGTTCCGCAAAGCGATGCGCa < 2:1035895/76‑1 (MQ=255) agagCGTCTGGAAGATGGCTCTTATCGGCTGGGGGGGGATAAATGGTTTTTTTTGGTTCCGCAAAGCGATGCGCAt > 1:1651166/1‑76 (MQ=255) cGTCTGGAAGATGGCTCTTATCGGCTGGTGGGGCATAAATGGTTTTTCTCGGTTCCGCAAAGCGATGCGCATCTgg < 1:736714/76‑1 (MQ=255) tGGAAGATGGCTCTTATCGGCTGGGGGGGCCTAAATGGGTTTTCTCGGTTCCGCAAAGCGATGCGCCTCTGGTGCt > 2:1646332/1‑76 (MQ=255) aaGATGGCTCTTATCGGCTGGGGGGGCATAAATGGTTTTTCTCGGTTCCGCAAAGCGATGCGCATCTGGTGCTGgc > 2:1006298/1‑76 (MQ=255) tGGCTCTTATCGGCTGGTGGGGCATAAATGGTTTTTCTCGGTTCCGCAAAGCGATGCGCATCTGGTGCTGGCGCAg < 1:768406/76‑1 (MQ=255) tGGCTCTTATCGGCTGGTGGGGCATAAATGGTTTTTCTCGGTTCCGCAAAGCGATGCGCATCTGGTGCTGGCGCAg < 2:2016120/76‑1 (MQ=255) ggCTCTTATCGGCTGGTGGGGCATAAATGGTTTTTCTCGGTTCCGCAAAGCGATGCGCATCTGGTGCTGGCGCAGa < 1:272618/76‑1 (MQ=255) ctctTATCGGCTGGGGGGACATAAATGGTTTTTCTCGGTTCCGCAAAGCGATGCGCATCTGGTGCTGGCGCAGAcc > 2:1589875/1‑76 (MQ=255) tATCGGCTGGGGGGGCATAAATGGTTTTTTTCGGTTCCGCAAAAAGGTGCGCCTCTGGTGCTGGGGCCGAACGCgg > 2:925530/1‑76 (MQ=255) tCGGCTGGTGGGGCATAAATGGTTTTTCTCGGTTCCGCAAAGCGATGCGCATCTGGTGCTGGCGCAGACCGCGggt < 1:967620/76‑1 (MQ=255) tggtggGGCATAAATGGTTTTTCTCGGTTCCGCAAAGCGATGCGCATCTGGTGCTGGCGCAGACCGCGGGTGGTCt < 2:2010981/76‑1 (MQ=255) tggtggGGCATAAATGGTTTTTCTCGGTTCCGCAAAGCGATGCGCATCTGGTGCTGGCGCAGACCGCGGGTGGTCt < 2:1146580/76‑1 (MQ=255) tggtggGGCATAAATGGTTTTTCTCGGTTCCGCAAAGCGATGCGCATCTGGTGCTGGCGCAGACCGCGGGTGGTCt < 2:1160899/76‑1 (MQ=255) tggtggGGCATAAATGGTTTTTCTCGGTTCCGCAAAGCGATGCGCATCTGGTGCTGGCGCAGACCGCGGGTGGTCt < 2:1323228/76‑1 (MQ=255) gtggGGCATAAATGGTTTTTCTCGGTTCCGCAAAGCGATGCGCATCTGGTGCTGGCGCAGACCGCGGGTGGTCTGt < 2:890363/76‑1 (MQ=255) tggGGCATAAATGGTTTTTCTCGGTTCCGCAAAGCGATGCGCATCTGGTGCTGGCGCAGACCGCGGGTGGTCTGTc > 2:1990126/1‑76 (MQ=255) | AAAGCAGGGCGGTTCCGATGTTATGAGCAACACCACCCGTGCAGAGCGTCTGGAAGATGGCTCTTATCGGCTGGTGGGGCATAAATGGTTTTTCTCGGTTCCGCAAAGCGATGCGCATCTGGTGCTGGCGCAGACCGCGGGTGGTCTGTC > REL606/4392987‑4393136 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 20 ≤ ATCG/ATCG < 30 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |