Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | AE015451 | 4,583,571 | 0 | G | . | 30.8% | 24.1 / 10.2 | 13 | intergenic (‑144/+27) | treSB/glgE | fused trehalose synthase B/maltokinase/Alpha‑1,4‑glucan:maltose‑1‑phosphate maltosyltransferase |
Reads supporting (aligned to +/- strand): ref base G (4/5); new base . (3/1); total (7/6) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 5.59e-01 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
GGGTCGTCGATGAAGGCTGCCGGGCGGGAACGCTTGGCCATGTGCGACTCCTTTCAATTGCATGCACAGACATCCGAGCTTGCTCGATCCCCCGTAGGAGCGGGTTTACCCGCGAAGATGACACCGCGGTGCCGAGGTGCATGGCACCGGCTTCGCCGGTGTTCGCGGATAAACCCGCTCCTACGGGGGGGGGGGCACAGTCTGAACCTTGTTAACGGGCTTTTTCGATGCGCCAGATGCCGAACGGCAGGTGCCAGGGCTCTATGCGCATCCATTGGGTCTTGCCGTACCACGACCAGCGGTGGCCGTTCATCAGGTCTTCACCCTGGGTTTCGGCGTTGTCGTCCAGCCCCAACTCCCACAGGGGCAACTCGAAGTGCGCCTCC > AE015451/4583387‑4583772 | gggTCGTCGATGAAGGCTGCCGGGCGGGAACGCTTGGCCATGTGCGACTCCTTTCAATTGCATGCACAGACATCCGAGCTTGCTCGATCCCCCGTAGGAGCGGGTTTACCCGCGAAGATGACACCGCGGTGCCGAGGTGCATGGCACCGGCTTCGCCGGTGTTCGCGGATAAACCCGCCCCTAC‑GGGGGGGGGGCACAGTCTGa > 1:32278/1‑204 (MQ=255) gcaCGCTTGGCCATGTGCGACTCCTTTCAATTGCATGCACAGACATCCGAGCTTGCTCGATCCCCCGTAGGAGCGGGTTTACCCGCGAAGATGACACCGCGGTGCCGAGGTGCATGGCACCGGCTTCGCCGGTGGTCGCGGATAAACCCGCTCCTACGGGGGGGGGGGcc > 2:139877/3‑169 (MQ=255) aCTCCTTTCAATTGCATGCACAGACATCCGAGCTTGCTCGATCCCCCGTAGGAGCGGGTTTACCCGCGAAGATGACACCGCGGTGCCGAGGTGCATGGCACCGGCTTCGCCGGTGTTCGCGGATAAACCCGCTCCTACGGGGGGGGGGGCACAGTCTGAACCTTGTTAACGGGCTTTTTCGATGCGCCAGATGCCGAACGGCAGGTGCCAGGGCTCTATGCGCATCCATTGGGTCtt < 1:171104/237‑1 (MQ=255) ggTTTACCCGCGAAGATGACACCGCGGTGCCGAGGTGCATGGCACCGGCTTCGCCGGTGTTCGCGGATAAACCCGCTCCTACGGGGGGGGGGGCACAGTCTGAACCTTGTTAACGGGCTTTTTCGATGCGCCAGATGCCGAACGGCAGGTGCCAGGGCTCTATGCGCATCCATTGGGTCTTGCCGTACCACGACCAGCGGTgg > 1:354339/1‑203 (MQ=255) ttACCCGCGAAGATGACACCGCGGTGCCGAGGTGCATGGCACCGGCTTCGCCGGTGTTCGCGGATAAACCCGCTCCTACGGGGGGGGGGGCACAGTCTGAACCTTGTTAACGGGCTTTTTCGATg > 1:92667/1‑125 (MQ=255) ggCACCGTCTTTGCCGGTGTTCGCGGATAAACCCGCTCCTAC‑GGGGGGGGGGCACAGTCTGaa > 2:345592/1‑63 (MQ=25) ggCTTCACCGGTGTTCGCGGGTAAACCCGCTCTTACGGGGGGGGGGGCACAGTCTGAACCTTGTTAACGGGCTTTTTCGATGCGCCAGATGCCGAACGGCAGGTGCCAGGGCTCTATg < 2:342801/118‑1 (MQ=255) gCTTCGCCGGTGTTCGCGGATAAACCCGCTCCTAC‑GGGGGGGGGGTACAGTCTGAACCtt > 2:173403/1‑60 (MQ=21) cTTCGCCGGTGTTCGCGGATAAACCCGCTCCTACGGGGGGGGGGGCACAGTCTGAACCTTGTTAACGGGCTTTTTCGATGCGCCAGATGCCGAACGGCAGGTGCCAGGGCTCTATGCGCATCCATTGGGTCTTGCCGTACCACGACCAGCGGTGGCCGTTCATCAGGTCTTCACCCTGGGTTTCGGCGTTGTCGTCCAGCCCCAACTCCCACAGGGGCAACTCGAAGTGCGcctcc > 1:248437/1‑236 (MQ=255) tgtTCGCGGATAACCCCGCTCCTAC‑GGGGGGGGGGCACAGTCTGAACCTTGTTAACGGGCTTTTTCGATGCGCCAGATGACGAACGTCAGGTGCCAGGGCTCTATGCGCATCCATTGGGTCTTGCCGTACCACGACCAGCGGt < 1:345594/143‑1 (MQ=255) cgcgGATAAACCCGCTCCTACGGGGGGGGGGGCACAGTCTGAACCTTGTTAACGGGCTTTTTCGATGCGCCAGATGCCGAACGGCAGGTGCCAGGGCTCTATGCGCATCCATTGGGTCTTGCCGTACCACGACCAGCGGTGGCCGTTCATCAGGTCTTCAc < 2:354337/161‑1 (MQ=255) aCGGGGGGGGGGGCACAGTCTGAACCTTGTTAACGGGCTTTTTCGATGCGCCAGATGCCGAACGGc < 2:275412/66‑1 (MQ=255) aCGGGGGGGGGGGCACAGTCTGAACCTTGTTAACGGGCTTTTTCGATGCGCCAGATGCCGAACGGCAGGTGCCAGGGCTCTATGCGCATCCATTGGGTCTTGCCGTACCACGACCAGCGGTGGCCGTTCATCAGGTCTTCACCCTGGGTTTCGGCGTTGTCGTCCAGCCCCAACTCCCACAgg < 1:209928/183‑1 (MQ=255) ggggggggggggCACAGTCTGAACCTTGTTAACGGGCTTTTTCGATGCGCCAGATGCCGAACGGCAGGTGCCAGGGCTCTATGCGCATCCATTGGGTCTTGCCGTACCACGACCAGCGGTGGCCGTtca < 2:248709/128‑1 (MQ=255) gggggggggggCACAGTCTGAACCTTGTTAACGGGCTTTTTCGATGCGCCAGATGCCGAACg < 1:70872/62‑1 (MQ=255) | GGGTCGTCGATGAAGGCTGCCGGGCGGGAACGCTTGGCCATGTGCGACTCCTTTCAATTGCATGCACAGACATCCGAGCTTGCTCGATCCCCCGTAGGAGCGGGTTTACCCGCGAAGATGACACCGCGGTGCCGAGGTGCATGGCACCGGCTTCGCCGGTGTTCGCGGATAAACCCGCTCCTACGGGGGGGGGGGCACAGTCTGAACCTTGTTAACGGGCTTTTTCGATGCGCCAGATGCCGAACGGCAGGTGCCAGGGCTCTATGCGCATCCATTGGGTCTTGCCGTACCACGACCAGCGGTGGCCGTTCATCAGGTCTTCACCCTGGGTTTCGGCGTTGTCGTCCAGCCCCAACTCCCACAGGGGCAACTCGAAGTGCGCCTCC > AE015451/4583387‑4583772 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 10 ≤ ATCG/ATCG < 23 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |