Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | NC_002947 | 4,980,590 | T→A | 8.2% | intergenic (‑29/+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,590 | 0 | T | A | 8.2% | 143.3 / 3.3 | 49 | intergenic (‑29/+39) | PP_4387/flgE | hypothetical protein/flagellar hook protein FlgE |
Reads supporting (aligned to +/- strand): ref base T (25/20); new base A (2/2); total (27/22) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 1.00e+00 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.68e-01 |
AGGTAGCGGCCGGTGACCAGGCTTTGCCATTGATAGAACGGGGCGGGGGCAGCGGGGGCAGGCAGAACCAGCACCAGCAAGGGGAACAAGGCCAGTAACGACAGCCGCTTCATTCGCGTGAACCTTCACATTCTTGG‑‑‑CGCTTTGCGGCCCCAATTACAAACACCAAAGCCAATACAATCATCAGGTCATCTGAATGATGGTCTGCATGATGGTGCTTTCAGTGGAAATGGTCTTGGCATTGGCCTGGTAGTTGCTCTGCGCCTTGATCAGGTTTACCAGCTCGCC > NC_002947/4980449‑4980733 | aggtagcggccggtgaccaggctttgccattgatagaacggggcgggggcagcgggggcaggcagaaccagcaccagcaaggggaacaaggccagtaacgacagccgcttcattcgcgtgaaccttcacattcttggggcCGCTTTGCg > 3:231502‑M2/141‑149 (MQ=255) gaccaggctttgccattgatagaacggggcgggggcagcgggggcaggcagaaccagcaccagcaaggggaacaaggccagtaacgacagccgcttcattcgcgtgaaccttcacattcttggggcCGCTTTGCGGCCCCAATTACaaa > 6:318817‑M2/127‑149 (MQ=255) ccaggctttgccattgatagaacggggcgggggcagcgggggcaggcagaaccagcaccagcaaggggaacaaggccagtaacgacagccgcttcattcgcgtgaaccttcacattcttggggcCGCTTTGCGGCCCCAATTACAAaca > 7:20278‑M2/125‑149 (MQ=255) caggctttgccattgatagaacggggcgggggcagcgggggcaggcagaaccagcaccagcaaggggaacaaggccagtaacgacagccgcttcattcgcgtgaaccttcacattcttggggcCGCTTTGCGGCCCCAATTACAAacac > 7:54932‑M2/124‑149 (MQ=255) attgatagaacggggcgggggcagcgggggcaggcagaaccagcaccagcaaggggaacaaggccagtaacgacagccgcttcattcgcgtgaaccttcacattcttggggcCGCTTTGCGGCCCCAATTACAAACACCAAAGCCAATa < 8:77961‑M2/37‑1 (MQ=255) atagaacggggcgggggcagcgggggcaggcagcaccagcaccagcaaggggaacaaggccagtaacgacagccgcttcattcgcgtgaaccttcacattcttggggcCGCTTTGCGGCCCCAATTACAAACCCCAAAGCCAATACaa > 2:316215‑M2/109‑148 (MQ=255) cggggcgggggcagcgggggcaggcagaaccagcaccagcaaggggaacaaggccagtaacgacagccgcttcattcgcgtgaaccttcacattattggggcCGCAAAGCGGccccaagaat < 8:181780‑M2/20‑5 (MQ=255) cggggcgggggcagcgggggcaggcagaaccagcaccagcaaggggaacaaggccagtaacgacagccgcttcattcgcgtgaaccttcacattattggggcCGCAAAGCGGccccaagaat > 7:181780‑M2/103‑118 (MQ=255) ggggcgggggcagcgggggcaggcagaaccagcaccagcaaggggaacaaggccagtaacgacagccgcttcattcgcgtgaaccttcacattcttggggcCGCTTTGCGGCCCCAATTAc > 3:321858‑M2/102‑121 (MQ=255) gggggcagcgggggcaggcagaaccagcaccagcaaggggaacaaggccagtaacgacagccgcttcattcgcgtgaaccttcacattcttggggcCGCTTTGCGGCCCCAATTACAAACACCAAAGCCAATACAATCATCAGGTCATc < 1:1451‑M2/53‑1 (MQ=255) gggggcagcgggggcaggcagaaccagcaccagcaaggggaacaaggccagtaacgacagccgcttcattcgcgtgaaccttcacattcttggggcCGCTTTGCGGCCCCAATTACAAACACCAAAGCCAATACAATCATCAGGTCATc < 4:231502‑M2/53‑1 (MQ=255) gcagcgggggcaggcagaaccagcaccagcaaggggaacaaggccagtaacgacagccgcttcattcgcgtgaaccttcacattcttggggcCGCTTTGCGGCCCCAATTACAAACACCAAAGCCAATACAATCATCAGGTCATCTGaa < 8:20278‑M2/57‑1 (MQ=255) gcagcgggggcaggcagaaccagcaccagcaaggggaacaaggccagtaacgacagccgcttcattcgcgtgaaccttcacattcttggggcCGCTTTGCGGCCCCAATTACAAACACCAAAGCCAATACAATCATCAGGTCATCTGaa > 4:135952‑M2/93‑149 (MQ=255) gggcaggcagaaccagcaccagcaaggggaacaaggccagtaacgacagccgcttcattcgcgtgaaccttcacattcttggggcCGCTTTGCGGCCCCAATTACAAACACCAAAGCCAATACAATCATCAGGTCATCTGAatgat < 5:226457‑M2/61‑1 (MQ=255) gggcaggcagaaccagcaccagcaaggggaacaaggccagtaacgacagccgcttcattcgcgtgaaccttcacattcttggggcCGCTTTGCGGCCCCAATTACAAACACCAAAGCCAATACAATCATCAGGTCATCTGAatgat > 6:226457‑M2/86‑146 (MQ=255) ggcagaaccagcaccagcaaggggaacaaggccagtaacgacagccgcttcattcgcgtgaaccttcacattcttggggcCGCTTTGCGGCCCCAATTACAAACACCAAAGCCAATACAATCATCAGGTCATCTGAATGATGGTCTGCa < 7:239482‑M2/69‑1 (MQ=255) caaggggaacaaggccagtaacgacagccgcttcattcgcgtgaaccttcacattcttggggcCGCTTTGCGGCCCCAATTACAAACACCAAAGCCAATACAATCATCAGGTCATCTGAATGATGGTCTGCATGATGGTGCTTTCAGTg < 8:54932‑M2/86‑1 (MQ=255) aaggggaacaaggccagtaacgacagccgcttcattcgcgtgaaccttcacattcttggggcCGCTTTGCGGCCCCAATTACAAACACCAAAGCCAATACAATCATCAGGTCATCTGAATGATGGTCTGCATGATGGTGCTTTCAGTgg < 1:316215‑M2/87‑1 (MQ=255) ggggaacaaggccagtaacgacagccgcttcattcgcgtgaaccttcacattcttggggcCGCTTTGCGGCCCCAATTACAAACACCAAAGCCAATACAATCATCAGTTCATCTGAatgatg < 7:304671‑M2/62‑1 (MQ=255) ggggaacaaggccagtaacgacagccgcttcattcgcgtgaaccttcacattcttggggcCGCTTTGCGGCCCCAATTACAAACACCAAAGCCAAATCAATCATCAGGTCATCTGAatgatg > 8:304671‑M2/61‑122 (MQ=255) ggaacaaggccagtaacgacagccgcttcattcgcgtgaaccttcacattcttggggcCGCTTTGCGGCCCCAATTACAAACACCAAAGCCAATACAATCATCAGGTCATCTGAATGATGGTCTGCatgatg < 8:229226‑M2/74‑1 (MQ=255) ggaacaaggccagtaacgacagccgcttcattcgcgtgaaccttcacattcttggggcCGCTTTGCGGCCCCAATTACAAACACCAAAGCCAATACAATCATCAGGTCATCAGAATGATGGTCTGCatgatg > 7:229226‑M2/59‑132 (MQ=255) ggccagtaacgacagccgcttcattcgcgtgaaccttcacattcttggggcCGCTTTGCGGCCCCAATTACAAACACCAAAGCCAATACAATCATCAGGTCATCTGAATGATGGTCTGCATGATGGTGCTTTCAGTGGAAATGGTCTTg > 2:36981‑M2/52‑149 (MQ=255) cgacagccgcttcattcgcgtgaaccttcacattcttggggcCGCTTTGCGGCCCCAATTACAAACACCAAAGCCAATACAATCATCAGGTCATCTGAATGATGGTCTGCATGATGGTGCTTTCAGTGGAAATGGTCTTGGCATTGGcc < 5:318817‑M2/107‑1 (MQ=255) gacagccgcttcattcgcgtgaaccttcacattcttggggcCGCTTTGCGGCCCCAATTACAAACACCAAAGCCAATACAATAATCAGGTCATCTGAATGATGGTCTGCATGATGGTGCTTTCAGTGGAAATGGTCTTGGCATTGGCCt > 4:179692‑M2/42‑149 (MQ=255) cgcttcattcgcgtgaaccttcacattcttggggcCGCTTTGCGGCCCCAATTACAAACACCAAAGCCAATACAATCATCAGGTCATCTGAATGATGGTCTGCATGATGGTGCTTTCAGTGGAAATGGTCTTGGCATTGGCCTgg < 5:57880‑M2/110‑1 (MQ=255) cGCTTCATTCGCGTGAACCTTCACATTCTTGGGGCCGCGTTGGGGCCCCCAATAAAAAAAACCAACCCAAAAAAAAAAAAAGGTAAAATGAATGATGGGATTGATTATTGTGGTTTGAGTGGAAATTGTGTTGGGATTTGGCTGGTatg > 6:57880/1‑147 (MQ=255) cttcattcgcgtgaaccttcacattcttggggcCGCTTTGCGGCCCCAATTACAAACCCCAAACCCAAAACAATCAACAGGTCATCTGAAAGATGGTGTGCATGATGGGGCAATCACTGGAAATGGACTTGGCATTGGCCTGGTAGTTg > 3:315071‑M2/34‑149 (MQ=255) cattcgcgtgaaccttcacattcttggggcCGCTTTGCGGCCCCAATTACAAACACCAAAGCCAATACAATCATCAGGTCATCTGAATGATGGTCTGCATGATGGTGCTTTCAGTGGAAATGGTCTTGGCATTGGCCTGGTAGTTGct < 1:36981‑M2/118‑1 (MQ=255) gcgtgaaccttcacattcttggggcCGCTTTGCGGCCCCAATTACAAACACCAAAGCCAATACAATCATCAGGTCATCt < 3:266824‑M2/54‑1 (MQ=255) gcgtgaaccttcacattcttggggcCGCTTTGCGGCCCCAATTACAAACACCAAAGCCAATACAATCATCAGGTCATCt > 4:266824‑M2/26‑79 (MQ=255) cgtgaaccttcacattcttggggcCGCTTTGCGGCCCCAATTACAAACACCAAAGCCAATACAATCATCAGGTCATCTGAATGATGGTCTGCATGATGGTGCTTTCAGTGGAAATGGTCTTGGCATTGGCCTGGTAGTTGCTCTGCGcc > 5:366133‑M2/25‑149 (MQ=255) gtgaaccttcacattcttggggcCGCTTTGCGGCCCCAATTACAAACACCAAAGCCAATACAATCATCAGGTCATCTGAATGATGGTCTGCATGATGGTGCTTTCAGt > 7:221715‑M2/24‑108 (MQ=255) gtgaaccttcacattcttggggcCGCTTTGCGGCCCCAATTACAAACACCAAAGCCAATACAATCATCAGGTCATCTGAATGATGGTCTGCATGATGGTGCTTTCAGt < 8:221715‑M2/85‑1 (MQ=255) gtgaaccttcacattcttggggcCGCTTTGCGGCCCCAATTACAAACACCAAAGCCAATACAATCATCAGGTCATCTGAATGATGGTCTGCATGATGGTGCTTTCAGTGGAAATGGTCTTGGCATTGGCCTGGTAGTTGCTCTGCGCCt > 2:42089‑M2/24‑149 (MQ=255) aaccttcacattcttggggcCGCTTTGCGGCCCCAATTACAAACACCa < 6:106138‑M2/28‑1 (MQ=255) aaccttcacattcttggggcCGCTTTGCGGCCCCAATTACAAACACCa > 5:106138‑M2/21‑48 (MQ=255) cactttcttggggcCGCTTTGCGGCCCCAATTTCAAACACCATAGCCAATACAATCATCAGGTCATCTGAATGTTGGTCTGCATGATGGTGCTTTCAGTGGAAATGGGCTTGGCATTGGCCTGGTAGTTGCTCTGCGCCTTGATCAgg < 4:315071‑M2/134‑1 (MQ=255) cacattcttggggcCGCTTTGCGGCCCCAATTACAAACACCAAAGCCAATACAATCATCAGGTCATCTGAATGATGGTCTGCATGATGGTGCTTTCAGTGGAAATGGTCTTGGCATTGGCCTGGTAGTTGCTCTGCGCCTTGATCAg < 3:135952‑M2/133‑1 (MQ=255) cattcttggggcCGCTTTGCGGCCCCAATTACAAACACCAAAGCCAATACAATCATAAGGTGATGGGAATGATGGTGTGCATGATGGTGCTTTAAGTGGAAATGGGCTTGGCATTGGCCGGGTAGTTGCTGTGCGCCTTGATCAGGttt > 8:210602‑M2/13‑149 (MQ=255) attcttggggcCGCTTTGCGGCCCCAATTACAAACACCAAAGCCAATACAATCATCAGGTCATCTGAATGATGGTCTGCATGATGGTGCTTTCAGTGGAAATGGTCTTGGCATTGGCCTGGTAGTTGCTCTGCGcc > 7:251414‑M2/12‑136 (MQ=255) attcttggggcCGCTTTGCGGCCCCAATTACAAACACCAAAGCCAATACAATCATCAGGTCATCTGAATGATGGTCTGCATGATGGTGCTTTCAGTGGAAATGGTCTTGGCATTGGCCTGGTAGTTGCTCTGCGcc < 8:251414‑M2/125‑1 (MQ=255) ttcttggggcCGCTTTGCGGCCCCAATTACAAACACCAAAGCCAATACAATCATCAGGTCATCTGAATGATGGTCTGCATGATGGTGCTTTCAGTGGAAATGGTCTTg > 6:312159‑M2/11‑108 (MQ=255) atcttggggcCGCTTTGCGGCCCCAATTACAAACACCAAAGCCAATACAATCATCAGGTCATCTGAATGATGGTCTGCATGATGGTGCTTTCAGTGGAAATGGTCTTg < 5:312159‑M2/98‑1 (MQ=255) ttggggcCGCAAAGCGGCCCCAATTACAAACACCAAAGCCAATACAATCATCAGGTCATCTGAATGATGGTCTGCATGATGGTGc < 7:260126‑M2/78‑1 (MQ=255) ttggggcCGCAAAGCGGCCCCAATTACAAACACCAAAGCCAATACAATCATCAGGTCATCTGAATGATGGTCTGCATGATGGTGc > 8:260126‑M2/8‑85 (MQ=255) ggggcCGCTTTGCGGCCCCAATTACAAACACCAAAGCCAATACAATCATCAGGTCATCTGAATGATGGTCTGCATGATGGTGCTTTCAGTGGAAATGGTCTTGGCATTGGCCTGGTAGTTGCTCTGCGCCTTGATCAGGTTTACCAGCt > 2:26077‑M2/6‑149 (MQ=255) ggcCGCTTTGCGGCCCCAATTACAAACACCAAAGCCAATACAATCATCAGGTCATCTGAATGATGGTCTGCATGATGGTGATGTCAGTGGAAATGGTCTTGGCATTGGCCTGGTAGTTGCTCTTCGCCTTGATCAGGTTTACAAGCTCg > 4:280520‑M2/4‑149 (MQ=255) cCGCTTTGCGGCCCCAATTACAAACACCAAAGCCAATACAATCATCAGGTCATCTGAATGATGGTCTGCATGATGGTGCTTTCAGTGGAAATGGTCTTGGCATTGGCCTGGTAGTTGCTCTGCGCCTTGATCAGGTTTACCAGCTCGcc > 7:68991‑M2/2‑149 (MQ=255) | AGGTAGCGGCCGGTGACCAGGCTTTGCCATTGATAGAACGGGGCGGGGGCAGCGGGGGCAGGCAGAACCAGCACCAGCAAGGGGAACAAGGCCAGTAACGACAGCCGCTTCATTCGCGTGAACCTTCACATTCTTGG‑‑‑CGCTTTGCGGCCCCAATTACAAACACCAAAGCCAATACAATCATCAGGTCATCTGAATGATGGTCTGCATGATGGTGCTTTCAGTGGAAATGGTCTTGGCATTGGCCTGGTAGTTGCTCTGCGCCTTGATCAGGTTTACCAGCTCGCC > NC_002947/4980449‑4980733 |
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: ‑ |