| Predicted mutation | |||||||
|---|---|---|---|---|---|---|---|
| evidence | seq id | position | mutation | freq | annotation | gene | description |
| RA | NC_002947 | 4,049,951 | A→G | 22.3% | intergenic (+40/‑61) | quiA → / → oprB‑III | quinate dehydrogenase/carbohydrate‑selective porin |
| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | NC_002947 | 4,049,951 | 0 | A | G | 22.3% | 24.4 / 4.7 | 18 | intergenic (+40/‑61) | quiA/oprB‑III | quinate dehydrogenase/carbohydrate‑selective porin |
| Reads supporting (aligned to +/- strand): ref base A (11/3); new base G (2/2); total (13/5) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 5.33e-01 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.68e-01 | |||||||||||
GCCGAAGACCGGCAAGCAGTACGTGGTGATCACTGCTGGCGGCGCCCGCCAGTCGACTGACCGTGGCGACTATGTGATTTCTTACGCTCTGCCGTAAAACCGCGTCGCCTGCTTCGCGGGCAAGCCCGCTCCTACAAGGC‑CCCTGTAGGAGCGGGCTTGCCCGCGAATGCCCCTCAAAGACAACCCGAGATTCTGTAATGCCATCCGCAATTCGCATCACCCCCACCCTGCTGCTGGCCCTGGCCAGCACCACTGCCCTTGCCAGCGACGATCTCATGACCCG > NC_002947/4049815‑4050097 | gCCGAAGACCGGCAAGCAGTACGTGGTGATCACTGCTGGCGGCGCCCGCCAGTCGACTGACCGTGGCGACTATGTGATTTCTTACGCTCTGCCGTAAAACCGCGTCGCCTGCTTCGCGGGCAAGCCCGCTCCTACAACCC‑CCCTGTAgaa > 2:152670/1‑148 (MQ=255) aTCACTGCTGGCGGCGCCCGCCAGTCGACTGACCGTGGCGACTATGTGATTTCTTACGCTCTGCCGTAAAACCGCGTCGCCTGCTTCGCGGGCAAGCCCGCTCCTACAAGGC‑CCCTGTAGGAGCGGGCTTGCCCGCGAATGCCCCTCaaa < 1:778285/150‑1 (MQ=255) aCTGCTGGCGGCGCCCGCCAGTCGACTGACCGTGGCGACTATGTGATTTCTTACGCTCTGCCGTAAAACCCACTACCCTGCTCCGCCCGCCCGCCCGCCGATAGGAGGCTGGCTG‑AGTGGCGGGCTTGCCCGCGACTGCCCCTccacgcc > 1:377070/1‑144 (MQ=255) ctgGCGGCGCCCGCCAGTCGACTGACCGTGGCGACTATGTGATTTCTTACGCTCTGCCGTAAAACCGCGTCGCCTGCTTCGCGGGCAAGCCCGCTCCTACAGGGG‑CCTTGTAGGAGc > 1:337670/1‑117 (MQ=38) ctgGCGGCGCCCGCCAGTCGACTGACCGTGGCGACTATGTGATTTCTTACGCTCTGCCGTAAAACCGCGTCGCCTGCTTCGCGGGCAAGCCCGCTCCTACAGGGG‑CCTTGTAGGAGc < 2:337670/117‑1 (MQ=38) gcgggcgccgccAGTCGACTGACCGTGGCGACTATGTGATTTCTTACGCTCTGCCGTAAAACCGCGTCGCCTACTGCGCGGGCAAGCACGCTCCTACAAAGC‑CCCTGTAGGAGCGGGCTTGCCCGCGAATGCCCCTCAAAGACAACCCga > 1:377619/8‑150 (MQ=39) ccAGTCGACTGACCGTGGCGACTATGTGATTTCTTACGCTCTGCCGTAAAACCGCGTCGCCTGCTTCGCGGGCAAGCCCGCTCCTACAAGGC‑CCCTGTAGGAGCGGGCTTGCCCGCGAATGCCCCTCaaa > 1:451192/1‑130 (MQ=255) ccAGTCGACTGACCGTGGCGACTATGTGATTTCTTACGCTCTGCCGTAAAACCGCGTCGCCTGCTTCGCGGGCAAGCCCGCTCCTACAAGGC‑CCCTGTAGGAGCGGGCTTGCCCGCGAATGCCCCTCaaa < 2:451192/130‑1 (MQ=255) gTGGCGACTATGTGATTTCTTACGCTCTGCCGTAAAACCGCGTCGCCTGCTTCGCGGGCAAGCCCGCTCCTACAAGGC‑CCCTGTAGGAGCGGGCTTGCCCGCGAATGCCCCCCAAAGACAACCCGAGATTCTGTAATGCCATCCCCACtt > 2:224061/1‑150 (MQ=255) aTGTGATTTCTTACGCTCTGCCGTAAAACCGCGTCGCCTGCTTCGCGGGCAAGCCCGCTCCTACAAGGC‑CCCTGTAGGAGCGGGCTTGCCCGCGAATGCCCCTCAAAGACAACCCGAGATTCTGTAATGCCATCCGCAATTCGCATCAcc > 1:688952/1‑150 (MQ=255) aaGCCCGCTCCTACAAGGC‑CCCTGTAGGAGCGGGCTTGCCCGCGAATGCCCCTCAAAGACAACCCGAGATTCTGTAATGCCATCCGCAATTCGCATCACCCGCACCCTGCTGCTGGCCCTGGCCAGCACCACTGCCCTTGCCAGcgacga > 2:534192/1‑150 (MQ=255) tCCTACAGGGG‑CCTTGTAGGAGCGGGCTTGCCCGCGAATGCCCCTCAAAGACAACCCGAGATTCTGTAATGCCATCCGCAAtt < 1:455328/83‑1 (MQ=25) tCCTACAGGGG‑CCTTGTAGGAGCGGGCTTGCCCGCGAATGCCCCTCAAAGACAACCCGAGATTCTGTAATGCCATCCGCAAtt > 2:455328/1‑83 (MQ=25) cTACAAGGC‑CCCTGTAGGAGCTGGCTTGCCCGCGAATGCCCCTCAAAGACAACCCGAGATTCTGTAATGCCATCCGCAATTCGCATCACCCCCACCCTGCTGCTGGCCCTGGCCAGCACCACTGCCCTTGCCAGCGACGATCTCATGAcc > 1:197668/1‑150 (MQ=255) aCAAGGC‑CCCTGTAGGAGCGGGCTTGCCCGCGAATGCCCCTCAAAGACAACCCGAGATTCTGTAATGCCATCCGCAATTCGCATCACCCCCACCCTGCTGCTGGCCCTGGCCAg > 2:138172/1‑114 (MQ=255) aCAAGGC‑CCCTGTAGGAGCGGGCTTGCCCGCGAATGCCCCTCAAAGACAACCCGAGATTCTGTAATGCCATCCGCAATTCGCATCACCCCCACCCTGCTGCTGGCCCTGGCCAGCACCACTGCCCTTGCCAGCGACGATCTCATGACCCg > 1:735306/1‑150 (MQ=255) aCAAGGC‑CCCTGTAGGAGCGGGCTTGCCCGCGAATGCCCCTCAAAGACAACCCGAGATTCTGTAATGCCATCCGCAATTCGCATCACCCCCACCCTGCTGCTGGCCCTGGCCAGCACCACTGCCCTTGCCAGCGACGATCTCATGACCCg > 2:643565/1‑150 (MQ=255) aCAAGGC‑CCCTGTAGGAGCGGGCTTGCCCGCGAAGGCCCCTCAAAGACAACCCGAGATTCTGTAATGCCATCCGCAATTCGCATCACCCCCACCCTGCTGCTGGCCCTGGCCAg < 1:138172/114‑1 (MQ=255) aaGGC‑CCCTGTAGGAGCGGGCTTGCCCGCGAATGCCCCTCAAAGACAACCCGAGATTCTGTAATGCCATCCGCAATTCGCATCACCCCCACCCTGCTGCTGGCCCTGGCCAGCACCACTGCCCTTGCCAGCGACGAtctc > 1:715462/1‑140 (MQ=255) aaGGC‑CCCTGTAGGAGCGGGCTTGCCCGCGAATGCCCCTCAAAGACAACCCGAGATTCTGTAATGCCATCCGCAATTCGCATCACCCCCACCCTGCTGCTGGCCCTGGCCAGCACCACTGCCCTTGCCAGCGACGAtctc < 2:715462/140‑1 (MQ=255) aGGC‑CCCTGTAGGAGCGGGCTTGCCCGCGAATGCCCCTCAAAGACAACCCGAGATTCTGTAATGCCATc > 1:494849/1‑69 (MQ=255) aGGC‑CCCTGTAGGAGCGGGCTTGCCCGCGAATGCCCCTCAAAGACAACCCGAGATTCTGTAATGCCATc < 2:494849/69‑1 (MQ=255) | GCCGAAGACCGGCAAGCAGTACGTGGTGATCACTGCTGGCGGCGCCCGCCAGTCGACTGACCGTGGCGACTATGTGATTTCTTACGCTCTGCCGTAAAACCGCGTCGCCTGCTTCGCGGGCAAGCCCGCTCCTACAAGGC‑CCCTGTAGGAGCGGGCTTGCCCGCGAATGCCCCTCAAAGACAACCCGAGATTCTGTAATGCCATCCGCAATTCGCATCACCCCCACCCTGCTGCTGGCCCTGGCCAGCACCACTGCCCTTGCCAGCGACGATCTCATGACCCG > NC_002947/4049815‑4050097 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 22 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |