Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | REL606 | 4,393,061 | 0 | T | G | 28.9% | 89.6 / 14.2 | 45 | V211G (GTG→GGG) | aidB | isovaleryl CoA dehydrogenase |
Reads supporting (aligned to +/- strand): ref base T (5/27); new base G (13/0); total (18/27) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 1.17e-07 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
AAAGCAGGGCGGTTCCGATGTTATGAGCAACACCACCCGTGCAGAGCGTCTGGAAGATGGCTCTTATCGGCTGGTGGGGCATAAATGGTTTTTCTCGGTTCCGCAAAGCGATGCGCATCTGGTGCTGGCGCAGACCGCGGGTGGTCTG > REL606/4392987‑4393134 | aaaGCAGGGCGGTTCCGATGTTATGAGCAACACCACCCGTGCAGAGCGTCTGGAAGATGGCTCTTATCGGCtggtg < 2:1690145/76‑1 (MQ=255) gCAGGGCGGTTCCGATGTTATGAGCAACACCACCCGTGCAGAGCGTCTGGAAGATGGCTCTTATCGGCTGGTgggg < 2:1625322/76‑1 (MQ=255) cAGGGCGGTTCCGATGTTATGAGCAACACCACCCGTGCAGAGCGTCTGGAAGATGGCTCTTATCGGCTGGTGGGGc < 1:1246996/76‑1 (MQ=255) aGGGCGGTTCCGATGTTATGAGCAACACCACCCGTGCAGAGCGTCTGGAAGATGGCTCTTATCGGCTGGTGGGGCa < 2:76129/76‑1 (MQ=255) gggCGGTTCCGATGTTATGAGCAACACCACCCGTGCAGAGCGTCTGGAAGATGGCTCTTATCGGCTGGTGGGGCan < 1:1559453/76‑2 (MQ=255) gTTCCGATGTTATGAGCAACACCACCCGTGCAGAGCGTCTGGAAGATGGCTCTTATCGGCTGGGGGGGGATAAATg > 2:1101164/1‑76 (MQ=255) ttCCGATGTTATGAGCAACACCACCCGTGCAGAGCGTCTGGAAGATGGCTCTTATCGGCTGGTGGGGCATAAATgg < 1:1574671/76‑1 (MQ=255) ttCCGATGTTATGAGCAACACCACCCGTGCAGAGCGTCTGGAAGATGGCTCTTATCGGCTGGGGGGGGGTAAATgg > 2:46208/1‑76 (MQ=255) tCCGATGTTATGAGCAACACCACCCGTGCAGAGCGTCTGGAAGATGGCTCTTATCGGCTGGTGGGGCATAAATGGt < 1:1753842/76‑1 (MQ=255) tCCGATGTTATGAGCAACACCACCCGTGCAGAGCGTCTGGAAGATGGCTCTTATCGGCTGGGGGGGGATAAATggg > 1:118028/1‑75 (MQ=255) cGATGTTATGAGCAACACCACCCGTGCAGAGCGTCTGGAAGATGGCTCTTATCGGCTGGTGGGGCATAAATGGttt < 2:1156780/76‑1 (MQ=255) aTGTTATGAGCAACACCACCCGTGCAGAGCGTCTGGAAGATGGCTCTTATCGGCTGGTGGGGCATAAATGGttttt < 2:1055399/76‑1 (MQ=255) tGTTATGAGCAACACCACCCGTGCAGAGCGTCTGGAAGATGGCTCTTATCGGCTGGGGGGGGGTAAATGGgttttt > 1:1010477/1‑75 (MQ=255) gTTATGAGCAACACCACCCGTGCAGAGCGTCTGGAAGATGGCTCTTATCGGCTGGTGGGGCATAAATGGTTTTtct < 1:311051/76‑1 (MQ=255) tATGAGCAACACCACCCGTGCAGAGCGTCTGGAAGATGGCTCTTATCGGCTGGTGGGGCATAAATGGTTTTTCTCg < 1:1113686/76‑1 (MQ=255) gAGCAACACCACCCGTGCAGAGCGTCTGGAAGATGGCTCTTATCGGCTGGTGGGGCATAAATGGTTTTTCTCGGtt < 2:437213/76‑1 (MQ=255) aGCAACACCACCCGTGCAGAGCGTCTGGAAGATGGCTCTTATCGGCTGGGGGGGGATAAATGGTTTTTTTCGGTTc > 2:1114671/1‑76 (MQ=255) cAACACCACCCGTGCAGAGCGTCTGGAAGATGGCTCTTATCGGCTGGTGGGGCATAAATGGTTTTTCTCGGTTCCg < 2:1391515/76‑1 (MQ=255) acacCACCCGTGCAGAGCGTCTGGAAGATGGCTCTTATCGGCTGGTGGGGCATAAATGGTTTTTCTCGGTTCCGCa < 2:903544/76‑1 (MQ=255) caccacCCGTGCAGAGCGTCTGGAAGATGGCTCTTATCGGCTGGTGGGGCATAAATGGTTTTTCTCGGTTCCGCaa < 1:836296/76‑1 (MQ=255) ccaccCGTGCAGAGCGTCTGGAAGATGGCTCTTATCGGCTGGTGGGGCATAAATGGTTTTTCTCGGTTCCGCAAAg < 2:895964/76‑1 (MQ=255) accCGTGCAGAGCGTCTGGAAGATGGCTCTTATCGGCTGGTGGGGCATAAATGGTTTTTCTCGGTTCCGCAAAGCg < 1:1421262/76‑1 (MQ=255) cGTGCAGAGCGTCTGGAAGATGGCTCTTATCGGCTGGTGGGGCATAAATGGTTTTTCTCGGTTCCGCAAAGCGATg < 2:1358150/76‑1 (MQ=255) cGTGCAGAGCGTCTGGAAGATGGCTCTTATCGGCTGGGGGGGCATAAATGGTTTTTCTCGGTTCCGCCAAGCGATg > 2:1615863/1‑76 (MQ=255) gTGCAGAGCGTCTGGAAGATGGCTCTTATCGGCTGGGGGGGGATAAATGGTTTTTCTCGGTTCCGCAAAGCGATgc > 1:1607150/1‑76 (MQ=255) agCGTCTGGAAGATGGCTCTTATCGGCTGGTGGGGCATAAATGGTTTTTCTCGGTTCCGCAAAGCGATGCGCATCt < 2:822152/76‑1 (MQ=255) agCGTCTGGAAGATGGCTCTTATCGGCTGGTGGGGCATAAATGGTTTTTCTCGGTTCCGCAAAGCGATGCGCATCt < 2:280905/76‑1 (MQ=255) gCGTCTGGAAGATGGCTCTTATCGGCTGGTGGGGCATAAATGGTTTTTCTCGGTTCCGCAAAGCGATGCGCATCTg < 2:830574/76‑1 (MQ=255) cGTCTGGAAGATGGCTCTTATCGGCTGGTGGGGCATAAATGGTTTTTCTCGGTTCCGCAAAGCGATGCGCATCTgg < 1:1138064/76‑1 (MQ=255) cGTCTGGAAGATGGCTCTTATCGGCTGGTGGGGCATAAATGGTTTTTCTCGGTTCCGCAAAGCGATGCGCATCTgg < 1:1545823/76‑1 (MQ=255) cGTCTGGAAGATGGCTCTTATCGGCTGGGGGGGCATAAATGGTTTTTCTCGGTTCCGCAAAGCGATGCGCATCTgg > 1:105919/1‑76 (MQ=255) gTCTGGAAGATGGCTCTTATCGGCTTGTGGGGCATAAATGGTTTTTCTCGGTTCCGCAAAGCGATGCGCATCTGGt < 2:1371456/76‑1 (MQ=255) tCTGGAAGATGGCTCTTATCGGCTGGTGGGGCATAAATGGTTTTTCTCGGTTCCGCAAAGCGATGCGCATCTGGTg < 2:1249374/76‑1 (MQ=255) cTGGAAGATGGCTCTTATCGGCTGGGGGGGGATAAATGGGTTTTCTCGGGTCCGCAAAGCGATGCGCATCTGGTGc > 1:324361/1‑76 (MQ=255) tGGAAGATGGCTCTTATCGGCTGGGGGGGCATAAATGGTTTTTCTCGGTTCCGCCAAACGATGCGCCTCTGGTGCt > 1:792139/1‑76 (MQ=255) gAAGATGGCTCTTATCGGCTGGGGGGGCATAAATGGGTTTTCTCGGTTCCGCAAAGCGATGCGCCTCTGGTGCTgg > 1:965598/1‑76 (MQ=255) aaGATGGCTCTTATCGGCTGGGGGGGCATAAATGGTTTTTCTCGGTTCCGCAAAGCGATGCGCATCTGGTGCTgac > 1:856303/1‑74 (MQ=255) gATGGCTCTTATCGGCTGGTGGGGCATAAATGGTTTTTCTCGGTTCCGCAAAGCGATGCGCATCTGGTGCTGgcgc < 2:1376149/76‑1 (MQ=255) aTGGCTCTTATCGGCTGGTGGGGCATAAATGGTTTTTCTCGGTTCCGCAAAGCGATGCGCATCTGGTGCTGGCGCa < 1:538186/76‑1 (MQ=255) tGGCTCTTATCGGCTGGTGGGGCATAAATGGTTTTTCTCGGTTCCGCAAAGCGATGCGCATCTGGTGCTGGCGCAg < 1:1671050/76‑1 (MQ=255) tGGCTCTTATCGGCTGGGGGGGGCTAAATGGGTTTTTTCGGCTCCGCAAAGCGATGCGCAACTGGTGCCGGCGCAg > 2:1610656/1‑76 (MQ=255) cGGCTGGTGGGGCATAAATGGTTTTTCTCGGTTCCGCAAAGCGATGCGCATCTGGTGCTGGCGCAGACCGCGggtg > 1:855862/1‑76 (MQ=255) cGGCTGGTGGGGCATAAATGGTTTTTCTCGGTTCCGCAAAGCGATGCGCATCTGGTGCTGGCGCAGACCGCGggtg > 2:1451875/1‑76 (MQ=255) ggCTGGTGGGGCATAAATGGTTTTTCTCGGTTCCGCAAAGCGATGCGCATCTGGTGCTGGCGCAGACCGCGggtgg < 2:377330/76‑1 (MQ=255) gCTGGTGGGGCATAAATGGTTTTTCTCGGTTCCGCAAAGCGATGCGCATCTGGTGCTGGCGCAGACCGCGggtggt > 2:1197348/1‑76 (MQ=255) cTGGTGGGGCATAAATGGTTTTTCTCGGTTCCGCAAAGCGATGCGCATCTGGTGCTGGCGCAGACCGCGGGTGGTc > 1:1289749/1‑76 (MQ=255) cTGGTCGGGCATAAATGGTTTTTCTCGGTTCCGCAAAGCGATGCGCATCTGGTGCTGGCGCAGACCGCGGGTGGTc > 1:164459/1‑76 (MQ=255) tggtggGGCATAAATGGTTTTTCTCGGTTCCGCAAAGCGATGCGCATCTGGTGCTGGCGCAGACCGCGGGTGGTCt < 2:716863/76‑1 (MQ=255) ggtggGGCATAAATGGTTTTTCTCGGTTCCGCAAAGCGATGCGCATCTGGTGCTGGCGCAGACCGCGGGTGGTCTg > 1:1158140/1‑76 (MQ=255) | AAAGCAGGGCGGTTCCGATGTTATGAGCAACACCACCCGTGCAGAGCGTCTGGAAGATGGCTCTTATCGGCTGGTGGGGCATAAATGGTTTTTCTCGGTTCCGCAAAGCGATGCGCATCTGGTGCTGGCGCAGACCGCGGGTGGTCTG > REL606/4392987‑4393134 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 23 ≤ ATCG/ATCG < 30 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |