Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_000913 | 4,492,578 | 0 | G | . | 62.5% | 29.1 / 82.9 | 48 | intergenic (‑53/+9) | idnT/idnO | L‑idonate/5‑ketogluconate/gluconate transporter/5‑keto‑D‑gluconate 5‑reductase |
Reads supporting (aligned to +/- strand): ref base G (7/11); new base . (16/14); total (23/25) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 3.83e-01 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.57e-01 | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
AACGAGGGCAATAAAGCCGTTAACTTTAAAGCCGATCATCAGGATAAGAAGCAGCGCGACGCCTGCCGCAATAATGATTAATGGCATAGTCTTTATCTCTTATTAGTGCCAGAAAAGCAGCGCTTTCGCGCTGCACGTAGGGGGGCTGTTAAACAGCCACTAACATGCCGCCATCAACAAACAACAGGTGGCCGTTTACGAAATCAGAGGCTTTTGAAGAAAGGAACACCGCAGCACCAATCAGTTCCTGCGGATCTCCCCAGCGTGCGGCGGGGG > NC_000913/4492439‑4492714 | aaCGAGGGCAATAAAGCCGTTAACTTTAAAGCCGATCATCAGGATAAGAAGCAGCGCGACGCCTGCCGCAATAATGATTAATGGCATAGTCTTTATCTCTTATTAGTGCCAGAAAAGCAGCGCTTTCGCGCTGCACGTAg > 1:834451/1‑140 (MQ=255) gAGGGCAATAAAGCCGTTAACTTTAAAGCCGATCATCAGGATAAGAAGCAGCGCGACGCCTGCCGCAATAATGATTAATGGCATAGTCTTTATCTCTTATTAGTGCCAGAAAAGCAGCGCTTTCGCGCTGCACGTAgggg < 1:1457781/140‑1 (MQ=255) aGGGCAATAAAGCCGTTAACTTTAAAGCCGATCATCAGGATAAGAAGCAGCGCGACGCCTGCCGCAATAATGATTAATGGCATAGTCTTTATCTCTTATTAGTGCCAGAAAAGCAGCGCTTTCGCGCTGCACGTAggggg > 1:764399/1‑140 (MQ=255) gggCAATAAAGCCGTTAACTTTAAAGCCGATCATCAGGATAAGAAGCAGCGCGACGCCTGCCGCAATAATGATTAATGGCATAGTCTTTATCTCTTATTAGTGCCAGAAAAGCAGCGCTTTCGCGCTGCACGTAgggggc > 2:464480/1‑139 (MQ=255) ggCAATAAAGCCGTTAACTTTAAAGCCGATCATCAGGATAAGAAGCAGCGCGACGCCTGCCGCAATAATGATTAATGGCATAGTCTTTATCTCTTATTAGTGCCAGAAAAGCAGCGCTTTCGCGCTGCACGTAggg > 1:1087480/1‑136 (MQ=255) ggCAATAAAGCCGTTAACTTTAAAGCCGATCATCAGGATAAGAAGCAGCGCGACGCCTGCCGCAATAATGATTAATGGCATAGTCTTTATCTCTTATTAGTGCCAGAAAAGCAGCGCTTTCGCGCTGCACGTAggg < 2:1087480/136‑1 (MQ=255) aTAAAGCCGTTAACTTTAAAGCCGATCATCAGGATAAGAAGCAGCGCGACGCCTGCCGCAATAATGATTAATGGCATAGTCTTTATCTCTTATTAGTGCCAGAAAAGCAGCGCTTTCGCGCTGCACGTA‑GGGGGCTGTTa > 2:989362/1‑140 (MQ=255) aaaGCCGTTAACTTTAAAGCCGATCATCAGGATAAGAAGCAGCGCGACGCCTGCCGCAATAATGATTAATGGCATAGTCTTTATCTCTTATTAGTGCCAGAAAAGCAGCGCTTTCGCGCTGCACGTA‑GGGGGCTGTTaaa < 1:989362/140‑1 (MQ=255) aaGCCGTTAACTTTAAAGCCGATCATCAGGATAAGAAGCAGCGCGACGCCTGCCGCAATAATGATTAATGGCATAGTCTTTATCTCTTATTAGTGCCAGAAAAGCAGCGCTTTCGCGCTGCACGTAGGGGGGCTGTTaaa < 1:543779/140‑1 (MQ=255) aaGCCGTTAACTTTAAAGCCGATCATCAGGATAAGAAGCAGCGCGACGCCTGCCGCAATAATGATTAATGGCATAGTCTTTATCTCTTATTAGTGCCAGAAAAGCAGCGCTTTCGCGCTGCACGTAGGGGGGCTGTTaaa < 2:1360367/140‑1 (MQ=255) aaGCCGTTAACTTTAAAGCCGATCATCAGGATAAGAAGCAGCGCGACGCCTGCCGCAATAATGATTAATGGCATAGTCTTTATCTCTTATTAGTGCCAGAAAAGCAGCGCTTTCGCGCTGCACGTA‑GGGGGCTGTTAAAc < 2:455602/140‑1 (MQ=255) aaaGCCGATCATCAGGATAAGAAGCAGCGCGACGCCTGCCGCAATAATGATTAATGGCATAGTCTTTATCTCTTATTAGTGCCAGAAAAGCAGCGCTTTCGCGCTGCACGTAGGGGGGCTGTTAAACAGCCACTAACATg < 1:1129768/140‑1 (MQ=255) atcatcAGGATAAGAAGCAGCGCGACGCCTGCCGCAATAATGATTAATGGCATAGTCTTTATCTCTTATTAGTGCCAGAAAAGCAGCGCTTTCGCGCTGCACGTA‑GGGGGCTGTTAAACAGCCACTAACATGCCGCCATc > 2:443234/1‑140 (MQ=255) atcatcAGGATAAGAAGCAGCGCGACGCCTGCCGCAATAATGATTAATGGCATAGTCTTTATCTCTTATTAGTGCCAGAAAAGCAGCGCTTTCGCGCTGCACGTA‑GGGGGCTGTTAAACAGCCACTAACATGCCGCCATc > 1:1056921/1‑140 (MQ=255) atcatcAGGATAAGAAGCAGCGCGACGCCTGCCGCAATAATGATTAATGGCATAGTCTTTATCTCTTATTAGTGCCAGAAAAGCAGCGCTTTCGCGCTGCACGTA‑GGGGGCTGTTAAACAGCCACTAACATGCCGCCATc > 1:617300/1‑140 (MQ=255) aGGATAAGAAGCAGCGCGACGCCTGCCGCAATAATGATTAATGGCATAGTCTTTATCTCTTATTAGTGCCAGAAAAGCAGCGCTTTCGCGCTGCACGTA‑GGGGGCTGTTAAACAGCCACTAACATGCCACCATCaacaaa < 1:1259374/140‑1 (MQ=255) gaagCAGCGCGACGCCTGCCGCAATAATGATTAATGGCATAGTCTTTATCTCTTATTAGTGCCAGAAAAGCAGCGCTTTCGCGCTGCACGTAGGGGGGCTGTTAAACAGCCACTAACATGCCGCCATCAACAAACAACAg > 1:867454/1‑140 (MQ=255) gcagcGCGACGCCTGCCGCAATAATGATTAATGGCATAGTCTTTATCTCTTATTAGTGCCAGAAAAGCAGCGCTTTCGCGCTGCACGTA‑GGGGGCTGTTAAACAGCCACTAACATGCCGCCATCAACAAACAACAGGTgg < 2:617300/140‑1 (MQ=255) agcGCGACGCCTGCCGCAATAATGATTAATGGCATAGTCTTTATCTCTTATTAGTGCCAGAAAAGCAGCGCTTTCGCGCTGCACGTAGGGGGGCTGTTAAACAGCCACTAACATGCCGCCATCaacaaa > 1:1167253/1‑129 (MQ=255) agcGCGACGCCTGCCGCAATAATGATTAATGGCATAGTCTTTATCTCTTATTAGTGCCAGAAAAGCAGCGCTTTCGCGCTGCACGTAGGGGGGCTGTTAAACAGCCACTAACATGCCGCCATCaacaaa < 1:1254164/129‑1 (MQ=255) agcGCGACGCCTGCCGCAATAATGATTAATGGCATAGTCTTTATCTCTTATTAGTGCCAGAAAAGCAGCGCTTTCGCGCTGCACGTAGGGGGGCTGTTAAACAGCCACTAACATGCCGCCATCaacaaa > 2:1254164/1‑129 (MQ=255) agcGCGACGCCTGCCGCAATAATGATTAATGGCATAGTCTTTATCTCTTATTAGTGCCAGAAAAGCAGCGCTTTCGCGCTGCACGTAGGGGGGCTGTTAAACAGCCACTAACATGCCGCCATCaacaaa < 2:1167253/129‑1 (MQ=255) agcGCGACGCCTGCCGCAATAATGATTAATGGCATAGTCTTTATCTCTTATTAGTGCCAGAAAAGCAGCGCTTTCGCGCTGCACGTA‑GGGGGCTGTTAAACAGCCACTAACATGCCGCCATCAACAAACAACAGGTGGcc > 1:704414/1‑140 (MQ=255) cgACGCCTGCCGCAATAATGATTAATGGCATAGTCTTTATCTCTTATTAGTGCCAGAAAAGCAGCGCTTTCGCGCTGCACGTA‑GGGGGCTGTTAAACAGCCACTAACATGCCGCCATCAACAAACAACAGGTGGCCGttt > 2:697110/1‑140 (MQ=255) cTGCCGCAATAATGATTAATGGCATAGTCTTTATCTCTTATTAGTGCCAGAAAAGCAGCGCTTTCGCGCTGCACGTAGGGGGGCTGTTAAACAGCCACTAACATGCCGCCATCAACAAACAACAGGTGGCCGTTTACGaa < 2:490330/140‑1 (MQ=255) ccGCAATAATGATTAATGGCATAGTCTTTATCTCTTATTAGTGCCAGAAAAGCAGCGCTTTCGCGCTGCACGTA‑GGGGGCTGTTAAACAGCCACTAACATGCCGCCATCAACAAACAACAGGTgg > 1:798451/1‑125 (MQ=255) ccGCAATAATGATTAATGGCATAGTCTTTATCTCTTATTAGTGCCAGAAAAGCAGCGCTTTCGCGCTGCACGTA‑GGGGGCTGTTAAACAGCCACTAACATGCCGCCATCAACAAACAACAGGTgg < 2:798451/125‑1 (MQ=255) cAATAATGATTAATGGCATAGTCTTTATCTCTTATTAGTGCCAGAAAAGCAGCGCTTTCGCGCTGCACGTAGGGGGGCTGTTAAACAGCCACTAACATGCCGCCATCAACAAACAACAGGTGGCCGTTTACGAAATCaga < 2:867454/140‑1 (MQ=255) tAATGGCATAGTCTTTATCTCTTATTAGTGCCAGAAAAGCAGCGCTTTCGCGCTGCACGTAGGGGGGCTGTTAAACAGCCACTAACATGCCGCCATCAACAAACAACAGGTGGCCGTTTACGAAATCAGAGGCTTTTgaa < 1:930713/140‑1 (MQ=255) aaTGGCATAGTCTTTATCTCTTATTAGTGCCAGAAAAGCAGCGCTTTCGCGCTGCACGTA‑GGGGGCTGTTAAACAGCCACTAACATGCCGCCATCAACAAACAACAGGTGGCCGTTTACGAAATCAGAGGCTTTTgaaga > 1:534080/1‑140 (MQ=255) aTGGCATAGTCTTTATCTCTTATTAGTGCCAGAAAAGCAGCGCTTTCGCGCTGCACGTA‑GGGGGCTGTTAAACAGCCACTAACATGCCGCCATCAACAAACAACAGGTGGCCGTTTACGAAATCAGAGGCTTTTgaagaa < 1:1001971/140‑1 (MQ=255) aTAGTCTTTATCTCTTATTAGTGCCAGAAAAGCAGCGCTTTCGCGCTGCACGTA‑GGGGGCTGTTAAACAGCCACTAACATGCCGCCATCAACAAACAACAGGTGGCCGTTTACGAAATCAGAGGCTTTTGAAGAAAGGaa < 2:704414/140‑1 (MQ=255) aTAGTCTTTATCTCTTATTAGTGCCAGAAAAGCAGCGCTTTCGCGCTGCACGTA‑GGGGGCTGTTAAACAGCCACTAACATGCCGCCATCAACAAACAACAGGTGGCCGTTTACGAAATCAGAGGCTTTTGAAGAAAGGaa < 2:534080/140‑1 (MQ=255) tAGTCTTTATCTCTTATTAGTGCCAGAAAAGCAGCGCTTTCGCGCTGCACGTAGGGGGGCTGTTAAACAGCCACTAACATGCCGCCATCAACAAACAACAGGTGGCCGTTTACGAAATCAGAGGCTTTTGAAGAAAGGaa < 1:955408/140‑1 (MQ=255) tAGTCTTTATCTCTTATTAGTGCCAGAAAAGCAGCGCTTTCGCGCTGCACGTA‑GGGGGCTGTTAAACAGCCACTAACATGCCGCCATCAACAAACAACAGGTGGCCGTTTACGAAATCAGAGGCTTTTGAAGAAAGGAac > 1:727698/1‑140 (MQ=255) tATCTCTTATTAGTGCCAGAAAAGCAGCGCTTTCGCGCTGCACGTA‑GGGGGCTGTTAAACAGCCACTAACATGCCGCCATCAACAAACAACAGGTGGCCGTTTACGAAATCAGAGGCTTTTGAAGAAAGGAACACCgcag > 1:789352/1‑140 (MQ=255) tATCTCTTATTAGTGCCAGAAAAGCAGCGCTTTCGCGCTGCACGTA‑GGGGGCTGTTAAACAGCCACTAACATGCCGCCATCAACAAACAACAGGTGGCCGTTTACGAAATCAGAGGCTTTTGAAGAAAGGAACACCgcag < 1:464480/140‑1 (MQ=255) aTCTCTTATTAGTGCCAGAAAAGCAGCGCTTTCGCGCTGCACGTAGGGGGGCTGTTAAACAGCCACTAACATGCCGCCATCAACAAACAACAGGTGGCCGTTTACGAAATCAGAGGCTTTTGAAGAAAGGAACACCgcag < 1:1062656/140‑1 (MQ=255) ttattaGTGCCAGAAAAGCAGCGCTTTCGCGCTGCACGTAGGGGGGCTGTTAAACAGCCACTAACATGCCGCCATCAACAAACAACAGGTGGCCGTTTACGAAATCAGAGGCTTTTGAAGAAAGGAACACCGCAGCACCa > 2:1182539/1‑140 (MQ=255) ttattaGTGCCAGAAAAGCAGCGCTTTCGCGCTGCACGTAGGGGGGCTGTTAAACAGCCACTAACATGCCGCCATCAACAAACAACAGGTGGCCGTTTACGAAATCAGAGGCTTTTGAAGAAAGGAACACCGCAGCACCa > 1:556854/1‑140 (MQ=255) ttattaGTGCCAGAAAAGCAGCGCTTTCGCGCTGCACGTAGGGGGGCTGTTAAACAGCCACTAACATGCCGCCATCAACAAACAACAGGTGGCCGTTTACGAAATCAGAGGCTTTTGAAGAAAGGAACACCGCAGCACCa > 1:1244774/1‑140 (MQ=255) ttattaGTGCCAGAAAAGCAGCGCTTTCGCGCTGCACGTA‑GGGGGCTGTTAAACAGCCACTAACATGCCGCCATCAACAAACAACAGGTGGCCGTTTACGAAATCAGAGGCTTTTGAAGAAAGGAACACCGCAGCACCaa > 1:1021667/1‑140 (MQ=255) attaGTGCCAGAAAAGCAGCGCTTTCGCGCTGCACGTA‑GGGGGCTGTTAAACAGCCACTAACATGCCGCCATCAACAAACAACAGGTGGCCGTTTACGAAATCAGAGGCTTTTTAAGAAAGGAACACCGCAGCACCAATc > 2:1404284/1‑140 (MQ=255) attaGTGCCAGAAAAGCAGCGCTTTCGCGCTGCACGTA‑GGGGGCTGTTAAACAGCCACTAACATGCCGCCATCAACAAACAACAGGTGGCCGTTTACGAAATCAGAGGCTTTTGAAGAAAGGAACACCg > 2:362341/1‑129 (MQ=255) attaGTGCCAGAAAAGCAGCGCTTTCGCGCTGCACGTA‑GGGGGCTGTTAAACAGCCACTAACATGCCGCCATCAACAAACAACAGGTGGCCGTTTACGAAATCAGAGGCTTTTGAAGAAAGGAACACCg < 1:362341/129‑1 (MQ=255) aGTGCCAGAAAAGCAGCGCTTTCGCGCTGCACGTA‑GGGGGCTGTTAAACAGCCACTAACATGCCGCCATCAACAAACAACAGGTGGCCGTTTACGAAATCAGAGGCTTTTGAAGAAAGGAACACCGCAGCACCAATCAGt < 2:789352/140‑1 (MQ=255) cAGAAAAGCAGCGCTTTCGCGCTGCACGTAGGGGGGCTGTTAAACAGCCACTAACATGCCGCCATCAACAAACAACAGGTGGCCGTTTACGAAATCAGAGGCTTTTGAAGAAAGGAACACCGCAGCACCAATCAGTTCCt > 1:263959/1‑140 (MQ=255) aGAAAAGCAGCGCTTTCGCGCTGCACGTAGGGGGGCTGTTAAACAGCCACTAACATGCCGCCATCAACAAACAACAGGTGGCCGTTTACGAAATCAGAGGCTTTTGAAGAAAGGAACACCGCAGCACCAATCAGTTCCTg < 2:263959/140‑1 (MQ=255) gAAAAGCAGTGCTTTCGCGCTGCACGTA‑GGGGGCTGTTAAACAGCCACTAACATGCCGCCATCAACAAACAACAGGTGGCCGTTTACGAAATCAGAGGCTTTTGAAGAAAGGAACACCGCAGCACCAATCAGTTCCTGCg > 2:667958/1‑140 (MQ=255) gAAAAGCAGTGCTTTCGCGCTGCACGTA‑GGGGGCTGTTAAACAGCCACTAACATGCCGCCATCAACAAACAACAGGTGGCCGTTTACGAAATCAGAGGCTTTTGAAGAAAGGAACACCGCAGCACCAATCAGTTCCTGCg < 1:667958/140‑1 (MQ=255) aaaaGCAGCGCTTTCGCGCTGCACGTA‑GGGGGCTGTTAAACAGCCACTAACATGCCGCCATCAACAAACAACAGGTGGCCGTTTACGAAATCAGAGGCTTTTGAAGAAAGGAACACCGCAGCACCAATCAGTTCCTGCgg > 1:677342/1‑140 (MQ=255) gcgcTTTCGCGCTGCACGTA‑GGGGGCTGTTAAACAGCCACTAACATGCCGCCATCAACAAACAACAGGTGGCCGTTTACGAAATCAGAGGCTTTTGAAGAAAGGAACACCGCAGCACCAATCAGTTCCTGCGGATCTccc < 1:443234/140‑1 (MQ=255) gcTTTCGCGCTGCACGTA‑GGGGGCTGTTAAACAGCCACTAACATGCCGCCATCAACAAACAACAGGTGGCCGTTTACGAAATCAGAGGCTTTTGAAGAAAGGAACACCGCAGCACCAATCAGTTCCTGCGGATCTCCCCa > 2:643398/1‑140 (MQ=255) gcTTTCGCGCTGCACGTA‑GGGGGCTGTTAAACAGCCACTAACATGCCGCCATCAACAAACAACAGGTGGCCGTTTACGAAATCAGAGGCTTTTGAAGAAAGGAACACCGCAGCACCAATCAGTTCCTGCGGATCTCCCCa < 2:727698/140‑1 (MQ=255) agaGGGGGGTGGTTAACCAGCCACTAACATGCCGCCATCAACAAACAACAGGTGGCCGTTTACGAAATCAGAGGCTTTTGAAGAAAGGAACACCGCAGCACCAACCAGTTCCTGCGGATCTCCCCAGCGTGCGGCggggg < 2:417431/138‑1 (MQ=255) | AACGAGGGCAATAAAGCCGTTAACTTTAAAGCCGATCATCAGGATAAGAAGCAGCGCGACGCCTGCCGCAATAATGATTAATGGCATAGTCTTTATCTCTTATTAGTGCCAGAAAAGCAGCGCTTTCGCGCTGCACGTAGGGGGGCTGTTAAACAGCCACTAACATGCCGCCATCAACAAACAACAGGTGGCCGTTTACGAAATCAGAGGCTTTTGAAGAAAGGAACACCGCAGCACCAATCAGTTCCTGCGGATCTCCCCAGCGTGCGGCGGGGG > NC_000913/4492439‑4492714 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |