Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_000913 | 4,429,984 | 0 | T | G | 36.8% | 50.7 / 23.2 | 38 | L41V (TTG→GTG) | cycA | serine/alanine/glycine/cycloserine:H(+)symporter |
Reads supporting (aligned to +/- strand): ref base T (3/21); new base G (14/0); total (17/21) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 7.03e-08 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
TACAGGAAGAAAAAAACATGGTAGATCAGGTAAAAGTCGTTGCCGATGATCAGGCTCCGGCTGAACAGTCGCTACGGCGCAATCTCACAAACCGACATATTCAGCTTATTGCCATTGGCGGTGCCATTGGTACGGGGTTGTTTATGGGGTCTGGCAAAACGATTAGCCTTGCCGGGCCGTCGATCATTTTCGTTTATATGATCATTGGTTTTATGCTCTTTTTCGTGATGCGGGCAATGGGGGAATTGCTGCTTTCGAATCTGGAATACAAATC > NC_000913/4429847‑4430120 | tACAGGAAGAAAAAAACATGGTAGATCAGGTAAAAGTCGTTGCCGATGATCAGGCTCCGGCTGAACAGTCGCTACGGCGCAATCTCACAAACCGACATATTCAGCTTATTGCCATTGGCGGTGGCATTGGTACGGGgttg > 1:465755/1‑140 (MQ=255) tACAGGAAGAAAAAAACATGGTAGATCAGGTAAAAGTCGTTGCCGATGATCAGGCTCCGGCTGAACAGTCGCTACGGCGCAATCTCACAAACCGACATATTCAGCTTATTGCCATTGGCGGTGCCATTGGTACGGGgt < 1:223538/138‑1 (MQ=255) tACAGGAAGAAAAAAACATGGTAGATCAGGTAAAAGTCGTTGCCGATGATCAGGCTCCGGCTGAACAGTCGCTACGGCGCAATCTCACAAACCGACATATTCAGCTTATTGCCATTGGCGGTGCCATTGGTACGGGgt > 2:223538/1‑138 (MQ=255) aaaaaaCATGGTAGATCAGGTAAAAGTCGTTGCCGATGATCAGGCTCCGGCTGAACAGTCGCTACGGCGCAATCTCACAAACCGACATATTCAGCTTATTGCCATTGGCGGTGCCATTGGTACGGGGTTGTTTATGGGGt < 1:542241/140‑1 (MQ=255) aaaCATGGTAGATCAGGTAAAAGTCGTTGCCGATGATCAGGCTCCGGCTGAACAGTCGCTACGGCGCAATCTCACAAACCGACATATTCAGCTTATTGCCATTGGCGGTGCCCTTGGTACGGGGGTGTTTTTGGGGTCTg > 2:32978/1‑140 (MQ=255) aCATGGTAGATCAGGTAAAAGTCGTTGCCGATGATCAGGCTCCGGCTGAACAGTCGCTACGGCGCAATCTCACAAACCGACATATTCAGCTTATTGCCATTGGCGGTGCCATTGGTACGGGGGTGTTTTTGGGGTCTGGc > 2:474261/1‑140 (MQ=255) ggTAGATCAGGTAAAAGTCGTTGCCGATGATCAGGCTCCGGCTGAACAGTCGCTACGGCGCAATCTCACAAACCGACATATTCAGCTTATTGCCATTGGCGGTGCCATTTGGACGGGGCTGATTATGGGGTCTGGCaaaa > 2:558031/1‑140 (MQ=255) aTCAGGTAAAAGTCGTTGCCGATGATCAGGCTCCGGCTGAACAGTCGCTACGGCGCAATCTCACAAACCGACATATTCAGCTTATTGCCATTGGCGGTGCCATTGGTACGGGGTTGTTTATGGGGTCTGGCaa < 2:272147/133‑1 (MQ=255) aTCAGGTAAAAGTCGTTGCCGATGATCAGGCTCCGGCTGAACAGTCGCTACGGCGCAATCTCACAAACCGACATATTCAGCTTATTGCCATTGGCGGTGCCATTGGTACGGGGGTGTTTATGGGGTCTGGCaa > 1:272147/1‑133 (MQ=255) cAGGTAAAAGTCGTTGCCGATGATCAGGCTCCGGCTGAACAGTCGCTACGGCGCAATCTCACAAACCGACATATTCAGCTTATTGCCATTGGCGGTGCCATTGGTACGGGGTTGTTTATGGGGTCTGGCAAAACGATTAg < 2:295175/140‑1 (MQ=255) gTAAAAGTCGTTGCCGATGATCAGGCTCCGGCTGAACAGTCGCTACGGCGCAATCTCACAAACCGACATATTCAGCTTATTGCCATTGGCGGTGCCATTGGTACGGGGGTGTTTATGGGGTCTGGCAAAACGATTAGCCt > 1:468157/1‑140 (MQ=255) tCGTTGCCGATGATCAGGCTCCGGCTGAACAGTCGCTACGGCGCAATCTCACAAACCGACATATTCAGCTTATTGCCATTGGCGGTGCCATTGGTACGGGGTTGTTTATGGGGTCTGGCAAAACGATTAGCCTTGCCggg < 1:430807/140‑1 (MQ=255) cGATGATCAGGCTCCGGCTGAACAGTCGCTACGGCGCAATCTCACAAACCGACATATTCAGCTTATTGCCATTGGCGGTGCCATTGGTACGGGGTTGTTTATGGGGTCTGGCAAAACGATTAGCCTTGCCGGGCCGTCGa < 1:315848/140‑1 (MQ=255) atCAGGCTCCGGCTGAACAGTCGCTACGGCGCAATCTCACAAACCGACATATTCAGCTTATTGCCATTGGCGGTGCCATTGGTACGGGGTTGTTTATGGGGTCTGGCAAAACGATTAGCCTTGCCGGGCCGTCGATCAtt < 1:429754/140‑1 (MQ=255) tCAGGCTCCGGCTGAACAGTCGCTACGGCGCAATCTCACAAACCGACATATTCAGCTTATTGCCACTGGCGGTGCCATTGGTACGGGGTTGTTTATGGGGTCTGGCAAAACGATTAGCCTTGCCGGGCCGTCGATCAttt < 1:45316/140‑1 (MQ=255) ggCTCCGGCTGAACAGTCGCTACGGCGCAATCTCACAAACCGACATATTCAGCTTATTGCCATTGGCGGTGCCATTGGTACGGGGTTGTTTATGGGGTCTGGCAAAACGATTAGCCTTGCCGGGCCGTCGATCATTTTCg < 2:758425/140‑1 (MQ=255) ggCTCCGGCTGAACAGTCGCTACGGCGCAATCTCACAAACCGACATATTCAGCTTATTGCCATTGGCGGTGCCATTGGTACGGGGTTGTTTATGGGGTCTGGCAAAACGATTAGCCTTGCCGGGCCGTCGATCATTTTCg < 1:378706/140‑1 (MQ=255) gCTGAACAGTCGCTACGGCGCAATCTCACAAACCGACATATTCAGCTTATTGCCATTGGCGGTGCCATTGGGACGGGGGTGTTTTTTGGGTCTTGCAAAAAGATTAGCCTTGCCGGGGCGTCGATCCTTTTTGTTTttat > 2:240067/1‑140 (MQ=255) gCTGAACAGTCGCTACGGCGCAATCTCACAAACCGACATATTCAGCTTATTGCCATTGGCGGTGCCATTGGGACGGGGGTGTTTATTGGGTCTGGCAAAAAGATTAGCCTTTCCGGGCCGTCGATCATTTTTGTTTatat > 1:397203/1‑140 (MQ=255) aaCAGTCGCTACGGCGCAATCTCACAAACCGACATATTCAGCTTATTGCCATTGGCGGTGCCATTGGTACGGGGGTGTTTATGGGGTCTGGCAAAACGATTAGCCTTGCCGGGCCGTCGATCATTTTCGTTTATATGATc > 1:409696/1‑140 (MQ=255) cAGTCGCTACGGCGCAATCTCACAAACCGACATATTCAGCTTATTGCCATTGGCGGTGCCCTTGGTACGGGGGTGTTTTTGGGGGCTGGCAAAACGATTAGCCTTGCCGGGCCGTCGATCATTTTCGTTTATATGATCAt > 2:461851/1‑140 (MQ=255) tCGCTACGGCGCAATCTCACAAACCGACATATTCAGCTTATTGCCATTGGCGGTGCCATTGGGACGGGGTTGTTTTTGGGGTCTGGCAAAACGATTAGCCTTGCCGGGCCGTCGATCATTTTCGTTTATATGATCATTgg > 1:350622/1‑140 (MQ=255) tACGGCGCAATCTCACAAACCGACATATTCAGCTTATTGCCATTGGCGGTGCCATTGGTACGGGGTTGTTTATGGGGTCTGGCAAAACGATTAGCCTTGCCGGGCCGTCGATCATTTTCGTTTATATGATCATTGGtttt < 2:409696/140‑1 (MQ=255) cGGCGCAATCTCACAAACCGACATATTCAGCTTATTGCCATTGGCGGTGCCATTGGTACGGGGTTGTTTATGGGGTCTGGCAAAACGATTAGCCTTGCCGGGCCGTCGATCATTTTCGTTTATATGATCATTGGTTTTAt < 2:587594/140‑1 (MQ=255) aCCGACATATTCAGCTTATTGCCATTGGCGGTGCCCTTGGTACCGGGGTGTTTATGGGGGCTTGCAAAACGATTTTGCTTTCCGGGCCGTCGATCATTTTCGGTTATATGATCATTGGTTTTATGCTCTTTTTCGTGATg > 2:73698/1‑140 (MQ=255) aCCGACATATTCAGCTTATTGCCATTGGCGGTGCCATTGGTACGGGGTTGTTTATGGGGTCTGGCAAAACGATTAGCCTTGCCGGGCCGTCGATCATTTTCGTTTATATGATCATTGGTTTTATGCTCTTTTTCGTGATg < 2:468157/140‑1 (MQ=255) aCCGACATATTCAGCTTATTGCCATTGGCGGTGCCATTGGTACGGGGTTGTTTATGGGGTCTGGCAAAACGATTAGCCTTGCCGGGCCGTCGATCATTTTCGTTTATATGATCATTGGTTTTATGCTCTTTTTCGTGATg < 2:465755/140‑1 (MQ=255) ccGACATATTCAGCTTATTGCCATTGGCGGTGCCATTGGTAC‑GGGTTGTTTATGGGGTCTGGCAAAACGATTAGCCTTGCCGGGCCGTCGATCATTTTCGTTTATATGATCATTGGTTTTATGCTCTTTTTCGTGATGCg < 2:557927/140‑1 (MQ=255) cGACATATTCAGCTTATTGCCATTGGCGGTGCCATTGGTACGGGGTTGTTTATGGGGTCTGGCAAAACGATTAGCCTTGCCGGGCCGTCGATCATTTTCGTTTATATGATCATTGGTTTTATGCTCTTTTTCGTGATGCg < 1:474261/140‑1 (MQ=255) cGACATATTCAGCTTATTGCCATTGGCGGTGCCATTGGTACGGGGTTGTTTATGGGGTCTGGCAAAACGATTAGCCTTGCCGGGCCGTCGATCATTTTCGTTTATATGATCATTGGTTTTATGCTCTTTTTCGTGATGCg < 1:48709/140‑1 (MQ=255) cGACATATTCAGCTTATTGCCATTGGCGGTGCCATTGGTACGGGGGTGTTTATGGGGTCTGGCAAAACGATTAGCCTTGCCGGGCCGTCGATCATTTTCGTTTATATGATCATTGGTTTTATGCTCTTTTTCGTGATGCg > 1:482464/1‑140 (MQ=255) cGACATATTCAGCTTATTGCCATTGGCGGTGCCATTGGGACGGGGGTGTTTATGGGGGCTGGCAAAACGATTAGCCTTTCCGGGCCGTCGGTCATTTTCGTTTATATTATCATTGGTTTTTTGCTCTTTTTCGTGATGCg > 2:585234/1‑140 (MQ=255) atTCAGCTTATTGCCATTGGCGGTGCCATTGGTACGGGGTTGTTTATGGGGTCTGGCAAAACGATTAGCCTTGCCGGGCCGTCGATCATTTTCGTTTATATGATCATTGGTTTTATGCTCTTTTTCGTGATGCGGGCAAt < 2:287681/140‑1 (MQ=255) atTCAGCTTATTGCCATTGGCGGTGCCATTGGTACGGGGTTGTTTATGGGGTCTGGCAAAACGATTAGCCTTGCCGGGCCGTCGATCATTTTCGTTTATATGATCATTGGTTTTATGCTCTTTTTCGTGATGCGGGCAAt < 2:11050/140‑1 (MQ=255) cTTATTGCCATTGGCGGTGCCATTTGTACGGGGGTGTTTTTGGGGGCTGGCAAAACGAATAGCCTTGCCGGGGCGTCGATCATTTTTGTTTATATGATCATTGGTTTTTTGCTCTTTTTTGTGATGCGGGCAATgggggg > 2:594550/1‑139 (MQ=255) ttATTGCCATTGGCGGTGCCATTGGTACGGGGTTGTTTATGGGGTCTGGCAAAACGATTAGCCTTGCCGGGCCGTCGATCATTTTCGTTTATATGATCATTGGTTTTATGCTCTTTTTCGTGATGCGGGCAATGGGGGaa < 1:415014/140‑1 (MQ=255) aTTGCCATTGGCGGTGCCATTGGTTCGGGGGTGTTTATGGGGTCTGGCAAAACGATTAGCCTTGCCGGGCCGTCGATCATTTTCGTTTATATGATCATTGGTTTTATGCTCTTTTTCGTGATGCGGGCAATGGGGGAAtt > 1:329644/1‑140 (MQ=255) aTTGCCATTGGCGGTGCCATTGGTACGGGGGTGTTTTTGGGGGCTGGCAAAACGATTAGCCCTCCCCGGCCGTCGATCATTTTCGTTTATATGATCATTGGTTTTTTGCTCTTTTTTGTGATGCGGGCAATGGGGGGAtt > 2:253022/1‑140 (MQ=255) ttGGCGGTGCCATTGGTACGGGGGTGTTTTTGGGGGCTGGCAAAACGATTAGCCTTTCCGGGCCATCGATCATTTTTGTTTATATGATCATTTGTTTTTTGCTCTTTTTTGTGATGCGGGCAATGGGGGGAATGCTGCtt > 2:758228/1‑140 (MQ=255) ttGGCGGTGCCATTGGTACGGGGGTGTTTATGGGGTCTGGCAAAACGATTAGCCTTGCCGGGCCGTCGATCATTTTCGTTTATATGATCATTGGTTTTATGCTCTTTTTCGTGATGCGGGCAATGGGGGAATTGCTGCtt > 2:223777/1‑140 (MQ=255) tGCCATTGGTACGGGGTTGTTTATGGGGTCTGGCAAAACGATTAGCCTTGCCGGGCCGTCGATCATTTTCGTTTATATGATCATTGGTTTTATGCTCTTTTTCGTGATGCGGGCAATGGGGGAATTGCTGCTTTCGAATc < 1:253022/140‑1 (MQ=255) ggggTTGTTTATGGGGTCTGGCAAAACGATTAGCCTTGCCGGGCCGTCGATCATTTTCGTTTATATGATCATTGGTTTTATGCTCTTTTTCGTGATGCGGGCAATGGGGGAATTGCTGCTTTCGAATCTGGAATACaa < 2:302118/138‑1 (MQ=255) ggggTTGTTTATGGGGTCTGGCAAAACGATTAGCCTTGCCGGGCCGTCGATCATTTTCGTTTATATGATCATTGGTTTTATGCTCTTTTTCGTGATGCGGGCAATGGGGGAATTGCTGCTTTCGAATCTGGAATACaa > 1:302118/1‑138 (MQ=255) gggTTGTTTATGGGGTCTGGCAAAACGATTAGCCTTGCCGGGCCGTCGATCATTTTCGTTTATATGATCATTGGTTTTATGCTCTTTTTCGTGATGCGGGCAATGGGGGAATTGCTGCTTTCGAATCTGGAATACAAATc > 2:22087/1‑140 (MQ=255) | TACAGGAAGAAAAAAACATGGTAGATCAGGTAAAAGTCGTTGCCGATGATCAGGCTCCGGCTGAACAGTCGCTACGGCGCAATCTCACAAACCGACATATTCAGCTTATTGCCATTGGCGGTGCCATTGGTACGGGGTTGTTTATGGGGTCTGGCAAAACGATTAGCCTTGCCGGGCCGTCGATCATTTTCGTTTATATGATCATTGGTTTTATGCTCTTTTTCGTGATGCGGGCAATGGGGGAATTGCTGCTTTCGAATCTGGAATACAAATC > NC_000913/4429847‑4430120 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 6 ≤ ATCG/ATCG < 22 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |