| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | NC_000913 | 3,742,599 | 0 | C | . | 65.7% | 31.8 / 51.9 | 35 | intergenic (‑67/‑134) | yiaJ/yiaK | transcriptional repressor for the yiaKLMNO‑lyxK‑sgbHUE operon/2,3‑diketo‑L‑gulonate reductase, NADH‑dependent |
| Reads supporting (aligned to +/- strand): ref base C (7/5); new base . (16/7); total (23/12) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 7.09e-01 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.75e-01 | |||||||||||
| Rejected as consensus: Frequency below/above cutoff threshold. | |||||||||||
CTTCCGGCTGGACGCTCTTTTTCCTGCGCCATCTCGTTCTCTTTTTTTCCCATCACTTCTTTCCCCATTTTGTCGCGTCCTGATGGTAGCGCAAAGTGTGCCGTAGTTCACGATCTCGACAGATAATTTATAACCAATTGATTTTTATGTCTTTTGAAATTCATCAATCAGATTGCCTTGTTAAAAAGTGATCGATATA > NC_000913/3742465‑3742663 | cTTCCGGCTGGACGCTCTTTTTCCTGCGCCATCTCGTTCTCTTTTTTTCCCATCACTTCTTTCCCCATTTTGTCGCGTCCTGATGGTAGCGCAAAGTGTGCCGTAGTTCACGATCTCGACAGATAATTTATAA‑CAATTg > 2:79668/1‑139 (MQ=255) ccGGCTGGACGCTCTTTTTCCTGCGCCATCTCGTTCTCTTTTTTTCCCATCACTTCTTTCCCCATTTTGTCGCGTCCTGATGGTAGCGCAAAGTGTGCCGTAGTTCACGATCTCGACAGATAATTTATAA‑CAATTGAtt > 2:450950/1‑139 (MQ=255) tctcGTTCTCTTTTTTTCCCATCACTTCTTTCCCCATTTTGTCGCGTCCTGATGGTAGCGCAAAGTGTGCCGTAGTTCACGATCTCGACAGATAATTTATAA‑CAATTGATTTTTAt > 2:586793/1‑116 (MQ=255) tctcGTTCTCTTTTTTTCCCATCACTTCTTTCCCCATTTTGTCGCGTCCTGATGGTAGCGCAAAGTGTGCCGTAGTTCACGATCTCGACAGATAATTTATAA‑CAATTGATTTTTAt < 1:586793/116‑1 (MQ=255) cGTTCTCTTTTTTTCCCATCACTTCTTTCCCCATTTTGTCGCGTCCTGATGGTAGCGCAAAGTGTGCCGTAGTTCACGATCTCGACAGATAATTTATAAC‑AATTGATTTTTATGTCTTTTGAAATTCATCAATCctgtt > 2:768072‑M2/1‑134 (MQ=255) cGTTCTCTTTTTTTCCCATCACTTCTTTCCCCATTTTGTCGCGTCCTGATGGTAGCGCAAAGTGTGCCGTAGTTCACGATCTCGACAGATAATTTATAAC‑AATTGATTTTTATGTCTTTTGAAATTCATCAATCctgtt > 2:680367‑M2/1‑134 (MQ=255) tttttttCCCATCACTTCTTTCCCCATTTTGTCGCGTCCTGATGGTAGCGCAAAGTGTGCCGTAGTTCACGATCTCGACAGATAATTTATAA‑CAATTGAttttt > 1:193484/1‑104 (MQ=255) tttttttCCCATCACTTCTTTCCCCATTTTGTCGCGTCCTGATGGTAGCGCAAAGTGTGCCGTAGTTCACGATCTCGACAGATAATTTATAA‑CAATTGAttttt < 2:193484/104‑1 (MQ=255) ttttttCCCATCACTTCTTTCCCCATTTTGTCGCGTCCTGATGGTAGCGCAAAGTGTGCCGTAGTTCACGATCTCGACAGATAATTTATAAC‑AATTGATTTTTATGTCTTTTGAAATTCATCAATCctgttgacaatta > 1:774510‑M2/1‑126 (MQ=255) ttttCCCATCACTTCTTTCCCCATTTTGTCGCGTCCTGATGGTAGCGCAAAGTGTGCCGTAGTTCACGATCTCGACAGATAATTTATAAC‑AATTGATTTTTATGTCTTTTGAAATTCATCAATCctgttgacaattaat > 2:423383‑M2/1‑124 (MQ=255) tCCCATCACTTCTTTCCCCATTTTGTCGCGTCCTGATGGTAGCGCAAAGTGAGCCGTAGTTCACGATCTCGACAGATAATTTATAAC‑AATTTATTTTTATGTCTTTTGAAATTCATCAATCctgttgacaattaatcat > 1:532747‑M2/1‑121 (MQ=255) aTCACTTCTTTCCCCATTTTGTCGCGTCCTGATGGTAGCGCAAAGTGTGCCGTAGTTCACGATCTCGACAGATAATTTATAA‑CAATTGATTTTTATGTCTTTTGaaa > 1:124112/1‑107 (MQ=255) aTCACTTCTTTCCCCATTTTGTCGCGTCCTGATGGTAGCGCAAAGTGTGCCGTAGTTCACGATCTCGACAGATAATTTATAA‑CAATTGATTTTTATGTCTTTTGaaa < 2:124112/107‑1 (MQ=255) ttactTCTTTCCCCATTTTGTCGCGTCCTTATGGTAGCGCAAAGTGTTCCGTAGTTCACGATCTCGACAGATAATTTATAAC‑AATTGATTTTTATGTCTGTTGAAATTCATCAATCctgttgacaattaatcatcggca > 2:441600‑M2/3‑116 (MQ=255) tCACTTCTTTCCCCATTTTGTCGCGTCCTGATGGTAGCGCAAAGTGTGCCGTAGTTCACGATCTCGACAGATAATTTATAAC‑AATTGATTTTTATGTCTTTTGAAATTCATCAATCctgttgacaattaatcatcggca > 2:753889‑M2/1‑116 (MQ=255) cACTTCTTTCCCCATTTTGTCGCGTCCTGATGGTAGCGCAAAGTGTGCCGTAGTTCACGATCTCGACAGATAATTTATAAC‑AATTGATTTTTATGTCTTTTGAAATTCATCAATCctgttgacaattaatcatcggcat > 1:152640‑M2/1‑115 (MQ=255) cttcttTCCCCATTTTGTCGCGTCCTGATGGTAGCGCAAAGTGTGCCGTAGTTCACGATCTCGACATATAATTTATAAC‑AATTGATTTTTATGTCTTTTGAAATTCATCAATCctgttgacaattaatcatcggcatag > 1:328117‑M2/1‑113 (MQ=255) tCCCCATTTTGTCGCGTCCTGATGGTAGCGCAAAGTGTGCCGTAGTTCACGATCTCGACAGATAATTTATAA‑CAATTGATTTTTATg > 1:103824/1‑87 (MQ=255) tCCCCATTTTGTCGCGTCCTGATGGTAGCGCAAAGTGTGCCGTAGTTCACGATCTCGACAGATAATTTATAA‑CAATTGATTTTTATg < 2:103824/87‑1 (MQ=255) gTCGCGTCCTGATGGTAGCGCAAAGTGTGCCGTAGTTCACGATCTCGACAGATAATTTATAAC‑AATTGATTTTTATGTCTTTTGAAATTCATCAATCctgtt < 1:678536‑M2/102‑6 (MQ=255) gTCGCGTCCTGATGGTAGCGCAAAGTGTGCCGTAGTTCACGATCTCGACAGATAATTTATAAC‑AATTGATTTTTATGTCTTTTGAAATTCATCAATCctgtt > 2:678536‑M2/1‑97 (MQ=255) tCCTGATGGTAGCGCAAAGTGTGCCGTAGTTCACGATCTCGACAGATAATTTATAAC‑AATTGATTTTTATGTCTTTTGAAATTCATCAATCctgttgacaattaatcat > 2:232038‑M2/1‑91 (MQ=255) tCCTGATGGTAGCGCAAAGTGTGCCGTAGTTCACGATCTCGACAGATAATTTATAAC‑AATTGATTTTTATGTCTTTTGAAATTCATCAATCctgttgacaattaatcat < 1:232038‑M2/109‑19 (MQ=255) aTGGTAGCGCAAAGTGTGCCGTAGTTCACGATCTCGACAGATAATTTATAAC‑AATTGATTTTTATGTCTTTTGAAATTCATCAATCctgttgacaattaatcatcggcatag < 1:53146‑M2/112‑27 (MQ=255) aTGGTAGCGCAAAGTGTGCCGTAGTTCACGATCTCGACAGATAATTTATAAC‑AATTGATTTTTATGTCTTTTGAAATTCATCAATCctgttgacaattaatcatcggcatag > 2:244549‑M2/1‑86 (MQ=255) aTGGTAGCGCAAAGTGTGCCGTAGTTCACGATCTCGACAGATAATTTATAAC‑AATTGATTTTTATGTCTTTTGAAATTCATCAATCctgttgacaattaatcatcggcatag < 1:244549‑M2/112‑27 (MQ=255) aTGGTAGCGCAAAGTGTGCCGTAGTTCACGATCTCGACAGATAATTTATAAC‑AATTGATTTTTATGTCTTTTGAAATTCATCAATCctgttgacaattaatcatcggcatag > 2:53146‑M2/1‑86 (MQ=255) cgcAAAGTGTGCCGTAGTTCACGATCTCGACAGATAATTTATAAC‑AATTGATTTTTATGTCTTTTGAAATTCATCAATCctgttgacaattaatcatcggc < 1:580438‑M2/101‑23 (MQ=255) cgcAAAGTGTGCCGTAGTTCACGATCTCGACAGATAATTTATAAC‑AATTGATTTTTATGTCTTTTGAAATTCATCAATCctgttgacaattaatcatcggc > 2:580438‑M2/1‑79 (MQ=255) aGTGTGCCGTAGTTCACGATCTCGACAGATAATTTATAAC‑AATTGATTTTTATGTCTTTTGAAATTCATCAATCctgttgacaattaatcatcggcatagtata < 1:22380‑M2/104‑31 (MQ=255) aGTGTGCCGTAGTTCACGATCTCGACAGATAATTTATAAC‑AATTGATTTTTATGTCTTTTGAAATTCATCAATCctgttgacaattaatcatcggcatagtata > 2:22380‑M2/1‑74 (MQ=255) gtgtGCCGTAGTTCACGATCTCGACAGATAATTTATAA‑CAATTGATTTTTATGTCTTTTGAAATTCatcaatc < 2:548721/73‑1 (MQ=255) gtgtGCCGTAGTTCACGATCTCGACAGATAATTTATAA‑CAATTGATTTTTATGTCTTTTGAAATTCatcaatc > 1:548721/1‑73 (MQ=255) tgCCGTAGTTCACGATCTCGACAGATAATTTATAAC‑AATTGATTTTTATGTCTTTTGAAATTCATCAATCctgttgacaattaat < 1:835253‑M2/85‑16 (MQ=255) tgCCGTAGTTCACGATCTCGACAGATAATTTATAAC‑AATTGATTTTTATGTCTTTTGAAATTCATCAATCctgttgacaattaat > 2:835253‑M2/1‑70 (MQ=255) | CTTCCGGCTGGACGCTCTTTTTCCTGCGCCATCTCGTTCTCTTTTTTTCCCATCACTTCTTTCCCCATTTTGTCGCGTCCTGATGGTAGCGCAAAGTGTGCCGTAGTTCACGATCTCGACAGATAATTTATAACCAATTGATTTTTATGTCTTTTGAAATTCATCAATCAGATTGCCTTGTTAAAAAGTGATCGATATA > NC_000913/3742465‑3742663 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 16 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |