Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_000913 | 3,401,459 | 0 | T | C | 29.2% | 76.7 / 10.5 | 48 | D625G (GAT→GGT) | csrD | regulator of CsrB and CsrC decay |
Reads supporting (aligned to +/- strand): ref base T (28/6); new base C (0/14); total (28/20) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 8.04e-08 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
GACAAGCACACTCCCATGAGACGCAGCGCGCATTATTCTACGTGAAAACGGATTAAACGGCAGGTTAAACCGAGTATCTTTGTGAATATTTTTTCACGTTAGTATCAAGTGGCTGTGAGGACGCGAAAAAATCCCCTTGCCCGCCTGTAACACCGCGCTGAATCAGGGTCTGCCACTCGCTTCGCGAACGCACGCCGGTGGCGTAAACCTGGGTGCTGGTCCCGGAGCAGGCTTCCACCAGGCTTTGAACCAGCAGCTGGTTCTCCGTTC > NC_000913/3401328‑3401597 | gACAAGCACACTCCCATGAGACGCAGCGCGCATTATTCTACGTGAAAACGGATTAAACGGCAGGTTAAACCGAGTATCTTTGTGAATATTTTTTCACGTTAGTATCAAGTGGCTGTGAGGACGCGAAAAAATCCCCTTGc < 2:296772/140‑1 (MQ=255) cAAGCACACTCCCATGAGACGCAGCGCGCATTATTCTACGTGAAAACGGATTAAACGGCAGGTTAAACCGAGTATCTTTGTGAATATTTTTTCACGTTAGTATCAAGTGGCTGTGAGGACGCGAAAAAATCCCCTTGccc < 1:607129/140‑1 (MQ=255) cAAGCACACTCCCATGAGACGCAGCGCGCATTATTCTACGTGAAAACGGATTAAACGGCAGGTTAAACCGAGTATCTTTGTGAATATTTTTTCACGTTAGTATCAAGTGGCTGTGAGGACGCGAAAAAACCCCCTTGccc < 2:377157/140‑1 (MQ=255) cAAGCACACTCCCATGAGACGCAGCGCGCATTATTCTACGTGAAAACGGATTAAACGGCAGGTTAAACCGAGTATCTTTGTGAATATTTTTTCACGTTAGTATCAAGTGGCTGTGAGGACGCGAAAAAACCCCCTTGccc < 2:781218/140‑1 (MQ=255) tCCCATGAGACGCAGCGCGCATTATTCTACGTGAAAACGGATTAAACGGCAGGTTAAACCGAGTATCTTTGTGAATATTTTTTCACGTTAGTATCAAGTGGCTGTGAGGACGCGAAAAAATCCCCTTGCCCGCCTGTAac > 1:336970/1‑140 (MQ=255) ccATGAAACGCAGCGCGCTTTTTTCTACGTGAAAACGGATTAAACGGCAGGTAAAACCGAGTATCTTTGTAAATATTTTTTCACGTTAGTATCAAGTGGCTGTGAGGCCGCGAAAAAATCCCCTTGCCCGCCTGTAacac < 2:347934/140‑1 (MQ=255) aTGAGACGCAGCGCGCATTATTCTACGTGAAAACGGATTAAACGGCAGGTTAAACCGAGTATCTTTGTGAATATTTTTTCACGTTAGTATCAAGTGGCTGTGAGGACGCGAAAAAATCCCCTTGCCCGCCTGTAACACcg > 1:347934/1‑140 (MQ=255) tGAGACGCAGCGCGCATTATTCTACGTGAAAACGGATTAAACGGCAGGTTAAACCGAGTATCTTTGTGAATATTTTTTCACGTTAGTATCAAGTGGCTGTGAGGACGCGAAAAAATCCCCTTGCCCGCCTGTAACACcgc > 1:695711/1‑140 (MQ=255) gagaCGCAGCGCGCATTATTCTACGAGAAAACGGATTAAAGGGCAGTTAAAACCGAGTTTCTTTGTGAATATTTTTTCACGTTAGTATAAAGTGGCTGTGAGGACGCAAAAAAATCCCCTTGCCCGCCTGTAACACcgcg < 2:336970/140‑1 (MQ=255) cGCAGCGCGCATTATTCTACGTGAAAACGGATTAAACGGCAGGTTAAACCGAGTATCTTTGTGAATATTTTTTCACGTTAGTATCAAGTGGCTGTGAGGACGCGAAAAAATCCCCTTGCCCGCCTGTAACACCGCGCTGa > 1:221532/1‑140 (MQ=255) gCAGCGCGCATTATTCTACGTGAAAACGGATTAAACGGCAGGTTAAACCGAGTATCTTTGTGAATATTTTTTCACGTTAGTATCAAGTGGCTGTGAGGACGCAAAAAAACCCCCTTGCCCGCCTGTAACACCGCGCTGaa < 1:238127/140‑1 (MQ=255) cAGCGCGCATTATTCTACGTGAAAACGGATTAAACGGCAGGTTAAACCGAGTATCTTTGTGAATATTTTTTCACGTTAGTATCAAGTGGCTGTGAGGACGCGAAAAAATCCCCTTGCCCGCCTGTAACACCGCGCTGAAt < 1:620968/140‑1 (MQ=255) cAGCGCGCATTATTCTACGTGAAAACGGATTAAACGGCAGGTTAAACCGAGTATCTTTGTGAATATTTTTTCACGTTAGTATCAAGTGGCTGTGAGGACGCGAAAAAATCCCCTTGCCCGCCTGTAACACCGCGCTGAAt > 1:253558/1‑140 (MQ=255) gcgcTTTATTCTACGTAAAAACGGATTAACCGGCAGTTTAAACCGAGTATCTTTGTGAATATTTTTTCTCTTTAGTTTCAAGGGGGTGTGGGGACGCGAAAAAACCCCCTTGCCCGCCTGTAACACCGCGCTGAATCAgg < 2:685353/140‑1 (MQ=255) gcATTATTCTACGTGAAAACGGATTAAACGGCAGGTTAAACCGAGTATCTTTGTGAATATTTTTTCACGTTAGTATCAAGTGGCTGTGAGGACGCGAAAAAACCCCCTTGCCCGCCTGTAACACCGCGCTGAATCAGGGt < 1:135056/140‑1 (MQ=255) cTACGTGAAAACGGATTAAACGGCAGGTTAAACCGAGTATCTTTGTGAATATTTTTTCACGTTAGTATCAAGTGGCTGTGAGGCCGCAAAAAAACCCCCTTGCCCGCCTGTAACACCGCGCTGAATCAGGGTCTGCCACt < 2:366187/140‑1 (MQ=255) gTGAAAACGGATTAAACGGCAGGTTAAACCGAGTATCTTTGTGAATATTTTTTCACGTTAGTATCAAGTGGCTGTGAGGACGCGAAAAAATCCCCTTGCCCGCCTGTAACACCGCGCTGAATCAGGGTCTGCCACt > 1:629552/1‑136 (MQ=255) gTGAAAACGGATTAAACGGCAGGTTAAACCGAGTATCTTTGTGAATATTTTTTCACGTTAGTATCAAGTGGCTGTGAGGACGCGAAAAAATCCCCTTGCCCGCCTGTAACACCGCGCTGAATCAGGGTCTGCCACTCGCt > 2:527115/1‑140 (MQ=255) gTGAAAACGGATTAAACGGCAGGTTAAACCGAGTATCTTTGTGAATATTTTTTCACGTTAGTATCAAGTGGCTGTGAGGACGCGAAAAAACCCCCTTGCCCGCCTGTAACACCGCGCTGAATCAGGGTCTGCCACt < 2:629552/136‑1 (MQ=255) aaaaAGGTTTAAACGCAAGGTAAAACCGAGTATTTTTGTAAATTTTTTTCCACGTTAGTATAAAGTGGCTGTGAGGCCGCAAAAAAACCCCCTGCCCCGCCTGTAACACCGCGCTGAATCAGGGTCTGCCACTCGCTTcg < 2:140895/140‑1 (MQ=255) aaaCGGATTAAACGGCAGGTTAAACCGAGTATCTTTGTGAATATTTTTTCACGTTAGTATCAAGTGGCTGTGAGGACGCGAAAAAATCCCCTTGCCCGCCTGTAACACCGCGCTGAATCAGGGTCTGCCACTCGCTTcgc > 1:568841/1‑140 (MQ=255) cGGATAAAACGGCAGGTAAAACCGAGTATCTTTGTAAATTTTTTTTCACGTTAGTATCAAGTGGTTGTGAGGACGCAAAAAAACCCCCTTGCCCGCCTGTAACACCGCGCTGAATCAGGGTCTGCCACTCGCTTCGCGaa < 2:221532/140‑1 (MQ=255) aGGTTAAACCGAGTATCTTTGTGAATATTTTTTCACGTTAGTATCAAGTGGCTGTGAGGACGCGAAAAAATCCCCTTGCCCGCCTGTAACACCGCGCTGAATCAGGGTCTGCCACTCGCTTCGCGAACGCACGCCTGtct > 1:140895/1‑138 (MQ=255) taaccggAGTTTTTTTGTGAATTTTTTTTCACGTTAGTTTCAGGGGGCTGGGAGGACGCAAAAAAACCCCCTTGCCCGCCTGTAACACCGCGCTGAATCAGGGTCTGCCACTCGCTTCGCGAACGCACGCCGGTGGCGTa < 2:568841/134‑1 (MQ=255) gTATCTTTGTGAATATTTTTTCACGTTAGTATCAAGTGGCTGTGAGGACGCGAAAAAATCCCCTTGCCCGCCTGTAACACCGCGCTGAATCAGGGTCTGCCACTCGCTTCGCGAACGCACGCCGGTGGCGTAAACCTggg > 2:195553/1‑140 (MQ=255) aaTCTTTGTGAATATTTTTTCACGTTAGTATCAAGTGGCTGTGAGGACGCGAAAAAATCCCCTTGCCCGCCTGTAACACCGCGCTGAATCAGGGTCTGCCACTCGCTTCGCGAACGCACGCCGGTGGCGTAAACCTGGGt > 2:746492/2‑140 (MQ=255) gtgAATATTTTTTCACGTTAGTATCAAGTGGCTGTGAGGACGCGAAAAAATCCCCTTGCCCGCCTGTAACACCGCGCTGAATCAGGGTCTGCCACTCGCTTCGCGAACGCACGCCGGTGGCGTAAACCTGGGTGCTGGTc > 1:617009/1‑140 (MQ=255) tttttttCAGTTTAGTATAAGGGGGTTGGGAGGACGAAAAAAACCCCCTTTCCCCGCCTGTAACACCGCGCTGAATCAGGGTCTGCCACTCGCTTCGCGAACGCACGCCGGTGGCGTAAACCTGGGTGCTGGTCCCGGAg < 1:999299/139‑1 (MQ=255) tttttttCACGTTAGTATCAAGTGGCTGGGAGGACGCGAAAAAACCCCCTTGCCCGCCTGTAACACCGCGCTGAATCAGGGTCTGCCACTCGCTTCGCGAACGCACGCCGGTGGCGTAAACCTGGGTGCTGGTCCCGGAg < 1:946111/139‑1 (MQ=255) atTTTTTCACGTTAGTATCAAGTGGCTGTGAGGACGCGAAAAAACCCCCTTGCCCGCCTGTAACACCGCGCTGAATCAGGGTCTGCCACTCGCTTCGCGAACGCACGCCGGTGGCGTAAACCTGGGTGCTGGTCCCGGAg < 1:385874/140‑1 (MQ=255) cGTTAGTATCAAGTGGCTGTGAGGACGCGAAAAAATCCCCTTGCCCGCCTGTAACACCGCGCTGAATCAGGGTCTGCCACTCGCTTCGCGAACGCACGCCGGTGGCGTAAACCTGGGTGCTGGTCCCGGAGCAGGCTTcc > 1:51050/1‑140 (MQ=255) aGTGGCTGTGAGGACGCGAAAAAATCCCCTTGCCCGCCTGTAACACCGCGCTGAATCAGGGTCTGCCACTCGCTTCGCGAACGCACGCCGGTGGCGTAAACCTGGGTGCTGGTCCCGGAGCAGGCTTCCACCAGGCTTTg > 2:667771/1‑140 (MQ=255) aGTGGCTGTGAGGACGCGAAAAAATCCCCTTGCCCGCCTGTAACACCGCGCTGAATCAGGGTCTGCCACTCGCTTCGCGAACGCACGCCGGTGGCGTAAACCTGGGTGCTGGTCCCGGAGCAGGCTTCCACCAGGCTTTg > 1:20372/1‑140 (MQ=255) aGTGGCTGTGAGGACGCGAAAAAATCCCCTTGCCCGCCTGTAACACCGCGCTGAATCAGGGTCTGCCACTCGCTTCGCGAACGCACGCCGGTGGCGTAAACCTGGGTGCTGGTCCCGGAGCAGGCTTCCACCAGGCTTTg > 2:460259/1‑140 (MQ=255) gggCTGTGGGGACGCAAAAAAACCCCCTTGCCCGCCTGTAACACCGCGCTGAATCAGGGTCTGCCACTCGCTTCGCGAACGCACGCCGGTGGCGTAAACCTGGGTGCTGGTCCCGGAGCAGGCTTCCACCAGGCTTTGaa < 2:509474/139‑1 (MQ=255) gCTGTGAGGACGCGAAAAAATCCCCTTGCCCGCCTGTAACACCGCGCTGAATCAGGGTCTGCCACTCGCTTCGCGAACGCACGCCGGTGGCGTAAACCTGGGTGCTGGTCCCGGAGCAGGCTTCCACCAGGCTTTg > 1:69132/1‑136 (MQ=255) gCTGTGAGGACGCGAAAAAATCCCCTTGCCCGCCTGTAACACCGCGCTGAATCAGGGTCTGCCACTCGCTTCGCGAACGCACGCCGGTGGCGTAAACCTGGGTGCTGGTCCCGGAGCAGGCTTCCACCAGGCTTTGAAcc > 2:732501/1‑140 (MQ=255) gCTGTGAGGACGCGAAAAAATCCCCTTGCCCGCCTGTAACACCGCGCTGAATCAGGGTCTGCCACTCGCTTCGCGAACGCACGCCGGTGGCGTAAACCTGGGTGCTGGTCCCGGAGCAGGCTTCCACCAGGCTTTGAAcc > 2:700813/1‑140 (MQ=255) gCTGTGAGGACGCGAAAAAATCCCCTTGCCCGCCTGTAACACCGCGCTGAATCAGGGTCTGCCACTCGCTTCGCGAACGCACGCCGGTGGCGTAAACCTGGGTGCTGGTCCCGGAGCAGGCTTCCACCAGGCTTTGAAcc > 2:116056/1‑140 (MQ=255) gCTGTGAGGACGCGAAAAAATCCCCTTGCCCGCCTGTAACACCGCGCTGAATCAGGGTCTGCCACTCGCTTCGCGAACGCACGCCGGTGGCGTAAACCTGGGTGCTGGTCCCGGAGCAGGCTTCCACCAGGCTTTGAAcc > 1:548500/1‑140 (MQ=255) gCTGTGAGGACGCGAAAAAATCCCCTTGCCCGCCTGTAACACCGCGCTGAATCAGGGTCTGCCACTCGCTTCGCGAACGCACGCCGGTGGCGTAAACCTGGGTGCTGGTCCCGGAGCAGGCTTCCACCAGGCTTTGAAcc > 1:509474/1‑140 (MQ=255) gCTGTGAGGACGCGAAAAAATCCCCTTGCCCGCCTGTAACACCGCGCTGAATCAGGGTCTGCCACTCGCTTCGCGAACGCACGCCGGTGGCGTAAACCTGGGTGCTGGTCCCGGAGCAGGCTTCCACCAGGCTTTGAAcc > 1:492624/1‑140 (MQ=255) gCTGTGAGGACGCGAAAAAACCCCCTTGCCCGCCTGTAACACCGCGCTGAATCAGGGTCTGCCACTCGCTTCGCGAACGCACGCCGGTGGCGTAAACCTGGGTGCTGGTCCCGGAGCAGGCTTCCACCAGGCTTTg < 2:69132/136‑1 (MQ=255) cTGTGAGGACGCGAAAAAATCCCCTTGCCCGCCTGTAACACCGCGCTGAATCAGGGTCTGCCACTCGCTTCGCGAACGCACGCCGGTGGCGTAAACCTGGGTGCTGGTCCCGGAGCAGGCTTCCACCAGGCTTTGAAcag > 2:949013/1‑138 (MQ=255) cTGTGAGGACGCGAAAAAATCCCCTTGCCCGCCTGTAACACCGCGCTGAATCAGGGTCTGCCACTCGCTTCGCGAACGCACGCCGGTGGCGTAAACCTGGGTGCTGGTCCCGGAGCAGGCTTCCACCAGGCTTTGAACca > 1:296647/1‑140 (MQ=255) gtgAGGACGCGAAAAAATCCCCTTGCCCGCCTGTAACACCGCGCTGAATCAGGGTCTGCCACTCGCTTCGCGAACGCACGCCGGTGGCGTAAACCTGGGTGCTGGTCCCGGAGCAGGCTTCCACCAGGCTTTGAACcagc > 2:982773/1‑140 (MQ=255) gAGGACGCGAAAAAATCCCCTTGCCCGCCTGTAACACCGCGCTGAATCAGGGTCTGCCACTCGCTTCGCGAACGCACGCCGGTGGCGTAAACCTGGGTGCTGGTCCCGGAGCAGGCTTCCACCAGGCTTTGAACcagcag > 2:894089/1‑140 (MQ=255) gcgAAAAAATCCCCTTGCCCGCCTGTAACACCGCGCTGAATCAGGGTCTGCCACTCGCTTCGCGAACGCACGCCGGTGGCGTAAACCTGGGTGCTGGTCCCGGAGCAGGCTTCCACCAGGCTTTGAACCAGCAGCTGGtt < 1:527115/140‑1 (MQ=255) aTCCCCTTGCCCGCCTGTAACACCGCGCTGAATCAGGGTCTGCCACTCGCTTCGCGAACGCACGCCGGTGGCGTAAACCTGGGTGCTGGTCCCGGAGCAGGCTTCCACCAGGCTTTGAACCAGCAGCTGGTTCTCCGTTc > 2:295284/1‑140 (MQ=255) | GACAAGCACACTCCCATGAGACGCAGCGCGCATTATTCTACGTGAAAACGGATTAAACGGCAGGTTAAACCGAGTATCTTTGTGAATATTTTTTCACGTTAGTATCAAGTGGCTGTGAGGACGCGAAAAAATCCCCTTGCCCGCCTGTAACACCGCGCTGAATCAGGGTCTGCCACTCGCTTCGCGAACGCACGCCGGTGGCGTAAACCTGGGTGCTGGTCCCGGAGCAGGCTTCCACCAGGCTTTGAACCAGCAGCTGGTTCTCCGTTC > NC_000913/3401328‑3401597 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 22 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |