Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | REL606 | 4,393,061 | 0 | T | G | 46.2% | 36.3 / 22.7 | 39 | V211G (GTG→GGG) | aidB | isovaleryl CoA dehydrogenase |
Reads supporting (aligned to +/- strand): ref base T (3/18); new base G (18/0); total (21/18) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 2.74e-08 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
CAGGGCGGTTCCGATGTTATGAGCAACACCACCCGTGCAGAGCGTCTGGAAGATGGCTCTTATCGGCTGGTGGGGCATAAATGGTTTTTCTCGGTTCCGCAAAGCGATGCGCATCTGGTGCTGGCGCAGACCGCGGGTGGTCTGTC > REL606/4392991‑4393136 | cAGGGCGGTTCCGATGTTATGAGCAACACCACCCGTGCAGAGCGTCTGGAAGATGGCTCTTATCGGCTGggggggg > 1:1270705/1‑75 (MQ=255) aGGGCGGTTCCGATGTTATGAGCAACACCACCCGTGCAGAGCGTCTGGAAGATGGCTCTTATCGGCTGGTGGGGCa < 1:246059/76‑1 (MQ=255) gggCGGTTCCGATGTTATGAGCAACACCACCCGTGCAGAGCGTCTGGAAGATGGCTCTTATCGGCTGggggggggg > 1:1403848/1‑73 (MQ=255) gCGGTTCCGATGTTATGAGCAACACCACCCGTGCAGAGCGTCTGGAAGATGGCTCTTATCGGCTGGTGGGGCATaa < 2:1429371/76‑1 (MQ=255) ggTTCCGATGTTATGAGCAACACCACCCGTGCAGAGCGTCTGGAAGATGGCTCTTATCGGCTGGTGGGGCATAAAt < 2:676709/76‑1 (MQ=255) gTTCCGATGTTATGAGCAACACCACCCGTGCAGAGCGTCTGGAAGATGGCTCTTATCGGCTGGGGGGGGGTAAATg > 1:1843843/1‑76 (MQ=255) cGATGTTATGAGCAACACCACCCGTGCAGAGCGTCTGGAAGATGGCTCTTATCGGCTGGTGGGGCATAAATGGttt < 2:1781563/76‑1 (MQ=255) cGATGTTATGAGCAACACCACCCGTGCAGAGCGTCTGGAAGATGGCTCTTATCGGCTGGGGGGGGGTAAAAGGttt > 1:255640/1‑76 (MQ=255) aTGTTATGAGCAACACCACCCGTGCAGAGCGTCTGGAAGATGGCTCTTATCGGCTGGTGGGGCATAAATGGttttt < 2:1158691/76‑1 (MQ=255) aTGTTATGAGCAACACCACCCGTGCAGAGCGTCTGGAAGATGGCTCTTATCGGCTGGGGGGGGATAAATGGttttt > 2:616059/1‑76 (MQ=255) gTTATGAGCAACACCACCCGTGCAGAGCGTCTGGAAGATGGCTCTTATCGGCTGGTGGGGCATAAATGGTTTTtct < 2:1412401/76‑1 (MQ=255) aTGAGCAACACCACCCGTGCAGAGCGTCTGGAAGATGGCTCTTATCGGCTGGGGGGGGGTAAATGGGTTTTTTCgg > 2:491875/1‑76 (MQ=255) tGAGCAACACCACCCGTGCAGAGCGTCTGGAAGATGGCTCTTATCGGCTGGGGGGGGATAAATGGGTTTTTTCGGt > 2:1108036/1‑76 (MQ=255) aGCAACACCACCCGTGCAGAGCGTCTGGAAGATGGCTCTTATCGGCTGGTGGGGCATAAATGGTTTTTCTCGGTTc < 2:1609980/76‑1 (MQ=255) cAACACCACCCGTGCAGAGCGTCTGGAAGATGGCTCTTATCGGCTGGTGGGGCATAAATGGTTTTTCTCGGTTCCg < 1:1452790/76‑1 (MQ=255) acacCACCCGTGCAGAGCGTCTGGAAGATGGCTCTTATCGGCTGGGGGGGGGTAAATGGGTTTTTTCGGGTTCGcc > 1:733124/1‑75 (MQ=255) caccacCCGTGCAGAGCGTCTGGAAGATGGCTCTTATCGGCTGGGGGGGGATAAATGGGTTTTTTCGGTTCCGcca > 2:898878/1‑74 (MQ=255) accaccCGTGCAGAGCGTCTGGAAGATGGCTCTTATCGGCTGGTGGGGCATAAATGGTTTTTCTCGGTTCCGCaaa < 2:1403848/76‑1 (MQ=255) tGCAGAGCGTCTGGAAGATGGCTCTTATCGGCTGGGGGGGGATAAATGGGTTTTTTCGGTTCCGCAAAGCGATgcg > 2:1214577/1‑76 (MQ=255) agagCGTCTGGAAGATGGCTCTTATCGGCTGGGGGGGGATAAATGGGTTTTTTCGGTTCCGCACAGCGATGCGCAt > 2:658971/1‑76 (MQ=255) agCGTCTGGAAGATGGCTCTTATCGGCTGGGGGGGCATAAATGGGTTTTTTTTGTTCCGCCAAGCGATGCGCATCt > 2:903515/1‑76 (MQ=255) cTGGAAGATGGCTCTTATCGGCTGGTGGGGCATAAATGGTTTTTCTCGGTTCCGCAAAGCGATGCGCATCTGGTGc < 1:1315174/76‑1 (MQ=255) cTGGAAGATGGCTCTTATCGGCTGGTGGGGCATAAATGGTTTTTCTCGGTTCCGCAAAGCGATGCGCATCTGGTGc < 1:1221440/76‑1 (MQ=255) tGGAAGATGGCTCTTATCGGCTGGGGGGGGATAAATGGTTTTTCTCGGTTCCGCAAAGCGATGCGCATCTGGTGCt > 1:1475379/1‑76 (MQ=255) tGGAAGATGGCTCTTATCGGCTGGGGGGGCATAAATGGGTTTTCTCGGTTCCGCAAAGCGATGCGCATCTGGTGCt > 2:1197232/1‑76 (MQ=255) ggAAGATGGCTCTTATCGGCTGGTGGGGCATAAATGGTTTTTCTCGGTTCCGCAAAGCGATGCGCATCTGGTGCTg < 1:1158223/76‑1 (MQ=255) aaGATGGCTCTTATCGGCTGGTGGGGGATAAATGGTTTTTCTCGGTTCCGCAAAGCGATGCGCATCTGGTGCTGgc > 2:1221913/1‑76 (MQ=255) aaGATGGCTCTTATCGGCTGGTGGGGCATAAATGGTTTTTCTCGGTTCCGCAAAGCGATGCGCATCTGGTGCTggn < 1:387133/76‑2 (MQ=255) aaGATGGCTCTTATCGGCTGGTGGGGCATAAAAGGTTTTTTTTGGTTCCGCAAAACGATGCGCATCTGGTGCTGgc > 2:150287/1‑76 (MQ=255) aaGATGGCTCTTATCGGCTGGGGGGGCATAAATGGTTTTTCTCGGTTCCGCAAAGCGATGCGCATCTGGGGCTggg > 2:258859/1‑75 (MQ=255) gATGGCTCTTATCGGCTGGGGGGGCATAAATGGTTTTTCTCGGGTCCGCAAAGCGATGCGCATCTGGTGCTGgcgc > 2:426456/1‑76 (MQ=255) tGGCTCTTATCGGCTGGTGGGGCATAAATGGTTTTTCTCGGTTCCGCAAAGCGATGCGCATCTGGTGCTGGCGCAg < 1:67450/76‑1 (MQ=255) gCTCTTATCGGCTGGTGGGGCATAAATGGTTTTTCTCGGTTCCGCAAAGCGATGCGCATCTGGTGCTGGCGCAGAc < 1:388108/76‑1 (MQ=255) gCTCTTATCGGCTGGGGGGGGATAAATGGTTTTTTTCGGTTCCGCAAAGCCATGCGCATCTGGGTCTGGCGCAGAc > 2:1344283/1‑76 (MQ=255) ttATCGGCTGGGGGGGCATAAATGGTTTTTCTCGGTTCCGCAAAGCGATGCGCATCTGGTGCTGGCGCAGACcgcg > 2:137471/1‑76 (MQ=255) ttATCGGCTGGGGGGGCATAAATGGGTTTTCTCGGTTCCGCAAAGCGATGCGCATCTGGTGCTGGCGCGGAGcgcg > 2:831707/1‑76 (MQ=255) tATCGGCTGGGGGGGGATAAAAGGTTTTTCTCGGTTCCGCAAAACGATGCGCATCTGGGGCTGGCGCAGACCGCgg > 2:1419651/1‑76 (MQ=255) tATCGGCTGGGGGGGCATAAATGGTTTTTCTCGGTTCCGCAAAGCGATGCGCATCTGGTGCTGGCGCAGACCGCgg > 2:1460437/1‑76 (MQ=255) tATCGGCTGGGGGGGCATAAATGGGTTTTCTCGGTTCCGCCAAGCGATGCGCATCTGGTGCTGGCGCCGACCGCgg > 2:1839386/1‑76 (MQ=255) tCGGCTGGTGGGGCATAAATGGTTTTTCTCGGTTCCGCAAAGCGATGCGCATCTGGTGCTGGCGCAGACCGCGggt < 2:1315055/76‑1 (MQ=255) tCGGCTGGTGGGGCATAAATGGTTTTTCTCGGTTCCGCAAAGCGATGCGCATCTGGTGCTGGCGCAGACCGCGggt < 1:1174622/76‑1 (MQ=255) cGGCTGGTGGGGCATAAATGGTTTTTCTCGGTTCCGCAAAGCGATGCGCATCTGGTGCTGGCGCAGACCGCGggtg < 2:551954/76‑1 (MQ=255) ggCTGGTGGGGCATAAATGGTTTTTCTCGGTTCCGCAAAGCGATGCGCATCTGGTGCTGGCGCAGACCGCGggtgg > 1:1119292/1‑76 (MQ=255) tggGGCATAAATGGTTTTTCTCGGTTCCGCAAAGCGATGCGCATCTGGTGCTGGCGCAGACCGCGGGTGGTCTGTc < 2:524692/76‑1 (MQ=255) | CAGGGCGGTTCCGATGTTATGAGCAACACCACCCGTGCAGAGCGTCTGGAAGATGGCTCTTATCGGCTGGTGGGGCATAAATGGTTTTTCTCGGTTCCGCAAAGCGATGCGCATCTGGTGCTGGCGCAGACCGCGGGTGGTCTGTC > REL606/4392991‑4393136 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 17 ≤ ATCG/ATCG < 29 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |