Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | NC_002947 | 1,058,264:1 | +T | 13.5% | intergenic (‑190/+30) | PP_0914 ← / ← sodB | GGDEF domain‑containing protein/superoxide dismutase |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_002947 | 1,058,264 | 1 | . | T | 13.5% | 114.6 / 16.1 | 37 | intergenic (‑190/+30) | PP_0914/sodB | GGDEF domain‑containing protein/superoxide dismutase |
Reads supporting (aligned to +/- strand): ref base . (15/17); new base T (3/2); total (18/19) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 6.60e-01 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.96e-01 |
CAGTGTATCGGTTGCCTGACCCCGCAACTTGAATGCGCTGTATCAATCTTCCAGAGGTTTGATGTTAGCGTGGTATGCGGCCTGAGCAAGCAAAACCGTGCCAGCCTGTTTGCGGCAACACAGGCACAAAAAAACCCGGCACATG‑GCCGGGTTTCTCGTTCAAGGCCGAATGACTCAGGCCTTGAAGGTCTTGCCTTCGAACTGCTCGGCAACGAAGGCCCAGTTGACCAGGTTCCAGAACGCCTCGACGTACTTCGGACGCAGGTTGCGGTAGTCGATG > NC_002947/1058120‑1058399 | aaGTGTATCGGTTGCCTGACCACCCAACTTGAGTGCGCGGAATCACCGTGCCGGATGTTGGATGGTATCGTGGTATGCGGCCTGAGCAAGCAAAACCGTGCCAGCCTGTTTGCGGCAACACAGGCACAAAAAAACCCGGCACATG‑GCCgg < 2:571139/149‑1 (MQ=255) tCNGTTGCCTGACCCCGCAACTTGAATNCGCTGTATCAATNTTCCAGNGNTTTGATGTTAGCGTGGTAGGCGGCCTGAGCAAGCAAAACCGTGCCAGCCTGTTTGCGGCAACACAGGCACAAAAAAACCCGGCACATG‑GCCGGGTTtctc > 1:496237/1‑150 (MQ=255) tCGGTTGCCTGACCCCGCAACTTGAATGCGCTGTATCAATCTTCCAGAGGTTTGATGTTAGCGTGGTATGCGGCCTGAGCAAGCAAAACCGTGCCAGCCTGTTTGCGGCAACACAGGCACAAAAAAACCCGGCACATG‑GCCGGGTTtctc > 2:377915/1‑150 (MQ=255) gTTGCCTGACCCCGCAACTTGAATGCGCTGTATCAATCTTCCAGAGGTTTGATGTTAGCGTGGTATGCGGCCTGAGCAAGCAAAACCGTACCAGCCTGTTTGCGGCAACACAGGCACAAAAAAACCCGGCACATG‑GCCGGGTTTCTCGtt < 1:151631/150‑1 (MQ=255) gCCTGACCCCGCAACTTGAATGCGCTGTATCAATCTTCCAGAGGTATGATGTTAGCGTAGAATGCGGCTTAAGCAATCAAAACTGTGCCAGCCTGTTTGCGGCAAAACAAAAAAAAAAAAACCCGGCACATG‑GCCGGGTTTCTCGTTCaa < 2:180367/150‑1 (MQ=255) ccTGACCCCGCAACTTGAATGCGCTGTATCAATCTTCCAGAGGTTTGATGTTAGCGTGGTATGCGGCCTGAGCAAGCAAAACCGTGCCAGCCTGTTTGCGGCAACACAGGCACAAAAAAACCCGGC‑CATGTGCCGGGtt > 1:325946/1‑139 (MQ=255) ccTGACCCCGCAACTTGAATGCGCTGTATCAATCTTCCAGAGGTTTGATGTTAGCGTGGTATGCGGCCTGAGCAAGCAAAACCGTGCCAGCCTGTTTGCGGCAACACAGGCACAAAAAAACCCGGC‑CATGTGCCGGGtt < 2:325946/139‑1 (MQ=255) cTGACCCCGCAACTTGAATGCGCTGTATCAATCTTCCAGAGGTTTGATGTTAGCGTGGTATGCGGCCTGAGCAAGCAAAACCGTGCCAGCCTGTTTGCGGCAACACAGGCACAAAAAAACCCGGCACATG‑GCCGGGTTTCTCGTTCAAg < 1:541881/149‑1 (MQ=255) cTGACCCCGCAACTTGAATGCGCTGTATCAATCTTCCAGAGGTTTGATGTTAGCGTGGTATGCGGCCTGAGCAAGCAAAACCGTGCCAGCCTGTTTGCGGCAACACAGGCACAAAAAAACCCGGCACATG‑GCCGGGGTTCTCGTTCAAg > 2:541881/1‑149 (MQ=255) gCAACTTGAACGCGCTGTATCAATCTTCCAGAGGTTTGATGTTAGCGTGGTATGCGGCCTGAGCAAGCAAAACCGTGCCAGCCTGTTTGCGGCAACACAGGCACAAAAAAACCCGGCACATG‑GCCGGGTTTCTCGTTCAAGTCCGAATGa < 1:358455/150‑1 (MQ=255) cTTGAATGCGCGGTATCAATCTTCCAGAGGTTTGATGTTAGCGTGGTATGCGGCCTGAGCAAGCAAAACCGTGCCAGCCTGTTTGCGGCAACACAGGCACAAAAAAACCCGGCACATG‑GCCGGGTTTCTCGTTCAAGGCCGAATGACTCa < 2:496237/150‑1 (MQ=255) ttGAATGCGCTGTATCAATCTTCCAGAGGTTTGATGTTAGCGTGGTATGCGGCCTGAGCAAGCAAAACCGTGCCAGCCTGTTTGCGGCAACACAGGCACAAAAAAACCCGGCACATG‑GCCGGGTTTCTCGTTCAAGGCCGAATGACTCAg < 2:430808/150‑1 (MQ=255) gTATCAATCTTCCAGAGGTTTGATGTTAGCGTGGTATGCGGCCTGAGCAAGCAAAACCGTGCCAGCCTGTTTGCGGCAACACAGGCACAAAAAAACCCGGCACATG‑GCCGGGTTTCTCGTTCAAGGCCGAATGACTCAGGCCTTGAAGGt < 1:271685/150‑1 (MQ=255) tCAATCTTCCAGAGGTTTGATGTTAGCGTGGTATGCGGCCTGAGCAAGCAAAACCGTGCCAGCCTGTTTGCGGCAACACAGGCACAAAAAAACCCGGCACATG‑GCCGGGTTTCTCGTTCAAGGCCGAATGACTCAGGCCTTGAAGGTCtt < 1:377915/150‑1 (MQ=255) ccAGAGGTTTGATGTTAGCGTGGTATGCGGCCTGAGCAGGCCAAACCGTGACCGCATGTTTACGGAAACCCACAAACGCGAGAACCCGGCACATG‑GCCGGGTTTCTCGTTCAAGGCCGAATGACTCAGGCCTTGAAGGTCTTGCCTTCGa > 1:96870/1‑150 (MQ=255) agagGTTTGATGTTAGCGTGGTATGCGGCCTGAGCAAGCAAAACCGTGCCAGCCTGTTTGCGGCAACACAGGCACAAAAAAACCCGGCACATG‑GCCGGGttt > 1:594978/1‑102 (MQ=255) agagGTTTGATGTTAGCGTGGTATGCGGCCTGAGCAAGCAAAACCGTGCCAGCCTGTTTGCGGCAACACAGGCACAAAAAAACCCGGCACATG‑GCCGGGttt < 2:594978/102‑1 (MQ=255) agagGTTTGATGTTAGCGTGGTATGCGGCCTGAGCAAGCAAAACCGTGCCAGCCTGTTTGCGGCAACACAGGCACAAAAAAACCCGGCACATG‑GCCGGGTTTCTCGTTCAAGGCCGAATGACTCAGGCCTTGAAGGTCTTGCCTTCgttg > 2:91005/1‑147 (MQ=255) gagGTTTGATGTTAGCGTGGTATGCGGCCTGAGCAAGCAAAACCGTGCCAGCCTGTTTGCGGCAACACAGGCACAAAAAAACCCGGCACATG‑GCCGGGTTTCTCGTTCAAGGCCGAATGACTCAGGCCTTGAAGGTCTTg > 1:596545/1‑140 (MQ=255) gagGTTTGATGTTAGCGTGGTATGCGGCCTGAGCAAGCAAAACCGTGCCAGCCTGTTTGCGGCAACACAGGCACAAAAAAACCCGGCACATG‑GCCGGGTTTCTCGTTCAAGGCCGAATGACTCAGGCCTTGAAGGTCTTg < 2:596545/140‑1 (MQ=255) tttGATGTTAGCGTGGTATGCGGCCTAAGCAAGCAAAAACGTGCCAGCCTGTTTGCGGCAACACAGGCACAAAAAAACCCGGCACATG‑GCCGGGTTTCTCGTTCAAGGCCGAATGACTCAGGCCTTGAAGGTCTTGCCTTCGAACTGCTc > 1:36957/1‑150 (MQ=255) aTGTTAGCGTGGTATGCGGCCTGAGCAAGCAAAACCGTGCCAGCCTGTTTGCGGCAACACAGGCACAAAAAAACCCGGCACATG‑GCCGGGTTTCTCGTTCAAGGCCGAATGACTCAGGCCTTGAAGGTCTTGCCTTCGAACTGCTCGGCa > 1:368152/1‑150 (MQ=255) tGCGGCCTGAGCAAGCAAAACCGTGCCAGCCTGTTTGCGGCAACACAGGCACAAAAAAACCCGGCACATG‑GCCGGGTTTCTCGTTCAAGGCCGAATGACTCAGGCCTTGAAGGTCTTGCCTTCGAACTGCTCGGCAACGAAGGCCCAGtt > 1:720274/1‑150 (MQ=255) gCGGCCTGAGCAAGCAAAACCGTGCCAGCCTGTTTGCGGCAACACAGGCACAAAAAAACCCGGCACATG‑GCCGGGTTTCTCGTTCAAGGCCGAATGACTCAGGCCTTGAAGGTCTTGCCTTCGAACTGCTCGGCAACGAAGGCCCAGTTg > 2:94341/1‑150 (MQ=255) gcaagcaaAACCGTGCCAGCCTGTTTGCGGCAACACAGTCACAAAAGAACCCGGCACATT‑GCCGGGTTTTTCGTTAAAGACCGCAGGACTCAGGCCTTGAAGGTCTTGCCTTCGAACTGCTCGGCAACGAAGGACCAGGTGATCAGGTTc > 1:578175/1‑150 (MQ=255) caagcaaAACCGTGCCAGCCTGTTTGCGGCAACACAGGCACAAAAAAACCCGGCACATG‑GCCGGGTTTCTCGTTCAAGGCCGAATGACTCAGGCCTTGAAGGTCTTGCCTTCGAACTGCTCGGCAACGAAGGCCCAGTTGACCAGGTTcc < 2:368152/150‑1 (MQ=255) aagcaaAACCGTGCCAGCCTGTTTGCGGCAACACAGGCACAAAAAAACCCGGCACATG‑GCCGGGTTTCTCGTTCAAGGCCGAATGACTCAGGCCTTGAAGGTCTTGCCTTCGAACTGCTCGGCAACGAAGGCCCAGTTGACCAGGTTCCa > 1:40128/1‑150 (MQ=255) ggCAACACAGGCACAAAATAACACGGCACATG‑GCCGGGTTTCTCGTTCAAGGCCGAATGACTCAGGCCTGGAAGGTCTTGCCTTCGAACTGCTCGGCAACGAAGGCCCAGTTGACCAGGTTCCAGAACGCCTCGACGTACTTCGGACGCa < 2:96870/150‑1 (MQ=255) cAACACAGGCACAAAAAAACCCGGCACATG‑GCCGGGTTTCTCGTTCAAGGCCGAATGACTCAGGCCTTGAAGGTCTTGCCTTCGAACTGCTCGGCAACGAAGGCCCAGTTGACCAGGTTCCAGAACGCCTCGACGTACTTCGGACGCAgg < 2:40128/150‑1 (MQ=255) acaAAAAAACCCGGCACATG‑GCCGGGTTTCTCGTTCAAGGCCGAATGACTCAGGCCTTGAAGGTCtt < 1:602239/67‑1 (MQ=255) acaAAAAAACCCGGCACATG‑GCCGGGTTTCTCGTTCAAGGCCGAATGACTCAGGCCTTGAAGGTCtt > 2:602239/1‑67 (MQ=255) cgagaaaCCCGGC‑CATGTGCCGGGTTTCTCGTTCAAGGCCGAATGACTCAGGCCTTGAAGGTCTTGCCTTCGAACTGCTCGGCAACGAAGGCCCAGTTGACCAGGTTCCAGAACGCCTCGACGTACTTCGGACGCAGGTTGCGGTAGTCg > 2:598232/5‑150 (MQ=255) gagaaaCCCGGC‑CATGTGCCGGGTTTCTCGTTCAAGGCCGAATGACTCAGGCCTTGAAGGTCTTGCCTTCGAACTGCTCGGCAACGAAGGCCCAGTTGACCAGGTTCCAGAACGCCTCGACGTACTTCGGACGCAGGTTGCGGTAGTCGa > 1:390959/4‑150 (MQ=255) aaaaaCCCGGCACATG‑GCCGGGTTTCTCGTTCAAGGCCGAATGACTCAGGCCTTGAAGGTCTTGCCTTCGAACTGCTCGGCAACGAAGGCCCAGTTGACCAGGTTCCAGAACGCCTCGACGTACTTCGGACGCAGGTTGCGGTAGTCGAt < 2:36957/150‑1 (MQ=255) gaaaCCCGGC‑CATGTGCCGGGTTTCTCGTTCAAGGCCGAATGACTCAGGCCTTGAAGGTCTTGCCTTCGAACTGCTCGGCAACGAAGGCCCAGTTGACCAGGTTCCAGAACGCCTCGACGTACTTCGGACGCAGGTTGCGGTAGTCGATg < 1:598232/149‑1 (MQ=255) acaTG‑GCCGGGTTTCTCGTTCAAGGCCGAATGACTCAGGCCTTGAAGGTCTTGCCTTCGAACTGCTCGGCAACGAAGGCCCa > 2:621732/1‑82 (MQ=255) acaTG‑GCCGGGTTTCTCGTTCAAGGCCGAATGACTCAGGCCTTGAAGGTCTTGCCTTCGAACTGCTCGGCAACGAAGGCCCa < 1:621732/82‑1 (MQ=255) | CAGTGTATCGGTTGCCTGACCCCGCAACTTGAATGCGCTGTATCAATCTTCCAGAGGTTTGATGTTAGCGTGGTATGCGGCCTGAGCAAGCAAAACCGTGCCAGCCTGTTTGCGGCAACACAGGCACAAAAAAACCCGGCACATG‑GCCGGGTTTCTCGTTCAAGGCCGAATGACTCAGGCCTTGAAGGTCTTGCCTTCGAACTGCTCGGCAACGAAGGCCCAGTTGACCAGGTTCCAGAACGCCTCGACGTACTTCGGACGCAGGTTGCGGTAGTCGATG > NC_002947/1058120‑1058399 |
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: ‑ |