Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | NC_002947 | 4,586,056:1 | +C | 100% | intergenic (+139/+76) | PP_4061 → / ← PP_4063 | hypothetical protein/long‑chain fatty acid‑‑CoA ligase |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_002947 | 4,586,056 | 1 | . | C | 96.4% | 87.4 / ‑2.9 | 28 | intergenic (+139/+76) | PP_4061/PP_4063 | hypothetical protein/long‑chain fatty acid‑‑CoA ligase |
Reads supporting (aligned to +/- strand): ref base . (0/1); new base C (14/13); total (14/14) | |||||||||||
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.78e-01 | |||||||||||
Rejected as polymorphism: E-value score below prediction cutoff. | |||||||||||
Rejected as polymorphism: Frequency below/above cutoff threshold. | |||||||||||
Rejected as polymorphism: Variant not supported by required number of reads on each strand. |
AGGCTTAAGGGTAGACCACGGCACCGGCGTGAACAGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTCTGC‑GGTAAACCCGCTCCTACCAGGGA‑C‑GTGCCGGCCTTGAATG‑CGGCGCTGCCCCTGGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCCGCGCCAGCACATCAACCCGCAGATGCAG‑CCGCAATCTCCGCCACACTGATCTCACGCATGCGGAACTTCTGCACCTTGC > NC_002947/4585910‑4586199 | aGGCTTAAGGGTAGACCACGGCACCGGCGTGAACAGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACCAGGGa < 4:174316/148‑1 (MQ=255) ggTAGACCACGGCACCGGCGTGAACAGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCGCGGGTAAACACGCTCCTACCAGGGACC‑GTGCCgg > 4:23980/1‑148 (MQ=39) cACGGCACCGGCGTGAACAGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACCAGGGACC‑GTGCCGGCCTTGAAt > 2:100731/1‑149 (MQ=39) cACGGCACCGGCGTGAACAGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTAATCGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACCAGGGACC‑GTGCCGGCCTTGAAt > 5:73495/1‑149 (MQ=38) cGGCACCGGCGTGAACAGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACCAGGGACC‑GTGCCGGCCTTGAATg‑c < 1:50803/149‑1 (MQ=39) ggCACCGGCGTGAACAGCTGATCTGCGCTGCTCCATTGAAAGGGTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTTATCGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACCAGGGACC‑GTGCCGGCCTTGAATg‑cg < 6:73495/149‑1 (MQ=38) aaCAGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACCAGGGACC‑GTGCCGGCCTTGAATg‑cg > 6:321739/1‑137 (MQ=38) aaCAGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACCAGGGACC‑GTGCCGGCCTTGAATg‑cg < 5:321739/137‑1 (MQ=38) cAGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACCAGGGACC‑GTGCCGGCCTTGAATG‑CGGCGCTGCCCCTGGt < 4:344126/149‑1 (MQ=39) aGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACCAGGGACC‑GTGCCGGCCTTGAATG‑CGGCGCTGCCCCTGGTa < 1:100731/149‑1 (MQ=39) cTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACCAGGGACC‑GTGCCGGCCTTGAATG‑CGGCGCTGCCCCTGGTAgg < 3:23980/149‑1 (MQ=39) cgcTGCTCCATTGAAGGGTTCAGTATCACCGGGGGCTTGCCACGGTGACTATCTTGCATTCAGCGGGCCGGCCTCTGCGCGGGTAACCCCGCTCCTACCAGGGACC‑GTGCCGGCCTTGAATG‑CGGCGCTGCCCCTGGTAGGAGCgg < 1:277574/146‑1 (MQ=255) cgcTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGAAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACCAGGGG‑CAGCGCC‑GCATTCAAGGCCGGCACGGTCCCTGGTAGGAGCg < 1:307536/145‑1 (MQ=9) ctCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACCAGGGACC‑GTGCCGGCCTTGAATG‑CGGCGCTGCCCCTGGTAGGAGCGGGTTTAc < 4:158241/147‑1 (MQ=38) ttGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACCAGGGACC‑GTGCCGGCCTTGAATG‑CGGCGCTGCCCCTgg > 4:335159/1‑119 (MQ=37) ttGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACCAGGGACC‑GTGCCGGCCTTGAATG‑CGGCGCTGCCCCTgg < 3:335159/119‑1 (MQ=37) ttGCATTGATCGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACCAGGGACC‑GTGCCGGCCTTGAATG‑CGGCGCTGCCCCTGGTAGGa < 1:13341/89‑1 (MQ=21) ttGCATTGATCGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACCAGGGACC‑GTGCCGGCCTTGAATG‑CGGCGCTGCCCCTGGTAGGa > 2:13341/1‑89 (MQ=21) gcgggTAAACCCGCTCCTACCAGGGACC‑GTGCCGGCCTTGAATG‑CGGCGCTGCCCCTGGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCCGCGCCAGCACATCAACCCGCAGATGCAG‑CCGCAATCTCCGCCACACTGATCTCACGCa > 2:249248/4‑149 (MQ=255) tAAACCCGCTCCTACCAGGGACC‑GTGCCGGCCTTGAATG‑CGGCGCTGCCCCTGGTAGGAGCGGGTGTACCCGCGAAGAGGCCGGCCGCGCCAGCACATCAACCCGCAGATGCAG‑CCGCAATCTCCGCCACACTGATCTCACGCATGCg > 1:342344/1‑148 (MQ=255) ccGCTCCTACCAGGGACC‑GTGCCGGCCTTGAATG‑CGGCGCTGCCCCTGGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCCGCGCCAGCACATCAACCCGCAg < 5:164089/103‑1 (MQ=255) ccGCTCCTACCAGGGACC‑GTGCCGGCCTTGAATG‑CGGCGCTGCCCCTGGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCCGCGCCAGCACATCAACCCGCAg > 6:164089/1‑103 (MQ=255) tCCTACCAGGGACC‑GTGCCGGCCTTGAATG‑CGGCGCTGCCCCTGGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCCGCGCCAGCACATCAACCCGCAGATGCAG‑CCGCAAt < 5:323892/112‑1 (MQ=255) tCCTACCAGGGACC‑GTGCCGGCCTTGAATG‑CGGCGCTGCCCCTGGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCCGCGCCAGCACATCAACCCGCAGATGCAG‑CCGCAAt > 6:323892/1‑112 (MQ=255) cTACCAGGGACC‑GTGCCGGCCTTGAATG‑CGGCGCTGCCCCTGGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCCGCGCCAGCACATCAACCCGCAGATGCAGCCCGCAATCTCCGCCACACTGATCTCACGCATGCGGAACTTCTg < 1:249248/147‑1 (MQ=255) tACCAGGGACC‑GTGCCGGCCTTGAATG‑CGGCGCTGCCCCTGGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCCGCGCCAGCACATCAACCCGCAGATGCAG‑CCGCAATCTCCGCCACACTGATCTCACGCATGCGGAACTTCTGCAcc > 3:205337/1‑149 (MQ=255) tACCAGGGACC‑GTGCCGGCCTTGAATG‑CGGCGCTGCCCCTGGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCCGCGCCAGCACATCAACCCGCAGATGCAG‑CCGCAATCTCCGCCACACTGATCTCACGCATGCGGAACTTCTGCAcc > 1:156956/1‑149 (MQ=255) ccAGGGACC‑GTGCCGGCCTTGAATG‑CGGCGCTGCCCCTGGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCCGCGCCAGCACATCAACCCGCAGATGCAG‑CCGCAATCTCCGCCACACTGATCTCACGCATGCGGAACTTCTGCACCtt > 6:103511/1‑149 (MQ=255) aGGGACC‑GTGCCGGCCTTGAATG‑CGGCGCTGCCCCTGGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCCGCGCCAGCACATCAACCCGCAGATGCAG‑CCGCAATCTCCGCCACACTGATCTCACGCATGCGGAACTTCTGCACCTTGc > 8:72279/1‑149 (MQ=255) | AGGCTTAAGGGTAGACCACGGCACCGGCGTGAACAGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTCTGC‑GGTAAACCCGCTCCTACCAGGGA‑C‑GTGCCGGCCTTGAATG‑CGGCGCTGCCCCTGGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCCGCGCCAGCACATCAACCCGCAGATGCAG‑CCGCAATCTCCGCCACACTGATCTCACGCATGCGGAACTTCTGCACCTTGC > NC_002947/4585910‑4586199 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |