Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | NC_002947 | 739,393 | T→A | 9.7% | E9V (GAA→GTA) ‡ | pilD ← | leader peptidase/N‑methyltransferase |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_002947 | 739,393 | 0 | T | A | 9.7% | 109.0 / 4.2 | 41 | E9V (GAA→GTA) ‡ | pilD | leader peptidase/N‑methyltransferase |
Reads supporting (aligned to +/- strand): ref base T (18/19); new base A (2/2); total (20/21) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 1.00e+00 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.27e-01 |
GTCGGCAGCCCCAATACCTCTTGCGCCTCACGCTGCCATTGCCGCTCCAGCATGATCGGCAGGCGATACACCAGCACGTTGATAAAGCTGCCCACCAGCAGGCCGAGCACCGTGGCAAGGGTGAGGAAATACGCTGGCTGTTCAGCCAGCAAAGCCCATAAAGTCATGTTCAGATCAAACTACCCAACTGGAAGATCGGCAGGTACATCGCCACCACCAGGCCACCCACCAGCAAGCCCAGGATCAGCACGATGGCCGGCTCCAGCAGGCTGGTCAACTGATCCAGCGCC > NC_002947/739252‑739541 | ggcggttgCCTCAATACCGCTGTCGCCTCACGCTGCCATTGCCGCTCCAGCATGATCGGCAGGCGATACACCAGCACGTTGATAAAGCTGCCCACCAGCAGGCCGAGCACCGTGGCAAGGGTGAGGAAATACGCTGGCTGTTcagccagc < 2:606445/143‑1 (MQ=255) aGCCCCAATACCTCTTGCGCCTCACGCTGCCATTGCCGCTCCAGCATGATCGGCAGGCGATACACCAGCACGTTGATAAAGCTGCCCACCAGCAGGTCGAGCACCGTGGCAAGGGTGAGGAAATACGCTGGCTGTTCAGCCAGCAAAGcc < 1:233610/150‑1 (MQ=255) cccAATACCTCTTGCGCCTCACGCTGCCATTGCCGCTCCAGCATGATCGGCAGGCGATACACCAGCACGTTGATAAAGCTGCCCACCAGCAGGCCGAGCACCGTGGCAAGGGTGAGGAAATACGCTGGCTGAAcagccagcgtattt < 1:218653/147‑7 (MQ=255) cccAATACCTCTTGCGCCTCACGCTGCCATTGCCGCTCCAGCATGATCGGCAGGCGATACACCAGCACGTTGATAAAGCTGCCCACCAGCAGGCCGAGCACCGTGGCAAGGGTGAGGAAATACGCTGGCTGAAcagccagcgtattt > 2:218653/1‑141 (MQ=255) cacAATACCTCTTGCGCCTCACGCTGCCATTGCCGCTCCAGCATGATCAGTCGTCTATACACCAGCACGTTGATAAAGCTGCCCACCAGCAGGCCGAGCACCGTGGCAAGGGTGAGGAAATACGCTGGCTGTTCAGCCAGCAAAGCCCAt > 1:617297/3‑150 (MQ=255) aaTACCTCTTGCGCCTCACGCTGCCATTGCCGCTCCAGCATGATCGGCAGGCGATACACCAGCACGTTGATAAAGCTGCCCACCAGCAGGCCGAGCACCGTGGCAAGGGTGAGGAAATACGCTGGCTGTTCAGCCAGCAAAGCCCATaaa < 2:297562/150‑1 (MQ=255) tctTGCGCCTCACGCTGCCATTGCCGCTCCAGCATGATCGGCAGGCGATACACCAGCACGTTGATAAAGCTGCCCACCAGCAGGCCGAGCACCGTGGCAAGGGTGAGGAAATACGCTGGCTGAAcagccag < 1:145902/131‑1 (MQ=255) tctTGCGCCTCACGCTGCCATTGCCGCTCCAGCATGATCGGCAGGCGATACACCAGCACGTTGATAAAGCTGCCCACCAGCAGGCCGAGCACCGTGGCAAGGGTGAGGAAATACGCTGGCTGAAcagccag > 2:145902/1‑131 (MQ=255) ttGCGCCTCACGCTGCCATTGCCGCTCCAGCATGATCGGCAGGCGATACACCAGCACGTTGATAAAGCTGCCCACCAGCAGGCCGAGCACCGTGGCAAGGGTGAGGAAATACGCTGGCTGTTCAGCCAGCAAAGCCCATAAAGTCATGtt < 2:750804/150‑1 (MQ=255) gcCTCACGCTGCCATTGCCGCTCCAGCATGATCGGCAGGCGATACACCAGCACGTTGATAAAGCTGCCCACAAGCAGGCCGAGCACCGTGGCAAGGGTGAGGAAATACGCTGGCTGTTCAGCCAGCAAAGCCCATAAAGTCATGTTCAGa < 1:420148/150‑1 (MQ=255) gCTGCCATTGCCGCTCCAGCATGATCGGCAGGCGATACACCAGCACGTTGATAAAGCTGCCCACCAGCAGGCCGAGCACCGTGGCAAGGGTGAGGAAATACGCTGGCTGTTCAGCCAGCAAAGCCCATAAAGTCATGTTCAGATCAAACt > 1:365729/1‑150 (MQ=255) aTGATCGGCAGGCGATACACCAGCACGTTGATAAAGCTGCCCACCAGCAGGCCGAGCACCGTGGCAAGGGTGAGGAAATACGCTGGCTGTTCAGCCAGCAAAGCCCATAAAGTCATGTTCAGATCAAACTACCCAACTGGAAGATCGGCa < 1:472305/150‑1 (MQ=255) tGATCGGCAGGCGATACACCAGCACGTTGATAAAGCTGCCCACCAGCAGGCCGAGCACCGTGGCAAGGGTGAGGAAATACGCTGGCTGTTCAGCCAGCAAAGCCCATAAAGTCATGTTCAGATCAAACTACCCAAc < 2:424350/136‑1 (MQ=255) tGATCGGCAGGCGATACACCAGCACGTTGATAAAGCTGCCCACCAGCAGGCCGAGCACCGTGGCAAGGGTGAGGAAATACGCTGGCTGTTCAGCCAGCAAAGCCCATAAAGTCATGTTCAGATCAAACTACCCAAc > 1:424350/1‑136 (MQ=255) aTCGGCAGGCGATACACCAGCACGTTGATAAAGCTGCCACACAGCAGGCCGAGCACCGTGGCACGGGTGAGTAATTTTTCTGGCTGTTCTGCCAGCAAAGCCCATAAAGTCATGTTCAGATCAAACTACCCAACTGGAAGATCGGCAGGt < 2:84031/150‑1 (MQ=255) ggCAGGCGATACACCAGCACGTTGATAAAGCTGCCCACCAGCAGGCCGAGCACCGTGGCAAGGGTGAGGAAATACGCTGGCTGTTCAGCCAGCAAAGCCCATAAAGTCATGTTCAGATCaaa < 1:556298/122‑1 (MQ=255) ggCAGGCGATACACCAGCACGTTGATAAAGCTGCCCACCAGCAGGCCGAGCACCGTGGCAAGGGTGAGGAAATACGCTGGCTGTTCAGCCAGCAAAGCCCATAAAGTCATGTTCAGATCaaa > 2:556298/1‑122 (MQ=255) gCGATACACCAGCACGTTGATAAAGCTGCCCACCAGCAGGCCGAGCACCGTGGCAAGGGTGAGGAAATACGCTGGCTGTTCAGCCATCAAAGCCCATAAAGTCATGTTCAGATCAAACTACCCAACTGGAAGATCGGCAGGTACATCGcc > 2:24068/1‑150 (MQ=255) acCAGCACGTTGATAAAGCTGCCCACCAGCAGGCCGAGCACCGTGGCAAGGGTGAGGAAATACGCTGGCTGTTCAGCCAGCAAAGCCCATAAAGTCATGTTCAGATCAAACTACCCAACTGGAAGATCGGCAGGTACATCGcca < 1:189398/144‑1 (MQ=255) acCAGCACGTTGATAAAGCTGCCCACCAGCAGGCCGAGCACCGTGGCAAGGGTGAGGAAATACGCTGGCTGTTCAGCCAGCAAAGCCCATAAAGTCATGTTCAGATCAAACTACCCAACTGGAAGATCGGCAGGTACATCGcca > 2:189398/1‑144 (MQ=255) ccAGCACGTTGATAAAGCTGCCCACCAGCAGGCCGAGCACCGTGGCAAGGGTGAGGAAATACGCTGGCTGTTCAGCCAGCAAAGCCCATAAAGTCATGTTCAGATCAAACTACCCAACTGGAAGATCGGCAGGTACATCGCCACCACCAg < 1:250845/150‑1 (MQ=255) gTTGATAAAGCTGCCCACCAGCAGGCCGAGCACCGTGGCAAGGGTGAGGAAATACGCTGGCTGTTCAGCCAGCAAAGCCCATAAAGTCATGTTCAGATCAAACTACCCAACTGGAAGATCGGCAGGTACATCGCCACCACCAGGccaccc > 1:180570/1‑150 (MQ=255) gATAAAGCTGCCCACCAGCAGGCCGAGCACCGTGGCAAGGGTGAGGAAATACGCTGGCTGTTCAGCCAGCAAAGCCCATAAAGTCATGTTCAGATCAAACTACCCAACTg < 2:449015/110‑1 (MQ=255) gATAAAGCTGCCCACCAGCAGGCCGAGCACCGTGGCAAGGGTGAGGAAATACGCTGGCTGTTCAGCCAGCAAAGCCCATAAAGTCATGTTCAGATCAAACTACCCAACTg > 1:449015/1‑110 (MQ=255) ccaccaGCAGGCCGAGCACCGTGGCAAGGGTGAGGAAATACGCTGGCTGTTCAGCCAGCAAAGCCCATAAAGTCATGTTCAGATCAAACTACCCAACTGGAAGATCGGCAGGTACATCGCCACCACCAGGCCACCCACCAGCAAGccc < 2:574399/148‑1 (MQ=255) ccaccaGCAGGCCGAGCACCGTGGCAAGGGTGAGGAAATACGCTGGCTGTTCAGCCAGCAAAGCCCATAAAGTCATGTTCAGATCAAACTACCCAACTGGAAGATCGGCAGGTACATCGCCACCACCAGGCCACCCACCAGCAAGccc > 1:574399/1‑148 (MQ=255) accaGCAGGCCGAGCACCGCGGCAAGGGTGAGGAAATACGCTGGCTGTTCAGCCAGCAAAGCCCAAAAAGTCATGTTCACATCAAACTACCCAACTAGAAGATCGGCAGGTACATCGCCACCACCAGGCCACCCACCAGCAAGCCCAGGa > 2:471283/1‑150 (MQ=255) ccaGCAGGCCGAGCACCGTGGCAAGGGTGAGGCAATAAGGTAAAAAATATCCTAGCAAAGGCCGTAAAGTGATGCTCAGCAAAAACTACACAACTACAAGATCGGCAGGTTCATCGCCACCACCAGGCCACCCACCAGCAAGCCCAGGAt > 1:578589/1‑150 (MQ=255) gcagGCCGAGCACCGTGGCAAGGGTGAGGAAATACGCTGGCTGTTCAGCCAGCAAAGCCCATAAAGTCATGTTCAGATCAAACTACCCAACTGGAAGATGGGCAGGTACATCGCCACCACCAGGCCAGCCAACAGCAAGCCCAGGATAAg > 1:369506/1‑150 (MQ=255) agGCCGAGCACCGTGGCAAGGGTGAGGAAATACGCTGGCTGTTCAGCCAGCAAAGCCCATAAAGTCATGTTCAGATCAAACTACCCAACTGGAAGATCGGCAGGTACATCGCCACCACCAGGCCACCCACCAGCAAGCCCAGGATCAGCa < 1:668853/150‑1 (MQ=255) cGAGCACCGTGGCAAGGGTGAGGAAATACGCTGGCTGTTCAGCCAGCAAAGCCCATAAAGTCATGTTCAGATCAAACTACCCAACTGGAAGATCGGCAGGTACATCGCCACCACCAGGCCACCCACCAGCAAGCCCAGGATCAGCACGAt > 1:465905/1‑150 (MQ=255) gAGCACCGTGGCAAGGGTGAGGAAATACGCTGGCTGTTCAGCCAGCAAAGCCCATAAAGTCATGTTCAGATCAAACTACCCAACTGGAAGATCGGCAGGTACATCGCCACCACCAGGccacccac < 1:489828/125‑1 (MQ=255) gAGCACCGTGGCAAGGGTGAGGAAATACGCTGGCTGTTCAGCCAGCAAAGCCCATAAAGTCATGTTCAGATCAAACTACCCAACTGGAAGATCGGCAGGTACATCGCCACCACCAGGccacccac > 2:489828/1‑125 (MQ=255) gCACCGTGGCAAGGGTGAGGAAATACGCTGGCTGTTCAGCCAGCAAAGCCCATAAAGTCATGTTCAGATCAAACTACCCAACTGGAAGATCGGCAGGTACATCGCCACCACCAGGCCACCCACCAGCAAGCCCAGGATCAGCACGATgga > 1:48881/1‑149 (MQ=255) cACCGTGGCAAGGGTGAGGAAATACGCTGGCTGTTCAGCCAGCAAAGCCCATAAAGTCATGTTCAGATCAAACTACCCAACTGGAAGATCGGCAGGTACATCGCCACCACCAGGCCACCCACCAGCAAGCCCAGGATCAGCACGATGGcc < 2:465905/150‑1 (MQ=255) gCAAGGGTGCGGAAGTAGGCTGGCTGTTCAGCCAGCAAAGCCCATAAAGTCATGTTCAGATCAAACTACCCAACTGGAAGATCGGCAGGTACATCGCCACCACCAGGCCACCCACCAGCAAGCCCAGGATCAGCACGATGGCCGGCTCca < 2:578589/150‑1 (MQ=255) ggTGAGGAAATACGCTGGCTGTTCAGCCAGCAAAGCCCATAAAGTCATGTTCAGATCAAACTACCCAACTGGAAGATCGGCAGGTACATCGCCACCACCAGGCCACCCACCAGCAAGCCCAGGATCAGCACGATGGCCGGCTCCAGCAgg > 1:434497/1‑150 (MQ=255) gAGGAAATACGCTGGCTGTTCAGCCAGCAAAGCCCATAAAGTCATGTTCAGATCAAACTACCCAACTGGAAGATCGGCAGGTACATCGCCACCACCAGGCCACCCACCAGCAAGCCCAGGATCAGCACGATGGCCGGCTCCAGCAGGCTg < 2:533427/150‑1 (MQ=255) aTACGCTGGCTGTTCAGCCAGCAAAGCCCATAAAGTCATGTTCAGATCAAACTACCCAACTGGAAGATCGGCAGGTACATCGCCACCACCAGGCCACCCACCAGCAAGCCCAGGATCAGCACGATGGCCGGCTCCAGCAGGCTGGTCAAc > 1:143313/1‑150 (MQ=255) aCGCTGGCTGTTCAGCCAGCAAAGCCCATAAAGTCATGTTCAGATCAAACTACCCAACTGGAAGATCGGCAGGTACATCGCCACCACCAGGCCACCCACCAGCAAGCCCAGGATCAGCACGATGGCCGGCTCCAGCAGGCTGGTCAACTg > 2:320296/1‑150 (MQ=255) ctggctgTTCAGCCAGCAAAGCCCATAAAGTCATGTTCAGATCAAACTACCCAACTGGAAGATCGGCAGGTACATCGCCACCACCAGGCCACCCACCAGCAAGCCCAGGATCAGCACGATGGCCGGCTCCAGCAGGCTGGTCAACTGATc < 2:180570/150‑1 (MQ=255) ttCAGCCAGCAAAGCCCATAAAGTCATGTTCAGATCAAACTACCCAACTGGAAGATCGGCAGGTACATCGCCACCACCAGGCCACCCACCAGCAAGCCCAGGATCAGCACGATGGCCGGCTCCAGCAGGCTGGTCAACTGATCCAGCGcc > 1:728897/1‑150 (MQ=255) | GTCGGCAGCCCCAATACCTCTTGCGCCTCACGCTGCCATTGCCGCTCCAGCATGATCGGCAGGCGATACACCAGCACGTTGATAAAGCTGCCCACCAGCAGGCCGAGCACCGTGGCAAGGGTGAGGAAATACGCTGGCTGTTCAGCCAGCAAAGCCCATAAAGTCATGTTCAGATCAAACTACCCAACTGGAAGATCGGCAGGTACATCGCCACCACCAGGCCACCCACCAGCAAGCCCAGGATCAGCACGATGGCCGGCTCCAGCAGGCTGGTCAACTGATCCAGCGCC > NC_002947/739252‑739541 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |