Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | REL606 | 4,393,061 | 0 | T | G | 35.5% | 47.2 / 15.1 | 31 | V211G (GTG→GGG) | aidB | isovaleryl CoA dehydrogenase |
Reads supporting (aligned to +/- strand): ref base T (5/15); new base G (11/0); total (16/15) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 6.77e-05 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
AAGCAGGGCGGTTCCGATGTTATGAGCAACACCACCCGTGCAGAGCGTCTGGAAGATGGCTCTTATCGGCT‑GGTGGGGCATAAATGGTTTTTCTCGGTTCCGCAAAGCGATGCGCATCTGGTGCTGGCGCAGACCGCGGGTGGTCTGTC > REL606/4392988‑4393136 | aaGCAGGGCGGTTCCGATGTTATGAGCAACACCACCCGTGCAGAGCGTCTGGAAGATGGCTCTTATCGGCt‑ggtgg < 2:1193278/76‑1 (MQ=255) aaGCAGGGCGGTTCCGATGTTATGAGCAACACCACCCGTGCAGAGCGTCTGGAAGATGGCTCTTATCGGCt‑ggtgg < 2:398789/76‑1 (MQ=255) aaGCAGGGCGGTTCCGATGTTATGAGCAACACCACCCGTGCAGAGCGTCTGGAAGATGGCTCTTATCGGCt‑ggggg > 2:800921/1‑76 (MQ=255) cAGGGCGGTTCCGATGTTATGAGCAACACCACCCGTGCAGAGCGTCTGGAAGATGGCTCTTATCGGCT‑GGTGGGGc < 1:1037752/76‑1 (MQ=255) aGGGCGGTTCCGATGTTATGAGCAACACCACCCGTGCAGAGCGTCTGGAAGATGGCTCTTATCGGCT‑GGTggggga > 2:58126/1‑74 (MQ=255) gCGGTTCCGATGTTATGAGCAACACCACCCGTGCAGAGCGTCTGGAAGATGGCTCTTATCGGCT‑GGTGGGGCATaa < 2:1326482/76‑1 (MQ=255) cGGTTCCGATGTTATGAGCAACACCACCCGTGCAGAGCGTCTGGAAGATGGCTCTTATTGGCT‑GGTGGGGCATaaa < 1:485274/76‑1 (MQ=255) ggTTCCGATGTTATGAGCAACACCACCCGTGCAGAGCGTCTGGAAGATGGCTCTTATCGGCT‑GGGGGGGGGTAAAt > 1:1153666/1‑76 (MQ=255) gTTCCGATGTTATGAGCAACACCACCCGTGCAGAGCGTCTGGAAGATGGCTCTTATCGGCT‑GGTGGGGCATAAATg < 1:1341086/76‑1 (MQ=255) gTTCCGATGTTATGAGCAACACCACCCGTGCAGAGCGTCTGGAAGATGGCTCTTATCGGCT‑GGGGGGGGGTaaaag > 1:1088720/1‑74 (MQ=255) gTTCCGATGTTATGAGCAACACCACCCGTGCAGAGCGTCTGGAAGATGGCTCTTATCGGCT‑GGGGGGGGATAAATg > 1:880396/1‑76 (MQ=255) ccGATGTTATGAGCAACACCACCCGTGCAGAGCGTCTGGAAGATGGCTCTTATCGGCT‑GGTGGGGGATAAATgggt > 1:367290/1‑74 (MQ=255) ccGATGTTATGAGCAACACCACCCGTGCAGAGCGTCTGGAAGATGGCTCTTATCGGCT‑GGTGGGGCATAAATGGtt < 1:159958/76‑1 (MQ=255) gATGTTATGAGCAACACCACCCGTGCAGAGCGTCTGGAAGATGGCTCTTATCGGCT‑GGGGGGGGGTAAAAGGgttt > 2:1400225/1‑76 (MQ=255) aTGTTATGAGCAACACCACCCGTGCAGAGCGTCTGGAAGATGGCTCTTATCGGCT‑GGTGGGGCATAAATGGttttt < 2:622601/76‑1 (MQ=255) tGTTATGAGCAACACCACCCGTGCAGAGCGTCTGGACGATGGCTCTTATCGGCTGGGGGGGGGATAATAGGttttt > 1:507321/1‑76 (MQ=255) gTTATGAGCAACACCACCCGTGCAGAGCGTCTGGAAGATGGCTCTTATCGGCT‑GGTGGGGCATAAATGGTTTTtct < 2:909107/76‑1 (MQ=255) gTTATGAGCAACACCACCCGTGCAGAGCGTCTGGAAGATGGCTCTTATCGGCT‑GGTGGGGCATAAATGGTTTTtct < 2:1109376/76‑1 (MQ=255) cAACACCACCCGTGCAGAGCGTCTGGAAGATGGCTCTTATCGGCT‑GGTGGGGCATAAATGGTTTTTCTCGGTTCCg < 2:1033160/76‑1 (MQ=255) cAACACCACCCGTGCAGAGCGTCTGGAAGATGGCTCTTATCGGCT‑GGGGGGGGATAAATGGTTTTTTTTGGTTCCg > 2:487411/1‑76 (MQ=255) accaccCGTGCAGAGCGTCTGGAAGATGGCTCTTATCGGCT‑GGTGGGGCATAAATGGTTTTTCTCGGTTCCGCaaa < 1:120067/76‑1 (MQ=255) caccCGTGCAGAGCGTCTGGAAGATGGCTCTTATCGGCT‑GGTGGGGCATAAATGGTTTTTCTCGGTTCCGCAAAGc < 1:991683/76‑1 (MQ=255) caccCGTGCAGAGCGTCTGGAAGATGGCTCTTATCGGCT‑GGTGGGGCATAAATGGTTTTTCTCGGTTCCGCAAAGc > 2:129818/1‑76 (MQ=255) accCGTGCAGAGCGTCTGGAAGATGGCTCTTATCGGCT‑GGTGGGGCATAAATGGTTTTTCTCGGTTGCGCAAAGCg < 2:1384747/76‑1 (MQ=255) accCGTGCAGAGCGTCTGGAAGATGGCTCTTATCGGCT‑GGGGGGGGGTAAAAGGGTTTTCTCGGCTCCCCAAAGCg > 2:88350/1‑76 (MQ=255) ccGTGCAGAGCGTCTGGAAGATGGCTCTTATCGGCT‑GGGGGGGGATAAATGGTTTTTTTCGGTTCCGCCAAGCGAt > 1:829844/1‑76 (MQ=255) tGCAGAGCGTCTGGAAGATGGCTCTTATCGGCT‑GGTGGGGCATAAATGGTTTTTCTCGGTTCCGCAAAGCGATgcg < 1:122472/76‑1 (MQ=255) tGCAGAGCGTCTGGAAGATGGCTCTTATCGGCT‑GGGGGGGCATAAATGGTTTTTTTCGGGTCCGCAAAGCGCTgcg > 2:341048/1‑76 (MQ=255) agagCGTCTGGAAGATGGCTCTTATCGGCT‑GGGGGGGGCTAAAAGGGTTTTTTCGGTTTCGCCAAACGATgcgcct > 1:363814/1‑74 (MQ=255) cGTCTGGAAGATGGCTCTTATCGGCT‑GGGGGGGCATAAATGGTTTTTTTCGGTTCCGCCAAGCGATGCGCATCTgg > 2:1500070/1‑76 (MQ=255) tGGAAGATGGCTCTTATCGGCT‑GGTGGGGCATAAATGGTTTTTCTCGGTTCCGCAAAGCGATGCGCATCTGGTGCt < 1:892469/76‑1 (MQ=255) aGATGGCTCTTATCGGCT‑GGTGGGGCATAAATGGTTTTTCTCGGTTCCGCAAAGCGATGCGCATCTGGTGCTGgcg < 1:58229/76‑1 (MQ=255) ctTATCGGCT‑GGTGGGGGGTAAATGGTTTTTCTCGGTTCCTCAAAGCGATGCGCATTTGGTGCAGGGCCAGACcgc > 2:150741/1‑76 (MQ=255) ctTATCGGCT‑GGGGGGGCATAAATGGGTTTTCTCGGTTCCGCAAAGCGATGCGCATCTGGTGCTGGCGCAGACcgc > 2:985671/1‑76 (MQ=255) cT‑GGTGGGGCATAAATGGTTTTTCTCGGTTCCGCAAAGCGATGCGCATCTGGTGCTGGCGCAGACCGCGGGTGGTc > 1:1442142/1‑76 (MQ=255) t‑ggtggGGCATAAATGGTTTTTCTCGGTTCCGCAAAGCGATGCGCATCTGGTGCTGGCGCAGACCGCGGGTGGTCt < 1:732176/76‑1 (MQ=255) ggtggGGCATAAATGGTTTTTCTCGGTTCCGCAAAGCGATGCGCATCTGGTGCTGGCGCAGACCGCGGGTGGTCTg > 2:354915/1‑76 (MQ=255) gtggGGCATAAATGGTTTTTCTCGGTTCCGCAAAGCGATGCGCATCTGATGCTGGCGCAGACCGCGGGTGGTCTGt > 1:541932/1‑76 (MQ=255) tggGGCATAAATGGTTTTTCTCGGTTCCGCAAAGCGATGCGCATCTGGTGCTGGCGCAGACCGCGGGTGGTCTGTc > 2:1129973/1‑76 (MQ=255) | AAGCAGGGCGGTTCCGATGTTATGAGCAACACCACCCGTGCAGAGCGTCTGGAAGATGGCTCTTATCGGCT‑GGTGGGGCATAAATGGTTTTTCTCGGTTCCGCAAAGCGATGCGCATCTGGTGCTGGCGCAGACCGCGGGTGGTCTGTC > REL606/4392988‑4393136 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 25 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |