| Predicted mutation | |||||||
|---|---|---|---|---|---|---|---|
| evidence | seq id | position | mutation | freq | annotation | gene | description |
| RA | NC_002947 | 739,205 | G→C | 13.7% | P72A (CCA→GCA) ‡ | pilD ← | leader peptidase/N‑methyltransferase |
| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | NC_002947 | 739,205 | 0 | G | C | 13.7% | 136.6 / 15.3 | 51 | P72A (CCA→GCA) ‡ | pilD | leader peptidase/N‑methyltransferase |
| Reads supporting (aligned to +/- strand): ref base G (23/21); new base C (4/3); total (27/24) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 1.00e+00 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.07e-01 | |||||||||||
AGCAGCGCGCTGGCCACTTCCACTACCGGGTAACGCAGGCTGATGCGGTTTTTGCAGGACGAACAGCGCCCGCCCAAGGCCAGGTAGCTGATGACCGGAATGTTCTCCCACGCACGAATCCGGTGCGCGCAGTGTGGGCACCGGGAGGCGGGCAGGCACAGGTCGAAACGCGCGTGCTGCGTGGTCGGCAGCCCCAATACCTCTTGCGCCTCACGCTGCCATTGCCGCTCCAGCATGATCGGCAGGCGATACACCAGCACGTTGATAAAGCTGCCCACCA > NC_002947/739069‑739348 | agcagcGCGCTGGCCACTTCCACTACCGGGTAACGCAGGCTGATGCGGTTTTTGCAGGACGAACAGCGCCCGCCCAAGGCCAGGTAGCTGATGACCGGAATGTTCTCCCACGCACGAATCCGGTGCGCGCAGTGTGGGCACCGGGAGGCg < 1:292054/150‑1 (MQ=255)agcagcGCGCTGGCCACTTCCACTACCGGGTAACGCAGGCTGATGCGGTTTTTGCAGGACGAACAGCGCCCGCCCAAGGCCAGGTAGCTGATGACCGGAATGTTCTCCCACGCACGAATCCGGTGCCCACACTGCGCgcaccggattcgt > 2:708520‑M2/1‑137 (MQ=255) gcgcgcTGGCCACGTCCACTACCGGGTAACGCAGGCTGATGCGGTTTTTGCAGGACGAACAGCGGCCGGCCAAGGCCAGGTAGCTGATGACCGGAATGTTCTCCCACGCACGAACCCGGTGCGCGCAGTGTGGGCACCGGGAGGCGggca < 2:37808/150‑1 (MQ=255) tGGCCACTTCCACTACCGGGTAACGCAGGCTGATGCGGTTTTTGCAGTACGAACAGCGCCCGCCCAAGGCCAGGTAGCTGATTACCGGAATGTTCTCCAACGCACGAATCCGGTGCGCGCAGTGTGGGCACCGGGAGGCGGGCAGGcaca > 2:449151/1‑150 (MQ=255) tGGCCACTTCCACTACCGGGTAACGCAGGCTGATGCGGTTTTTGCAGGACGAACAGCGCCCGCCCAAGGCCAGGTAGCTGATGACCGGAATGTTCTCCCACGCACGAATCCGGTGCGCGCAGTGTGGGCACCGGGAGGCGGGCAGGcaca > 2:77123/1‑150 (MQ=255) ccACTTCCACTACCGGGTAACGCAGGCTGATGCGGTTTTTGCAGGACGAACAGCGCCCGCCCAAGGCCAGGTAGCTGATGACCGGAATGTTCTCCCACGCACGAATCCGGTGCCCACACTGCGCgcaccggattcgtgcgt < 1:1099569‑M2/141‑18 (MQ=255) ccACTTCCACTACCGGGTAACGCAGGCTGATGCGGTTTTTGCAGGACGAACAGCGCCCGCCCAAGGCCAGGTAGCTGATGACCGGAATGTTCTCCCACGCACGAATCCGGTGCCCACACTGCGCgcaccggattcgtgcgt > 2:1099569‑M2/1‑124 (MQ=255) cACTTCCACTACCGGGTAACGCAGGCTGATGCGGTTTTTGCAGGACGAACAGCGCCCGCCCAAGGCCAGGTAGCTGATGACCGGAATGTTCTCCCACGCACGAATCCGGTGCGCGCAGTGTGGGCACCGGGAGGCgg < 1:1083252/137‑1 (MQ=255) cACTTCCACTACCGGGTAACGCAGGCTGATGCGGTTTTTGCAGGACGAACAGCGCCCGCCCAAGGCCAGGTAGCTGATGACCGGAATGTTCTCCCACGCACGAATCCGGTGCGCGCAGTGTGGGCACCGGGAGGCgg > 2:1083252/1‑137 (MQ=255) cACTACCGGGTAACGCAGGCTGATGCGGTTTTTGCAGGACGAACAGCGCCCGCCCAAGGCCAGGTAGCTGATGACCGGAATGTTCTCCCACGCACGAATCCGGTGCGCGCAGTGTGGGCACCGGGAGGCGGGCAGGCACAGGTCGAAAcg > 2:935635/1‑150 (MQ=255) ccGGGTAACGCAGGCTGATGCGGTTTTTGCAGGACGAACAGCGCCCGCCCAAGGCCAGGTAGCTGATGACCGGAATGTTCTCCCACGCACGAATCCGGTGCGCGCAGTGTGGGCACCGGGAGGCGGGCAGGCACAGGTCGAAACGCGCGt > 1:1158728/1‑150 (MQ=255) tGATGCGGTTTTTGCAGGACGAACAGCGCCCGCCCAAGGCCAGGTAGCTGATGACCGGAATGTTCTCCCACGCACGAATCCGGTGCGCGCAGTGTGGGCACCGGGAGGCGGGCAGGCACAGGTCGAAAc > 2:838865/1‑129 (MQ=255) tGATGCGGTTTTTGCAGGACGAACAGCGCCCGCCCAAGGCCAGGTAGCTGATGACCGGAATGTTCTCCCACGCACGAATCCGGTGCGCGCAGTGTGGGCACCGGGAGGCGGGCAGGCACAGGTCGAAAc < 1:838865/129‑1 (MQ=255) tGCGGTTTTTGCAGGACGAACAGCGCCCGCCCAAGGCCAGGTAGCTGATGACCGGAATGTTCTCCCACGCACGAATCCGGTGCGCGCAGTGTGGGCACCGGGAGGCGGGCAGGCACAGGTCGaa < 1:426206/124‑1 (MQ=255) tGCGGTTTTTGCAGGACGAACAGCGCCCGCCCAAGGCCAGGTAGCTGATGACCGGAATGTTCTCCCACGCACGAATCCGGTGCGCGCAGTGTGGGCACCGGGAGGCGGGCAGGCACAGGTCGaa > 2:426206/1‑124 (MQ=255) cGGTTTTTGCAGGATGAACAGCGCCCGCCCAAGGCCAGGTAGCTGATGACCAGAATGTTCTCCCACGCACGAATCCGGTGCGCGCAGTGTGGGCACCGGGAGGCGGGCAGGCACAGGTCGAAACGCGCGTGCTGCGTGGTCGGCAGcccc < 2:372080/150‑1 (MQ=255) ggACGAACAGCGCCCGCCCAAGGCCAGGTAGCTGATGACCGGAATGTTCTCCCACGCACGAATCCGGTGCGCGCAGTGTGGGCACCGGGAGGCGGGCAGGCACAGGTCGAAAcgcgc > 2:841000/1‑117 (MQ=255) ggACGAACAGCGCCCGCCCAAGGCCAGGTAGCTGATGACCGGAATGTTCTCCCACGCACGAATCCGGTGCGCGCAGTGTGGGCACCGGGAGGCGGGCAGGCACAGGTCGAAAcgcgc < 1:841000/117‑1 (MQ=255) gACGAACAGCGCCCGCCCAAGGCCAGGTAGCTGATGACCGGAATGTTCTCCCACGCACGAATCCGGTGCGCGCAGTGTGGGCACCGGGAGGCGGGCAGGCACAGGTCGAAACGCGCGTGCTGCGTGGTCGGCAGcc > 2:946872/1‑136 (MQ=255) gACGAACAGCGCCCGCCCAAGGCCAGGTAGCTGATGACCGGAATGTTCTCCCACGCACGAATCCGGTGCGCGCAGTGTGGGCACCGGGAGGCGGGCAGGCACAGGTCGAAACGCGCGTGCTGCGTGGTCGGCAGcc < 1:946872/136‑1 (MQ=255) gACGAACAGCGCCCGCCCAAGGCCAGGTAGCTGATGACCGGAATGTTCTCCCACGCACGAATCCGGTGCGCGCAGTGTGGGCACCGGGAGGCGGGCAGGCACAGGTCGAAACGCGCGTGCTGCGTGGTCGGCAGCCCCAATACCTCTTgc > 2:1184642/1‑150 (MQ=255) gACGAACAGCGCCCGCCCAAGGCCAGGTAGCTGATGACCGGAATGTTCTCCCACGCACGAATCCGGTGCCCACACTGCGCgcaccggattcgt < 1:208935‑M2/93‑14 (MQ=255) gACGAACAGCGCCCGCCCAAGGCCAGGTAGCTGATGACCGGAATGTTCTCCCACGCACGAATCCGGTGCCCACACTGCGCgcaccggattcgt > 2:208935‑M2/1‑80 (MQ=255) gAACAGCGCCCGCCCAAGGCCAGGTAGCTGATGACCGGAATGTTCTCCCACGCACGAATCCGGTGCGCGCAGTGTGGGCACCGGGAGGCGGGCAGGCACAGGTCGAAACGCGCGTGCTGCGTGGTCGGCAGCCCCAATACCTCTTGCGcc < 2:517556/150‑1 (MQ=255) aCAGCGCCCGCCCAAGGCCAGGTAGCTGATGACCGGAATGTTCTCCCACGCACGAATCCGGTGCGCGAAGTGTGGGCACTGGGAGGCGGGCAGGCACAGGTCGAAACGCGCGTGCTGCGTGGTCGGCAGCCCCAATACCTCTTGCGCCTc < 2:887059/150‑1 (MQ=255) gcgcCCGCCCAAGGCCAGGTAGCTGATGACCGGAATGTTCTCCCACGCACGAATCCGGTGCGCGCAGTGTGGGCACCGGGAGGCGGGCAGGCACAGGTCGAAACGCGCGTGCTGCGTGGTCGGCAGCCCCAATACCTCTTGCGCCTCACg > 2:729659/1‑150 (MQ=255) cccgcccAAGGCCAGGTAGCTGATGACCGGAATGTTCTCCCACGCACGAATCCGGTGCGCGCAGTGTGGGCACCGGGAGGCGGGCAGGCACAGGTCGAAACGCGCGTGCTGCGTGGTCGGCAGCCCCAATACCTCTTGCGCCTCACGCTg > 1:60964/1‑150 (MQ=255) ccAAGGCCAGGTAGCTGATGACCGGAATGTTCTCCCACGCACGAATCCGGTGCGCGCAGTGTGGGCACCGGGAGGCGGGCAGGCACAGGTCGAAAcgcgc < 1:18973/100‑1 (MQ=255) ccAAGGCCAGGTAGCTGATGACCGGAATGTTCTCCCACGCACGAATCCGGTGCGCGCAGTGTGGGCACCGGGAGGCGGGCAGGCACAGGTCGAAAcgcgc > 2:18973/1‑100 (MQ=255) ccAAGGCCAGGTAGCTGATGACCGGAATGTTCTCCCACGCACGAATCCGGTGCCCACACTGCGCgcaccggattcgtgcgtggga > 2:1041301‑M2/1‑64 (MQ=255) ccAAGGCCAGGTAGCTGATGACCGGAATGTTCTCCCACGCACGAATCCGGTGCCCACACTGCGCgcaccggattcgtgcgtggga < 1:1041301‑M2/85‑22 (MQ=255) cTGATGACCGGAATGTTCTCCCACGCACGAATCCGGTGCGCGCAGTGTGGGCACCGGGAGGCGGGCAGGCACAGGTCGAAACGCGCGTGCTGCGTGGTCGGCAGCCCCAATACCTCTTGCGCCTCACGCTGCCATTGCCGCTCCAGCATg > 1:601522/1‑150 (MQ=255) atgaCCGGAATGTTCTCCCACGCACGAATCCGGTGCGCGCAGTGTGGGCACCGGGAGGCGGGCAGGCACAGGTCGAAACGCGCGTGCTGCGTGGTCGGCAGCCCCAATACCTCTTGCGCCTCACGCTGCCATTGCCGCTCCAGCATGATc < 2:29778/150‑1 (MQ=255) ccGGAATGTTCTCCCACGCACGAATCCGGTGCGCGCAGTGTGGGCACCGGGAGGCGGGCAGGCACAGGTCGAAACGCGCGTGCTGCGTGGTCGGCAGCCCCAATACCTCTTGCGCCTCACGCTGCCATTGCCGCTCCAGCATGATCGGCa < 2:601522/150‑1 (MQ=255) ccGGAATGTTCTCCCACGCACGAATCCGGTGCGCGCAGTGTGGGCACCGGGAGGCGGGCAGGCACAGGTCGAAACGCGCGTGCTGCGTGGTCGGCAGCCCCAATACCTCTTGCGCCTCACGCTGCCATTGCCGCTCCAGCATGATCGGCa > 1:211110/1‑150 (MQ=255) gAATGTTCTCCCACGCACGAATCCGGTGCGCGCAGTGTGGGCACCGGGAGGCGGGCAGGCACAGGTCGAAACGCGCGTGCTGCGTGGTCGGCAGCCCCAATACCTCTTGCGCCTCACGCTg < 2:39766/121‑1 (MQ=255) gAATGTTCTCCCACGCACGAATCCGGTGCGCGCAGTGTGGGCACCGGGAGGCGGGCAGGCACAGGTCGAAACGCGCGTGCTGCGTGGTCGGCAGCCCCAATACCTCTTGCGCCTCACGCTg > 1:39766/1‑121 (MQ=255) gAATGTTCTCCCACGCACGAATCCGGTGCGCGCAGTGTGGGCACCGGGAGGCGGGCAGGCACAGGTCGAAACGCGCGTGCTGCGTGGTCGGCAGCCCCAATACCTCTTGCGCCTCACGCTGCCATTGCCGCTCCAGCATGATCGGCAGGc > 1:95129/1‑150 (MQ=255) gTTCTCCCACGCACGAATCCGGTGCGCGCAGTGTGGGCACCGGGAGGCGGGCAGGCACAGGTCGAAACGCGCGTGCTGCGTGGTCGGCAGCCCCAATACCTCTTGCGCCTCACGCTGCCATTGCCGCTCCAGCATGATCGGCAGGCGATa < 1:1184642/150‑1 (MQ=255) ttCTCCCACGCACGAATCCGGTGCGCGCAGTGTGGGCACCGGGAGGCGGGCAGGCACAGGTCGAAACGCGCGTGCTGCGTGGTCGGCAGCCCCAATACCTCTTGCGCCTCACGCTGc > 2:164997/1‑117 (MQ=255) ttCTCCCACGCACGAATCCGGTGCGCGCAGTGTGGGCACCGGGAGGCGGGCAGGCACAGGTCGAAACGCGCGTGCTGCGTGGTCGGCAGCCCCAATACCTCTTGCGCCTCACGCTGc < 1:164997/117‑1 (MQ=255) ctcCCACGCACGAATCCGGTGCGCGCAGTGTGGGCACCGGGAGGCGGGCAGGCACAGGTCGAAACGCGCGTGCTGCGTGGTCGGCAGCCCCAATACCTCTTGCGCCTCACGCTGCCATTGCCGCTCCAGCATGATCGGCAGGCGATacac < 1:77123/150‑1 (MQ=255) gAATCCGGTGCGCGCAGTGTGGGCACCGGGAGGCGGGCAGGCACAGGTCGAAACGCGCGTGCTGCGTGGTCGGCAGCCCCAATACCTCTTGCGCCTCACGCTGCCATTGCCGCTCCAGCATGATCGGCAGGCGATACACCAGCACGTTGa > 1:882818/1‑150 (MQ=255) ggTGCGCGCAGTGTGGGCACCGGGAGGCGGGCAGGCACAGGTCGAAACGCGCGTGCTGCGTGGTCGGCAGCCCCAATACCTCTTGCGCCTCACGCTGCCATTGCCGCTCCAGCATGATCGGCAGGCGATacac > 2:548115/1‑133 (MQ=255) ggTGCGCGCAGTGTGGGCACCGGGAGGCGGGCAGGCACAGGTCGAAACGCGCGTGCTGCGTGGTCGGCAGCCCCAATACCTCTTGCGCCTCACGCTGCCATTGCCGCTCCAGCATGATCGGCAGGCGATacac < 1:548115/133‑1 (MQ=255) cgcAGTGTGGGCACCGGGAGGCGGGCAGGCACAGGTCGAAACGCGCGTGCTGCGTGGTCGGCAGCCCCAATACCTCTTGCGCCTCACGCTGCCATTGCCGCTc > 1:158911/1‑103 (MQ=255) cgcAGTGTGGGCACCGGGAGGCGGGCAGGCACAGGTCGAAACGCGCGTGCTGCGTGGTCGGCAGCCCCAATACCTCTTGCGCCTCACGCTGCCATTGCCGCTc < 2:158911/103‑1 (MQ=255) gcAGTGTGGGCACCGGGAGGCGGGCAGGCACAGGTCGAAACGCGCGTGCTGCGTGGTCGGCAGCCCCAATACCTCTTGCGCCTCACGCTGCCATTGCCGCTCCAGCATGATCGGCAGGCGa > 2:377184/1‑121 (MQ=255) gcAGTGTGGGCACCGGGAGGCGGGCAGGCACAGGTCGAAACGCGCGTGCTGCGTGGTCGGCAGCCCCAATACCTCTTGCGCCTCACGCTGCCATTGCCGCTCCAGCATGATCGGCAGGCGa < 1:377184/121‑1 (MQ=255) gcAGTGTGGGCACCGGGAGGCGGGCAGGCACAGGTCGAAACGCGCGTGCTGCGTGGTCGGCAGCCCCAATACCTCTTGCGCCTCACGCTGCCATTGCCGCTCCAGCATGATCGGCAGGCGATACACCAGCACGTTGATAAAGCTGCccac > 2:696450/1‑150 (MQ=255) aGTGTGGGCACCGGGAGGCGGGCAGGCACAGGTCGAAACGCGCGTGCTGCGTGGTCGGCAGCCCCAATACCTCTTGCGCCTCACGCTGCCATTGCCGCTCCAGCATGATCGGCAGGCGATACACCAGCACGTTGATAAAGCTGCccacca < 2:1158728/150‑1 (MQ=255) ggCACCGGGAGGCGGGCAGGCACAGGTCGAAACGCGCGTGCTGCGTGGTCGGCAGCCCCAATACCTCTTGCGCCTCACGCTGCCa > 2:782954/1‑85 (MQ=255) ggCACCGGGAGGCGGGCAGGCACAGGTCGAAACGCGCGTGCTGCGTGGTCGGCAGCCCCAATACCTCTTGCGCCTCACGCTGCCa < 1:782954/85‑1 (MQ=255) | AGCAGCGCGCTGGCCACTTCCACTACCGGGTAACGCAGGCTGATGCGGTTTTTGCAGGACGAACAGCGCCCGCCCAAGGCCAGGTAGCTGATGACCGGAATGTTCTCCCACGCACGAATCCGGTGCGCGCAGTGTGGGCACCGGGAGGCGGGCAGGCACAGGTCGAAACGCGCGTGCTGCGTGGTCGGCAGCCCCAATACCTCTTGCGCCTCACGCTGCCATTGCCGCTCCAGCATGATCGGCAGGCGATACACCAGCACGTTGATAAAGCTGCCCACCA > NC_002947/739069‑739348 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 22 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |