Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | NC_002947 | 5,458,946 | Δ1 bp | 17.2% | intergenic (+23/‑39) | leuS → / → PP_4795 | leucine‑‑tRNA ligase/lipoprotein |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_002947 | 5,458,946 | 0 | T | . | 17.2% | 75.2 / 13.0 | 29 | intergenic (+23/‑39) | leuS/PP_4795 | leucine‑‑tRNA ligase/lipoprotein |
Reads supporting (aligned to +/- strand): ref base T (14/10); new base . (2/3); total (16/13) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 6.32e-01 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.76e-01 |
CGAAGAAGTCGAAGCCGCTGCCCGCAGCAACGAGAATGTCCTGCGCTTCATCGATGGCCTGACCATCCGCAAGGTCATCGTGGTACCGGGCAAGCTGGTCAACATCGTCGCCAACTGATGAAAACGCCCACCCGC‑ATCGATTGCGGGTGGGCGGAACAGGCCCAAAAGGGAGCAACAACATGATCAAACGCAATTTGCTGGTAATCGGCCTGGCGGTCATGCTCAGCGCCTGCGGTTTCCAGCTGCGCGGCACCGGCGCCACCGAGCTGAGCGTGAAGGAAATGGAC > NC_002947/5458806‑5459092 | cgAAGAAGTCGAAGCCGCTGCCCGCAGCAACGAGAATGTCCTGCGCTTCATCGATGGCCTGACCATCCGCAAGGTCATCGTGGTACCGCGCAAGCTGGTCAACATCGTCGCCAACTGATGAAAACGCCCACCCGCAATCGA‑TGCGGGTgg < 2:310843/150‑1 (MQ=255) cgAAGAAGTCGAAGCCGCCGCCCGCAGCAACGAGAATGTCCTGCGCCAATTCCAGGGCCTGACCATCCGCAAGGTCATCGTGGTACCGGGCAAGCTGGTCAACATCGTCGCCGACTGATGAAAACGCCCACCCGC‑ATCGATTGCGGGTgg < 1:499732/150‑1 (MQ=255) gAAGCCGCTGCCCGCAGCAACGAGAATGTCCTGCGCTTCATCGATGGCCTGACCATCCGCAAGGTCATCGTGGTACCGGGCAAGCTGGTCAACATCGTCGCCAACTGATGAAAACGCCCACCCGC‑ATCGATTGCGGGTGGGCGGAACAgg > 2:113843/1‑150 (MQ=255) gAAGCCGCTGCCCGCAGCAACGAGAATGTCCTGCGCTTCATCGATGGCCTGACCATCCGCAAGGTCATCGTGGTACCGGGCAAGCTGGTCAACATCGTCGCCAACTGATGAAAACGCCCACCCGC‑ATCGATTGCGGGTGGGCGGAACAgg < 1:593065/150‑1 (MQ=255) cTGCCCGCAGCAACGAGAATGTCCTGCGCTTCATCGATGGCCTGACCATCCGCAAAGTCATCGTGGTACCGGGCAAGCTGGTCAACATCGTCGCCAACTGATGAAAACGCCCACCCGC‑ATCGATTGCGGGTGGGCGGAACAGGCCCaaaa < 1:40208/150‑1 (MQ=255) gagaATGTCCTGCGCTTCATCGATGGCCTGACCATCCGCAAGGTCATCGTGGTACCGGGCAAGCTGGTCAACATCGTCGCCAACTGATGAAAACGCCCACCCGC‑ATCGATTGCGGGTGGGCGGAACAGGCCCAAAAGGGAGCAACAACAt > 1:499148/1‑150 (MQ=255) ttCATCGATGGCCTGACCATCCGCAAGGTCATCGTGGTACCGGGCAAGCTGGTCAACATCGTCGCCAACTGATGAAAACGCCCACCCGCAATCGA‑TGCGGGTggg > 2:63658/1‑105 (MQ=39) ttCATCGATGGCCTGACCATCCGCAAGGTCATCGTGGTACCGGGCAAGCTGGTCAACATCGTCGCCAACTGATGAAAACGCCCACCCGCAATCGA‑TGCGGGTggg < 1:63658/105‑1 (MQ=39) tCGATGGCCTGACCATCCGCAAGGTCATCGTGGTACCGGGCAAGCTGGTCAACATCGTCGCCAACTGATGAAAACGCCCACCCGC‑ATCGATTGCGGGTGGGCGGAACAGGCCCAAAAGGGAGCAACAACATGATCAAACGCAATTTGCTg < 2:499148/150‑1 (MQ=255) ccTGACCATCCGCAAGGTCATCGTGGTACCGGGCAAGCTGGTCAACATCGTCGCCAACTGATGAAAACGCCCACCCGC‑ATCGATTGCGGGTGGGCGGAACAGGCCCAAAAGGGAGCAACAACATGATCAAACGCAATTTGCTGGTAATCg > 2:208207/1‑150 (MQ=255) cTGACCATCCGCAAGGTCATCGTGGTACAGGGCAAGCTGGTCAACATCGTCGCCAACTGATGAAAACGCCCACCCGC‑ATCGATTGCGGGTGGGCGGAACAGGCCCAAAAGGGAGCAACAACATGATCAAACGCAATTTGCTGGTAATCgg > 2:198677/1‑150 (MQ=255) tGACCATCCGCAAGGTCATCGTGGTACCGGGCAAGCTGGTCAACATCGTCGCCAACTGATGAAAACGCCCACCCGCAATCGA‑TGCGGGTgg > 2:479859/1‑91 (MQ=39) tGACCATCCGCAAGGTCATCGTGGTACCGGGCAAGCTGGTCAACATCGTCGCCAACTGATGAAAACGCCCACCCGCAATCGA‑TGCGGGTgg < 1:479859/91‑1 (MQ=39) aaGGTCATCGTGGTACCGGGCAAGCTGGTCAACATCGTCGCCAACTGATGAAAACGCCCACCCGC‑ATCGATTGCGGGTGGGCg < 2:280008/83‑1 (MQ=255) aaGGTCATCGTGGTACCGGGCAAGCTGGTCAACATCGTCGCCAACTGATGAAAACGCCCACCCGC‑ATCGATTGCGGGTGGGCg > 1:280008/1‑83 (MQ=255) cATCGTGGTACCGGGCAAGCTGGTCAACATCGTCGCCAACTGATGAAAACGCCCACCCGC‑ATCGATTGCGGGTGGGCGGAACAGGCCCAAAAGGGAGCAACAACATGATCAAACGCAATTTGCTGGTAATCGGCCTGGCGGTCATGCTCa < 1:113843/150‑1 (MQ=255) aTCGTGGTACCGGGCAAGCTGGTCAACATCGTCGCCAACTGATGAAAACGCCCACCCGC‑ATCGATTGCGGGTGGGCGGaa < 2:260856/80‑1 (MQ=255) aTCGTGGTACCGGGCAAGCTGGTCAACATCGTCGCCAACTGATGAAAACGCCCACCCGC‑ATCGATTGCGGGTGGGCGGaa > 1:260856/1‑80 (MQ=255) tGGTACCGGGCAAGCTGGTCAACATCGTCGCCAACTGATGAAAACGCCCACCCGC‑ATCGATTGCGGGTGGGCGGAACAGGCCCAAAAGGGAGCAACAACATGATCAAACGCAATTTGCTGGTAATCGGCCTGGCGGTCATGCTCAGCGcc > 1:348930/1‑150 (MQ=255) ccGGGCAAGCTGGTCAACATCGTCGCCAACTGATGAAAACGCCCACCCGC‑ATCGATTGCGGGTGGGCGGAACAGGCCCAAAAGGGAGCAACAACATGATCAAACGCAATTTGCTGGTAATCGGCCTGGCGGTCATGCTCAGCGCCTGCgg > 1:599707/1‑150 (MQ=255) ccGGGCAAGCTGGTCAACATCGTCGCCAACTGATGAAAACGCCCACCCGC‑ATCGATTGCGGGTGGGCGGAACAGGCCCAAAAGGGAGCAACAACATGATCAAACGCAATTTGCTGGTAATCGGCCTGGCGGTCATGCTCAGCGCCTGCgg > 2:505625/1‑150 (MQ=255) gggCAAGCTGGTCAACATCGTCGCCAACTGATGAAAACGCCCACCCGC‑ATCGATTGCGGGTGGGCGGAACAGGCCCAAAAGGGAGCAACAACATGATCAAACGCAATTTGCTGGTAATCGGCCTGGCGGTCATGCTCAGCGCCTGCGGtt > 2:166544/1‑150 (MQ=255) gCTGGTCAACATCGTCGCCAACTGATGAAAACGCCCACCCGC‑ATCGATTGCGGGTGGGCGGAACAGGCCCAAAAGGGAGCAACAACATGATCAAACGCAATTTGCTGGTAATCGGCCTGGCGGTc < 1:337968/125‑1 (MQ=255) gCTGGTCAACATCGTCGCCAACTGATGAAAACGCCCACCCGC‑ATCGATTGCGGGTGGGCGGAACAGGCCCAAAAGGGAGCAACAACATGATCAAACGCAATTTGCTGGTAATCGGCCTGGCGGTc > 2:337968/1‑125 (MQ=255) tCAACATCGTCGCCAACTGATGAAAACGCCCACCCGC‑ATCGATTGCGGGTGGGCGGAACAGGCCCAAAAGGGAGCAACAACATGATCAAACGCAATTTGCTGGTAATCGGCCTGGCGGTCATGCTCAGCGCCTGCGGTTTCCAGCTgcgc > 1:81698/1‑150 (MQ=255) aaCATCGTCGCCAACTGATGAAAACGCCCACCCGC‑ATCGATTGCGGGTGGGCGGAACAGGCCCAAAAGGGAGCAACAACATGATCAAACGCAATTTGCTGGTAATCGGCCTGGCGGTCATGCTCAGCTCCTGCGGTTTCCAGCTGCGCgg > 2:320987/1‑150 (MQ=255) gaAAACGCCCACCCGC‑ATCGATTGCGGGTGGGCGGAACAGGCCCAAAAGGGAGCAACAACATGATCAAACGCAATTTGCTGGTAATCGGCCTGGCGGTCATGCTCAGCGCCTGCGGTTTCCAGCTGCGCGGCACCGGCGCCACCGAGCTg > 1:468237/1‑150 (MQ=255) cACCCGC‑ATCGATTGCGGGTGGGCGGAACAGGCCCAAAAGGGAGCAACAACATGATCAAACGCAATTTGCTGGTAATCGGCCTGGCGGTCATGCTCAGCGCCTGCGGTTTCCAGCTGCGCGGCACCGGCGCCACCGAGCTGAGCGTGAAg < 1:208207/150‑1 (MQ=255) cGATTGCGGGTGGGCGGAACAGGCCCAAAAGGGAGCAACAACATGATCAAACGCAATTTGCTGGTAATCGGCCTGGCGGTCATGCTCAGCGCCTGCGGTTTCCAGCTGCGCGGCACCGGCGCCACCGAGCTGAGCGTGAAGGAAATGGAc < 2:599707/150‑1 (MQ=255) | CGAAGAAGTCGAAGCCGCTGCCCGCAGCAACGAGAATGTCCTGCGCTTCATCGATGGCCTGACCATCCGCAAGGTCATCGTGGTACCGGGCAAGCTGGTCAACATCGTCGCCAACTGATGAAAACGCCCACCCGC‑ATCGATTGCGGGTGGGCGGAACAGGCCCAAAAGGGAGCAACAACATGATCAAACGCAATTTGCTGGTAATCGGCCTGGCGGTCATGCTCAGCGCCTGCGGTTTCCAGCTGCGCGGCACCGGCGCCACCGAGCTGAGCGTGAAGGAAATGGAC > NC_002947/5458806‑5459092 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |