Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | NC_002947 | 1,728,698 | A→G | 20.9% | intergenic (+50/+73) | 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,698 | 0 | A | G | 20.9% | 39.5 / 7.1 | 24 | intergenic (+50/+73) | cspA‑I/PP_1523 | major cold shock protein/hypothetical protein |
Reads supporting (aligned to +/- strand): ref base A (10/9); new base G (2/3); total (12/12) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 1.00e+00 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.21e-01 |
GCATCGTACCCTGGTCGAAGGGCAGCGGGTGGAGTACGCCTGCGTGCAGGGCCAGAAAGGCCTGCAAGCCGAAGACGTAGTGGGCCTCTGAGCCTCTGGTTTCAAGCCGTAGGAGCGGGCTTGCCCGCGAAGCAGCCACCACGGTGGCTGCCACGG‑GCATAGCCCGTATTCGCGGGTAAACCCATTCCTACAGAAACGAGCAAGCCTGCCGTTTCATGATGTACGCCAGGTGATTTCCTCTTCACCGTCGGCGCTGATGCGGACCCAGCGGTCAGCATCTTCTTCCC > NC_002947/1728558‑1728844 | gCATCGTACCCTGGTCGAAGGGCAGCGGGTGGAGTACGCCTGCGTGCAGGGCCAGAAAGGCCTGCAAGCCGAAGACGTAGTGGGCCTCTGAGCCTCTGGTTTCAAGCCGTAGGAGCGGGCTTGCCCGCGAAGCAGCCACCACGGTGGCTg > 2:268057/1‑150 (MQ=255) gCATCGTACCCTGGTCGAAGGGCAGCGGGTGGAGTACGCCTGCGTGCAGGGCCAGAAAGGCCTGCAAGCCGAAGACGTAGTGGGCCTCTGAGCCTCTGGTTTCAAGCCGTAGGAGCGGGCTTGCCCGCGAAGCAGCCACCACGGTGGCTg < 2:342616/150‑1 (MQ=255) aTCGTACCCTGGTCGAAGGGCAGCGGGTGGAGTACGCCTGCGTGCAGGGCCAGAAAGGCCTGCAAGCCGAAGACGTAGTGGGCCTCTGAGCCTCTGGTTTCAAGCCGTAGGAGCGGGCTTGCCCGCGAAGCAGCCACCACGGTGGCTGcc > 1:42169/1‑150 (MQ=255) tACCCTGGTCGAAGGGCAGCGGGTGGAGTACGCCTGCGTGCAGGGCCAGAAAGGCCTGCAAGCCGAAGACGTAGTGGGCCTCTGAGCCTCTGGTTTCAAGCCGTAGGAGCGGGCTTGCCCGTG‑‑GCAGCCACCGTGGTGGCTGcttcgcgg < 2:379467/150‑8 (MQ=38) cccTGGTCGAAGGGCAGCGGGTGGAGTACGCCTGCGTGCAGGGCGAGAAAGGCCTGCAAGCCGAAGACGTAGTGGGCCTCTGAGCCTCTGGTTTCAAGCCGTAGGAGCGGGCTTGCCCGAGAAGCCGCCACCGTGGTGGCTGcttctcgg < 2:30271/150‑8 (MQ=255) ccTGGTCGAAGGGCAGCGGGTGGAGTACGCCTGCGTGCAGGGCCAGAAAGGCCTGCAAGCCGAAGACGTAGTGGGCCTCTGAGCCTCTGGTTTCAAGCCGTAGGAGCGGGCTTGCCCGCGAAGCAGCCACCACGGTGGCTGCCACGG‑GCa > 2:378875/1‑150 (MQ=255) ccTGGTCGAAGGGCAGCGGGTGGAGTACGCCTGCGTGCAGGGCCAGAAAGGCCTGCAAGCCGAAGACGTAGTGGGCCTCTGAGCCTCTGGTTTCAACTCCCATGGCTTGGCTTGCCCGCGAAGCAGCCACCACGGTGGCTGCCACGG‑GCa < 1:378875/150‑1 (MQ=255) ggAGTACGCCTGCGTGCAGGGCCAGAAAGGCCTGCAAGCCGAAGACGTAGTGGGCCTCTGAGCCTCTGGTTTCAAGCCGTAGGAGCGGGCTTGCCCGCGAAGCAGCCACCACGGTGGCTGCCACGG‑GCATAGCCCGTATTCGCGGGTaaa < 1:77913/150‑1 (MQ=255) gCCTGCGTGCAGGGCCAGAAAGGCCTGCAAGCCGAAGACGTAGTGGGCCTCTGAGCCTCTGGTTTCAAGCCGTAGGAGCGGGCTTGCCCGCGAAGCAGCCACCACGGTGGCTGCCACGG‑GCATAGCCCGTATTCGCGGGTaaa > 1:462192/1‑143 (MQ=255) gCCTGCGTGCAGGGCCAGAAAGGCCTGCAAGCCGAAGACGTAGTGGGCCTCTGAGCCTCTGGTTTCAAGCCGTAGGAGCGGGCTTGCCCGCGAAGCAGCCACCACGGTGGCTGCCACGG‑GCATAGCCCGTATTCGCGGGTaaa < 2:462192/143‑1 (MQ=255) gTGCAGGGCCAGAAAGGCCTGCAAGCCGAAGACGTAGTGGGCCTCTGAGCCTCTGGTTTCAAGCCGTAGGAGCGGGCTTGCCCGCGAAACAGCCACCACGGTGGCTGCCACGG‑GCATAGCCCGTATTCGCGGGTAAACCCATTCCTACaa > 2:478423/1‑149 (MQ=255) gCCAGAAAGGCCTGCAAGCCGAAGACGTAGTGGGCCTCTGAGCCTCTGGTTTCAAGCCGTAGGAGCGGGCTTGCCCGCGAAGCAGCCACCACGGTGGCTGCCACGG‑GCATAGCCCGTATTCGCGGGTAAACCCATTCCTACAGAAACGaa > 1:232830/1‑149 (MQ=255) aaaGGCCTGCAAGCCGAAGACGTAGTGGGCCTCTGAGCCTCTGGTTTCAAGCCGTAGGAGCGGGCTTGCCCGCGAGGCAGCCACCGTGGTGGCTGcttcgc > 2:260271/1‑96 (MQ=38) aaaGGCCTGCAAGCCGAAGACGTAGTGGGCCTCTGAGCCTCTGGTTTCAAGCCGTAGGAGCGGGCTTGCCCGCGAGGCAGCCACCGTGGTGGCTGcttcgc < 1:260271/101‑6 (MQ=38) ccGAAGACGTAGTGGGCCTCTGAGCCTCTGGTTTCAAGCCGTAGGAGCGGGCTTGCCCGCGAAGCAGCCACCACGGTGGCTGCCACGG‑GCATAGCCCGTATTCGCGGGTAAACCCATTCCTACAGAAACGAGCAAGCCTGCCGCTTcgtc < 2:192084/150‑4 (MQ=255) tctGAGCCTCTGGTTTCAAGCCGTAGGAGCGGGCTTGCCCGCGAAGCAGCCACCACGGTGGCTGCCACGG‑GCATAGCCCGTATTCGCGGGTAAACCCATTCCTACAGAAACGAGCAAGCCTGCCGTTTCATGATGTACGCCAGGTGAttt > 1:35875/1‑150 (MQ=255) gAGCCTCTGGTTTCAAGCCGTAGGAGCGGGCTTGCCCGCGAAGCAGCCACCACGGTGGCTGCCACGG‑GCATAGCCCGTATTCGCGGGTAAACCCATTCCTACAGAAACGAGCAAGCCTGCCGTTTCATGATGTACGCCAGGTGATTTCct < 2:218055/150‑1 (MQ=255) aGCCTCTGGTTTCAAGCCGTAGGAGCGGGCTTGCCCGCGAAGCAGCCACCACGGTGGCTGCCACGG‑GCATAGCCCGTATTCGCGGGTAAACCCATTCCTACAGAAACGAGCAAGCCTGCCGTTTCATGATGTACGCCAGGTGATTTCctc < 1:250078/150‑1 (MQ=255) tttCAAGCCGTAGGAGCGGGCTTGCCCGCGAAGCAGACACCACGGTGGATCCAACGGCGCA‑AGCCCGTATTCGCGGGTAAACCCATTCCTACAGAAACGAGCAAGCCTGCCGTTTCATGATGTACGCCAGGTGATTTCCTCTTCACCGTc > 1:60626/1‑150 (MQ=255) cAAGCCGTAGGAGCGGGCTTGCCCGCGAAGCAGCCAACACGGTGGCTGCCAGGG‑GCATAGCCCGTATTCGCGGGTAGGCCGATTCGTACAGAAACGAGCAAGCCTGCCGTTTCATGATGTACGCCAGGTGATTTCCTCTTCACCGTCGgc < 1:478423/150‑1 (MQ=255) gTAGGAGCGGGCTTGCCCGCGAAGCAGCCACCACGGTGGCTGCCACGG‑GCATAGCCCGTATTCGCGGGTAAACCCATTCCTACAGAAACGAGCAAGCCTGCCGTTTCATGATGTACGCCAGGTGATTTCCTCTTCACCGTCGGCGCTGAt > 1:395006/1‑150 (MQ=255) gAGCGGGCTTGCCCGCGAAGCAGCCACCACGGTGGCTGCCACGG‑GCATAGCCCGTATTCGCGGGTAAACCCATTCCTACAGAAACGAGCAAGCCTGCCGTTTCATGATGTACGCCAGGTGATTTCCTCTTCACCGTCGGCGCTGATGCgg < 2:42169/150‑1 (MQ=255) gCTTGCCCGCGAAGCAGCCACCACGGTGGCTGCCACGG‑GCATAGCCCGTATTCGCGGGTAAACCCATTCCTACAGAAACGAGCAAGCCTGCCGTTTCATGATGTACGCAAGAGGATTTCCTCTTCACCGTCGGCGCTGATGCGGACCCAg > 1:374397/1‑150 (MQ=255) gcagcCACCGTGGTGGCTGCCACGG‑GCATAGCCCGTATTCGCGGGTAAACCCATTCCTACAGAAACGAGCAAGCCTGCCGTTTCATGATGTACGCCAGGTGATTTCCTCTTCACCGTCGGCGCTGATGCGGACCCAGCGGTCAGCAtctt > 2:356452/1‑150 (MQ=255) ccaccaCGGTGGCTGCCACGG‑GACTACCTTGTACCGGCGGGTAAACCCATTCCTACAGAAACGAGCAAGCCTGCCGTTTCATGATGTACGCCAGGTGATTTCCTCTTCACCGTCGGCGCTGATGCGGACCCAGCGGTCAGCATCttcttc < 2:374397/150‑1 (MQ=255) accacGGTGGCTGCCACGG‑GCATAGCCCGTATTCGCGGGTAAACCCATTCCTACAGAAACGAGCAAGCCTGCCGTTTCATGATGTACGCCAGGTGATTTCCTCTTCACCGTCGGCGCTGATGCGGACCCAGCGGTCAGCATCTTCTTccc > 2:293222/1‑150 (MQ=255) | GCATCGTACCCTGGTCGAAGGGCAGCGGGTGGAGTACGCCTGCGTGCAGGGCCAGAAAGGCCTGCAAGCCGAAGACGTAGTGGGCCTCTGAGCCTCTGGTTTCAAGCCGTAGGAGCGGGCTTGCCCGCGAAGCAGCCACCACGGTGGCTGCCACGG‑GCATAGCCCGTATTCGCGGGTAAACCCATTCCTACAGAAACGAGCAAGCCTGCCGTTTCATGATGTACGCCAGGTGATTTCCTCTTCACCGTCGGCGCTGATGCGGACCCAGCGGTCAGCATCTTCTTCCC > NC_002947/1728558‑1728844 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 22 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |