| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | REL606 | 4,393,061 | 0 | T | G | 31.8% | 87.9 / 12.2 | 44 | V211G (GTG→GGG) | aidB | isovaleryl CoA dehydrogenase |
| Reads supporting (aligned to +/- strand): ref base T (2/28); new base G (14/0); total (16/28) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 1.04e-09 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 | |||||||||||
| Rejected as consensus: Frequency below/above cutoff threshold. | |||||||||||
AGCAGGGCGGTTCCGATGTTATGAGCAACACCACCCGTGCAGAGCGTCTGGAAGATGGCTCTTATCGGCTGGTGGGGCATAAATGGTTTTTCTC‑GGTTCCGCAAAGCGATGCGCATCTGGTGCTGGCGCAGACCGCGGGTGGTCTG > REL606/4392989‑4393134 | aGCAGGGCGGTTCCGATTTTATGAGCAACACCCCCCGTGCAGAGCGTCTGGAAGATGGCTCTTATCGGCTGGTggg < 1:1024925/76‑1 (MQ=255)aGCAGGGCGGTTCCGATGTTATGAGCAACACCACCCGTGCAGAGCGTCTGGAAGATGGCTCTTATCGGCTGGTggg < 1:727069/76‑1 (MQ=255)aGCAGGGCGGTTCCGATGTTATGAGCAACACCACCCGTGCAGAGCGTCTGGAAGATGGCTCTTATCGGCTGGTggg < 1:134639/76‑1 (MQ=255) gCAGGGCGGTTCCGATGTTATGAGCAACACCACCCGTGCAGAGCGTCTGGAAGATGGCTCTTATCGGCTGGTgggg < 2:490431/76‑1 (MQ=255) gCAGGGCGGTTCCGATGTTATGAGCAACACCACCCGTGCAGAGCGTCTGGAAGATGGCTCTTATCGGCTGGTgggg < 2:1629112/76‑1 (MQ=255) cAGGGCGGTTCCGATGTTATGAGCAACACCACCCGTGCAGAGCGTCTGGAAGATGGCTCTTATCGGCTGggggggg > 2:1551251/1‑75 (MQ=255) gggCGGTTCCGATGTTATGAGCAACACCACCCGTGCAGAGCGTCTGGAAGATGGCTCTTATCGGCTGggggggggt > 1:1449026/1‑73 (MQ=255) ggCGGTTCCGATGTTATGAGCAACACCACCCGTGCAGAGCGTCTGGAAGATGGCTCTTATCGGCTGGTGGGGCATa < 1:273521/76‑1 (MQ=255) gCGGTTCCGATGTTATGAGCAACACCACCCGTGCAGAGCGTCTGGAAGATGGCTCTTATCGGCTGGGGGGGGATaa > 1:218131/1‑76 (MQ=255) ggTTCCGATGTTATGAGCAACACCACCCGTGCAGAGCGTCTGGAAGATGGCTCTTATCGGCTGGTGGGGCATAAAt < 2:1469278/76‑1 (MQ=255) gTTCCGATGTTATGAGCAACACCACCCGTGCAGAGCGTCTGGAAGATGGCTCTTATCGGCTGGGGGGGCATAAATg > 2:1563828/1‑76 (MQ=255) gTTCCGATGTTATGAGCAACACCACCCGTGCAGAGCGTCTGGAAGATGGCTCTTATCGGCTGGGGGGGCATAAATg > 1:1621794/1‑76 (MQ=255) ccGATGTTATGAGCAACACCACCCGTGCAGAGCGTCTGGAAGATGGCTCTTATCGGCTGGGGGGGGGTAAATgggt > 1:1810038/1‑74 (MQ=255) gATGTTATGAGCAACACCACCCGTGCAGAGCGTCTGGAAGATGGCTCTTATCGGCTGGTGGGGCATAAATGGtttt < 2:987824/76‑1 (MQ=255) aTGTTATGAGCAACACCACCCGTGCAGAGCGTCTGGAAGATGGCTCTTATCGGCTGGTGGGGCATAAATGGttttt < 2:518771/76‑1 (MQ=255) aTGTTATGAGCAACACCACCCGTGCAGAGCGTCTGGAAGATGGCTCTTATCGGCTGGGGGGGGGTAAATGGgtttt > 1:1563173/1‑76 (MQ=255) tGTTATGAGCAACACCACCCGTGCAGAGCGTCTGGAAGATGGCTCTTATCGGCTGGTGGGGCATAAATGGTTTTtc < 2:905627/76‑1 (MQ=255) tGTTATGAGCAACACCACCCGTGCAGAGCGTCTGGAAGATGGCTCTTATCGGCTGGGGGGGGGTAAATGGtttttt > 1:715670/1‑75 (MQ=255) gTTATGAGCAACACCACCCGTGCAGAGCGTCTGGAAGATGGCTCTTATCGGCTGGTGGGGCATAAATGGTTTTtct < 1:996527/76‑1 (MQ=255) tATGAGCAACACCACCCGTGCAGAGCGTCTGGAAGATGGCTCTTATCGGCTGGGGGGGCCTAAATGGGTTTTCTC‑g > 2:1115682/1‑76 (MQ=255) aGCAACACCACCCGTGCAGAGCGTCTGGAAGATGGCTCTTATCGGCTGGTGGGGCATAAATGGTTTTTCTC‑GGTTc < 2:1708824/76‑1 (MQ=255) gCAACACCACCCGTGCAGAGCGTCTGGAAGATGGCTCTTATCGGCTGGGGGGGGATAAATGGTTTTTTTC‑GGTTcc > 1:1379286/1‑76 (MQ=255) cAACACCACCCGTGCAGAGCGTCTGGAAGATGGCTCTTATCGGCTGGTGGGGCATAAATGGTTTTTCTC‑GGTTCCg < 2:1430710/76‑1 (MQ=255) cAACACCACCCGTGCAGAGCGTCTGGAAGATGGCTCTTATCGGCTGGTGGGGCATAAATGGTTTTTCTC‑GGTTCCg < 2:359230/76‑1 (MQ=255) aaCACCACCCGTGCAGAGCGTCTGGAAGATGGCTCTTATCGGCTGGGGGGGGATAAATGGTTTTTTTT‑GGTTCCGc > 2:156940/1‑76 (MQ=255) acacCACCCGTGCAGAGCGTCTGGAAGATGGCTCTTATCGGCTGGGGGGGGATAAATGGGTTTTTTC‑GGTTCCCCa > 2:1778796/1‑76 (MQ=255) accaccCGTGCAGAGCGTCTGGAAGATGGCTCTTATCGGCTGGTGGGGCATAAATGGTTTTTCTC‑GGTTCCGCaaa < 1:1459951/76‑1 (MQ=255) accCGTGCAGAGCGTCTGGAAGATGGCTCTTATCGGCTGGGGGGGGATAAATGGGTTTTTTCGGGTTCCCCaaacc > 1:395551/1‑74 (MQ=255) ccGTGCAGAGCGTCTGGAAGATGGCTCTTATCGGCTGGTGGGGCATAAATGGTTTTTCTC‑GGTTCCGCAAAGCGAt < 2:1139141/76‑1 (MQ=255) ccGTGCAGAGCGTCTGGAAGATGGCTCTTATCGGCTGGTGGGGCATAAATGGTTTTTCTC‑GGTTCCGCAAAGCGAt < 2:334600/76‑1 (MQ=255) gTGCAGAGCGTCTGGAAGATGGCTCTTATCGGCTGGGGGGGCATAAATGGTTTTTCTC‑GGTTCCGCAAAGCGCTgc > 1:1729243/1‑76 (MQ=255) agagCGTCTGAAAGATGGCTCTTATCGGCTGGGGGGGCATAAATGGGTTTTCTC‑GGTTCCGCCAAGCGATGCGCAt > 2:1617795/1‑76 (MQ=255) gagCGTCTGGAAGATGGCTCTTATCGGCTGGTGGGGCATAAATGGTTTTTCTC‑GGTTCCGCAAAGCGATGCGCATc < 1:824454/76‑1 (MQ=255) agCGTCTGGAAGATGGCTCTTATCGGCTGGTGGGGCATAAATGGTTTTTCTC‑GGTTCCGCAAAGCGATGCGCATCt < 1:1091163/76‑1 (MQ=255) agCGTCTGGAAGATGGCTCTTATCGGCTGGGGGGGGATAAATGGGTTTTTTC‑GGTTCCGCCAAGCGATGCGCATCt > 2:175009/1‑76 (MQ=255) gCGTCTGGAAGATGGCTCTTATCGGCTGGTGGGGCATAAATGGTTTTTCTC‑GGTTCCGCAAAGCGATGCGCATCTg < 1:1250138/76‑1 (MQ=255) gTCTGGAAGATGGCTCTTATCGGCTGGTGGGGCATAAATGGTTTTTCTC‑GGTTCCGCAAAGCGATGCGCATCTGGt < 1:430911/76‑1 (MQ=255) gTCTGGAAGATGGCTCTTATCGGCTGGGGGGGGGTTAATGGGTTTTTTC‑GGTGCCGGAACGCGATGCTAATCCGGt > 2:1904945/1‑76 (MQ=255) tCTGGAAGAGCGCTCTTATCGGCTGGTTGGGCATAAATGGTTTTTCTC‑GGTTCCGCAAAGCGATGCGCATCTGGTg < 2:1546733/76‑1 (MQ=255) ggAAGATGGCTCTTATCGGCTGGTGGGGCATAAATGGTTTTTCTC‑GGTTCCGCAAAGCGATGCGCATCTGGTGCTg < 2:1556006/76‑1 (MQ=255) ggAAGATGGCTCTTATCGGCTGGGGGGGCATAAATGGGTTTTCTC‑GGTTCCGCAAAGCGATGCGCATCTGgtggtg > 2:689033/1‑76 (MQ=255) gAAGATGGCTCTTATCGGCTGGTGGGGCATAAATGGTTTTTCTC‑GGTTACGCAAAGCGATGCGCATCTGGTGCTgg < 1:1115682/76‑1 (MQ=255) tGGCTCTTATCGGCTGGTGGGGCATAAATGGTTTTTCTC‑GGTTCCGCAAAGCGATGCGCATCTGGTGCTGGCGCAg < 1:808843/76‑1 (MQ=255) gCTCTTATCGGCTGGTGGGGCATAAATGGTTTTTCTC‑GGTTCCGCAAAGCGATGCGCATCTGGTGCTGGCGCAGAc < 1:914532/76‑1 (MQ=255) tATCGGCTGGTGGGGCATAAATGGTTTTTCTC‑GGTTCCGCAAAGCGATGCGCATCTGGTGCTGGCGCAGACCGCgg < 1:1030623/76‑1 (MQ=255) ggCTGGTGGGGCATAAATGGTTTTTCTC‑GGTTCCGCAAAGCGATGCGCATCTGGTGCTGGCGCAGACCGCGggtgg < 2:646262/76‑1 (MQ=255) cTGGTGGGGCATAAATGGTTTTTCTC‑GGTTCCGCAAAGCGATGCGCATCTGGTGCTGGCGCAGACCGCGGGTGGTc > 1:1831098/1‑76 (MQ=255) cTGGTGGGGCATAAATGGTTTTTCTC‑GGTTCCGCAAAGCGATGCGCATCTGGTGCTGGAGCAGACCGCGGGTGGTc > 1:1065905/1‑76 (MQ=255) tggtggGGCATAAATGGTTTTTCTC‑GGTTCCGCAAAGCGATGCGCATCTGGTGCTGGCGCAGACCGCGGGTGGTCt > 2:1009757/1‑76 (MQ=255) tggtggGGCATAAATGGTTTTTCTC‑GGTTCCGCAAAGCGATGCGCATCTGGTGCTGGCGCAGACCGCGGGTGGTCt > 1:635869/1‑76 (MQ=255) ggtggGGCATAAATGGTTTTTCTC‑GGTTCCGCAAAGCGATGCGCATCTGGTGCTGGCGCAGACCGCGGGTGGTCTg < 2:1469659/76‑1 (MQ=255) ggtggGGCATAAATGGTTTTTCTC‑GGTTCCGCAAAGCGATGCGCATCTGGTGCTGGCGCAGACCGCGGGTGGTCTg < 2:1449026/76‑1 (MQ=255) ggtggGGCATAAATGGTTTTTCTC‑GGTTCCGCAAAGCGATGCGCATCTGGTGCTGGCGCAGACCGCGGGTGGTCTg > 1:528641/1‑76 (MQ=255) | AGCAGGGCGGTTCCGATGTTATGAGCAACACCACCCGTGCAGAGCGTCTGGAAGATGGCTCTTATCGGCTGGTGGGGCATAAATGGTTTTTCTC‑GGTTCCGCAAAGCGATGCGCATCTGGTGCTGGCGCAGACCGCGGGTGGTCTG > REL606/4392989‑4393134 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 21 ≤ ATCG/ATCG < 30 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |