| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | REL606 | 4,393,061 | 0 | T | G | 34.0% | 87.1 / 18.9 | 47 | V211G (GTG→GGG) | aidB | isovaleryl CoA dehydrogenase |
| Reads supporting (aligned to +/- strand): ref base T (1/30); new base G (16/0); total (17/30) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 1.13e-11 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 | |||||||||||
| Rejected as consensus: Frequency below/above cutoff threshold. | |||||||||||
AAAGCAGGGCGGTTCCGATGTTATGAGCAACACCACCCGTGCAGAGCGTCTGGAAGATGGCTCTTATCGGCT‑GGTGGGGCATAAATGGTTTTTCTCGGTTCCGCAAAGCGATGCGCATCTGGTGCTGGCGCAGACCGCGGGTGGTCTGTC > REL606/4392987‑4393136 | aaaGCAGGGCGGTTCCGATGTTATGAGCAACACCACCCGTGCAGAGCGTCTGGAAGATGGCTCTTATCGGCt‑ggtg < 2:600517/76‑1 (MQ=255) aaGCAGGGCGGTTCCGATGTTATGAGCAACACCACCCGTGCAGAGCGTCTGGAAGATGGCTCTTATCGGCt‑ggtgg < 2:3260/76‑1 (MQ=255) aGCAGGGCGGTTCCGATGTTATGAGCAACACCACCCGTGCAGAGCGTCTGGAAGATGGCTCTTATCGGCT‑GGTggg < 1:1286168/76‑1 (MQ=255) gCAGGGCGGTTCCGATGTTATGAGCAACACCACCCGTGCAGAGCGTCTGGAAGATTGCTCTTATCGGCT‑GGTgggg < 2:730385/76‑1 (MQ=255) cAGGGCGGTTCCGATGTTATGAGCAACACCACCCGTGCAGAGCGTCTGGAAGATGGCTCTTATCGGCT‑GGTGGGGc < 2:429278/76‑1 (MQ=255) cAGGGCGGTTCCGATGTTATGAGCAACACCACCCGTGCAGAGCGTCTGGAAGATGGCTCTTATCGGCT‑GGTGGGGc < 2:1033115/76‑1 (MQ=255) cAGGGCGGTTCCGATGTTATGAGCAACACCACCCGTGCAGAGCGTCTGGAAGATGGCTCTTATCGGCT‑GGTGGGGc < 1:887255/76‑1 (MQ=255) ggCGGTTCCGATGTTATGAGCAACACCACCCGTGCAGAGCGTCTGGAAGATGGCTCTTATCGGCT‑GGTGGGGCATa < 1:1240255/76‑1 (MQ=255) ggCGGTTCCGATGTTATGAGCAACACCACCCGTGCAGAGCGTCTGGAAGATGGCTCTTATCGGCT‑GGTGGGGCATa < 1:1671616/76‑1 (MQ=255) cGGTTCCGATGTTATGAGCAACACCACCCGTGCAGAGCGTCTGGAAGATGGCTCTTATCGGCT‑GGGGGGGGGTaaa > 2:805720/1‑76 (MQ=255) ggTTCCGATGTTATGAGCAACACCACCCGTGCAGAGCGTCTGGAAGATGGCTCTTATCGGCT‑GGTGGGGCATAAAt < 2:1818069/76‑1 (MQ=255) ggTTCCGATGTTATGAGCAACACCACCCGTGCAGAGCGTCTGGAAGATGGCTCTTATCGGCT‑GGTGGGGCATAAAt < 1:883431/76‑1 (MQ=255) ttCCGATGTTATGAGCAACACCACCCGTGCAGAGCGTCTGGAAGATGGCTCTTATCGGCT‑GGTGGGGCATAAATgg < 2:670659/76‑1 (MQ=255) ccGATGTTATGAGCAACACCACCCGTGCAGAGCGTCTGGAAGATGGCTCTTATCGGCT‑GGTGGGGCATAAATGGtt < 1:1867122/76‑1 (MQ=255) gATGTTATGAGCAACACCACCCGTGCAGAGCGTCTGGAAGATGGCTCTTATCGGCT‑GGTGGGGCATAAATGGtttt < 2:804020/76‑1 (MQ=255) aTGTTATGAGCAACACCACCCGTGCAGAGCGTCTGGAAGATGGCTCTTATCGGCT‑GGGGAGGCGTAAAAGGctttt > 2:563448/1‑76 (MQ=255) tGTTATGAGCAACACCACCCGTGCAGAGCGTCTGGAAGATGGCTCTTATCGGCT‑GGGGGGGGGTAAATGGgttttt > 1:1250503/1‑75 (MQ=255) ttATGAGCAACACCACCCGTGCAGAGCGTCTGGAAGATGGCTCTTATCGGCT‑GGGGGGGGATAAATGGgttttttt > 2:1014169/1‑73 (MQ=255) aGCAACACCACCCGTGCAGAGCGTCTGGAAGATGGCTCTTATCGGCT‑GGTGGGGCATAAATGGTTTTTCTCGGTTc < 2:621418/76‑1 (MQ=255) aGCAACACCACCCGTGCAGAGCGTCTGGAAGATGGCTCTTATCGGCT‑GGTGGGGCATAAATGGTTTTTCTCGGTTc < 1:984398/76‑1 (MQ=255) aGCAACACCACCCGTGCAGAGCGTCTGGAAGATGGCTCTTATCGGCT‑GGGGGGGCATAAAAGGGTTTTTTCGGGTc > 1:92898/1‑76 (MQ=255) cAACACCACCCGTGCAGAGCGTCTGGAAGATGGCTCTTATCGGCT‑GGTGGGGCATAAATGGTTTTTCTCGGTTCCg < 2:723385/76‑1 (MQ=255) acacCACCCGTGCAGAGCGTCTGGAAGATGGCTCTTATCGGCTGGGGGGGGGGTAAAGGGTTTTTTTTGGTTCCGc > 2:932091/1‑76 (MQ=255) acacCACCCGTGCAGAGCGTCTGGAAGATGGCTCTTATCGGCT‑GGGGGGGGATAAATGGGTTTTTTCGGTTCCGcc > 1:1085375/1‑75 (MQ=255) ccaccCGTGCAGAGCGTCTGGAAGATGGCTCTTATCGGCT‑GGTGGGGCATAAATGGTTTTTCTCGGTTCCGCAAAg < 2:1145239/76‑1 (MQ=255) ccaccCGTGCAGAGCGTCTGGAAGATGGCTCTTATCGGCT‑GGGGGGGCATAAATGGTTTTTCTCGGTTCCGCAAAg > 2:1876061/1‑76 (MQ=255) caccCGTGCAGAGCGTCTGGAAGATGGCTCTTATCGGCT‑GGTGGGGCATAAATGGTTTTTCTCGGTTCCGCAAAGc < 2:616550/76‑1 (MQ=255) ccGTGCAGAGCGTCTGGAAGATGGCTCTTATCGGCT‑GGTGGGGCATAAATGGTTTTTCTCGGTTCCGCAAAGCGAt < 1:1679967/76‑1 (MQ=255) cGTGCAGAGCGTCTGGAAGATGGCTCTTATCGGCT‑GGGGGGGGATAAATGGTTTTTTTCGGTTTCGCAAAACGATg > 2:294329/1‑76 (MQ=255) gCAGAGCGTCTGGAAGATGGCTCTTATCGGCT‑GGTGGGGCATAAATGGTTTTTCTCGGTTCCGCAAAGCGATgcgc < 1:1515857/76‑1 (MQ=255) cAGAGCGTCTGGAAGATGGCTCTTATCGGCT‑GGGGGGGGGTAAATGGGTTTTTTCGGTTCCGCAAAGCGATGCGCa > 1:1632638/1‑76 (MQ=255) agagCGTCTGGAAGATGGCTCTTATCGGCT‑GGTGGGGCATAAATGGTTTTTCTCGGTTCCGCAAAGCGATGCGCAt < 1:461721/76‑1 (MQ=255) agagCGTCTGGAAGATGGCTCTTATCGGCT‑GGGGGGGGATAAATGGTTTTTTTCGGTTCCGCAAAGCGATGCGCAt > 2:1420511/1‑76 (MQ=255) agagCGTCTGGAAGATGGCTCTTATCGGCT‑GGGGGGGCATAAATGGGTTTTTTCGGGTCCGCCAAGCGATGCGCAt > 2:1882955/1‑76 (MQ=255) gCGTCTGGAAGATGGCTCTTATCGGCT‑GGTGGGGCATAAATGGTTTTTCTCGGTTCCGCAAAGCGATGCGCATCTg < 2:604814/76‑1 (MQ=255) gTCTGGAAGATGGCTCTTATCGGCT‑GGTGGGGCATAAATGGTTTTTCTCGGTTCCGCAAAGCGATGCGCATCTGGt < 1:204785/76‑1 (MQ=255) gTCTGGAAGATGGCTCTTATCGGCT‑GGTGGGGCATAAATGGTTTTTCTCGGTTCCGCAAAGCGATGCGCATCTGGt < 1:272139/76‑1 (MQ=255) tCTGGAAGATGGCTCTTATCGGCT‑GGTGGGGCATAAATGGTTTTTCTCGGTTCCGCAAAGCGATGCGCATCTGGTg < 2:352287/76‑1 (MQ=255) cTGGAAGATGGCTCTTATCGGCT‑GGTGGGGCATAAATGGTTTTTCTCGGTTCCGCAAAGCGATGCGCATCTGGTGc < 1:1541173/76‑1 (MQ=255) tGGAAGATGGCTCTTATCGGCT‑GGGGGGGCATAAATGGGTTTTTTCCGGTCCGCAAAGCGATGCGCCTCTGGTGCt > 2:238453/1‑76 (MQ=255) ggAAGATGGCTCTTATCGGCT‑GGGGGGGCGTAAATGGTTTTTCTCGGTTCCGCCAAGCGGTGCGCATCTGGCGCTg > 2:1963651/1‑76 (MQ=255) ggAAGATGGCTCTTATCGGCT‑GGGGGGGCATAAATGGTTTTTTTCGGTTCCGCAAAGCGATGCGCATCTGGTGCTg > 2:1896406/1‑76 (MQ=255) gATGGCTCTTATCGGCT‑GGTGGGGCATAAATGGTTTTTCTCGGTTCCGCAAAGCGATGCGCATCTGGTGCTGgcgc < 2:1787371/76‑1 (MQ=255) gtctTATCTGCT‑GGTGGGGCATAAATTGTTTTTCTCGGTTCCGCAAAGCGATGCGCATCTGGTGCTGGCGCAGAcc < 2:965705/75‑1 (MQ=255) tctTATCGGCT‑GGGGGGGCATAAATGGTTTTTCTCGGTTCCGCAAAGCGATGCGCATCTGGTGCTGGGGCAGACcg > 2:1609404/1‑76 (MQ=255) ttATCGGCT‑GGGGGGGGATAAATGGTTTTTCTCGGTTCCGCAAAGCGATGCGCATCTGGTGCTGGCGCCGACcggg > 2:1252488/1‑74 (MQ=255) aTCGGCT‑GGTGGGGCATAAATGGTTTTTCTCGGTTCCGCAAAGCGATGCGCATCTGGTGCTGGCGCAGACCGCggg < 2:1552493/76‑1 (MQ=255) cGGCT‑GGTGGGGCATAAATGGTTTTTCTCGGTTCCGCAAAGCGATGCGCATCTGGTGCTGGCGCAGACCGCGggtg < 2:223466/76‑1 (MQ=255) cGGCT‑GGTGGGGCATAAATGGTTTTTCTCGGTTCCGCAAAGCGATGCGCATCTGGTGCTGGCGCAGACCGCGggtg < 2:1071541/76‑1 (MQ=255) cGGCT‑GGTGGGGCATAAATGGTTTTTCTCGGTTCCGCAAAGCGATGCGCATCTGGTGCTGGCGCAGACCGCGggtg < 1:969460/76‑1 (MQ=255) ggCT‑GGTGGGGCATAAATGGTTTTTCTCGGTTCCGCAAAGCGATGCGCATCTGGTGCTGGCGCAGACCGCGggtgg > 2:1984539/1‑76 (MQ=255) t‑ggtggGGCATAAATGGTTTTTCTCGGTTCCGCAAAGCGATGCGCATCTGGTGCTGGCGCAGACCGCGGGTGGTCt < 2:524146/76‑1 (MQ=255) gtggGGCATAAATGGTTTTTCTCGGTTCCGCAAAGCGATGCGCATCTGGTGCTGGCGCAGACCGCGGGTGGTCTGt > 1:495880/1‑76 (MQ=255) tggGGCATAAATGGTTTTTCTCGGTTCCGCAAAGCGATGCGCATCTGGTGCTGGCGCAGACCGCGGGTGGTCTGTc > 2:873262/1‑76 (MQ=255) tggGGCATAAATGGTTTTTCTCGGTTCCGCAAAGCGATGCGCATCTGGTGCTGGCGCAGACCGCGGGTGGTCTGTc < 1:667617/76‑1 (MQ=255) | AAAGCAGGGCGGTTCCGATGTTATGAGCAACACCACCCGTGCAGAGCGTCTGGAAGATGGCTCTTATCGGCT‑GGTGGGGCATAAATGGTTTTTCTCGGTTCCGCAAAGCGATGCGCATCTGGTGCTGGCGCAGACCGCGGGTGGTCTGTC > 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: ‑ |