| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | NC_000913 | 3,742,598 | 0 | C | . | 44.4% | 119.8 / 102.0 | 63 | intergenic (‑66/‑135) | 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 (20/15); new base . (13/15); total (33/30) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 4.53e-01 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 | |||||||||||
| Rejected as consensus: Frequency below/above cutoff threshold. | |||||||||||
TGGCTTCCGGCTGGACGCTCTTTTTCCTGCGCCATCTCGTTCTCTTTTTTTCCCATCACTTCTTTCCCCATTTTGTCGCGTCCTGATGGTAGCGCAAAGTGTGCCGTAGTTCACGATCTCGACAGATAATTTATAACCAATTGATTTTTATGTCTTTTGAAATTCATCAATCAGATTGCCTTGTTAAAAAGTGATCGATATAT > NC_000913/3742462‑3742664 | tGGCTTCCGGCTGGACGCTCTTTTTCCTGCGCCATCTCGTTCTCTTTTTTTCCCATCACTTCTTTCCCCATTTTGTCGCGTCCTGATGGTAGCGCAAAGTGTGCCGTAGTTCACGATCTCGACAGATAATTTATAAcaa < 2:210565/139‑3 (MQ=255) gCTTCCGGCTGGACGCTCTTTTTCCTGCGCCATCTCGTTCTCTTTTTTTCCCATCACTTCTTTCCCCATTTTGTCGCGTCCTGATGGTAGCGCAAAGTGTGCCGTAGTTCACGATCTCGACAGATAATTTATAAcaatt < 1:1055544/139‑5 (MQ=255) gCTTCCGGCTGGACGCTCTTTTTCCTGCGCCATCTCGTTCTCTTTTTTTCCCATCACTTCTTTCCCCATTTTGTCGCGTCCTGATGGTAGCGCAAAGTGTGCCGTAGTTCACGATCTCGACAGATAATTTATAAcaatt < 2:1340998/139‑5 (MQ=255) ccGGCTGGACGCTCTTTTTCCTGCGCCATCTCGTTCTCTTTTTTTCCCATCACTTCTTTCCCCATTTTGTCGCGTCCTGATGGTAGCGCAAAGTGTGCCGTAGTTCACGATCTCGACAGATAATTTATAA‑CAATTGAtt > 2:242674/1‑139 (MQ=255) ggCTGGACGCTCTTTTTCCTGCGCCATCTCGTTCTCTTTTTTTCCCATCACTTCTTTCCCCATTTTGTCGCGTCCTGATGGTAGCGCAAAGTGTGCCGTAGTTCACGATCTCGACAGATAATTTATAA‑CAATTGAtttt > 2:1286896/1‑139 (MQ=255) tttCCTGCGCCATCTCGTTCTCTTTTTTTCCCATCACTTCTTTCCCCATTTTGTCGCGTCCTGATGGTAGCGCAAAGTGTGCCGTAGTTCACGATCTCGACAGATAATTTATAA‑CAATTGATTTTTATGTCTTTTGaaa < 1:89189/139‑1 (MQ=255) tttCCTGCGCCATCTCGTTCTCTTTTTTTCCCATCACTTCTTTCCCCATTTTGTCGCGTCCTGATGGTAGCGCAAAGTGTGCCGTAGTTCACGATCTCGACAGATAATTTATAA‑CAATTGATTTTTATGTCTTTTGaaa < 1:388759/139‑1 (MQ=255) tttCCTGCGCCATCTCGTTCTCTTTTTTTCCCATCACTTCTTTCCCCATTTTGTCGCGTCCTGATGGTAGCGCAAAGTGTGCCGTAGTTCACGATCTCGACAGATAATTTATAA‑CAATTGATTTTTATGTCTTTTGaaa < 1:1821848/139‑1 (MQ=255) tttCCTGCGCCATCTCGTTCTCTTTTTTTCCCATCACTTCTTTCCCCATTTTGTCGCGTCCTGATGGTAGCGCAAAGTGTGCCGTAGTTCACGATCTCGACAGATAATTTATAA‑CAATTGATTTTTATGTCTTTTGaaa < 2:282506/139‑1 (MQ=255) tCCTGCGCCATCTCGTTCTCTTTTTTTCCCATCACTTCTTTCCCCATTTTGTCGCGTCCTGATGGTAGCGCAAAGTGTGCCGTAGTTCACGATCTCGACAGATAATTTATAA‑CAATTGATTTTTATGTCTTTTGAAAtt < 2:540101/139‑1 (MQ=255) gcgcCATCTCGTTCTCTTTTTTTCCCATCACTTCTTTCCCCATTTTGTCGCGTCCTGATGGTAGCGCAAAGTGTGCCGTAGTTCACGATCTCGACAGATAATTTATAA‑CAATTGATTTTTATGTCTTTTGAAATtcatc > 2:1704347/1‑139 (MQ=255) cgcCATCTCGTTCTCTTTTTTTCCCATCACTTCTTTCCCCATTTTGTCGCGTCCTGATGGTAGCGCAAAGTGTGCCGTAGTTCACGATCTCGACAGATAATTTATAA‑CAATTGATTTTTATGTCTTTTGAAATtcatca < 1:242674/139‑1 (MQ=255) gcCATCTCGTTCTCTTTTTTTCCCATCACTTCTTTCCCCATTTTGTCGCGTCCTGATGGTAGCGCAAAGTGTGCCGTAGTTCACGATCTCGACAGATAATTTATAA‑CAATTGATTTTTATGTCTTTTGAAATTCatcaa > 1:717326/1‑139 (MQ=255) tctcGTTCTCTTTTTTTCCCATCACTTCTTTCCCCATTTTGTCGCGTCCTGATGGTAGCGCAAAGTGTGCCGTAGTTCACGATCTCGACAGATAATTTATAAC‑AATTGATTTTTATGTCTTTTGAAATTCATCAATCct < 2:391590‑M2/139‑3 (MQ=255) ttCTCTTTTTTTCCCATCACTTCTTTCCCCATTTTGTCGCGTCCTGATGGTAGCGCAAAGTGTGCCGTAGTTCACGATCTCGACAGATAATTTATAAC‑AATTGATTTTTATGTCTTTTGAAATTCATCAATCctgttga < 1:56350‑M2/139‑8 (MQ=255) ttCTCTTTTTTTCCCATCACTTCTTTCCCCATTTTGTCGCGTCCTGATGGTAGCGCAAAGTGTGCCGTAGTTCACGATCTCGACAGATAATTTATAAC‑AATTGATTTTTATGTCTTTTGAAATTCATCAATCctgttga < 2:1027905‑M2/139‑8 (MQ=255) ctctTTTTTTCCCATCACTTCTTTCCCCATTTTGTCGCGTCCTGATGGTAGCGCAAAGTGTGCCGTAGTTCACGATCTCGACAGATAATTTATAAC‑AATTGATTTTTATGTCTTTTGAAATTCATCAATCctgttgaca < 1:1390461‑M2/139‑10 (MQ=255) tttttttCCCATCACTTCTTTCCCCATTTTGTCGCGTCCTGATGGTAGCGCAAAGTGTGCCGTAGTTCACGATCTCGACAGATAATTTATAA‑CAATTGATTTTTAt > 2:1500824/1‑106 (MQ=255) tttttttCCCATCACTTCTTTCCCCATTTTGTCGCGTCCTGATGGTAGCGCAAAGTGTGCCGTAGTTCACGATCTCGACAGATAATTTATAA‑CAATTGATTTTTAt < 1:1500824/106‑1 (MQ=255) tttCCCATCACTTCTTTCCCCATTTTGTCGCGTCCTGATGGTAGCGCAAAGTGTGCCGTAGTTCACGATCTCGACAGATAATTTATAAC‑AATTGATTTTTATGTCTTTTGAAATTCATCAATCctgttgacaattaatc > 1:1729312‑M2/1‑123 (MQ=255) tCCCATCACTTCTTTCCCCATTTTGTCGCGTCCTGATGGTAGCGCAAAGTGTGCCGTAGTTCACGATCTCGACAGATAATTTATAAC‑AATTGATTTTTATGTCTTTTGAAATTCATCAATCctgttgacaattaatcat > 2:93472‑M2/1‑121 (MQ=255) tCCCATCACTTCTTTCCCCATTTTGTCGCGTCCTGATGGTAGCGCAAAGTGTGCCGTAGTTCACGATCTCGACAGATAATTTATAAC‑AATTGATTTTTATGTCTTTTGAAATTCATCAATCctgttgacaattaatcat > 1:391590‑M2/1‑121 (MQ=255) ccATCACTTCTTTCCCCATTTTGTCGCGTCCTGATGGTAGCGCAAAGTGTGCCGTAGTTCACGATCTCGACAGATAATTTATAAC‑AATTGATTTTTATGTCTTTTGAAATTCATCAATCctgttgacaattaatcatcg > 1:392326‑M2/1‑119 (MQ=255) aTCACTTCTTTCCCCATTTTGTCGCGTCCTGATGGTAGCGCAAAGTGTGCCGTAGTTCACGATCTCGACAGATAATTTATAA‑CAATTGATTTTTATGTCttt > 2:240056/1‑102 (MQ=255) aTCACTTCTTTCCCCATTTTGTCGCGTCCTGATGGTAGCGCAAAGTGTGCCGTAGTTCACGATCTCGACAGATAATTTATAA‑CAATTGATTTTTATGTCttt < 1:240056/102‑1 (MQ=255) tCACTTCTTTCCCCATTTTGTCGCGTCCTGATGGTAGCGCAAAGTGTGCCGTAGTTCACGATCTCGACAGATAATTTATAAC‑AATTGATTTTTATGTCTTTTGAAATTCATCAATCctgttgacaattaatcatcggca > 1:1704347‑M2/1‑116 (MQ=255) tCACTTCTTTCCCCATTTTGTCGCGTCCTGATGGTAGCGCAAAGTGTGCCGTAGTTCACGATCTCGACAGATAATTTATAAC‑AATTGATTTTTATGTCTTTTGAAATTCATCAATCctgttgacaattaatcatcggca > 1:791777‑M2/1‑116 (MQ=255) cACTTCTTTCCCCATTTTGTCGCGTCCTGATGGTAGCGCAAAGTGTGCCGTAGTTCACGATCTCGACAGATAATTTATAAC‑AATTGATTTTTATGTCTTTTGAAATTCATCAATCctgttgacaattaatcatcggcat > 1:1397213‑M2/1‑115 (MQ=255) cACTTCTTTCCCCATTTTGTCGCGTCCTGATGGTAGCGCAAAGTGTGCCGTAGTTCACGATCTCGACAGATAATTTATAAC‑AATTGATTTTTATGTCTTTTGAAATTCATCAATCctgttgacaattaatcatcggcat > 1:1286896‑M2/1‑115 (MQ=255) cACTTCTTTCCCCATTTTGTCGCGTCCTGATGGTAGCGCAAAGTGTGCCGTAGTTCACGATCTCGACAGATAATTTATAAC‑AATTGATTTTTATGTCTTTTGAAATTCATCAATCctgttgacaattaatcatcggcat > 2:56350‑M2/1‑115 (MQ=255) cttTCCCCATTTTGTCGCGTCCTGATGGTAGCGCAAAGTGTGCCGTAGTTCACGATCTCGACAGATAATTTATAAC‑AATTGATTTTTATGTCTTTTGAAATTCATCAATCctgttgacaattaatcatcggcatagtat < 1:1498205‑M2/139‑30 (MQ=255) ttCCCCATTTTGTCGCGTCCTGATGGTAGCGCAAAGTGTGCCGTAGTTCACGATCTCGACAGATAATTTATAAC‑AATTGATTTTTATGTCTTTTGAAATTCATCAATCctgttgacaattaatcatcggcatagtatat > 2:891854‑M2/1‑108 (MQ=255) aTTTTGTCGCGTCCTGATGGTAGCGCAAAGTGTGCCGTAGTTCACGATCTCGACAGATAATTTATAAC‑AATTGATTTTTATGTCTTTTGAAATTCATCAATCctgttgacaattaatcat > 2:530021‑M2/1‑102 (MQ=255) aTTTTGTCGCGTCCTGATGGTAGCGCAAAGTGTGCCGTAGTTCACGATCTCGACAGATAATTTATAAC‑AATTGATTTTTATGTCTTTTGAAATTCATCAATCctgttgacaattaatcat < 1:530021‑M2/120‑19 (MQ=255) aTTTTGTCGCGTCCTGATGGTAGCGCAAAGTGTGCCGTAGTTCACGATCTCGACAGATAATTTATAAC‑AATTGATTTTTATGTCTTTTGAAATTCATCAATCctgtt > 1:34517‑M2/1‑102 (MQ=255) aTTTTGTCGCGTCCTGATGGTAGCGCAAAGTGTGCCGTAGTTCACGATCTCGACAGATAATTTATAAC‑AATTGATTTTTATGTCTTTTGAAATTCATCAATCctgtt < 2:34517‑M2/107‑6 (MQ=255) aTTTTGTCGCGTCCTGATGGTAGCGCAAAGTGTGCCGTAGTTCACGATCTCGACAGATAATTTATAA‑CAATTGATTTTTAt > 1:373743/1‑81 (MQ=255) aTTTTGTCGCGTCCTGATGGTAGCGCAAAGTGTGCCGTAGTTCACGATCTCGACAGATAATTTATAA‑CAATTGATTTTTAt < 2:373743/81‑1 (MQ=255) aTTTTGTCGCGTCCTGATGGTAGCGCAAAGTGTGCCGTAGTTCACGATCTCGACAGATAATTTATAA‑CAATTGATTTTTATg < 1:1396155/82‑1 (MQ=255) aTTTTGTCGCGTCCTGATGGTAGCGCAAAGTGTGCCGTAGTTCACGATCTCGACAGATAATTTATAA‑CAATTGATTTTTATg > 2:1396155/1‑82 (MQ=255) tttGTCGCGTCCTGATGGTAGCGCAAAGTGTGCCGTAGTTCACGATCTCGACAGATAATTTATAAC‑AATTGATTTTTATGTCTTTTGAAATTCATCAATCctgttgacaattaatcat > 2:305302‑M2/1‑100 (MQ=255) tttGTCGCGTCCTGATGGTAGCGCAAAGTGTGCCGTAGTTCACGATCTCGACAGATAATTTATAAC‑AATTGATTTTTATGTCTTTTGAAATTCATCAATCctgttgacaattaatcat < 1:305302‑M2/118‑19 (MQ=255) tttGTCGCGTCCTGATGGTAGCGCAAAGTGTGCCGTAGTTCACGATCTCGACAGATAATTTATAAC‑AATTGATTTTTATGTCTTTTGAAATTCATCAATCctgttgacaattaat < 1:956084‑M2/115‑16 (MQ=255) tttGTCGCGTCCTGATGGTAGCGCAAAGTGTGCCGTAGTTCACGATCTCGACAGATAATTTATAAC‑AATTGATTTTTATGTCTTTTGAAATTCATCAATCctgttgacaattaat > 2:956084‑M2/1‑100 (MQ=255) tGTCGCGTCCTGATGGTAGCGCAAAGTGTGCCGTAGTTCACGATCTCGACAGATAATTTATAA‑CAATTGATTTTTATGTCTTTTGAAATtcat < 1:1583438/93‑1 (MQ=255) tGTCGCGTCCTGATGGTAGCGCAAAGTGTGCCGTAGTTCACGATCTCGACAGATAATTTATAA‑CAATTGATTTTTATGTCTTTTGAAATtcat > 2:1583438/1‑93 (MQ=255) cgTCCTGATGGTAGCGCAAAGTGTGCCGTAGTTCACGATCTCGACAGATAATTTATAA‑CAATTGATTTTTATGTCTTTTGAAATtcatca > 2:1267975/1‑90 (MQ=255) cgTCCTGATGGTAGCGCAAAGTGTGCCGTAGTTCACGATCTCGACAGATAATTTATAA‑CAATTGATTTTTATGTCTTTTGAAATtcatca < 1:1267975/90‑1 (MQ=255) tCCTGATGGTAGCGCAAAGTGTGCCGTAGTTCACGATCTCGACAGATAATTTATAAC‑AATTGATTTTTATGTCTTTTGAAATTCATCAATCctgttgacaattaatcat < 2:1474785‑M2/109‑19 (MQ=255) tCCTGATGGTAGCGCAAAGTGTGCCGTAGTTCACGATCTCGACAGATAATTTATAAC‑AATTGATTTTTATGTCTTTTGAAATTCATCAATCctgttgacaattaatcat > 1:1474785‑M2/1‑91 (MQ=255) aGCGCAAAGTGTGCCGTAGTTCACGATCTCGACAGATAATTTATAA‑CAATTGATTTTTATGTCTTTTGaaa < 2:586020/71‑1 (MQ=255) aGCGCAAAGTGTGCCGTAGTTCACGATCTCGACAGATAATTTATAA‑CAATTGATTTTTATGTCTTTTGaaa > 1:586020/1‑71 (MQ=255) gcgcAAAGTGTGCCGTAGTTCACGATCTCGACAGATAATTTATAA‑CAATTGATTTTTAt > 2:926984/1‑59 (MQ=255) gcgcAAAGTGTGCCGTAGTTCACGATCTCGACAGATAATTTATAA‑CAATTGATTTTTAt < 1:926984/59‑1 (MQ=255) tgCCGTAGTTCACGATCTCGACAGATAATTTATAAca > 2:1542922/1‑36 (MQ=255) tgCCGTAGTTCACGATCTCGACAGATAATTTATAAca < 1:1542922/37‑2 (MQ=255) ccGTAGTTCACGATCTCGACAGATAATTTATAAC‑AATTGATTTTTATGTCTTTTGAAATTCATCAATCctgttg > 2:152184‑M2/1‑68 (MQ=255) ccGTAGTTCACGATCTCGACAGATAATTTATAAC‑AATTGATTTTTATGTCTTTTGAAATTCATCAATCctgttg < 1:152184‑M2/74‑7 (MQ=255) ccGTAGTTCACGATCTCGACAGATAATTTATAAC‑AATTGATTTTTATGTCTTTTGAAATTCATCAATCctgtt < 2:1242417‑M2/73‑6 (MQ=255) ccGTAGTTCACGATCTCGACAGATAATTTATAAC‑AATTGATTTTTATGTCTTTTGAAATTCATCAATCctgtt > 1:1242417‑M2/1‑68 (MQ=255) ccGTAGTTCACGATCTCGACAGATAATTTATAAC‑AATTGATTTTTATGTCTTTTGAAATTCATCAATCctg > 2:1612171‑M2/1‑68 (MQ=255) ccGTAGTTCACGATCTCGACAGATAATTTATAAC‑AATTGATTTTTATGTCTTTTGAAATTCATCAATCctg < 1:1612171‑M2/71‑4 (MQ=255) ttCACGATCTCGACAGATAATTTATAAC‑AATTGATTTTTATGTCTTTTGAAATTCATCAATCctgttgacaattaatcat < 2:1807845‑M2/80‑19 (MQ=255) ttCACGATCTCGACAGATAATTTATAAC‑AATTGATTTTTATGTCTTTTGAAATTCATCAATCctgttgacaattaatcat > 1:1807845‑M2/1‑62 (MQ=255) ttCACGATCTCGACAGATAATTTATAA‑CAATTGATTTTTATGTCTTTTGAAATTCatcaatc < 2:1453897/62‑1 (MQ=255) ttCACGATCTCGACAGATAATTTATAA‑CAATTGATTTTTATGTCTTTTGAAATTCatcaatc > 1:1453897/1‑62 (MQ=255) aaTTTATAAC‑AATTGATTTTTATGTCTTTTGAAATTCATCAATCctgttgacaattaatcatcg < 1:608865‑M2/64‑21 (MQ=255) aaTTTATAAC‑AATTGATTTTTATGTCTTTTGAAATTCATCAATCctgttgacaattaatcatcg > 2:608865‑M2/1‑44 (MQ=255) | TGGCTTCCGGCTGGACGCTCTTTTTCCTGCGCCATCTCGTTCTCTTTTTTTCCCATCACTTCTTTCCCCATTTTGTCGCGTCCTGATGGTAGCGCAAAGTGTGCCGTAGTTCACGATCTCGACAGATAATTTATAACCAATTGATTTTTATGTCTTTTGAAATTCATCAATCAGATTGCCTTGTTAAAAAGTGATCGATATAT > NC_000913/3742462‑3742664 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 30 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |