Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | NC_002947 | 4,950,972 | T→G | 20.3% | intergenic (‑122/+35) | fliL ← / ← fliK | flagellar protein FliL/flagellar hook‑length control protein FliK |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_002947 | 4,950,972 | 0 | T | G | 20.3% | 58.9 / 10.6 | 36 | intergenic (‑122/+35) | fliL/fliK | flagellar protein FliL/flagellar hook‑length control protein FliK |
Reads supporting (aligned to +/- strand): ref base T (12/16); new base G (2/5); total (14/22) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 6.76e-01 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 |
CTTCGCTCTTCGCCATGCCAATAATCCGTCGTCCATCGGGGGTTCAAAGGAGAATGGCTTGGCTTGAGCAAGTGTTATGCCAGATTTTGCA‑‑GGGGATGGATTCCAGCAGGTTTTGCGGTGCTGGCTTCTTCGCG‑GGCTCGCCCGCGAAGAAGCCAGCCCAGGCCAACAGTTATCAGGCGTAGTAATCGACCGCACTGTCGCCGACAACCACCTGCTGCTCCAAAGGCCGGGCGCTATCGACCGCTTCGGCCCCACCCTGCTCGGCCCGGCGCGCCGCCACCC > NC_002947/4950835‑4951116 | cTTCGCTCTTCGCCATGCCAATAATCCGTCGTCCATCGGGGGTTCAAAGGAGAATGGCTTGGCTTGAGCAAGTGTTATGCCAGATTTTGCA‑‑GGGGATGGATTCCAGCAGGTTTTGCGGTGCTGGCTTCTTCGCG‑GGCTCGCCCGCga < 2:877/147‑1 (MQ=255) cGCTCTTTGCCATGCCAATAATCCGTGGGCCCTCGGGGGTTCAAATGATAATGGCTTGGCTTGAGCAAGTGTTATGCCAGATTTTGCATCGGGGATGGATTCCAGCAGGTTTTGCGGTGCTGGCTTCTTCGCG‑GGCGAGCCCGCgaag < 1:125085/148‑1 (MQ=255) cGCTCTTCGCCCTGCCAATAATCCGTCGTCCATCGGGGGTTCAAAGGAGAATGGCTTGGCTTGAGCAAGTGTTATGCCAGATTTTGCA‑‑GGGGATGGATTCCAGCAGGTTTTGCGGTGCTGGCTTCTTCGCG‑GGCTCGCCCGCgaagaa < 3:68761/148‑1 (MQ=255) cGCTCTTCGCCATGCCAATAATCCGTCGTCCATCGGGGGTTCAAAGGAGAATGGCTTTGCTTGAGCAAGTGTTATGCCAGATTTTGCA‑‑GGGGATGGATTCCAGCAGGTTTTGCGGTGCTGGCTTCTTCGCG‑GGCGAGCCCGCgaagaa < 2:415751/148‑1 (MQ=255) cGCTCTTCGCCATGCCAATAATCCGTCGTCCATCGGGGGTTCAAAGGAGAATGGCTTGGCTTGAGCAAGTGTTATGCCAGATTTTGCA‑‑GGGGATGGATTCCAGCAGGTTTTGCGGTGCTGGCTTCTTCGCG‑GGCTCGCCCGCgaagaa < 6:234739/148‑1 (MQ=255) cGCTCTTCGCCATGCCAATAATCCGGCGACCATCGGGGGTTCACAGGAGAATGGCTTGGCTTGAGCACGGGTTATGCCAGAGATTGCA‑‑GGCGAGGGATTCCAGCAGGTTTTGCGGTGCTGGCTTCTTCGCG‑GGCTCGCCCGCtaaga < 5:303637/147‑1 (MQ=255) gCTCATCGCCATGCCAATAATCCGTCGTCCATCGGGGGTTCAAAGGAGAATGGCTTGGCTTGAGCAAGTGTTATGCCAGATTTTGCA‑‑GGGGATGGATTCCAGCAGGTTTTGCGGTGCTGGCTTCTTCGCG‑GGCCCGCCCGCGaagaag < 7:368997/148‑1 (MQ=255) ttCGCCATGCCAATAATCCGTCGTCCATCGGGGGTTCAAAGGAGAATGGCTTGGCTTGAGCAAGTGTTATTCCAGCTCTTGCA‑‑GGGGATGGATTCCAGCAGGTTTTGCGGTGCTAGCTTCTTCGCG‑GGCGAGCCCGCGAAGAagcca < 1:354831/147‑1 (MQ=255) aataatCCGTCGTCCATCGGGGGTTCAAAGGAGAATGGCTTGGCTTGAGCAAGTGTTATGCCAGATTTTGCA‑‑GGGGATGGATTCCAGCAGGTTTTGCGGTGCTGGCTTCTTCGCG‑GGCTCGCCCGCGAAGAAGCCCGCCCAGGCCAAcc > 7:262398/1‑148 (MQ=255) aatCCGTCGTCCATCGGGGGTTCAAAGGAGAATGGCTTGGCTTGAGCAAGTGTTATGCCAGATTTTGCA‑‑GGGGATGGATTCCAGCAGGTTTTGCGGTGCTGGCTTCTGCGCG‑GGCGCGCCCGCGAAGAAGCCAGCCCAGGCCAACAGtt < 1:408079/149‑1 (MQ=255) cgtcgtCCATCGGGGGTTCAAAGGAGAATGGCTTGGCTTGAGCAAGTGTTATGCCAGATTTTGCA‑‑GGGGATGGATTCCAGCAGGTTTTGCGGTGCTGGCTTCTTCGCG‑GGCTCGCCCGCGAAGAAGCCAGCCCAGGCCAAc < 8:317165/141‑1 (MQ=255) cgtcgtCCATCGGGGGTTCAAAGGAGAATGGCTTGGCTTGAGCAAGTGTTATGCCAGATTTTGCA‑‑GGGGATGGATTCCAGCAGGTTTTGCGGTGCTGGCTTCTTCGCG‑GGCTCGCCCGCGAAGAAGCCAGCCCAGGCCAAc > 7:317165/1‑141 (MQ=255) gtcgtcCATCGGGGGTTCAAAGGAGAATGGCTTGGCTTGAGCAAGTGTTATGCCAGATTTTGCA‑‑GGGGATGGATTCCAGCAGGTTTTGCGGTGCTGGCTTCTTCGCG‑GGCTCGCCCGCGAAGAAGCCAGCCCAGGCCAACAGTTATCAg < 3:411941/149‑1 (MQ=255) cATCGGGGGTTCAAAGGAGAATGGCTTGGCTTGAGCAAGTGTTATGCCAGATTTTGCA‑‑GGGGATGGATTCCAGCAGGTTTTGCGGTGCTGGCTTCTTCGCG‑GGCTCGCCCGCGAAGAAGCCAGCCCAGGCCAACAGTTATCAGGCgta < 6:220913/148‑1 (MQ=255) ccgTTCTAAGGAGATTGGCTTGGCTTGAGCAAGTGTTATGCCAGATTTTGCA‑‑GGGGATGGATTCCAGCAGGTTTTGCGGTGCTGGCTTCTTCGCG‑GGCTCGCCCGCGAAGAAGCCAGCCCAGGCCAACAGTTATCAGGCGTAGTAATc < 3:231318/146‑1 (MQ=255) aaaTGAGTATGGCTTGGCTTGAGCATGTGTTAAGCCAGATTTTGCA‑‑GGGGATGGATTCCAGCAGGTTTTGCGGTGCTTGCTTCTTCGCGTCGCTCGCCCGCGAAGAAGCCAGCCCAGGCCAACAGTTATCAGGCGTAGTAATCGACCg < 1:119811/148‑1 (MQ=255) aaaGGAGAATGGCTTGGCTTGAGCAAGTGTTATGCCAGATTTTGCA‑‑GGGGATGGATTCCAGCAGGTTTTGCGGTGCTGGCTTCTTCGCG‑GGCTCGCCCGCGAAGAAGCCAGCCCAGGCCAACAGTTATCAGGCGTAGTAATCGACCGCa > 2:304440/1‑149 (MQ=255) aaaGGAGAATGGCTTGGCTTGAGCAAGTGTTATGCCAGATTTTGCA‑‑GGGGATGGATTCCAGCAGGTTTTGCGGTGCTGGCTTCTTCGCG‑GGCTCGCCCGCGAAGAAGCCAGCCCAGGCCAACAGTTATCAGGCGTAGTAATCGACCGCa < 1:304440/149‑1 (MQ=255) aGGAGAATGGCTTGGCTTGAGCAAGTGTTATGCCAGATTTTGCA‑‑GGGGATGGATTCCAGCAGGTTTTGCGGTGCTGGCTTCTTCGCG‑GGCTCGCCCGCGAAGAAGCCAGCCCAGGCCAACAGTTATCAgg > 8:314680/1‑130 (MQ=255) aGGAGAATGGCTTGGCTTGAGCAAGTGTTATGCCAGATTTTGCA‑‑GGGGATGGATTCCAGCAGGTTTTGCGGTGCTGGCTTCTTCGCG‑GGCTCGCCCGCGAAGAAGCCAGCCCAGGCCAACAGTTATCAgg < 7:314680/130‑1 (MQ=255) aGGAGAATGGCTTGGCTTGAGCAAGTGTTATGCCAGATTTTGCA‑‑GGGGATGGATTCCAGCAGGTTTTGCGGTGCTGGCTTCTTCGCG‑GGCTCGCCCGCGAAGAAGCCACCCCAGGCCAACACTTAAAAGGCGTAGTAATCGagagagg > 4:231318/1‑142 (MQ=255) gCAAGTGTTATGCCAGATTTTGCA‑‑GGGGATGGATTCCAGCAGGTTTTGCGGTGCTGGCTTCTTCGCG‑GGCTCGCCCGCGAAGAAGCCAGCCCAGGCCAACAGTTATCAGGCGTAGTAATCGACCGCACTGTCGCCGACAACCACctgct > 7:333075/1‑149 (MQ=255) gCCAGATTTTGCA‑‑GGGGATGGATTCCAGCAGGTTTTGCGGTGCTGGCTTCTTCGCG‑GGCTCGCCCGCGAAGAAGCCAGCCCAGGCCAACAGTTATCAGGCGTAGTAATCGACCGCACTGTCGCCGACAACCACCTGCTGCTCCAAAGGc > 5:131418/1‑149 (MQ=255) gCCAGATTTTGCA‑‑GGGGATGGATTCCAGCAGGTTTTGCGGTGCTGGCTTCTTCGCG‑GGCTCGCCCGCGAAGAAGCCAGCCCAGGCCAACAGTTATCAGGCGTAGTAATCGACCGCACTGTCGCCGACAACCACCTGCTGCTCCAAAGGc > 1:329461/1‑149 (MQ=255) gATTTTGCA‑‑GGGGATGGATTCCAGCTGGTTTTGCGGTGCTGGCTTCTTCGCG‑GGCGCGCCCGCGAAGAAGCCAGCCCAGGCCAACAGTTATCAGGCGTAGTAATCGACCGCACTGTCGCCGACAACCACCTGCTGCTCCAAAGGCCggg < 6:128716/149‑1 (MQ=255) ttGCA‑‑GGGGATGGATTCCAGCAGGTTTTGCGGTGCTGGCTTCTTCGCG‑GGCTCGCCCGCGAAGAAGCCAGCCCAGGCCAACAGTTATCAGGCGTAGTAATCGACCGCACTGTCGCCGACAACCACCTGCTGCTCCAAAGGCCGGgcgc < 8:333075/148‑1 (MQ=255) cagGTTTTGCGGTGCTGGCTTCTTCGCG‑GGCTCGCCCGCGAAGAAGCCAGCCCAGGCCAACAGTTATCAGGCGTAGTAATCGACCGCACTGTCGCCGACAACCACCTGCTGCTCCAAAGGCCGGGCGCTATCGACCGCTTCGGCCCCAc > 7:199248/1‑149 (MQ=255) ggTGCTGGCTTCTTCGCG‑GGCTCGCCCGCGAAGAAGCCAGCCCAGGCCAACAGTTATCAGGCGTAGTAATCGACCGCACTGTCGCCGACAACCACCTGCTGCTCCAAAGGCCGGGCGCTATCGACCGCTTCGGCCCCACCCTGCTCGGc > 4:145965/1‑149 (MQ=255) gTGCTGGCTTCTTCGCG‑GGCTCGCCCGCGAAGAAGCCAGCCCAGGCCAACAGTTATCAGGCGTAGTAATCGACCGCACTGTCGCCGACAACCACCTGCTGCTCCAAAGGCCGGGCGCTATCGACCGCTTCGGCCCCACCCTGCTCGGcc < 4:411859/149‑1 (MQ=255) ttcttcGCG‑GGCGAGCCCGCGAAGAAGCCAGCCCAGGCCAACAGTTATCAGGCGTAGTAATCGACCGCACTGTCGCCGACAACCACCTGCTGCTCCAAAGGCCGGGCGCTAACGACCGCTTCGGCCCCACACTCCTCGGCCCGgagcgc > 8:245453/1‑149 (MQ=255) tcttcGCG‑GGCGAGCCCGCGAAGAAGCCAGCCCAGGCCAACAGTTATCAGGCGTAGTAATCGACCGCACTGTCGCCGACAACCACCTGCTGCTCCAAAGGCCGGGCGCTATCGACCGCTTCGGCCCCACCCTGCTCGGCCCGGCGcgcc > 3:250124/1‑149 (MQ=255) cttcGCG‑GGCTCGCCCGCGAAGAAGCCAGCCCAGGCCAACAGTTATCAGGCGTAGTAATCGACCGCACTGTCGCCGTCAACCACCTGCTGCTCCAAAGGGCGGGCGCTATCGACCGCTTCGGCCCCACCCTGCTCGGCCCGgcgcgc < 4:350717/147‑1 (MQ=255) cttcGCG‑GGCTCGCCCGCGAAGAAGCCAGCCCAGGCCAACAGTTATCAGGCGTAGTAATCGACCGCACTGTCGCCGACAACCACCTGCTGCTCCAAAGGCCGGGCGCTATCGACCGCTTCGGCCCCACCCTGCTCGGCCCGgcgcgc < 5:272322/147‑1 (MQ=255) tcGCG‑GGCTCGCCCGCGAAGAAGCCAGCCCAGGCCAACAGTTATCAGGCGTAGTAATCGACCGCACTGTCGCCGACAACCACCTGCTGCTCCAAAGGCCGGGCGCTATCGACCGCTTCGGCCCCACCCTGCTCGGCCCGGCGcgccgc > 4:294548/1‑148 (MQ=255) tcGCG‑GGCTCGCCCGCGAAGAAGCCAGCCCAGGCCAACAGTTATCAGGCGTAGTAATCGACCGCACTGTCGCCGACAACCACCTGCTGCTCCAAAGGCCGGGCGCTATCGACCGCTTCGGCCCCACCCTGCTCGGCCCGGCGCgccgcc > 6:157544/1‑149 (MQ=255) g‑ggCTCGCCCGCGAAGAAGCCAGCCCAGGCCAACAGTTATCAGGCGTAGTAATCGACCGCACTGTCGCCGACAACCACCTGCTGCTCCAAAGGCCGGGCGCTATCGACCGCTTCGGCCCCACCCTGCTCGGCCCGGCGCGCCGCCAccc < 3:145965/149‑1 (MQ=255) | CTTCGCTCTTCGCCATGCCAATAATCCGTCGTCCATCGGGGGTTCAAAGGAGAATGGCTTGGCTTGAGCAAGTGTTATGCCAGATTTTGCA‑‑GGGGATGGATTCCAGCAGGTTTTGCGGTGCTGGCTTCTTCGCG‑GGCTCGCCCGCGAAGAAGCCAGCCCAGGCCAACAGTTATCAGGCGTAGTAATCGACCGCACTGTCGCCGACAACCACCTGCTGCTCCAAAGGCCGGGCGCTATCGACCGCTTCGGCCCCACCCTGCTCGGCCCGGCGCGCCGCCACCC > NC_002947/4950835‑4951116 |
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: ‑ |