| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | NC_000913 | 2,040,075 | 0 | C | G | 52.9% | 19.5 / 43.6 | 35 | L200V (CTG→GTG) | yedY | membrane‑anchored, periplasmic TMAO, DMSO reductase |
| Reads supporting (aligned to +/- strand): ref base C (5/11); new base G (18/0); total (23/12) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 1.53e-05 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.39e-05 | |||||||||||
GTGGTGCCGTGGATTGGTTTTCCGCTGCACAAATTGCTGGCGCTTGCCGAACCTACCAGCAATGCGAAGTATGTCGCTTTCGAAACAATTTATGCACCGGAACAGATGCCAGGCCAGCAGGACCGCTTTATCGGCGGTGGGCTGAAATATCCTTATGTCGAAGG‑ATTGCGTCTCGACGAAGCAATGCATCCGCTCACACTGATGACCGTAGGTGTTTATGGCAAGGCGTTACCGCCACAAAATGGCGCGCCGGTGCGACTGATTGTGCCGTGGAAATATGGCTTTAA > NC_000913/2039934‑2040220 | gtggtgCCGTGGATTGGATTTCCGCTGCACAAATTGCTGGCGCTTGCCGAACCTACCAGCAATGCGAAGTATGTCGCTTTCGAAACAATTTATGCGCCGGAACAGATGCCAGGCCAGCAGGACCGCTTTATCGGCGGTGGGCTGAAatag < 2:144703/150‑2 (MQ=255) ccGTGGATTGGTTTTCCGCTGCACAAATTGCTGGCGCTTGCCGAACCTACCAGCAATGCGAAGTATGTCGCTTTCGAAACAATTTATGCACCGGAACAGATGCCAGGCCAGCAGGACCGCTTTATCGGCGGTGGGGTGAAata > 1:2686/1‑143 (MQ=255) cGTGGATTGGTTTTCCGCTGCACAAATTGCTGGCGCTTGCCGAACCTACCAGCAATGCGAAGTATGTCGCTTTCGAAACAATTTATGCACCGGAACAGATGCCAGGCCAGCAGGACCGCTTTATCGGCGGCGGGGTGAAATATCCTTAtat > 2:470215/1‑149 (MQ=255) gATTGGTTTTCCGCTGCACAAATTGCTGGCGCTTGCCGAACCTACCAGCAATGCGAAGTATGTCGCTTTCGAAACAATTTATGCACCGGAACAGATGCCAGGCCAGCAGGACCGCTTTATCGGCGGTGGGCTGAAATATCCTTATGTCGaa < 1:269912/151‑1 (MQ=255) tGCACAAATTGCTGGCGCTTGCCGAACCTACCAGCAATGCGAAGTATGTCGCTTTCGAAACAATTTATGCACCGGAACAGATGCCAGGCCAGCAGGACCGCTTTATCGGCGGTGGGGTGGAATATCCTTTTTTCGAAAG‑ATTTCg > 2:562930/1‑145 (MQ=255) gCACAAATTGCTGGCGCTTGCCGAACCTACCAGCAATGCGAAGTATGTCGCTTTCGAAACAATTTATGCACCGGAACAGATGCCAGGGCAGCAGGACCGCTTTATCGGCTGTGGGCTGAAATATCCTTATGTCGAAGG‑ATTGCGTCTcgc < 2:560183/150‑2 (MQ=255) acaAATTGCTGGCGCTTGCCGAACCTACCAGCAATGCGAAGTATGTCGCTTTCGAAACAATTTATGCACCGGAACAGATGCCAGGCCAGCAGGACCGCTTTATCGGCGGTGGGCTGAAATATCCTTATGTCGAAGG‑ATTGCGTCTcgacga < 2:604701/151‑1 (MQ=255) aTTGCTGGCGCTTGCCGAACCTACCAGCAATGCGAAGTATGTCGCTTTCGAAACAATTTATGCACCGGAACAGATGCCAGGCCAGCAGGACCGCTTTATCGGCGGTGGGGTGGAAAATCCTTTTGTCGAAGG‑ATTTCGGCTCGACGaa > 1:217111/1‑148 (MQ=255) aTTGCTGGCGCTTGCCGAACCTACCAGCAATGCGAAGTATGTCGCTTTCGAAACAATTTATGCACCGGAACAGATGCCAGGCCAGCAGGACCGCTTTATCGGCGGTGGGGGGGAAAATTCCTTTTTCGGAAG‑ATTTCGTCTcgccgga > 1:700347/1‑146 (MQ=255) cgcTTGCCGAACCTACCAGCAATGCGAAGTATGTCGCTTTCGAAACAATTTATGCACCGGAACAGATGCCAGGCCAGCAGGACCGCTTTATCGGCGGTGGGCTGAAATATCCTTATGTCGAAGG‑ATTGCGTCTCGACGAAGCAATGCATcc < 1:482896/151‑1 (MQ=255) gcTTGCCGAACCTACCAGCAATGCGAAGTATGTCGCTTTCGAAACAATTTATGCACCGGAACAGATGCCAGGCCAGCAGGACCGCTTTATCGGCGGCGGGGTGAAATATCCTTATTTCCAAAG‑GATGCTGCTCGACGCAACAAAGCAAcc > 2:63246/1‑150 (MQ=255) ccGAACCTACCAGCAATGCGAAGTATGTCGCTTTCGAAACAATTTATGCACCGGAACAGATGCCAGGCCAGCAGGACCGCTTTATCGGCGGTGGGCTGAAATATCCTTATGTCGAAGG‑ATTGCGTCTCGACGAAGCAATGCATCCGCTcac < 2:721272/151‑1 (MQ=255) ccAGCAATGCGAAGTATGTCGCTTTCGAAACAATTTATGCACCGGAACAGATGCCAGGCCAGCAGGACCGCTTTATCGGCGGTGGGCTGAAATATCCTTATGTCGAAGG‑ATTGCGTCTCGACGAAGCAATGCATCCGCTCACACTGATGAc < 1:155703/151‑1 (MQ=255) cAATGCGAAGTATGTCGCTTTCGAAACAATTTATGCACCGGAACAGATGCCAGGCCAGCAGGACCGCTTTATCGGCGGTGGGGTGAAATATTCTTATTTCGGAGG‑GTTTTGTTTCGACGAAGCAAAGCATCCCCTCCCCCTGATTAcc > 1:430433/1‑148 (MQ=255) cAATGCGAAGTATGTCGCTTTCGAAACAATTTATGCACCGGAACAGATGCCAGGCCAGCAGGACCGCTTTATCGGCGGTGGGGGGAAATATACCTTTGTCGGAGG‑ATTTCGGCTCGACGGAGGAAAGCATCCGCCCACACtgctgaac > 2:358286/1‑146 (MQ=255) gCGAAGTATGTCGCTTTCGAAACAATTTATGCACCGGAACAGATGCCAGGCCAGCAGGACCGCTTTATCGGCGGTGGGGGGGAAAATTCTTTTGTTGAAGG‑AGTGGGTTTCGACGGAGCAAAGCATCCCCCCCCCCCCATGAACGTAgg > 2:585197/1‑149 (MQ=255) cGAAGTATGTCGCTTTCGAAACAATTTATGCACCGGAACAGATGCCAGGCCAGCAGGACCGCTTTATCGGCGGTGGGGTGGAATATCCTTTTGTCGAAGG‑AATGCGTCTCCACGGAGCAAAGCATCCCCCCACACTGATGACCGTAgggg > 1:320599/1‑148 (MQ=255) tGTCGCTTTCGAAACAATTTATGCACCGGAACAGATGCCAGGCCAGCAGGACCGCTTTATCGGCGGTGGGGTGAAATATCCTTTTTTCGAAGG‑ATTGCGTCTCGACGAAGCAAAGCATCCCCCCACACCGATGGCCGGAgggggtttttg > 2:708193/1‑141 (MQ=255) gCTTTCGAAACAATTTATGCACCGGAACAGATGCCAGGCCAGCAGGACCGCTTTATCGGCGGTGGGCTGAAATATCCTTATGTCGTAGG‑GATTTGTCTCGAAGAAGCAAAACATTCGCCCCCACTGATTACCGGAGGTGTTTATGGCa > 1:383692/1‑148 (MQ=255) aaCAATTTATGCACCGGAACAGATGCCAGGCCAGCAGGACCGCTTTATCGGCGGTGGGCTGAAATATCCTTATGTCGAAGG‑ATTGCGTCTCGACGAAGCAATGCATCCGCTCACACTGATGACCGTAGGTGTTTATGGCAAGGCGTTACCg < 1:479995/151‑1 (MQ=255) tGCACCGGAACAGATGCCAGGCCAGCAGGACCGCTTTATCGGCGGTGGGGTGGAATATTCTTTTAGCGAAGG‑AAATCGGCTTGGCGGAGGAAAGCACCCCCCCACACTGATGGCCGGAGGGGTTTTTGGCAAGGGGGTAACACCACcaaa > 2:274780/1‑150 (MQ=255) gCACCGGAACAGATGCCAGGCCAGCAGGACCGCTTTATCGGCGGTGGGCTGGAAAAACCTTATGTTGAAGG‑ATTGCGTCTCGGCGAAGCAAAGCAACCCCTCACACTTATGACCGTAGGGGTTTTTGGCAAGGCGTTAACCCCCCaaaa > 1:390376/1‑149 (MQ=255) cACCGGAACAGATGCCAGGCCAGCAGGACCGCTTTATCGGCGGTGGGGTGAAATATCCTTATGTCGAAGG‑ATTTCGTCTCCACGAAGCAATGCATCCCCTCACACTGATGACCGTAGGTGTTTTTGGCAAAGCGGTACCGCCACaaaa > 1:349464/1‑148 (MQ=255) cGGAACAGATGCCAGGCCAGCAGGACCGCTTTATCGGCGGTGGGGTGAAATATCCTTATTTTGAAGG‑ATTTCGTCTCGACGGAGCAATGCATCCGCTCCCACTGATGACCGTAGGTGTTTTTGGCAAGGCGTTTCCGCCCCAAAAAGgcg > 1:624131/1‑150 (MQ=255) gccagggccagAAGGCCCGTTTAAGCGGGGGGGGGCGAAAATTTCTTTATG‑CGAAGGAATTGCCCCTGGAGAAAGCAATCAACCCCCTCCCACTGACGACCGGGGTTTTTTATGGAAAGTCTTACCCGCCACAAAATGGCGCGCCGGTGc < 2:228266/145‑1 (MQ=255) caggccagCAGGACCGCTTTATCGGCGGTGGGGTGAAAAATCCTTATGTCGGAGG‑ATTGCGTCTCGACGAAGCAAAGCATCCGCTCACACTGATGGCCGGAGGGGGTTTTTGCAAGGGGGTACCGCCCCAAAAAAGCGCGCCggcg > 1:419598/1‑144 (MQ=255) cagcccagGGGGACCGCTTTGGCGGGGTTGGGAGAAAAATTCTTTGGGTCGAAGG‑ATTGCGTCTGGACGAAGCAATGCACCCGCCCACACTAATGACCGGAGGTTTTTATGGCAAGGGGTTACCGCCACAAAATGGCGCGCCGGTGc < 2:641371/147‑1 (MQ=255) cagcagGACCGCTTTATCGGCGGTGGGGTGGAAAATTCTTATTTCGGAGG‑ATTGCGTCTCCACGAAGCAAAGCATCCCCTCACACTGATGACCGTAGGGGTTTATTGCAAGGCGGTACCGCCACAAAAAAGCGCGCCGGGGCGGCTGAtt > 1:328040/1‑150 (MQ=255) agcagGACCGCTTTATCGGCGGTGGGGGGAAAAATCCTTTTTTCGAAGG‑AGTGCGGCTCGGCGGAGGAATGCATCCGCCCACACTGATGACCGGAGGGGTTTTTTGGAAAGCGGTAACCCCACAAAAAGGCGCGCCGGGGGGGCTGAtt > 2:167100/1‑149 (MQ=255) ggACCGCTTTATGGGCGGTGGGCTGAAATATCCTTATGGCGAAGG‑TTTGCGTCTCGACGAAGCAATGCATCCGCTCACACTGATGCCCGTAGGTGTTTAGGGCAGGGCGTTACCGCCACAAAATGGCGCGCCGGTGCGACTGAtt < 2:687164/145‑1 (MQ=255) gACCGCTTTATCGGCGGTGGGGTGAAAAATTCTTATTTCGCAGG‑AATGCGTCTCGCCGGAGCAATTCATTCGCTCACACCGATTACCGGAGGGGGTTTTGGCCAAGCGTTAACGCCCCAAAAAAGCGCGCCCGGGGGGATGATtgttgc > 2:535797/1‑146 (MQ=255) ttAGCGGCGGGGGGCTAAAATACCCTTATGTCGAAGG‑ATTGCGTCTCGACGAAGCAATGCATCCGCTCCCACTGAGGACCGTGGGTTTTTATGGCAAGGCGTAACCGCCACAAAATGGCGCGCCGGTGCGACTGATTGTGCCGTGGAAat < 2:519019/150‑1 (MQ=255) tGGGCTGAAATATCCTTATGTCGAAGG‑ATTGCGTCTCGACGAAGCAATGCATCCGCTCACACTGATGACCGTAGGTGTTTATGGCAAGGCGTTACCGCCACAAAATGGCGCGCCGGTGCGACTGATTGTGCCGTGGAAATATGGCTTTaa > 2:83803/1‑150 (MQ=255) gggCTGAAATATCCTTATGTCGAAGG‑ATTGCGTCTCGACGAAGCAATGCATCCGCTCACACTGATGACCGTAGGTGTTTATGGCAAGGCGTTACCGCCACAAAATGGCGCGCCGGTGCGACTGATTGTGCCGTGGAAATATGGCTTTaa > 1:356939/1‑149 (MQ=255) gggCTGAAATATCCTTATGTCGAAGG‑ATTGCGTCTCGACGAAGCAATGCATCCGCTCACACTGATGACCGTAGGTGTTTATGGCAAGGCGTTACCGCCACAAAATGGCGCGCCGGTGCGACTGATTGTGCCGTGGAAATATGGCTTTaa > 2:106035/1‑149 (MQ=255) | GTGGTGCCGTGGATTGGTTTTCCGCTGCACAAATTGCTGGCGCTTGCCGAACCTACCAGCAATGCGAAGTATGTCGCTTTCGAAACAATTTATGCACCGGAACAGATGCCAGGCCAGCAGGACCGCTTTATCGGCGGTGGGCTGAAATATCCTTATGTCGAAGG‑ATTGCGTCTCGACGAAGCAATGCATCCGCTCACACTGATGACCGTAGGTGTTTATGGCAAGGCGTTACCGCCACAAAATGGCGCGCCGGTGCGACTGATTGTGCCGTGGAAATATGGCTTTAA > NC_000913/2039934‑2040220 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 24 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |