| Predicted mutation | |||||||
|---|---|---|---|---|---|---|---|
| evidence | seq id | position | mutation | freq | annotation | gene | description |
| RA | NC_002947 | 4,980,590 | T→A | 16.1% | 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 | 16.1% | 77.3 / 5.6 | 36 | intergenic (‑29/+39) | PP_4387/flgE | hypothetical protein/flagellar hook protein FlgE |
| Reads supporting (aligned to +/- strand): ref base T (14/16); new base A (4/2); total (18/18) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 6.58e-01 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 | |||||||||||
CACAGGTAGCGGCCGGTGACCAGGCTTTGCCATTGATAGAACGGGGCGGGGGCAGCGGGGGCAGGCAGAACCAGCACCAGCAAGGGGAACAAGGCCAGTAACGACAGCCGCTTCATTCGCGTGAACCTTCACATTCTT‑‑‑GGCGCTTTGCGGCCCCAATTACAAACACCAAAGCCAATACAATCATCAGGTCATCTGAATGATGGTCTGCATGATGGTGCTTTCAGTGGAAATGGTCTTGGCATTGGCCTGGTAGTTGCTCTGCGCCTTGATCAGGTTTACCAGCTCGCCGGT > NC_002947/4980446‑4980736 | cacatgtagcggcgggtgaccaggctttggcattgatagaacggggcgggggcagcgggggcaggcagaaccagcaccagcaaggggaacaaggccagtaacgacagccgcttcattcgcgtgaaccttcacattcttggggcCGCttt < 6:304437‑M2/6‑1 (MQ=255) gcggccggtgaccaggctttgccattgatagaacggggcgggggcagcgggggcaggcagaaccagcaccagcaaggggaacaaggccagtaacgacagccgcttcattcgcgtgaaccttcacattcttggggcCGCAAAGCGGcccc < 1:139146‑M2/14‑1 (MQ=255) caggctttgccattgatagaacggggcgggggcagcgggggcaggcagaaccagcaccagcaaggggaacaaggccagtaacgacagccgcttcattcgcgtgaaccttcacattcttggggcCGCTTTGCGGCCCCAATTACAAacac > 2:263705‑M2/124‑149 (MQ=255) gccattgatagaacggggcgggggcagcgggggcaggcagaaccagcaccagcaaggggaacaaggccagtaacgacagccgcttcattcgcgtgaaccttcacattcttggggcCGCTTTGCGGCCCCAATTACAAACACCAAAGcc > 8:223935‑M2/116‑148 (MQ=255) ccatttatagaacgtggcgggggcagcgggggcaggcagaactagcaccagcaaggtgaacaaggccagtaacgacagccgcttcattcgcgtgaaccttcacattcttggggcCGCTTTGCGGCCCCAATTACAAACACCAAAGCCaa < 4:293659‑M2/35‑1 (MQ=255) ccattgatagaacggggcgggggcagcgggggcaggcagaaccagcaccagcaaggggaacaaggccagtaacgacagccgcttcattcgcgtgaaccttcacattcttggggcCGCTTTGCGGCCCCAATTACCAACACCAAAGCCa > 3:293659‑M2/115‑148 (MQ=255) ccattgatagaacggggcgggggcagcgggggcaggcagaaccagcaccagcaaggggaacaaggccagtaacgacagccgcttcattcgcgtgaaccttcacattcttggggcCGCTTTGCGGCCCCAATTACAAACACCAAAGCCa < 5:22633‑M2/34‑1 (MQ=255) ttgatagaacgggccgggggcagcgggggcaggcagaaccagcaccagcaaggggaacaaggccagtaacgacagccgcttcattcgcgtgaaccttcacattctctggggcCGCTTTGCGGCCCCAATTACAAACACCAAAGCCAATa < 1:263705‑M2/37‑1 (MQ=255) ggggcgggggcagcgggggcaggcagaaccagcaccagcaaggggaacaaggccagtaacgacagccgcttcattcgcgtgaaccttcacattcttggggcCGCTTTGCGGCCCCAATTACAAACACCAAAGCCAATACAATCATCAgg > 7:217313‑M2/102‑149 (MQ=255) gggcagcgggggcaggcagaaccaccaccagcaaggggaacaaggccagtaacgacagccgcttcattcgcgtgaaccttcacattcttggggcCGCTATGCGGCCCCAAATACAAACACCACAGCCAAAACAATCAACAGCTCATCt > 8:336779‑M2/95‑148 (MQ=255) gcagcgggggcaggcagaaccagcaccagcaaggggaacaaggccagtaacgacagccgcttcattcgcgtgaaccttcacattcttggggcCGCTTTGCGGCCCCAATTACAAACACCAAAGCCAATACAATCATCAGGTCATCTGaa < 8:204803‑M2/57‑1 (MQ=255) gcgggggcaggcagaaccagcaccagcaaggggaacaaggccagtaacgacagccgcttcattcgcgtgaaccttcacattcttggggcCGCTTTGCGGCCCCAATTACAAACACCAAAGCCAATACAATCATCAGGTCATCTGAatga < 3:357225‑M2/60‑1 (MQ=255) gggggcaggcagaaccagcaccagcaaggggaacaaggccagtaacgacagccgcttcattcgcgtgaaccttcacattcttggggcCGCTTTGCGGCCCCAATTACAAACACCAAAGCCAATACAATCATCAGGTCATCTGAatgatg < 3:326620‑M2/62‑1 (MQ=255) agcaccagcaaggggaacaaggccagtaacgacagccgcttcattcgcgtgaaccttcacattcttggggcCGCTTTGCGGCCCCAATTACAAACACCAAAGCCAATACAATCATCAGGTCATCTGAATGATGGTGTGCATGATGGTGc > 7:57932‑M2/72‑149 (MQ=255) ggggaacaaggccagtaacgacagccgcttcattcgcgtgaaccttcacattcttggggcCGCTTTGCGGCCCCCATTAAAAACACCAAAGCCAAAAAAATAATGAAGTAATCTGAATGATGGTCTGCATGATGGTGCTTTCAGTGGa > 2:171669‑M2/61‑148 (MQ=255) aacaaggccagtatcgacagcctcttatttcgcgtgaaccttcacattcttgggggCGCTTTGCGGCCCGAATTACAAACAGCAAAGCCAATACAATCTTCAGGTCATCTGTATGATGGTCTGCATGATGGTGCTTTCAGTGGAAATgg < 7:336779‑M2/93‑1 (MQ=255) aacaaggccagtaacgacagccgcttcattcgcgtgaaccttcacattcttggggcCGCTTTGCGGCCCCCATTACAAACACCAAAGCCAATACAATCATAAGGTCATATGAATGATGGGCTGCATGATGGTGCGGTAAGTGGAAATgg > 2:321029‑M2/57‑149 (MQ=255) cAGTAACGACAGCCGCTTCATTCGCGTGAACCTTCACATTCTTGGGGGCGCGTTGGGGCGCCCATTAGAAACACCAAAACCAAAAAAAAAAACAGGTCAAATGAATGATGGTGTGCAAGATGGTGGTTTAAGTAGTAATGGTGTTTGcg > 2:183860/1‑148 (MQ=255) aacgacagccgcttcattcgcgtgaaccttcacattcttggggcCGCTTTGCGGCCCCAATTACAAACACCAAAGCCAATACAATCATCAGGTCATCTGAATGATGGTCTGCATGATGGTGCTTTCAGTGGAAATGGTCTTGGCATTg < 8:217313‑M2/104‑1 (MQ=255) cgcgcgaaccttcacattcttggggcCGCGTTGCGGCCCCAATTACAAACACCAAAGCCAATACAATCATCAGGGCATCTGAATGATGGTCTGCATGATGGTGCTTTCAGTGGAAATGGTCTTGGCATTGGCCTGGTAGTTGCTCTgc < 1:321029‑M2/122‑1 (MQ=255) cgtgaaccttcacattcttggggcCGCTTTGCGGCCCCAATTACAAACACCAAAGCCAATACAATCATCAGGTCATCTGAATGATGGTCTGCATGATGGTGCTTTCAGTGGAAATGGTCTTGGCATTGGCCTGGTAGTTGCTCTGCGcc < 1:183860‑M2/125‑1 (MQ=255) gtgaaccttcacattcttggggcCGCTTTGCGGCCCCAATTACAAACACCAAAGCCAATACAATCATCAGGTCATCTGAATGATGGTCTGCATGATGGTGCTTTCAGTGGAAATGGTCTTGGCATTGGCCTGGTAGTTGCTCTGCGCCt > 8:357605‑M2/24‑149 (MQ=255) tgaaccttcacattcttggggcCGCTTTGCGGCCCCAATTACAAACACCAAAGCCAATACAATCATCAGGTCATCTGAATGATGGTCTGCATGATGGTGCTTTCAGTGGAAATGGTCTTGGCATTGGCCTGGTAGTTGCTCTGCGCCtt > 3:164589‑M2/23‑149 (MQ=255) ccttcacattcttggggcCGCTTTGCGGCCCCAATTACAAACACCAAAGCCAATACAATCATCAGGTCATCTGAATGATGGTCTGCATGATGGTGCTTTCAGTGGAAATGGTCTTGGCATTGGCCTGGTAGTTGCTCTGCGCCTTGATc < 4:164589‑M2/131‑1 (MQ=255) cttcacattcttggggcCGCTTTGCGGCCCCAATTACAAACACCAAAGCCAATACAATCATCAGGTCATCTGAATGATGGTCTGCATGATGGTGCTTTCAGTGGAAATGGTCTTGGCATTGGCCTGGTAGTTGCTCTGCGCCTTGAt > 3:107756‑M2/18‑147 (MQ=255) tcacattcttggggcCGCTTTGCGGCCCCAATTACAAACACCAAAGCCAATACAATCATCAGGTCATCTGAATGATGGTCTGCATGATGGTGCTTTCAGTGGAAATGGTCTTGGCATTGGCCTGGTAGTTGCTCTgc < 4:73322‑M2/122‑1 (MQ=255) tcacattcttggggcCGCTTTGCGGCCCCAATTACAAACACCAAAGCCAATACAATCATCAGGTCATCTGAATGATGGTCTGCATGATGGTGCTTTCAGTGGAAATGGTCTTGGCATTGGCCTGGTAGTTGCTCTgc > 3:73322‑M2/16‑137 (MQ=255) cacattcttggggcCGCTTTGCGGCCCCAATTACAAACACCAAAGCCAATACAATCATCAGGTCATCTGAATGATGGTCTGCATGATGGTGCTTTCAGTGGAAATGGTCTTGGCATTGGCCTGGTAGTTGCTCTGCGCCTTGATCAgg < 1:171669‑M2/134‑1 (MQ=255) cttggggcCGCTTTGCGGCCCCAATTACAAACACCAAAGCCAATACAATCATCAGGTCATCTGAATGATGGTCTGCATGATGGTGCTTTCAGTGGAAATGGTCTTGGCATTGGCCTGGTAGTTGCTCTGCGCCTTGATCAGGTTTAcc < 7:357605‑M2/140‑1 (MQ=255) attggggcCGCAAAGCGGCCCCAATTACAAACACCAAAGCCAATACAATCATCAGGTCATCTGAATGATGGTCTGCATGATGGTGCTTTCAGTGGAAATGGTCTTGGCATTGGc < 4:122983‑M2/106‑1 (MQ=255) attggggcCGCAAAGCGGCCCCAATTACAAACACCAAAGCCAATACAATCATCAGGTCATCTGAATGATGGTCTGCATGATGGTGCTTTCAGTGGAAATGGTCTTGGCATTGGc > 3:122983‑M2/9‑114 (MQ=255) ttggggcCGCTTTGCGGCCCCAATTACAAACACCAAAGCCAATACAATCATCAGGTCATCTGAATGATGGTCt < 7:253339‑M2/66‑1 (MQ=255) ttggggcCGCTTTGCGGCCCCAATTACAAACACCAAAGCCAATACAATCATCAGGTCATCTGAATGATGGTCt > 8:253339‑M2/8‑73 (MQ=255) ttggggcCGCAAAGCGGCCCCAATTACAAACACCAAAGCCAATACAATCATCAGGTCATCTGAATGATGGTCTGCATGATGGTGCTTTCAGTGGAAATGGTCTTGGCATTGGCCTGGTAGTTGCTCTGCGCCTTGATCAGGTTTACCAg > 1:353656‑M2/8‑149 (MQ=255) ttggggcCGCAAAGCGGCCCCAATTACAAACACCAAAGCCAATACAATCATCAGGTCATCTGAATGATGGTCTGCATGATGGTGCTTTCAGTGGAAATGGTCTTGGAATTGGCCTGGTAGTTGCTCTGAGCATTGATCATGTTtaa > 7:206757‑M2/8‑145 (MQ=255) gcCGCTTTGCGGCCCCAATTACAAACACCAAAGCCAATACAATCATCAGGTCATCTGAATGATGGTCTGCATGATGGTGCTTTCAGTGGAAATGGTCTTGGCATTGGCCTGGTAGTTGCTCTGCGCCTTGATCAGGTTTACCAGCTCGc > 6:289261‑M2/3‑149 (MQ=255) cTTTGCGGCCCCAATTACAAACACCAAAGCCAATACAATCATCAGGTCATCTGAATGATGGTCTGCATGATGGTGCTTTCAGTGCAAATGGTCTTGGCATTGGCCTGGTAGTTGCTCTGCGCCTTGATCAGGTTTACCAGCTCTCCGGt < 5:266555/149‑1 (MQ=255) tttGCGGCCCCAATTACAAACACCAAAGCCAATACAATCATCAGGTCATCTGAATGATGGTCTGCATGATGGTGCTTTCAGTGGAAATGGTCTTGGCATTGGCCTGGTAGTTGCTCTGCGCCTTGATCAGGTTTACCAGCTCGCCGGt > 3:240674/1‑148 (MQ=255) | CACAGGTAGCGGCCGGTGACCAGGCTTTGCCATTGATAGAACGGGGCGGGGGCAGCGGGGGCAGGCAGAACCAGCACCAGCAAGGGGAACAAGGCCAGTAACGACAGCCGCTTCATTCGCGTGAACCTTCACATTCTT‑‑‑GGCGCTTTGCGGCCCCAATTACAAACACCAAAGCCAATACAATCATCAGGTCATCTGAATGATGGTCTGCATGATGGTGCTTTCAGTGGAAATGGTCTTGGCATTGGCCTGGTAGTTGCTCTGCGCCTTGATCAGGTTTACCAGCTCGCCGGT > NC_002947/4980446‑4980736 |
| 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: ‑ |