Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | CP006881 | 924,066:1 | +A | 100% | intergenic (‑46/‑233) | U712_04490 ← / → U712_04495 | Vegetative catalase/Aliphatic sulfonates import ATP‑binding protein SsuB |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | CP006881 | 924,066 | 1 | . | A | 82.4% | 48.0 / 3.3 | 17 | intergenic (‑46/‑233) | U712_04490/U712_04495 | Vegetative catalase/Aliphatic sulfonates import ATP‑binding protein SsuB |
Reads supporting (aligned to +/- strand): ref base . (2/1); new base A (1/13); total (3/14) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 6.32e-02 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 | |||||||||||
Rejected as polymorphism: Variant not supported by required number of reads on each strand. | |||||||||||
Rejected as polymorphism: Polymorphic indel expands or contracts a homopolymer stretch. |
GTCCGCGAGAACCGGCAGTCATTGAGTTTTGATTATCTCCAACCGGAGCGCCCCAGCTAGTTGTCAGTTTATTTGAACTCATGTTATCACCTCTTGGAATTTATAAGAACATAATATCATCTCTAA‑GT‑AAAA‑AGTCAATATTATTTTATAATTATTATAAAATAGTTATCTCTATTAATCGATATATAGTATCGAAATTGTTGAAAAACCTTGTTGTAACAGCTTTTTATGACACTT > CP006881/923939‑924175 | gTCCGCGAGAACCGGCAGTCATTGAGTTTTGATTATCTCCAACCGGAGCGCCCCAGCTAGTTGTCAGTTTATTTGAACTCATGTTATCACCGCTTGGAATTTATAAGAACATAATTTCATCTCTAA‑GT‑aaaa‑aagtcac > 5:35452/1‑133 (MQ=255) ccGCGAGAACCGGCAGTCATTGAGTTTTGATTATCTCCAACCGGAGCGCCCCAGCTAGTTGTCAGTTTATTTGAACTCATGTTATCACCTCTTGGAATTTATAAGAACATAATATCATCTCTAA‑GTAAAAA‑AGTCAata < 7:83553/139‑1 (MQ=255) gCAGTCATTGAGTTTTGATTATCTCCAACCGGAGCGCCCCAGCTAGTTGTCAGTTTATTTGAACTCATGTTATCACCTCTTGGAATTTATAAGAACATAATATCATCTCTAA‑GTAAAAA‑AGTCAATATTATTTTATAat < 4:135819/139‑1 (MQ=255) gCAGTCATTGAGTTTTGATTATCTCCAACCGGAGCGCCCCAGCTAGTTGTCAGTTTATTTGAACTCATGTTATCACCTCTTGGAATTTATAAGAACATAATATCATCTCTAA‑GTAAAAA‑AGTCAATATTATTTTATAat < 8:133346/139‑1 (MQ=255) ttGAGTTTTGATTATCTCCAACCGGAGCGCCCCAGCTAGTTGTCAGTTTATTTGAACTCATGTTATCACCTCTTGGAATTTATAAGAACATAATATCATCTCTAA‑GTAAAAA‑AGTCAATATTATTTTATAATTATtata < 1:39571/139‑1 (MQ=255) aTTATCTCCAACCGGAGCGCCCCAGCTAGTTGTCAGTTTATTTGAACTCATGTTATCACCTCTTGGAATTTATAAGAACATAATATCATCTCTAA‑GTAAAAA‑AGTCAATATTATTTTATAATTATTATAAAATAGTTAt < 3:115809/139‑1 (MQ=255) aTTATCTCCAACCGGAGCGCCCCAGCTAGTTGTCAGTTTATTTGAACTCATGTTATCACCTCTTGGAATTTATAAGAACATAATATCATCTCTAA‑GTAAAAA‑AGTCAATATTATTTTATAATTATTATAAAATAGTTAt < 5:127626/139‑1 (MQ=255) aCCGGAGCGCCCCAGCTAGTTGTCAGTTTATTTGAACTCATGTTATCACCTCTTGGAATTTATAAGAACATAATATCATCTCTAA‑GTAAAAA‑AGTCAATATTATTTTATAATTATTATAAAATAGTTATCTCTATTAAt < 2:283818/139‑1 (MQ=255) aCCGGAGCGCCCCAGCTAGTTGTCAGTTTATTTGAACTCATGTTATCACCTCTTGGAATTTATAAGAACATAATATCATCTCTAA‑GTAAAAA‑AGTCAATATTATTTTATAATTATTATAAAATAGTTATCTCTATTAAt < 6:253875/139‑1 (MQ=255) cGGAGCGCCCCAGCTAGTTGTCAGTTTATTTGAACTCATGTTATCACCTCTTGGAATTTATAAGAACATAATATCATCTCTAA‑GTAAAAA‑AGTCAATATTATTTTATAATTATTATAAAATAGTTATCTCTATTAATCg < 3:143698/139‑1 (MQ=255) gcgcCCCAGCTAGTTGTCAGTTTATTTGAACTCATGTTATCACCTCTTGGAATTTATAAGAACATAATATCATCTCTAA‑GT‑AAAAGAGTCAATATTATTTTATAATTATTATAAAATAGTTATCTCTATTaa > 5:116099/1‑132 (MQ=255) gcgcCCCAGCTAGTTGTCAGTTTATTTGAACTCATGTTATCACCTCTTGGAATTTATAAGAACATAATATCATCTCTAA‑GT‑AAAACAGTCAATATTATTTTATAATTATTATAAAATAGTTATCTCTATTaa < 6:116099/132‑1 (MQ=255) cgcCCCAGCTAGTTGTCAGTTTATTTGAACTCATGTTATCACCTCTTGGAATTTATAAGAACATAATATCATCTCTAA‑GTAAAAA‑AGTCAATATTATTTTATAATTATTATAAAATAGTTATCTCTATTAATCGatata < 5:218372/139‑1 (MQ=255) cccAGCTAGTTGTCAGTTTATTTGAACTCATGTTATCACCTCTTGGAATTTATAAGAACATAATATCATCTCTAA‑GTAAAAA‑AGTCAATATTATTTTATAATTATTATAAAATAGTTATCTCTATTAATCGATATATAg < 4:233246/139‑1 (MQ=255) gTTGTAAGTTTATTTGAACTCATGTTATCACCTCTTGGAATTTATAAGAACATAATATCATCTCTAA‑GTAAAAA‑AGTCAATATTATTTTATAATTATTATAAAATAGTTATCTCTATTAATCGATATATAGTATCGa < 1:127289/137‑1 (MQ=255) ttGTCAGTTTATTTGAACTCATGTTATCACATCTTGGAATTTATAAGAACATAATATAAACTATAATAA‑AAAA‑AGTCAATATTATTTTATAATTATTATAAAATAGTTATCTCTAATAATCTATATATAGTATCTAAAt > 8:252635/1‑139 (MQ=255) tttATTTGAACTCATGTTATCTCCTCTTGTAATTTATAAGAACATAATATCATCTCTAA‑GTAAAAA‑AGTCAATATTATTTTATAATTATTATAAAATAGTTATCTCTATTAATCGATATATAGTATAGAAATTGTTGaa < 7:252635/139‑1 (MQ=255) ttATAAGAACATAATATCATCTCTAA‑GTAAAAA‑AGTCAATATTATTTTATAATTATTATAAAATAGTTATCTCTATTAATCGATATATAGTATCGAAATTGTTGAAAAACCTTGTTGTAACAGCTTTTTATGACACtt > 8:77410/1‑138 (MQ=255) | GTCCGCGAGAACCGGCAGTCATTGAGTTTTGATTATCTCCAACCGGAGCGCCCCAGCTAGTTGTCAGTTTATTTGAACTCATGTTATCACCTCTTGGAATTTATAAGAACATAATATCATCTCTAA‑GT‑AAAA‑AGTCAATATTATTTTATAATTATTATAAAATAGTTATCTCTATTAATCGATATATAGTATCGAAATTGTTGAAAAACCTTGTTGTAACAGCTTTTTATGACACTT > CP006881/923939‑924175 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |