Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | NC_002947 | 4,049,954 | C→G | 18.2% | intergenic (+43/‑58) | 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,954 | 0 | C | G | 18.2% | 43.4 / 6.3 | 22 | intergenic (+43/‑58) | quiA/oprB‑III | quinate dehydrogenase/carbohydrate‑selective porin |
Reads supporting (aligned to +/- strand): ref base C (13/5); new base G (2/2); total (15/7) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 5.65e-01 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.83e-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: ‑ |