Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | NC_002947 | 1,728,699 | C→T | 11.5% | intergenic (+51/+72) | cspA‑I → / ← PP_1523 | major cold shock protein/hypothetical protein |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_002947 | 1,728,699 | 0 | C | T | 11.5% | 130.1 / 8.6 | 52 | intergenic (+51/+72) | cspA‑I/PP_1523 | major cold shock protein/hypothetical protein |
Reads supporting (aligned to +/- strand): ref base C (23/23); new base T (3/3); total (26/26) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 1.00e+00 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.35e-01 |
GGGCATCGTACCCTGGTCGAAGGGCAGCGGGTGGAGTACGCCTGCGTGCAGGGCCAGAAAGGCCTGCAAGCCGAAGACGTAGTGGGCCTCTGAGCCTCTGGTTTCAAGCCGTAGGAGCGGGCTTGCCCGCGAAGCAGCCACCACGGTGGCTGC‑‑CACGGGCATAGCCCGTATTCGCGGGTAAACCCATTCCTACAGAAACGAGCAAGCCTGCCGTTTCATGATGTACGCCAGGTGATTTCCTCTTCACCGTCGGCGCTGATGCGGACCCAGCGGTCAGCATCTTCTTCCCC > NC_002947/1728556‑1728845 | gggCATCGTACCCTGGTCGAAGGGCAGCGGGTGGAGTACGCCTGCGTGCAGGGCCAGAAAGGCCTGCAAGCCGAAGACGTAGTGGGCCTCTGAGCCTCTGGTTTCAAGCCGTAGGAGCGGGCTTGCCCGCGAAGAAGCCACCACGGTggg > 2:334950/1‑149 (MQ=255) cGTACCCTGGTCGAAGGGCAGCGGGTGGAGTACGCCTGCGGGCAGTGCCAGAAAGGCCTGCAAGCCGAAGACGTAGTCGGCCTCTGAGCCTCTGGTTTCAAGCCGTAGGAGCGGGCTTGCCCGCGAAGCAGCCACCGTGGTGGCTGcttc < 1:33773/150‑4 (MQ=255) ccTGGTCGAAGGGCAGCGGGTGGAGTACGCCTGCGTGCAGGGCCAGAAAGGCCTGCAAGCCGAAGACGTAGTGGGCCTCTGAGCCTCTGGTTTCAAGCCGTAGGAGCGGGCTTGCCCGCGAAGCAGCCACCACGGTGGCTGC‑‑CACGGGCa > 1:176405/1‑150 (MQ=255) cGAAGGGCAGCGGGTGGAGTACGCCTGCGTGCAGGGCCAGAAAGTCCTGCAAGCCGAAGACGTAGTGGGCCTCTGAGCCTCTGGTTTCAAGCCGTAGGAGCGGGCTTGCCCGCGAAGCAGCCACCACGGTGGCTGC‑‑CACGGGCATAGccc < 2:270225/150‑1 (MQ=255) aGGGCAGCGGGTGGAGTACGCCTGCGTGCAGGGCCAGAAAGGCCTGCAAGCCGAAGACGTAGTGGGCCTCTGAGCCTCTGGTTTCAAGCCGTAGGAGCGGGCTTGCCCGCGAAGCAGCCACCACGGTGGCTGC‑‑CACGGGCATAGCCCGTa < 2:546291/150‑1 (MQ=255) cGGGTGGAGTACGCCTGCGTGCAGGGCCAGAAAGGCCTGCAAGCCGAAGACGTAGTGGGCCTCTGAGCCTCTGGTTTCAAGCCGTAGGAGCGGGCTTGCCCGCGAAGCAGCCACCACGGTGGCTGC‑‑CACGGGCATAGCCCGTATTCGCgg < 2:392386/150‑1 (MQ=255) aGTACGCCTGCGTGCAGGGCCAGAAAGGCCTGCAAGCCGAAGACGTAGTGGGCCTCTGAGCCTCTGGTTTCAAGCCGTAGGAGCGGGCTTGCCCGCGAAGCAGCCACCACGGTGGCTGC‑‑CACGGGCATAGCCCGTATTCGCGGGTAAAcc > 2:1137955/1‑150 (MQ=255) ccTGCGTGCAGGGCCAGAAAGGCCTGCAAGCCGAAGACGTAGTGGGCCTCTGAGCCTCTGGTTTCAAGCCGTAGGAGCGGGCTTGCCCGCGAAGCAGCCACCACGGTGGCTGC‑‑CACGGGCATAGCCCGTATTCGCGGGTAAACCCATTcc < 1:1137955/150‑1 (MQ=255) cTGCGTGCAGGGCCAGAAAGGCCTGCAAGCCGAAGACGTAGTGGGCCTCTGAGCCTCTGGTTTCAAGCCGTAGGAGCGGGCTTGCCCGCGAAGCAGCCACCACGGTGGCTGC‑‑CACGGGCATAGCCCGTATTcgcg < 2:129981/135‑1 (MQ=255) cTGCGTGCAGGGCCAGAAAGGCCTGCAAGCCGAAGACGTAGTGGGCCTCTGAGCCTCTGGTTTCAAGCCGTAGGAGCGGGCTTGCCCGCGAAGCAGCCACCACGGTGGCTGC‑‑CACGGGCATAGCCCGTATTcgcg > 1:129981/1‑135 (MQ=255) gCGTGCAGGGCCAGAAAGGCCTGCAAGCCGAAGACGTAGTGGGCCTCTGAGCCTCTGGTTTCAAGCCGTAGGAGCGGGCTTGCCCGCGAAGCAGCCACCACGGTGGCTGC‑‑CACGGGCATAGCCCGTATTCGCGGGTAAACCCATTCCTAc < 2:661604/150‑1 (MQ=255) gCGGGCAGGGCCAGAAAGGCCTGCAAGCCGAAGACGTAGTGGGCCTCTGAGCCTCTGGTTTCAAGCCGTAGGAGCGGGCTTGCCCGCGAAGCAGCCACCACGGTGGCTGC‑‑CACGGGCATAGCCCGTATTCGCGGGTAAACCCATTCCTAc < 1:482021/150‑1 (MQ=255) gCAAGCCGAAGACGTAGTGGGCCTCTGAGCCTCTGGTTTCAAGCCGTAGGAGCGGGCTTGCCCGCGAAGCAGCCACCACGGTGGCTGC‑‑CACGGGCATAGCCCGTATTCGCGGGTAAACCCATTCCTACAGAAACGAGCCAGCCTGCCGtt > 1:572332/1‑150 (MQ=255) cGAAGACGTAGTGGGCCTCTGAGCCTCTGGTTTCAAGCCGTAGGAGCGGGCTTGCCCGCGAAGCAGCCACCACGGTGGCTGC‑‑CACGGGCATAGCCCGTATGCGCGGGTAAACCCATTCCTACAGGAACGAGCAAGCCTGCCGTTTCatga < 2:176405/150‑1 (MQ=255) gAAGACGTAGTGGGCCTCTGAGCCTCTGGTTTCAAGCCGTAGGAGCGGGCTTGCCCGCGAAGCAGCCACCACGGTGGCTGC‑‑CACGGGCATAGCCCGTATTCGCGGGTAAACCCATTCCTACAGa < 2:1151425/124‑1 (MQ=255) gAAGACGTAGTGGGCCTCTGAGCCTCTGGTTTCAAGCCGTAGGAGCGGGCTTGCCCGCGAAGCAGCCACCACGGTGGCTGC‑‑CACGGGCATAGCCCGTATTCGCGGGTAAACCCATTCCTACAGa > 1:1151425/1‑124 (MQ=255) aGACGTAGTGGGCCTCTGAGCCTCTGGTTTCAAGCCGTAGGAGCGGGCTTGCCCGCGAAGCAGCCACCGTGGTGGCTGc < 1:772321/79‑1 (MQ=38) aGACGTAGTGGGCCTCTGAGCCTCTGGTTTCAAGCCGTAGGAGCGGGCTTGCCCGCGAAGCAGCCACCGTGGTGGCTGc > 2:772321/1‑79 (MQ=38) aGACGTAGTGGGCCTCTGAGCCTCTGGTTTCAAGCCGTAGGAGCGGGCTTGCCCGCGAAGCAGCCACCACGGTGGCTGC‑‑CACGGGCATAGCCCGTATTCGCGGGTAAACCCATTCCTACAGAAACGAGCAAGc > 1:133954/1‑133 (MQ=255) aGACGTAGTGGGCCTCTGAGCCTCTGGTTTCAAGCCGTAGGAGCGGGCTTGCCCGCGAAGCAGCCACCACGGTGGCTGC‑‑CACGGGCATAGCCCGTATTCGCGGGTAAACCCATTCCTACAGAAACGAGCAAGc < 2:133954/133‑1 (MQ=255) ccTCTGAGCCTCTGGTTTCAAGCCGTAGGAGCGGGCTTGCCCGCGAAGCAGCCACCACGGTGGCTGC‑‑CACGGGCATAGCCCGTATTCGCGGGTAAAc < 2:517269/97‑1 (MQ=255) ccTCTGAGCCTCTGGTTTCAAGCCGTAGGAGCGGGCTTGCCCGCGAAGCAGCCACCACGGTGGCTGC‑‑CACGGGCATAGCCCGTATTCGCGGGTAAAc > 1:517269/1‑97 (MQ=255) ccTCTGAGCCTCTGGTTTCAAGCCGTAGGAGCGGGCTTGCCCGCGAAGCAGCCACCACGGTGGCTGC‑‑CACGGGCATAGCCCGTATTCGCGGGTAAACCCATTCCTACAGAAACGAGCAAGCCTGCCGTTTCATGATGTACGCCAGGTGAt > 1:488334/1‑150 (MQ=255) ctctGAGCCTCTGGTTTCAAGCCGTAGGAGCGGGCTTGCCCGCGAAGCAGCCACCACGGTGGCTGC‑‑CACGGGCATAGCCCGTATTCGCGGGTAAACCCATTCCTACAGAAACGAGCAAGCCCGCCGTTTCATGATGTACGCCAGGTGAtt > 2:573856/1‑150 (MQ=255) tctGAGCCTCTGGTTTCAAGCCGTAGGAGCGGGCTTGCCCGCGAAGCAGCCACCACGGTGGCTGC‑‑CACGGGCATAGCCCGTATTCGCGGGTAAACCCATTCCTACAGAAACGAGCAAGCCTGCCGTTTCATGATGTACGCCAGGTGAttt > 2:890726/1‑150 (MQ=255) ccTCTGGTTTCAAGCCGTAGGAGCGGGCTTGCCCGCGAAGCAGCCACCACGGTGGCTGC‑‑CACGGGCATAGCCCGTATTCGCGGGTAAACCCATTCCTACAGAAACGAGCAAGCCTGCCGTTTCATGATGTACGCCAGGTGATTTCCTCtt < 2:684381/150‑1 (MQ=255) ccTCTGGTTTCAAGCCGTAGGAGCGGGCTTGCCCGCGAAGCAGCCACCACGGTGGCTGC‑‑CACGGGCATAGCCCGGATTCGCGGGTAAACCCATTCCTACAGAAACGGGCAAGCCTGCCGTTTCATGATGTACGCCAGGTGATTTCCTCtt > 1:684381/1‑150 (MQ=255) tGGTTTCAAGCCGTAGGAGCGGGCTTGCCCGCGAAGCAGCCACCACGGTGGCTGC‑‑CACGGGCATAGCCCGTATTCGCGGGTAAACCCATTCCTACAGAAACGAGCAAGCCTGc > 1:1041489/1‑113 (MQ=255) tGGTTTCAAGACGTAGGAGCGGGCTTGCCCGCGAAGCAGCCACCACGGTGGCTGC‑‑CACGGGCATAGCCCGTATTCGCGGGTAAACCCATTCCTACAGAAACGAGCAAGCCTGc < 2:1041489/113‑1 (MQ=255) tttCAAGCCGTAGGAGCGGGCTTGCCCGCGAAGCAGCCACCACGGTGGCTGC‑‑CACGGGCATAGCCCGTATTCGCGGGTAAACCCATTCCTACAGAAACGAGCAAGCCTGCCGTTTCATGATGTACGCCAGGTGATTTCCTCTTCACCGTc > 2:108489/1‑150 (MQ=255) tCAAGCCGTAGGAGCGGGCTTGCCCGCGAAGCAGCCACCACGGTGGCTGC‑‑CACGGGCATATCCCGTATTCGCGGGTAAACCCATTCCTACATAAACGAGCAAGCCTGCCGTTTCATGATGTACGCCAGGTGATTTCCTCTTCACCGTCgg < 1:193962/150‑1 (MQ=255) aaGCCGTAGGAGCGGGCTTGCCCGCGAAGCAGCCACCACGGTGGCTGC‑‑CACGGGCATAGCCCGTATTCGCGGGTAAACCCATTCCTACAGAAACGAGCAAGCCTGCCGTTTCATGATGTACGCCAGGTGATTTc < 1:704067/134‑1 (MQ=255) aaGCCGTAGGAGCGGGCTTGCCCGCGAAGCAGCCACCACGGTGGCTGC‑‑CACGGGCATAGCCCGTATTCGCGGGTAAACCCATTCCTACAGAAACGAGCAAGCCTGCCGTTTCATGATGTACGCCAGGTGATTTc > 2:704067/1‑134 (MQ=255) aaGCCGTAGGAGCGGGCTTGCCCGCGAAGCAGCCACCACGGTGGCTGC‑‑CACGGGCATAGCCCGTATTCGCGGGTAAACCCATTCCTACAGAAACGAGCAAGCCTGCCGTTTCATGATGTACGCCAGGTGATTTCCTCTTCACCGTCGgcg > 2:178494/1‑150 (MQ=255) aaGCCGTAGGAGCGGGCTTGCCCGCGAAGCAGCCACCACGGTGGCTGC‑‑CACGGGCATAGCCCGTATTCGCGGGTAAACCCATTCCTACAGAAACGAGCAAGCCTGCCGTTTCATGATGTACGCCAGGTGATTTCCTCTTCACCGTCGgcg < 1:890726/150‑1 (MQ=255) aaGCCGTAGGAGCGGGCTTGCCCGCGAAGCAGCCACCACGGTGGCTGC‑‑CACGGGCATAGCCCGTATTAGCGGGTAAACCCATTCCTACAGAAACGAGCAAGCCTGCCGTTTTATGACGTACGCCCGATGATTTTCTCTTAACCGTgggac > 1:249945/1‑148 (MQ=255) gCCGTAGGAGCGGGCTTGCCCGCGAAGCAGCCACCACGGTGGCTGC‑‑CACGGGCATAGCCCGTATTCGCGGGTAAACCCATTCCTACAGAAACGAGCAAGCCTGCCGTTTCATGATGTACGCCAGGt > 2:381962/1‑126 (MQ=255) gCCGTAGGAGCGGGCTTGCCCGCGAAGCAGCCACCACGGTGGCTGC‑‑CACGGGCATAGCCCGTATTCGCGGGTAAACCCATTCCTACAGAAACGAGCAAGCCTGCCGTTTCATGATGTACGCCAGGt < 1:381962/126‑1 (MQ=255) gCTTGCCCGCGAAGCAGCCACCACGGTGGCTGC‑‑CACGGGCATAGCCCGTATTCGCGGGTAAACCCATTCCTACAGAAACGAGCAAGCCTGCCGTTTCATGATGTACGCCAGGTGATTTCCTCTTCACCGTCGGCGCTGATGCGGACCCAg > 2:852612/1‑150 (MQ=255) ccGCGAAGCAGCCACCACGGTGGCTGC‑‑CACGGGCATAGCCCGTATTCGCGGGTAAACCCATTCCTACAGAAACGAGCAAGCCTGCCGTTTCATGATGTACGCCAGGTGATTTCCTCTTCACCGTCGGCGCTGATGCGGACCCAGCGGTCa > 2:456262/1‑150 (MQ=255) ccGCGAAGCAGCCACCACGGTGGCTGC‑‑CACGGGCATAGCCCGTATTCGCGGGTAAACCCATTCCTACAGAAACGAGCAAGCCTGCCGTTTCATGATGTACGCCAGGTGATTTCCTCTTCACCGTCGGCGCTGATGCGGACCCAGCGGTCa < 1:456262/150‑1 (MQ=255) cgcgAAGCAGCCACCACGGTGGCTGC‑‑CACGGGCATAGCCCGTATTCGCNGGAAAACCCATTCCTACAGAAACGAGCAAGCCTGCCGTTTCATGATGTACGCGAGGTGATTTACTCTTCACCGTCGGCGCTGATGCGGACCCAGCGGTCAg > 1:831008/1‑150 (MQ=255) gcgAAGCAGCCACCACGGTGGCTGC‑‑CACGGGCATAGCCCGTATTCGCGGGTAAACCCATTCCTACAGAAACGAGCAAGCCTGCCGTTTCATGATGTACGCCAGGTGATTTCCTCTTCACCGTCGGCGCTGATGCg < 2:25556/135‑1 (MQ=255) gcgAAGCAGCCACCACGGTGGCTGC‑‑CACGGGCATAGCCCGTATTCGCGGGTAAACCCATTCCTACAGAAACGAGCAAGCCTGCCGTTTCATGATGTACGCCAGGTGATTTCCTCTTCACCGTCGGCGCTGATGCg > 1:25556/1‑135 (MQ=255) ccgtggcagcCACCGTGGTGGCTGCTTCGCGGGCATAGCCCGTATTCGCGGGTAAACCCATTCCTACAGAAACGAGCAAGCCTGCCGTTTCATGATGTACGCCAGGTGATTTCCTCTTCACCGTCGGCGCTGATGCGGACCCAGCGGTCa > 1:409443/6‑150 (MQ=38) ccgtggcagcCACCGTGGTGGCTGC‑‑CACGGGCATAGCCCGTATTCGCGGGTAAACCCATTCCTACAGAAACGAGCAAGCCTGCCGTTTCATGATGTACGCCAGGTGATTTCCTCTTCACCGTCGGCGCTGATGCGGACCCAGc > 1:156590/6‑143 (MQ=255) ccgtggcagcCACCGTGGTGGCTGC‑‑CACGGGCATAGCCCGTATTCGCGGGTAAACCCATTCCTACAGAAACGAGCAAGCCTGCCGTTTCATGATGTACGCCAGGTGATTTCCTCTTCACCGTCGGCGCTGATGCGGACCCAGc < 2:156590/138‑1 (MQ=255) gAAGCAGCCACCACGGTGGCTGC‑‑CACGGGCATAGCCCGTATTCGCGGGTAAACCCATTCCTACAGAAACGAGCAAGCCTGCCGTTTCATGATGTACGCCAGGTGATTTCCTCTTCACCGTCg < 1:430004/122‑1 (MQ=255) gAAGCAGCCACCACGGTGGCTGC‑‑CACGGGCATAGCCCGTATTCGCGGGTAAACCCATTCCTACAGAAACGAGCAAGCCTGCCGTTTCATGATGTACGCCAGGTGATTTCCTCTTCACCGTCg > 2:430004/1‑122 (MQ=255) aaGCAGCCACCACGGTGGCTGC‑‑CACGGGCATAGCCCGTATTCGCGGGTAAACCCATTCCTACAGAAACGAGCAAGCCTGCCGTTTCATGATGAACGCCAGGTGATTTCCTCTTCACCGTCGGCGCTGATGCGGACCCAGCGGTCAGCAtc < 1:647202/150‑1 (MQ=255) ccaccaCGGTGGCTGT‑‑CACGGGCATATCCCGTATTCGCGGGTAAACCCATTCCTACAGAAACGAGCAAGCCTGCCGTTTCATGATGTACGCCAGGTGATTTCCTCTTCACCGTCGGCGCTGATGCGGACCCAGCGGTCAGCATCttcttc < 2:249945/150‑1 (MQ=255) ccaccaCGGTGGCTGC‑‑CACGGGCATAGCCCGTATTCGCGGGTAAACCCATTCCTACAGAAACGAGCAAGCCTGCCGTTTCATGATGTACGCCAGGTGATTTCCTCTTCACCGTCGGCGCTGATGCGGACCCAGCGGTCAGCATCttcttc < 1:108489/150‑1 (MQ=255) accacGGTGGCTGC‑‑CACGGGCATAGCCCGTATTCGCGGGTAAACCCATTCCTACAGAAACGAGCAAGCCTGCCGTTTCATGATGTACGCCAGGTGATTTCCTCTTCACCGTCGGCGCTGATGCGGACCCAGCGGTc > 2:422138/1‑136 (MQ=255) accacGGTGGCTGC‑‑CACGGGCATAGCCCGTATTCGCGGGTAAACCCATTCCTACAGAAACGAGCAAGCCTGCCGTTTCATGATGTACGCCAGGTGATTTCCTCTTCACCGTCGGCGCTGATGCGGACCCAGCGGTc < 1:422138/136‑1 (MQ=255) accacGGTGGCTGC‑‑CACGGGCATAGCCCGTATTCGCGGGTAAACCCATTCCTACAGAAACGAGCAAGCCTGCCGTTTCATGATGTACGCCAGGTGATTTCCTCTTCACCGTCGGCGCTGATGCGGACCCAGCGGTCAGCATCTTCTTccc < 2:488334/150‑1 (MQ=255) ccacGGTGGCTGC‑‑CACGGGCATAGCCCGTATTCGCGGGTAAACCCATTCCTACAGAAACGAGCAAGCCTGCCGTTTCATGATGTACGCCAGGTGATTTCCTCTTCACCGTCGGCGCTGATGCGGACCCAGCGGTCAGCATCTTCTTcccc < 1:191263/150‑1 (MQ=255) ccacGGTGGCTGC‑‑CACGGGCATAGCCCGTATTCGCGGGTAAACCCATTCCTACAGAAACGAGCAAGCCTGCCGTTTCATGATGTACGCCAGGTGATTTCCTCTTCACCGTCGGCGCTGATGCGGACCCAGCGGTCAGCATCTTCTTcccc > 2:5380/1‑150 (MQ=255) | GGGCATCGTACCCTGGTCGAAGGGCAGCGGGTGGAGTACGCCTGCGTGCAGGGCCAGAAAGGCCTGCAAGCCGAAGACGTAGTGGGCCTCTGAGCCTCTGGTTTCAAGCCGTAGGAGCGGGCTTGCCCGCGAAGCAGCCACCACGGTGGCTGC‑‑CACGGGCATAGCCCGTATTCGCGGGTAAACCCATTCCTACAGAAACGAGCAAGCCTGCCGTTTCATGATGTACGCCAGGTGATTTCCTCTTCACCGTCGGCGCTGATGCGGACCCAGCGGTCAGCATCTTCTTCCCC > NC_002947/1728556‑1728845 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |