| Predicted mutation | |||||||
|---|---|---|---|---|---|---|---|
| evidence | seq id | position | mutation | freq | annotation | gene | description |
| RA | NC_002947 | 1,728,699 | C→T | 20.1% | 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 | 20.1% | 36.5 / 7.3 | 25 | intergenic (+51/+72) | cspA‑I/PP_1523 | major cold shock protein/hypothetical protein |
| Reads supporting (aligned to +/- strand): ref base C (10/10); new base T (2/3); total (12/13) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 1.00e+00 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.98e-01 | |||||||||||
GCATCGTACCCTGGTCGAAGGGCAGCGGGTGGAGTACGCCTGCGTGCAGGGCCAGAAAGGCCTGCAAGCCGAAGACGTAGTGGGCCTCTGAGCCTCTGGTTTCAAGCCGTAGGAGCGGGCTTGCCCGCGAAGCAGCCACCACGGTGGCTGCCACGG‑GCATAGCCCGTATTCGCGGGTAAACCCATTCCTACAGAAACGAGCAAGCCTGCCGTTTCATGATGTACGCCAGGTGATTTCCTCTTCACCGTCGGCGCTGATGCGGACCCAGCGGTCAGCATCTTCTTCCCCGTC > NC_002947/1728558‑1728848 | gCATCGTACCCTGGTCGAAGGGCAGCGGGTGGAGTACGCCTGCGTGCAGGGCCAGAAAGGCCTGCAAGCCGAAGACGTAGTGGGCCTCTGAGCCTCTGGTTTCAAGCCGTAGGAGCGGGCTTGCCCGCGAAGCAGCCACCACGGTGGCTg < 2:342616/150‑1 (MQ=255)gCATCGTACCCTGGTCGAAGGGCAGCGGGTGGAGTACGCCTGCGTGCAGGGCCAGAAAGGCCTGCAAGCCGAAGACGTAGTGGGCCTCTGAGCCTCTGGTTTCAAGCCGTAGGAGCGGGCTTGCCCGCGAAGCAGCCACCACGGTGGCTg > 2:268057/1‑150 (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 < 2:462192/143‑1 (MQ=255) gCCTGCGTGCAGGGCCAGAAAGGCCTGCAAGCCGAAGACGTAGTGGGCCTCTGAGCCTCTGGTTTCAAGCCGTAGGAGCGGGCTTGCCCGCGAAGCAGCCACCACGGTGGCTGCCACGG‑GCATAGCCCGTATTCGCGGGTaaa > 1:462192/1‑143 (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) cGGTGGCTGCCACGG‑GCATAGCCCGTATTCGCGGGTAAACCCATTCCTACAGAAACGAGCAAGCCTGCCGTTTCATGATGTACGCCAGGTGATTTCCTCTTCACCGTCGGCGCTGATGCGGACCCAGCGGTCAGCATCTTCTTCCCCGTc < 2:60626/150‑1 (MQ=255) | GCATCGTACCCTGGTCGAAGGGCAGCGGGTGGAGTACGCCTGCGTGCAGGGCCAGAAAGGCCTGCAAGCCGAAGACGTAGTGGGCCTCTGAGCCTCTGGTTTCAAGCCGTAGGAGCGGGCTTGCCCGCGAAGCAGCCACCACGGTGGCTGCCACGG‑GCATAGCCCGTATTCGCGGGTAAACCCATTCCTACAGAAACGAGCAAGCCTGCCGTTTCATGATGTACGCCAGGTGATTTCCTCTTCACCGTCGGCGCTGATGCGGACCCAGCGGTCAGCATCTTCTTCCCCGTC > NC_002947/1728558‑1728848 |
| 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: ‑ |