Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_000913 | 4,429,984 | 0 | T | G | 29.7% | 62.9 / 10.2 | 37 | L41V (TTG→GTG) | cycA | serine/alanine/glycine/cycloserine:H(+)symporter |
Reads supporting (aligned to +/- strand): ref base T (4/22); new base G (11/0); total (15/22) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 1.60e-06 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
GGTACAGGAAGAAAAAAACATGGTAGATCAGGTAAAAGTCGTTGCCGATGATCAGGCTCCGGCTGAACAGTCGCTACGGCGCAATCTCACAAACCGACATATTCAGCTTATTGCCATTGGCGGTGCCATTGGTACGGGGTTGTTTATGGGGTCTGGCAAAACGATTAG‑CCTTGCCGGGCCGTCGATCATTTTCGTTTATATGATCATTGGTTTTATGCTCTTTTTCGTGATGCGGGCAATGGGGGAATTGCTGCTTTCGAATCTGGAATACAA > NC_000913/4429845‑4430117 | ggTACAGGAAGAAAAAAACATGGTAGATCAGGTAAAAGTCGTTGCCGATGATCAGGCTCCGGCTGAACAGTCGCTACGGCGCAATCTCACAAACCGACATATTCAGCTTATTGCCATTGGCGGTGCCATTGGTACGGGgt < 2:28535/140‑1 (MQ=255) ggTACAGGAAGAAAAAAACATGGTAGATCAGGTAAAAGTCGTTGCCGATGATCAGGCTCCGGCTGAACAGTCGCTACGGCGCAATCTCACAAACCGACATATTCAGCTTATTGCCATTGGCGGTGCCATTGGTACGGGgt < 2:159892/140‑1 (MQ=255) ggTACAGGAAGAAAAAAACATGGTAGATCAGGTAAAAGTCGTTGCCGATGATCAGGCTCCGGCTGAACAGTCGCTACGGCGCAATCTCACAAACCGACATATTCAGCTTATTGCCATTGGCGGTGCCATTGGTACGGGgt < 1:251196/140‑1 (MQ=255) tACAGGAAGAAAAAAACATGGTAGATCAGGTAAAAGTCGTTGCCGATGATCAGGCTCCGGCTGAACAGTCGCTACGGCGCAATCTCACAAACCGACATATTCAGCTTATTGCCATTGGCGGTGCCATTGGTACGGGgttg > 1:684027/1‑140 (MQ=255) tACAGGAAGAAAAAAACATGGTAGATCAGGTAAAAGTCGTTGCCGATGATCAGGCTCCGGCTGAACAGTCGCTACGGCGCAATCTCACAAACCGACATATTCAGCTTATTGCCATTGGCGGTGCCATTGGTACGGGgttg > 2:680392/1‑140 (MQ=255) aGGAAGAAAAAAACATGGTAGATCAGGTAAAAGTCGTTGCCGATGATCAGGCTCCGGCTGAACAGTCGCTACGGCGCAATCTCACAAACCGACATATTCAGCTTATTGCCATTGGCGGTGCCATTGGTACGGGGGTGttt > 1:294664/1‑140 (MQ=255) ggAAGAAAAAAACATGGTAGATCAGGTAAAAGTCGTTGCCGATGATCAGGCTCCGGCTGAACAGTCGCTACGGCGCAATCTCACAAACCGACATATTCAGCTTATTGCCATTGGCGGTGCCATTGGTACGGGGTTGTTTa < 1:680392/140‑1 (MQ=255) gaagaaAAAAACATGGTAGATCAGGTAAAAGTCGTTGCCGATGATCAGGCTCCGGCTGAACAGTCGCTACGGCGCAATCTCACAAACCGACATATTCAGCTTATTGCCATTGGCGGTGCCATTGGTACGGGGTTGTTTAt < 1:683569/140‑1 (MQ=255) gaagaaAAAAACATGGTAGATCAGGTAAAAGTCGTTGCCGATGATCAGGCTCCGGCTGAACAGTCGCTACGGCGCAATCTCACAAACCGACATATTCAGCTTATTGCCATTGGCGGTGCCATTGGTACGGGGTTGTTTAt > 1:66462/1‑140 (MQ=255) gaagaaAAAAACATGGTAGATCAGGTAAAAGTCGTTGCCGATGATCAGGCTCCGGCTGAACAGTCGCTACGGCGCAATCTCACAAACCGACATATTCAGCTTATTGCCATTGGCGGTGCCATTGGTACGGGGGTGTTTAt > 1:119901/1‑140 (MQ=255) aagaaAAAAACATGGTAGATCAGGTAAAAGTCGTTGCCGATGATCAGGCTCCGGCTGAACAGTCGCTACGGCGCAATCTCACAAACCGACATATTCAGCTTATTGCCATTGGCGGTGCCATTGGTACGGGGTTGTTTATg < 2:660619/140‑1 (MQ=255) agaaAAAAACATGGTAGATCAGGTAAAAGTCGTTGCCGATGATCAGGCTCCGGCTGAACAGTCGCTACGGCGCAATCTCACAAACCGACATATTCAGCTTATTGCCATTGGCGGTGCCATTGGTACGGGGTTGTTTATgg > 2:157754/1‑140 (MQ=255) gaaAAAAACATGGTAGATCAGGTAAAAGTCGTTGCCGATGATCAGGCTCCGGCTGAACAGTCGCTACGGCGCAATCTCACAAACCGACATATTCAGCTTATTGCCATTGGCGGTGCCATTGGTACGGGGTTGTTTATggg < 2:684027/140‑1 (MQ=255) aaaaaaaCATGGTAGATCAGGTAAAAGTCGTTGCCGATGATCAGGCTCCGGCTGAACAGTCGCTACGGCGCAATCTCACAAACCGACATATTCAGCTTATTGCCATTGGCGGTGCCATTGGTACGGGGGTGTTTATgggg > 1:176795/1‑140 (MQ=255) aaaaaaaCATGGTAGATCAGGTAAAAGTCGTTGCCGATGATCAGGCTCCGGCTGAACAGTCGCTACGGCGCAATCTCACAAACCGACATATTCAGCTTATTGCCATTGGCGGCGCCACTTGTACGGGGGTGTTTCTgtgg > 2:395406/1‑140 (MQ=255) aTGGTAGATCAGGTAAAAGTCGTTGCCGATGTTCAGGCTCCGGCTGAACAGTCGCTACGGCGCAATCTCACAAACCGACATATTCAGCTTATTGCCATTGGCGGTGCCATTGGTACGGGGTTGTTTATGGGGTCTGGCaa < 1:55449/140‑1 (MQ=255) ggTAGATCAGGTAAAAGTCGTTGCCGATGATCAGGCTCCGGCTGAACAGTCGCTACGGCGCAATCTCACAAACCGACATATTCAGCTTATTGCCATTGGCGGTGCCATTGGTACGGGGGTGTTTATGGGGTCTGGCaaaa > 2:461791/1‑140 (MQ=255) aGATCAGGTAAAAGTCGTTGCCGATGATCAGGCTCCGGCTGAACAGTCGCTACGGCGCAATCTCACAAACCGACATATTCAGCTTATTGCCATTGGCGGTGCCATTGGTACGGGGTTGTTTATGGGGTCTGGCAAAACGa < 1:270320/140‑1 (MQ=255) aGGTAAAAGTCGTTGCCGATGATCAGGCTCCGGCTGAACAGTCGCTACGGCGCAATCTCACAAACCGACATATTCAGCTTATTGCCATTGGCGGTGCCATTTGGTTGGGGGTGTTTTTGGGGGCTGGCCAAAAGaattac > 2:185108/1‑135 (MQ=255) gTAAAAGTCGTTGCCGATGATCAGGCTCCGGCTGAACAGTCGCTACGGCGCAATCTCACAAACCGACATATTCAGCTTATTGCCATTGGCGGTGCCCTTGGTACGGGGGTGTTTATTGGGGCTGGCAAAACGATTAA‑CCt > 2:541128/1‑140 (MQ=255) tCGTTGCCGATGATCAGGCTCCGGCTGAACAGTCGCTACGGCGCAATCTCACAAACCGACATATTCAGCTTATTGCCATTGGCGGTGCCATTGGTACGGGGTTGTTTATGGGGTCTGGCAAAACGATTAG‑CCTTGCCggg < 1:514944/140‑1 (MQ=255) ttGCCGATGATCAGGCTCCGGCTGAACAGTCGCTACGGCGCAATCTCACAAACCGACATATTCAGCTTATTGCCATTGGCGGTGCCATTGGTACGGGGTTGTTTATGGGGTCTGGCAAAACGATTAG‑CCTTGCCGGGCCg < 1:734188/140‑1 (MQ=255) tgatCAGGCTCCGGCTGAACAGTCGCTACGGCGCAATCTCACAAACCGACATATTCAGCTTATTGCCATTGGCGGTGCCATTTGTACGGGGGTGTTTTTGGGGGCTGGCAAAAAGATTTG‑CCTTGCCCGGGCGGCGGTcc > 2:399065/1‑139 (MQ=255) gatCAGGCTCCGGCTGAACAGTCGCTACGGCGCAATCTCACAAACCGACATATTCAGCTTATTGCCATTGGCGGTGCCATTGGTACGGGGTTGTTTATGGGGTCTGGCAAAACGATTAG‑CCTTGCCGGGCCGTCGATCAt < 2:338456/140‑1 (MQ=255) atCAGGCTCCGGCTGAACAGTCGCTACGGCGCAATCTCACAAACCGACATATTCAGCTTATTGCCATTGGCGGTGCCATTGGTACGGGGTTGTTTATGGGGTCTGGCAAAACGATTAG‑CCTTGCCGGGCCGTCGATCAtt < 2:176795/140‑1 (MQ=255) tCAGGCTCCGGCTGAACAGTCGCTACGGCGCAATCTCACAAACCGACATATTCAGCTTATTGCCATTGGCGGTGCCATTGGTACGGGGTTGTTTATGGGGTCTGGCAAAACGATTAG‑CCTTGCCGGGCCGTCGATCAttt < 1:541128/140‑1 (MQ=255) ggCTCCGGCTGAACAGTCGCTACGGCGCAATCTCACAAACCGACATATTCAGCTTATTGCCATTGGCGGTGCCATTGGTACGGGGTTGTTTATGGGGTCTGGCAAAACGATTAG‑CCTTGCCGGGCCGTCGATCATTTTCg < 2:66462/140‑1 (MQ=255) ggCTCCGGCTGAACAGTCGCTACGGCGCAATCTCACAAACCGACATATTCAGCTTATTGCCATTGGCGGTGCCATTGGTACGGGGTTGTTTATGGGGTCTGGCAAAACGATTAG‑CCTTGCCGGGCCGTCGATCATTTTCg < 2:663711/140‑1 (MQ=255) gCTCCGGCTGAACAGTCGCTACGGCGCAATCTCACAAACCGACATATTCAGCTTATTGCCATTGGCGGTGCCATTGGTACGGGGTTGTTTATGGGGTCTGGCAAAACGATTAG‑CCTTGCCGGGCCGTCGATCATTTTCGt < 1:157754/140‑1 (MQ=255) cGGCTGAACAGTCGCTACGGCGCAATCTCACAAACCGACATATTCAGCTTATTGCCATTGGCGGTGCCATTGGTACGGGGTTGTTTATGGGGTCTGGCAAAACGATTAG‑CCTTGCCGGGCCGTCGATCATTTTCGTTtat < 1:495360/140‑1 (MQ=255) cGGCTGAACAGTCGCTACGGCGCAATCTCACAAACCGACATATTCAGCTTATTGCCATTGGCGGTGCCATTGGTACGGGGGTGTTTATGGGGGCTGGCAAAACGATTAG‑CCTTGCCGGGCCGTCGATCATTTTCGTTtat > 2:767771/1‑140 (MQ=255) gCTGAACAGTCGCTACGGCGCAATCTCACAAACCGACATATTCAGCTTATTGCCATTGGCGGTGCCATTGGTACGGGGTTGTTTATGGGGTCTGGCAAAACGATTAG‑CCTTGCCGGGCCGTCGATCATTTTCGTTTatat < 1:780385/140‑1 (MQ=255) gCTGAACAGTCGCTACGGCGCAATCTCACAAACCGACATATTCAGCTTATTGCCATTGGCGGTGCCATTGGTACGGGGGTGTTTATGGGGTCTGGCAAAACGATTAG‑CCTTGCCGGGCCGTCGATCATTTTCGTTTatat > 1:700236/1‑140 (MQ=255) gCTGAACAGTCGCTACGGCGCAATCTCACAAACCGACATATTCAGCTTATTGCCATTGGCGGTGCCATTGGTACGGGGGTGTTTATGGGGTCTGGCAAAACGATTAG‑CCTTGCCGGGCCGTCGATCATTTTCGTTTatat > 2:42302/1‑140 (MQ=255) aaCAGTCGCTACGGCGCAATCTCACAAACCGACATATTCAGCTTATTGCCATTGGCGGTGCCATTGGTACGGGGTTGTTTATGGGGTCTGGCAAAACGATTAG‑CCTTGCCGGGCCGTCGATCATTTTCGTTTATATGATc > 1:220985/1‑140 (MQ=255) gTCGCTACGGCGCAATCTCACAAACCGACATATTCAGCTTATTGCCATTGGCGGTGGCCATGGGAAGGGGGTGTTTTTTGGGGCTTGCAAAAAGAGTAGCCCTTGCGGGGCGGT‑GATCCTTTTTGTTTTTATTATCAttt > 2:279390/1‑139 (MQ=255) tACGGCGCAATCTCACAAACCGACATATTCAGCTTATTGCCATTGGCGGTGCCATTGGTACGGGGTTGTTTATGGGGTCTGGCAAAACGATTAG‑CCTTGCCGGGCCGTCGATCATTTTCGTTTATATGATCATTGGtttt < 1:395406/140‑1 (MQ=255) ggCGCAATCTCACAAACCGACATATTCAGCTTATTGCCATTGGCGGTGCCATTGGTACGGGGTTGTTTATGGGGTCTGGCAAAACGATTAG‑CCTTGCCGGGCCGTCGATCATTTTCGTTTATATGATCATTGGTTTTATg < 2:700236/140‑1 (MQ=255) tctcACAAACCGACATATTCAGCTTATTGCCATTGGCGGTGCCCTTGGTACGGGGGTGTTTTTGGGGTCTGGCAAAACGATTAG‑CCTTGCCGGGCCGTCGATCATTTTCGTTTATATGATCATTGGTTTTATGCTCtttt > 2:687513/1‑140 (MQ=255) ctcACAAACCGACATATTCAGCTTATTGCCATTGGCGGGGCCCTTGGGACGGGGGTGTTTATGGGGTCTGTCAAAAAGAGTAG‑CCCTGCCGGGCCGTCGAACACTTTTGTTTATATGATCATTGGTTTTTTGCTCttttt > 2:477984/1‑140 (MQ=255) cATATTCAGCTTATTGCCATTGGCGGTGCCATTGGTACGGGGTTGTTTATGGGGTCTGGCAAAACGATTAG‑CCTTGCCGGGCCGTCGATCATTTTCGTTTATATGATCATTGGTTTTATGCTCTTTTTCGTGATGCGGGc < 1:461791/140‑1 (MQ=255) atTCAGCTTATTGCCATTGGCGGTGCCATTGGTACGGGGTTGTTTATGGGGTCTGGCAAAACGATTAG‑CCTTGCCGGGCCGTCGATCATTTTCGTTTATATGATCATTGGTTTTATGCTCTTTTTCGTGATGCGGGCAAt < 2:119901/140‑1 (MQ=255) aGCTTATTGCCATTGGCGGTGCCATTGGTACGGGGGTTTTTTTGGGGGCTGGCAAAAAGATTAG‑CCTTGCCGGGCCGGCGGTCATTTTCGTTTATATGATCATTGGTTTTATGCTCTTTTTCGTGGTGCGGGCAATgggg > 2:280991/1‑140 (MQ=255) aGCTTATTGCCATTGGCGGTGCCATTGGTACGGGGGTGTTTATGGGGGCTGGCAAAACGATTAG‑CCTTGCCGGGCCGTCGATCATTTTCGTTTATATGATCATTGGTTTTTTGCTCTTTTTCGTGATGCGGGCAATgggg > 2:712142/1‑140 (MQ=255) aTTGCCATTGGCGGTGCCATTGGTACGGGGGTGTTTTTGGGGTCTGGCAAAACGATTAG‑CCTTGCCGGGGCGTCGATCATTTTTGTTTATATGATCATTGGTTTTATGCTCTTTTTCGTGATGCGGGCAATGGGGGGAtt > 2:538498/1‑140 (MQ=255) ttGCCATTGGCGGTGCCATTGGTACGGGGTTGTTTATGGGGTCTGGCAAAACGATTAG‑CCTTGCCGGGCCGTCGATCATTTTCGTTTATATGATCATTGGTTTTATGCTCTTTTTCGTGATGCGGGCAATGGGGGAATtg < 1:423401/140‑1 (MQ=255) ttGCCATTGGCGGTGCCATTGGTACGGGGTTGTTTATGGGGTCTGGCAAAACGATTAG‑CCTTGCCGGGCCGTCGATCATTTTCGTTTATATGATCATTGGTTTTATGCTCTTTTTCGTGATGCGGGCAATGGGGGAATtg < 1:295256/140‑1 (MQ=255) cATTGGCGGTGCCATTGGTACGGGGGTGTTTATTGGGGCTGGCAAAAAGATTAG‑CCTTTCCGGGCCGTCGGTCATTTTTGTTTATATTATCATTGGTTTTTTGCTCTTTTTTGTTATGCGGGCAATGGGGGGATttctgc > 2:295256/1‑140 (MQ=255) cATTGGCGGTGCCATTGGTACGGGGGTGTTTATGGGGTCTGGCAAAACGATTAG‑CCTTGCCGGGGCGGCGATCATTTTCGTTTATATGATCATTTGGTTTATGCTCTTTTTTGTGATGCGGGCAATGGGGGGATtgctgt > 2:423401/1‑139 (MQ=255) aCGGGGTTGTTTATGGGGTCTGGCAAAACGATTAG‑CCTTGCCGGGCCGTCGATCATTTTCGTTTATATGATCATTGGTTTTATGCTCTTTTTCGTGATGCGGGCAATGGGGGAATTGCTGCTTTCGAATCTGGAATACaa > 2:267000/1‑140 (MQ=255) | GGTACAGGAAGAAAAAAACATGGTAGATCAGGTAAAAGTCGTTGCCGATGATCAGGCTCCGGCTGAACAGTCGCTACGGCGCAATCTCACAAACCGACATATTCAGCTTATTGCCATTGGCGGTGCCATTGGTACGGGGTTGTTTATGGGGTCTGGCAAAACGATTAG‑CCTTGCCGGGCCGTCGATCATTTTCGTTTATATGATCATTGGTTTTATGCTCTTTTTCGTGATGCGGGCAATGGGGGAATTGCTGCTTTCGAATCTGGAATACAA > NC_000913/4429845‑4430117 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |