Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | REL606 | 4,393,061 | 0 | T | G | 46.7% | 29.6 / 14.7 | 30 | V211G (GTG→GGG) | aidB | isovaleryl CoA dehydrogenase |
Reads supporting (aligned to +/- strand): ref base T (3/13); new base G (14/0); total (17/13) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 4.79e-06 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
AAAGCAGGGCGGTTCCGATGTTATGAGCAACACCACCCGTGCAGAGCGTCTGGAAGATGGCTCTTATCGGCTGGTGGGGCATAAAT‑GGTTTTTCTCGGTTCCGCAAAGCGATGCGCATCTGGTGCTGGCGCAGACCGCGGGTGGTCTGTC > REL606/4392987‑4393136 | aaaGCAGGGCGGTTCCGATGTTATGAGCAACACCACCCGTGCAGAGCGTCTGGAAGATGGCTCTTATCGGCtggtg < 1:1709552/76‑1 (MQ=255) cAGGGCGGTTCCGATGTTATGAGCAACACCACCCGTGCAGAGCGTCTGGAAGATGGCTCTTATCGGCTGggggggg > 2:702162/1‑75 (MQ=255) cAGGGCGGTTCCGATGTTATGAGCAACACCACCCGTGCAGAGCGTCTGGAAGATGGCTCTTATCGGCTGGTGGGGc < 1:1656053/76‑1 (MQ=255) aGGGCGGTTCCGATGTTATGAGCAACACCACCCGTGCAGAGCGTCTGGAAGATGGCTCTTATCGGCTGGGGGGGCa > 2:301006/1‑76 (MQ=255) ttCCGATGTTATGAGCAACACCACCCGTGCAGAGCGTCTGGAAGATGGCTCTTATCGGCTGGTGGGGCATAAAT‑gg < 1:41153/76‑1 (MQ=255) ccGATGTTATGAGCAACACCACCCGTGCAGAGCGTCTGGAAGATGGCTCTTATCGGCTGGTGGGGCATAAAT‑GGtt < 1:1491850/76‑1 (MQ=255) gATGTTATGAGCAACACCACCCGTGCAGAGCGTCTGGAAGATGGCTCTTATCGGCTGGTGGGGCATAAAT‑GGtttt < 1:926133/76‑1 (MQ=255) gATGTTATGAGCAACACCACCCGTGCAGAGCGTCTGGAAGATGGCTCTTATCGGCTGGGGGGGGGTAAAT‑GGgttt > 1:696464/1‑76 (MQ=255) aTGTTATGAGCAACACCCCCCGTGCAGAGCGCCTGGAAGATGGCTCTTATCGGCTGGTGGGGCATAAAT‑GGttttt < 2:229946/76‑1 (MQ=255) tGTTATGAGCAACACCACCCGTGCAGAGCGTCTGGAAGATGGCTCTTATCGGCTGGGGGGGGGTAAAT‑GGtttttt > 2:850670/1‑75 (MQ=255) gTTATGAGCAACACCACCCGTGCAGAGCGTCTGGAAGATGGCTCTTATCGGCTGGGGGGGGATAAAT‑GGgtttttt > 1:1588936/1‑74 (MQ=255) ttATGAGCAACACCACCCGTGCAGAGCGTCTGGAAGATGGCTCTTATCGGCTGGGGGGGGATAAATAGGttttttc > 1:622055/1‑74 (MQ=255) ttATGAGCAACACCACCCGTGCAGAGCGTCTGGAAGATGGCTCTTATCGGCTGGGGGGGCGTAAAT‑GGGTTTtttc > 2:323132/1‑76 (MQ=255) tGAGCAACACCACCCGTGCAGAGCGTCTGGAAGATGGCTCTTATCGGCTGGTGGGGCATAAAT‑GGTTTTTCTCGGt < 1:1320376/76‑1 (MQ=255) aaCACCACCCGTGCAGAGCGTCTGGAAGATGGCTCTTATCGGCTGGGGGGGGATAAAT‑GGTTTTTTTTGGGTcctc > 2:860751/1‑74 (MQ=255) accaccCGTGCAGAGCGTCTGGAAGATGGCTCTTATCGGCTGGGGGGGCATAAAT‑GGTTTTTCTCGGTTCCGCaaa > 1:1203841/1‑76 (MQ=255) ccaccCGTGCAGAGCGTCTGGAAGATGGCTCTTATCGGCTGGGGGGGGGTAAAT‑GGGTTTTTTTGGTTTCGCCAAg > 2:1091804/1‑76 (MQ=255) caccCGTGCAGAGCGTCTGGAAGATGGCTCTTATCGGCTGGGGGGGGATAAAT‑GGGTTTTTTCGGTTCCGCAAAGc > 1:1326580/1‑76 (MQ=255) accCGTGCAGAGCGTCTGGAAGATGGCTCTTATCGGCTGGTGGGGCGTAAAT‑GGGTTTTTTCGGGTCCGCCAAGCg > 2:515780/1‑76 (MQ=255) tGCAGAGCGTCTGGAAGATGGCTCTTATCGGCTGGTGGGGCATAAAT‑GGTTTTTCTCGGTTCCGCAAAGCGATgcg < 2:820426/76‑1 (MQ=255) cAGAGCGTCTGGAAGATGGCTCTTATCGGCTGGGGGGGGATAAAT‑GGTTTTTTTTGGTTCCGCCAAGCGATGCGCa > 1:807589/1‑76 (MQ=255) agagCGTCTGGAAGATGGCTCTTATCGGCTGGTGGGGCATAAAT‑GGTTTTTCTCGGTTCCGCAAAGCGATGCGCAt < 2:541369/76‑1 (MQ=255) agagCGTCTGGAAGATGGCTCTTATCGGCTGGGGGGGGATAAAA‑GGGTTTTCTCGGTGCCCCAAAGGGATGCGCAt > 2:1369802/1‑76 (MQ=255) gagCTTCTGGAAGATGGCTCTTATCGGCTGGTGGGGCATAAAT‑GGTTTTTCTCGGTTCCGCAAAGCGATGCGCATc < 1:432348/76‑1 (MQ=255) agCGTCTGGAAGATGGCTCTTATCGGCTGGGGGGGGATAAAT‑GGTTTTTTTCGGGTCCGCCAAGCGATGCGCATCt > 2:1381208/1‑76 (MQ=255) cGTCTGGAAGATGGCTCTTATCGGCTGGGGGGGGATAAAT‑GGGTTTTTTCGGTTCCGCAAAGCCATGCGCCTCTgg > 1:37663/1‑76 (MQ=255) gTCTGGAAGATGGCTCTTATCGGCTGGGGGGGCATAAAT‑GGTTTTTTTCGGTTTCGCCAAGCGATGCGCATCTggg > 1:820817/1‑75 (MQ=255) cTGGAAGATGGCTCTTATCGGCTGGGGGGGCATAAAT‑GGTTTTTCTCGGTTCCGCAAAGCGATGCGCATCTGGTGc > 2:299880/1‑76 (MQ=255) ggAAGATGGCTCTTATCGGCTGGTGGGGCATAAAT‑GGTTTTTCTCGGTTCCGCAAAGCGATGCGCATCTGGTGCTg < 1:523589/76‑1 (MQ=255) aaGATGGCTCTTATCGGCTGGTGGGGGATAAAT‑GGGTTTTTTCGGGTTCGCGAAGCGATGCGCATCTGGGGTTggg > 1:99346/1‑75 (MQ=255) aGATGGCTCTTATCGGCTGGTGGGGCATAAAT‑GGTTTTTCTCGGTTCCGCAAAGCGATGCGCATCTGGTGCTGgcg > 1:1252864/1‑76 (MQ=255) tGGCTCTTATCGGCTGGGGGGGCATAAAT‑GGTTTTTTTCGGTTCCGCAAAGCGATGCGCATCTGGGGCTGGCGCAg > 2:479461/1‑76 (MQ=255) ggCTCTTATCGGCTGGGGGGGGATAAAT‑GGGTTTTCTCGGTTCCGCAAAGCGATGCGCATCTGGTGCTGgcgccgg > 2:882402/1‑73 (MQ=255) ctctTATCGGCTGGTGGGGCATAAAT‑GGTTTTTCTCGGTTCCGCAAAGCGATGCGCATCTGGTGCTGGCGCAGAcc < 2:32251/76‑1 (MQ=255) ctctTATCGGCTGGTGGGGCATAAAT‑GGTTTTTCTCGGTTCCGCAAAGCGATGCGCATCTGGTGCTGGCGCAGAcc < 2:1846315/76‑1 (MQ=255) cGGCTGGTGGGGCATAAAT‑GGTTTTTCTCGGTTCCGCAAAGCGATGCGCATCTGGTGCTGGCGCAGACCGCGggtg < 2:1830608/76‑1 (MQ=255) cTGGTGGGGCATAAAT‑GGTTTTTCTCGGTTCCGCAAAGCGATGCGCATCTGGTGCTGGCGCAGACCGCGGGTGGTc > 1:1507336/1‑76 (MQ=255) gtggGGCATAAAT‑GGTTTTTCGCGGTTCCGCAAAGCGATGCGCATCTGGTGCTGGCGCAGACCGCGGGTGGTCTGt < 2:103875/76‑1 (MQ=255) tggGGCATAAAT‑GGTTTTTCTCGGTTCCGCAAAGCGATGCGCATCTGGTGCTGGCGCAGCCCGCGGGTGGTCTGTc > 2:1290590/1‑76 (MQ=255) | AAAGCAGGGCGGTTCCGATGTTATGAGCAACACCACCCGTGCAGAGCGTCTGGAAGATGGCTCTTATCGGCTGGTGGGGCATAAAT‑GGTTTTTCTCGGTTCCGCAAAGCGATGCGCATCTGGTGCTGGCGCAGACCGCGGGTGGTCTGTC > REL606/4392987‑4393136 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 20 ≤ ATCG/ATCG < 28 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |