Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | NC_002947 | 4,980,638 | T→A | 12.7% | intergenic (‑32/+39) | PP_4387 ← / ← flgE | hypothetical protein/flagellar hook protein FlgE |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_002947 | 4,980,638 | 0 | T | A | 12.7% | 150.9 / 13.6 | 55 | intergenic (‑32/+39) | PP_4387/flgE | hypothetical protein/flagellar hook protein FlgE |
Reads supporting (aligned to +/- strand): ref base T (25/23); new base A (3/4); total (28/27) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 7.05e-01 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.26e-01 |
AACACAGGTAGCGGCCGGTGACCAGGCTTTGCCATTGATAGAACGGGGCGGGGGCAGCGGGGGCAGGCAGAACCAGCACCAGCAAGGGGAACAAGGCCAGTAACGACAGCCGCTTCATTCGCGTGAACCTTCACATTCTTGGGGCCGCTTTGCGGCCCCAATTACAAACACCAAAGCCAATACAATCATCAGGTCATCTGAATGATGGTCTGCATGATGGTGCTTTCAGTGGAAATGGTCTTGGCATTGGCCTGGTAGTTGCTCTGCGCCTTGATCAGGTTTACCAGCTCGCCGGT > NC_002947/4980489‑4980784 | aaCACAGGTAGCGGCCGGTGACCAGGCTTTGCCATTGATAGAACGGGGCGGGGGCAGCGGGGGCAGGCAGAACCAGCACCAGCAAGGGGAACAAGGCCAGTAACGACAGCCGCTTCATTCGCGTGAACCTTCACATTCTTGGGGCCGCtt > 1:599416/1‑150 (MQ=255) cacaGGTAGCGGCCGGTGACCAGGCTTTGCCATTGATAGAACGGGGCGGGGGCAGCGGGGGCAGGCAGAACCAGCACCAGCAAGGGGAACAAGGCCAGTAACGACAGCCGCTTCATTCGCGTGAACCTTCACATTCTTGGGGCCGCTTTg > 1:584102/1‑150 (MQ=255) acaGGTAGCGGCCGGTGACCAGGCTTTGCCATTGATAGAACGGGGCGGGGGCAGCGGGGGCAGGCAGAACCAGCACCAGCAAGGGGAACAAGGCCAGTAACGACAGCCGCTTCATTCGCGTGAACCTTCACATTCTTGGGGCCGCTTTGc < 1:120397/150‑1 (MQ=255) aGGTAGCGGCCGGTGACCAGGCTTTGCCATTGATAGAACGGGGCGGGGGCAGCGGGGGCAGGCAGAACCAGCACCAGCAAGGGGAACAAGGCCAGTAACGACAGCCGATCCTGTCGCGTGAACCTTCACATTCTTGGGGCCGCTTTGCgg > 1:374056/1‑150 (MQ=255) gtgttcGGCCGGTGACCAGGCTTTGCCATTGATAGAACGGGGCGGGGGCAGCGGGGGCAGGCAGAACCAGCACCAGCAAGGGGAACAAGGCCAGTAACGACAGCCGCTTCATTCGCGTGAACCTTCACATTCTTGGGGCCGCTTTGCGGc < 1:439378/145‑1 (MQ=255) cGGCCGGTGACCAGGCTTTGCCATTGATAGAACGGGGCGGGGGCAGCGGGGGCAGGCAGAACCAGCACCAGCAAGGGGAACAAGGCCAGTAACGACAGCCGCTTCATTCGCGTGAACCTTCACATTCTTGGGGCCGCAAAGCGGCCCCaa < 1:77729/150‑1 (MQ=255) ggCCGGTGACCAGGCTTTGCCATTGATAGAACGGGGCGGGGGCAGCGGGGGCAGGCAGAACCAGCACCAGCAAGGGGAACAAGGCCAGTAACGACAGCCGCTTCATTCGCGTGAACCTTCACATTCTTGGGGCCGCTTTGCGGCCCCAAt < 2:1410839/150‑1 (MQ=255) cGGTGACCAGGCTTTGCCATTGATAGAACGGGGCGGGGGCAGCGGGGGCAGGCAGAACCAGCACCAGCAAGGGGAACAAGGCCAGTAACGACAGCCGCTTCATTCGCGTGAACCTTCACATTCTTGGGGCCGCTTTGCGGCCCCaagaat < 2:584102/150‑5 (MQ=255) cGGTGACCAGGCTTTGCCATTGATAGAACGGGGCGGGGGCAGCGGGGGCAGGCAGAACCAGCACCAGCAAGGGGAACAAGGCCAGTAACGACAGCCGCTTCATTCGCGTGAACCTTCACATTCTTGGGGCCGCTTTGCGGCCCCAATTAc > 2:1279707/1‑150 (MQ=255) ccAGGCTTTGCCATTGATAGAACGGGGCGGGGGCAGCGGGGGCAGGCAGAACCAGCACCAGCAAGGGGAACAAGGCCAGTAACGACAGCCGCTTCATTCGCGTGAACCTTCACATTCTTGGGGCCGCAAAGCGGcccc < 2:1306871/138‑1 (MQ=255) ccAGGCTTTGCCATTGATAGAACGGGGCGGGGGCAGCGGGGGCAGGCAGAACCAGCACCAGCAAGGGGAACAAGGCCAGTAACGACAGCCGCTTCATTCGCGTGAACCTTCACATTCTTGGGGCCGCAAAGCGGcccc > 1:1306871/1‑138 (MQ=255) ccAGGCTTTGCCATTGATAGAACGGGGCGGGGGCAGCGGGGGCAGGCAGAACCAGCACCAGCAAGGGGAACAAGGCCAGTAACGACAGCCGCTTCATTCGCGTGAACCTCCACATTCTTGGGGCCGCTTTGCGGCCCCAATTACAAacac < 1:946283/150‑1 (MQ=255) gCTTTGCCATTGATAGAACGGGGCGGGGGCAGCGGGGGCAGGCAGAACCAGCACCAGCAAGGGGAACAAGGCCAGTAACGACAGCCGCTTCATTCGCGTGAACCTTCACATTCTTGGGGCCGCTTTGCGGCCCCAATTACAAACACCaaa < 2:599416/150‑1 (MQ=255) gtcATTGATAGAACGGGGCGGGGGCCGCGGGGGCAGGCAGAACCAGCACCAGCAAGGGGAACAAGGCCAGTAACGACAGCCGCTTCATTCGCGTGAACCTTCACATTCTTGGGGCCGCTTTGCGGCCCCAATTACAAACACCAAAGCCaa > 1:490914/3‑150 (MQ=255) cATTGATAGAACGGGGCGGGGGCAGCGGGGGCAGGCAGAACCAGCACCAGCAAGGGGAACAAGGCCAGTAACGACAGCCGCTTCATTCGCGTGAACCTTCACATTCTTGGGGCCGCTTTGCGGCCCCAATTACAAACACCAAAGCCAATa > 2:1274274/1‑150 (MQ=255) ttGATAGAACGGGGCGGGGGCAGCGGGGGCAGGCAGAACCAGCACCAGCAAGGGGAACAAGGCCAGTAACGACAGCCGCTTCATTCGCGTGAACCTTCACATTCTTGGGGCCGCTTTGCGGCCCCAATTACAAACACCAAAGCCAATACa > 1:820228/1‑150 (MQ=255) gATAGAACGGGGCGGGGGCAGCGGGGGCAGGCAGAACCAGCACCAGCAAGGGGAACAAGGCCAGTAACGACAGCCGCTTCATTCGCGTGAACCTTCACATTCTTGGGGCCGCTTTGCGGCCCCAATTACAAACACCAACGCCAATACAat < 2:354714/150‑1 (MQ=255) tAGAACGGGGCGGGGGCAGCGGGGGCAGGCAGAACCAGCACCAGCAAGGGGAACAAGGCCAGTAACGACAGCCGCTTCATTCGCGTGAACCTTCACATTCTTGGGGCCGCAAAGCGGCCCCaa < 2:181973/123‑1 (MQ=255) tAGAACGGGGCGGGGGCAGCGGGGGCAGGCAGAACCAGCACCAGCAAGGGGAACAAGGCCAGTAACGACAGCCGCTTCATTCGCGTGAACCTTCACATTCTTGGGGCCGCAAAGCGGCCCCaa > 1:181973/1‑123 (MQ=255) gcagcgggggcagGCAGAACCAGCACCAGCAAGGGGAACAAGGCCAGTAACGACAGCCGCTTCATTCGCGTGAACCTTCACATTCTTGGGGCCGCTTTGCGGCCCCAATTACAAACACCAAAGCCAATACAATCATCAGGTCATCTGAat > 2:1044684/1‑150 (MQ=255) gcagAACCAGCACCAGCAAGGGGAACAAGGCCAGTAACGACAGCCGCTTCATTCGCGTGAACCTTCACATTCTTGGGGCCGCTTTGCGGCCCCAATTACAAACACCAAAGCCAATACAATCATCAGGTCATCTGAATGATGGTCTGCatg > 2:580590/1‑150 (MQ=255) agAACCAGCACCAGCAAGGGGAACAAGGCCAGTAACGACAGCCGCTTCATTCGCGTGAACCTTCACATTCTTGGGGCCGCTTTGCGGCCCCAATTACAAACACCAAAGCCAATACAATCATCAGGTCATCTGAATGATGGTCTGCatgat < 2:820228/150‑1 (MQ=255) gAACCAGCACCAGCAAGGGGAACAAGGCCAGTAACGACAGCCGCTTCATTCGCGTGAACCTTCACATTCTTGGGGCCGCTTTGCGGCCCCAATTACAAACACCAAAGCCAATACAATCATCAGGTCATCTGAATGATGGTCTGCatgatg < 2:1359869/150‑1 (MQ=255) gAACCAGCACCAGCAAGGGGAACAAGGCCAGTAACGACAGCCGCTTCATTCGCGTGAACCTTCACATTCTTGGGGCCGCTTTGCGGCCCCAATTACAAACACCAAAGCCAATACAATCATCAGGTCATCTGAATGATGGTCTGCatgatg < 1:1279707/150‑1 (MQ=255) accagcaccagcAAGGGGAACAAGGCCAGTAACGACAGCCGCTTCATTCGCGTGAACCTTCACATTCTTGGGGCCGCAAAGCGGCCCCaa > 1:971440/1‑90 (MQ=255) accagcaccagcAAGGGGAACAAGGCCAGTAACGACAGCCGCTTCATTCGCGTGAACCTTCACATTCTTGGGGCCGCAAAGCGGCCCCaa < 2:971440/90‑1 (MQ=255) ccagcaAGGGGAACAAGGCCAGTAACGACAGCCGCTTCATTCGCGTGAACCTTCACATTCTTGGGGCCGCTTTGCGGCCCCAATTACAAACACCAAAGCCAATACAATCATCAGGTCATCTGAATGATGGTCTGCa < 2:207685/136‑1 (MQ=255) ccagcaAGGGGAACAAGGCCAGTAACGACAGCCGCTTCATTCGCGTGAACCTTCACATTCTTGGGGCCGCTTTGCGGCCCCAATTACAAACACCAAAGCCAATACAATCATCAGGTCATCTGAATGATGGTCTGCa > 1:207685/1‑136 (MQ=255) cagcaAGGGGAACAAGGCCAGTAACGACAGCCGCTTCATTCGCGTGAACCTTCACATTCTTGGGGCCGCTTTGCGGCCCCAATTACAAACACCAAAGCCAATACAATCATCAGGTCATCTGAATGATGGTCTGCATGATGGTGCTTTCAg > 2:590635/1‑150 (MQ=255) agcaAGGGGAACAAGGCCAGTAACGACAGCCGCTTCTCTCGCGTGAACCTTCACATTCTTGGGGCCGCTTTGCGGCCCCAATTACAAACACCAAAGCCAATACAATCATCAGGTCATCTGAATGATGGTCTGCATGATGGTGCTTTCAGt < 1:1143138/150‑1 (MQ=255) gcaTGGGGAACAAGGCCAGTAACGACAGCCGCTTCATTCGCGTGAACCTTCACATTCTTGGGGCCGCTTTGCGGCCCCAATTACAAACACCAAAGCCAATACAATCATCAGGTCATCTGaa > 2:721363/1‑121 (MQ=255) gcaTGGGGAACAAGGCCAGTAACGACAGCCGCTTCATTCGCGTGAACCTTCACATTCTTGGGGCCGCTTTGCGGCCCCAATTACAAACACCAAAGCCAATACAATCATCAGGTCATCTGaa < 1:721363/121‑1 (MQ=255) gcaAGGGGAACAAGGCCAGTAACGACAGCCGCTTCATTCGCGTGAACCTTCACATTCTTGGGGCCGCTTTGCgg > 2:356387/1‑74 (MQ=255) gcaAGGGGAACAAGGCCAGTAACGACAGCCGCTTCATTCGCGTGAACCTTCACATTCTTGGGGCCGCTTTGCgg < 1:356387/74‑1 (MQ=255) aaGGGGAACAAGGCCAGTAACGACAGCCGCTTCATTCGCGTGAACCTTCACATTCTTGGGGCCGCTTTGCGGCCCCAATTACAAACACCAAAGCCAATACAATCATCAGGTCATCTGAatgatg < 1:1131037/124‑1 (MQ=255) aaGGGGAACAAGGCCAGTAACGACAGCCGCTTCATTCGCGTGAACCTTCACATTCTTGGGGCCGCTTTGCGGCCCCAATTACAAACACCAAAGCCAATACAATCATCAGGTCATCTGAatgatg > 2:1131037/1‑124 (MQ=255) aGGGGAACAAGGCCAGTAACGACAGCCGCTTCATTCGCGTGAACCTTCACATTCTTGGGGCCGCTTTGCGGCCCCAATTACAAACACCAAAGCCAATACAATCATCAGGTCATCTGAATGATGGTCTGCATGATGGTGCTTTCAGTGGaa > 2:491454/1‑150 (MQ=255) cAAGGCCAGTAACGACAGCCGCTTCATGCCCGTGAACCTTCACATTCTTGGGGCCGCTTTGCGGCCCCAATTACAAACACCAAAGCCAATACAATCATCAGGTCATCTGAATGATGGTCTGCATGATGGTGCTTTCAGTGGAAATGGTCt < 1:220027/150‑1 (MQ=255) aaGGCCAGTAACGACAGCCGCTTCATTCGCGTGAACCTTCACATTCTTGGGGCCGCTTTGCGGCCCCAATTACAAACACCAAAGCCAATACAATCATCAGGTCATCTGAATGATGGTCTGCATGATGGTGCTTTCAGTGGAAATGGTCtt > 1:1389080/1‑150 (MQ=255) gaaccccGCTTCATTCGAGTGAATCTTCACATTCTTGGGGCCGCTTTGCGGCCCCAATTACAAACACCAAAGCCAATACACACACCAGGTCATCTGAATGATGGTCTGAATGATGGTGCTTTCAGTGGAAATGGTCTTGGCATTGGCCTg < 2:490914/145‑1 (MQ=255) gCTTCATTCGCGTGAACCTTCACATTCTTGGGGCCGCTTTGCGGCCCCAATTACAAACACCAAAGCCAATACAATCATCATGTCATCTGAAAAATTGTCTGCATGATGGTGCTTTCAGTGGAAATGGTCTTGGCATTGGCCTGGTAGTTg > 2:550363/1‑150 (MQ=255) ttCGCGTGAACCTTCACATTCTTGGGGCCGCTTTGCGGCCCCAATTACAAACACCAAAGCCAATACAATCATCAGGTCATCTGAatgat < 2:722771/89‑1 (MQ=255) ttCGCGTGAACCTTCACATTCTTGGGGCCGCTTTGCGGCCCCAATTACAAACACCAAAGCCAATACAATCATCAGGTCATCTGAatgat > 1:722771/1‑89 (MQ=255) cgTGAACCTTCACATTCTTGGGGCCGCTTTGCGGCCCCAATTACAAACACCAAAGCCAATACAATCATCAGGTCATCTGAATGATGGTCTGCATGATGGTGCTTTCAGTGGAAATGGTCTTGGCATTGGCCTGGTAGTTGCTCTGCGCCt < 1:590635/150‑1 (MQ=255) cTTCACATTCTTGGGGCCGCTTTGCGGCCCCAATTACAAACACCAAAGCCAATACAATCATCAGGTCATCTGAATGATGGTCTGCATGATGGTGCTTTCAGTGGAAATGGTCTTGGCATTGGCCTGGTAGTTGCTCTGCGCCTTGATCAg < 1:491454/150‑1 (MQ=255) cTTCACATTCTTGGGGCCGCTTTGCGGCCCCAATTACAAACACCAAAGCCAATACAATCATCAGGTCATCTGAATGATGGTCTGCATGATGGTGCTTTCAGTGGAAATGGAACGGTCATTGGCCTGGTAGTTGCTCTGCGCCTTGATCAg > 1:435024/1‑150 (MQ=255) ttCTTGGGGCCGCTTTGCGGCCCCAATTACAAACACCAAAGCCAATACAATCATCAGGTCATCTGAATGATGGTCTGCATGATGGTGCTTTCAGTGGAAATGGTCTTGGCATTGGCCTGGTAGTTGCTCTGCGCCTTGATCAGGTTTAcc > 1:863031/1‑150 (MQ=255) ttGGGGCCGCTTTGCGGCCCCAATTACAAACACCAAAGCCAATACAATCATCAGGTCATCTGAATGATGGTCTGc > 2:466552/1‑75 (MQ=255) ttGGGGCCGCTTTGCGGCCCCAATTACAAACACCAAAGCCAATACAATCATCAGGTCATCTGAATGATGGTCTGc < 1:466552/75‑1 (MQ=255) ggggCCGCTTTGCGGCCCCAATTACAAACACCAAAGCCAATACAATCATCAGGTCATCTGAATGATGGTCTGCATGATGGTGCTTTCAGTGGAAATGGTCTTGGCATTGGCCTGGTAGTTGCTCTGCGCCTTGATCAGGTTTACCAGCTc > 1:752532/1‑150 (MQ=255) ccGCTTTGCGGCCCCAATTACAAACACCAAAGCCAATACAATCATCAGGTCATCTGAATGATGGTCTGCATGATGGTGCTTTCAGTGGAAATGGTCTTGGCATTGGCCTGGTAGTTGCTCTGCGCCTTGATCAGGTTTACCAGCTCGCCg > 1:1199735/1‑150 (MQ=255) gCTTTGCGGCCCCAATTACAAACACCAAAGCCAATACAATCATCAGGTCATGAGGATGATGGTCTGCATGATGGTGCTTTCAGTGGAAATGGTCTTGGCATTGGCCTGGTAGTTGCTCTGCGCCTTGATCAGGTTTACCAGCTCGCCGGt > 1:812649/1‑150 (MQ=255) gCTTTGCGGCCCCAATTACAAACACCAAAGCCAATACAATCATCAGGTCATCTGAATGATGGTCTGCATGATGGTGCTTTCAGTGGAAATGGTCTTGGCATTGGCCTGGTAGTTGCTCTGCGCCTTGATc < 1:980439/130‑1 (MQ=255) gCTTTGCGGCCCCAATTACAAACACCAAAGCCAATACAATCATCAGGTCATCTGAATGATGGTCTGCATGATGGTGCTTTCAGTGGAAATGGTCTTGGCATTGGCCTGGTAGTTGCTCTGCGCCTTGATc > 2:980439/1‑130 (MQ=255) cTTTGCGGCCCCAATTACAAACACCAAAGCCAATACAATCATCAGGTCATCTGAATGATGGTCTGCATGATGGTGCTTTCAGTGGAAATGGTCTTGGCATTGGCCTGGTAGTTGCTCTGCGcc > 1:649351/1‑123 (MQ=255) cTTTGCGGCCCCAATTACAAACACCAAAGCCAATACAATCATCAGGTCATCTGAATGATGGTCTGCATGATGGTGCTTTCAGTGGAAATGGTCTTGGCATTGGCCTGGTAGTTGCTCTGCGcc < 2:649351/123‑1 (MQ=255) | AACACAGGTAGCGGCCGGTGACCAGGCTTTGCCATTGATAGAACGGGGCGGGGGCAGCGGGGGCAGGCAGAACCAGCACCAGCAAGGGGAACAAGGCCAGTAACGACAGCCGCTTCATTCGCGTGAACCTTCACATTCTTGGGGCCGCTTTGCGGCCCCAATTACAAACACCAAAGCCAATACAATCATCAGGTCATCTGAATGATGGTCTGCATGATGGTGCTTTCAGTGGAAATGGTCTTGGCATTGGCCTGGTAGTTGCTCTGCGCCTTGATCAGGTTTACCAGCTCGCCGGT > NC_002947/4980489‑4980784 |
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: ‑ |