| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | REL606 | 4,393,061 | 0 | T | G | 43.8% | 38.6 / 16.0 | 33 | V211G (GTG→GGG) | aidB | isovaleryl CoA dehydrogenase |
| Reads supporting (aligned to +/- strand): ref base T (2/16); new base G (14/0); total (17/16) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 5.09e-07 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 | |||||||||||
| Rejected as consensus: Frequency below/above cutoff threshold. | |||||||||||
AGCAGGGCGGTTCCGATGTTATGAGCAACACCACCCGTGCAGAGCGTCTGGAAGATGGCTCTTATCGGCTGGTGGGGCATAAATGGTTTTTCTCGGTTCCGCAAAGCGATGCGCATCTGGTGCTGGCGCAGACCGCGGGTGGTCTGT > REL606/4392989‑4393135 | aGCAGGGCGGTTCCGATGTTATGAGCAACACCACCCGTGCAGAGCGTCTGGAAGATGGCTCTTATCGGCtggtgg < 2:183312/75‑1 (MQ=255)aGCAGGGCGGTTCCGATGTTATGAGCAACACCACCCGTGCAGAGCGTCTGGAAGATGGCTCTTATCGGCtggggg > 1:1225543/1‑75 (MQ=255) gCAGGGCGGTTCCGATGTTATGAGCAACACCACCCGTGCAGAGCGTCTGGAAGATGGCTCTTATCGGCTgggggg > 2:1288955/1‑75 (MQ=255) gCAGGGCGGTTCCGATGTTATGAGCAACACCACCCGTGCAGAGCGTCTGGAAGATGGCTCTTATCGGCTGGTggg < 1:153356/75‑1 (MQ=255) aGGGCGGTTCCGATGTTATGAGCAACACCACCCGTGCAGAGCGTCTGGAAGATGGCTCTTATCGGCTGGTGGGGc < 2:582071/75‑1 (MQ=255) gCGGTTCCGATGTTATGAGCAACACCACCCGTGCAGAGCGTCTGGAAGATGGCTCTTATCGGCTGGTGGGGCATa < 1:1395671/75‑1 (MQ=255) ggTTCCGATGTTATGAGCAACACCACCCGTGCAGAGCGTCTGGAAGATGGCTCTTATCGGCTGGGGGGGGGTaaa > 1:1208275/1‑75 (MQ=255) tttCCTATTTTATGAGCAACACCGCCCGTGCAGAGCGTCTGGAAGATGGCTCTTATCGGCTGGTGGGGCATAAAt < 2:493991/74‑1 (MQ=255) ttCCGATGTTATGAGCAACACCACCCGTGCAGAGCGTCTGGAAGATGGCTCTTATCGGCTGGGGGGGGATAAATg > 1:1598712/1‑75 (MQ=255) ttCCGATGTTATGAGCAACACCACCCGTGCAGAGCGTCTGGAAGATGGCTCTTATCGGCTGGGGGGGCATAAATg > 1:1595227/1‑75 (MQ=255) ccGATGTTATGAGCAACACCACCCGTGCAGAGCGTCTGGAAGATGGCTCTTATCGGCTGGTGGGGCATAAATGGt < 1:906305/75‑1 (MQ=255) gATGTTATGAGCAACACCACCCGTGCAGAGCGTCTGGAAGATGGCTCTTATCGGCTGGTGGGGCATAAATGGttt < 2:933206/75‑1 (MQ=255) gATGTTATGAGCAACACCACCCGTGCAGAGCGTCTGGAAGATGGCTCTTATCGGCTGGTGGGGCATAAATGGttt < 1:742697/75‑1 (MQ=255) tGTTATGAGCAACACCACCCGTGCAGAGCGTCTGGAAGATGGCTCTTATCGGCTGGTGGGGCATAAATGGttttt < 2:63985/75‑1 (MQ=255) tGTTATGAGCAACACCACCCGTGCAGAGCGTCTGGAAGATGGCTCTTATCGGCTGGGGGGGGATAAATGGttttt > 1:407886/1‑75 (MQ=255) tATGAGCAACACCACCCGTGCAGAGCGTCTGGAAGATGGCTCTTATCGGCTGGTGGGGCATAAATGGTTTTtctc < 1:1331441/75‑1 (MQ=255) tATGAGCAACACCACCCGTGCAGAGCGTCTGGAAGATGGCTCTTATCGGCTGGGGGGGGATAAATGGtttttttt > 2:475516/1‑72 (MQ=255) aaCACCACCCGTGCAGAGCGTCTGGAAGATGGCTCTTATCGGCTGGGGGGGCGTAAATGGTTTTTCTCCGTTCCg > 2:172710/1‑75 (MQ=255) acacCACCCGTGCAGAGCGTCTGGAAGATGGCTCTTATCGGCTGGTGGGGCATAAATGGTTTTTCTCGGTTCCGc < 1:174850/75‑1 (MQ=255) caccacCCGTGCAGAGCGTCTGGAAGATGGCTCTTATCGGCTGGGGGGGCATAAATGGTTTTTTTCGGTTCCGCa > 2:1547082/1‑75 (MQ=255) caccCGTGCAGAGCGTCTGGAAGATGGCTCTTATCGGCTGGTGGGGCATAAATGGTTTTTCTCGGTTCCGCAAAg < 1:689966/75‑1 (MQ=255) accCGTGCAGAGCGTCTGGAAGATGGCTCTTATCGGCTGGGGGGGCGTAAATGGTTTTTCTCTGGTCCGCaaaac > 2:310295/1‑73 (MQ=255) accCGTGCAGAGCGTCTGGAAGATGGCTCTTATCGGCTGGGGGGGCATAAATGGGTTTTCTCGGTTCCGCAAAGc > 1:1405989/1‑75 (MQ=255) cGTGCAGAGCGTCTGGAAGATGGCTCTTATCGGCTGGGGGGGGATAAATGGGTTTTTTCGGGTCCGCCAAGCGAt > 1:209849/1‑75 (MQ=255) gTGCAGAGCGTCTGGAAGATGGCTCTTATCGGCTGGTGGGGCATAAATGGCTTTTTTCTGGTCCCCCAACCGATg > 1:801331/1‑75 (MQ=255) agagCGTCTGGAAGATGGCTCTTATCGGCTGGTGGGGCATAAATGGTTTTTCTCGGTTCCGCAAAGCGATGCGCa < 1:1170399/75‑1 (MQ=255) agCGTCTGGAAGATGGCTCTTATCGGCTGGTGGGGCATAAATGGTTTTTCTCGGTTCCGCAAAGCGATGCGCATc < 1:572148/75‑1 (MQ=255) gTCTGGAAGATGGCTCTTATCGGCTGGTGGGGCATAAATGGTTTTTCTCGGTTCCGCAAAGCGATGCGCATCTgg < 2:1124811/75‑1 (MQ=255) gTCTGGAAGATGGCTCTTATCGGCTGGTGGGGCATAAATGGTTTTTCTCGGTTCCGCAAAGCGATGCGCATCTgg < 1:1288956/75‑1 (MQ=255) tCTGGAAGATGGCTCTTATCGGCTGGGGGGGCATAAATGGGTTTTCTCGGTTTCGCAAAGCGATGCGCATCCggg > 1:496491/1‑74 (MQ=255) cTGGAAGATGGCTCTTATCGGCTGGGGGGGGATAAATGGTTTTTCTCGGTTCCGCAAAGCGATGCGCATCTgggg > 2:1561017/1‑73 (MQ=255) gAAGATGGCTCTTATCGGCTGGTGGGGCATAAATGGTTTTTCTCGGTTCCGCAAAGCGATGCGCATCTGGTGCTg < 2:1546348/75‑1 (MQ=255) gAAGATGGCTCTTATCGGCTGGGGGGGCATAAATGGTTTTTCTCGGTTCCGCCAAGCGATGCGCATCTGGGGCTg > 1:22623/1‑75 (MQ=255) aGATGGCTCTTATCGGCTGGGGGGGCATAAATGGTTTTTCTCGGTTCCGCAAAGCGGTGCGCATCTGGTGCTGgc > 2:368777/1‑75 (MQ=255) aGATGGCTCTTATCGGCTGGCGGGGCATAAATGGCTTTTTTCGGTTCCCCAAAGCGGGGCGCATCTGGGGTTggg > 2:1423836/1‑74 (MQ=255) ttATCGGCTGGTGGGGGATAAATGGTTTTTCTCGGTTCCGCAAAGCGATGCGCATCTGGTGCTGGCGCAGACcgc > 2:1370421/1‑75 (MQ=255) aTCGGCTGGTGGGGCATAAATGGGTTTTCTTCGTTCCCCAAACCGATGCGGCTCTTGCGCTGGCGCGCACCGCgg > 2:1124474/1‑75 (MQ=255) tggtggGGCATAAATGGTTTTTCTCGGTTCCGCAAAGCGATGCGCATCTGGTGCTGGCGCAGACCGCGGGTGGTc < 1:964372/75‑1 (MQ=255) tggtggGGCATAAATGGTTTTTCTCGGTTCCGCAAAGCGATGCGCATCTGGTGCTGGCGCAGACCGCGGGTGGTc < 1:1033592/75‑1 (MQ=255) ggtggGGCATAAATGGTTTTTCTCGGTTCCGCAAAGCGATGCGCATCTGGTGCTGGCGCAGACCGCGGGTGGTCt > 1:1209450/1‑75 (MQ=255) tggGGCATAAATGGTTTTTCTCGGTTCCGCAAAGCGATGCGCATCTGGTGCTGGCGCAGACCGCGGGTGGTCTGt > 2:899938/1‑75 (MQ=255) | AGCAGGGCGGTTCCGATGTTATGAGCAACACCACCCGTGCAGAGCGTCTGGAAGATGGCTCTTATCGGCTGGTGGGGCATAAATGGTTTTTCTCGGTTCCGCAAAGCGATGCGCATCTGGTGCTGGCGCAGACCGCGGGTGGTCTGT > REL606/4392989‑4393135 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 16 ≤ ATCG/ATCG < 26 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |