Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | NC_002947 | 3,919,380 | T→C | 10.2% | intergenic (+33/‑29) | mexB → / → mobA | multidrug resistance protein MexB/molybdenum cofactor guanylyltransferase |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_002947 | 3,919,380 | 0 | T | C | 10.2% | 98.3 / 3.2 | 39 | intergenic (+33/‑29) | mexB/mobA | multidrug resistance protein MexB/molybdenum cofactor guanylyltransferase |
Reads supporting (aligned to +/- strand): ref base T (18/17); new base C (2/2); total (20/19) | |||||||||||
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.64e-01 |
TCAGCGCCACCCTGGCGGTGGTGTTCGTGCCGGTGTTCTTCGTGGTGGTGATGCGCCTGTCCGGGCGACGCCAGGCGCATGACAGCGACGGGCAACCCGTGCCCCGCGAAAGCTGACTGGACAAGCCCCGCCTGAAATGGGAGGGTTCTCGGCAATGACTGCCCGGAACCCTTTTTCATGCCTGACGCACTCCCCCCCTGCTCCATCCTCATTCTCGCCGGTGGCCGCGGCCAGCGCATGGGCGGCCGTGACAAGGGCCTGGTGGACTGGCAAGGCGAGCCGTTGATTGCAC > NC_002947/3919232‑3919523 | tCAGCGCCACCCTGGCGGTGGTGTTCGTGCCGGTGTTCTTCGTGGTGGTGATGCGCCTGTCCGGGCGACGCCAGGCGCATGACAGCGACGGGCAACCCGTGCCCCGCGAAAGCTGACTGGACAAGCCCCGCCTGAAATGGGAGGGTtct > 2:7635/1‑149 (MQ=255) ccAGCGCCACCCTGGCGGTGGTGTTCGTGCCGGTGTTCTTCGTGGTGGTGATGCGCCTGTCCGGGCGACGCCAGGCGCATGACAGCGACGGGCAACCCGTGCCCCGCGAAAGCTGACTGGACAAGCCCCGCCTGAAATGGGAGGGTtct > 5:142473/2‑149 (MQ=255) aGCGCCACCCTGGCGGTGGTGTTCGTGCCGGTGTTCTTCGTGGTGGTGATGCGCCTGTCCGGGCGACGCCAGGCGCATGACAGCGACGGGCAACCCGTGCCCCGCGAAAGCTGACTGGACAAGCCCCGCCTGAAATGGGAGGGTtctc < 4:79842/148‑1 (MQ=255) gCGGTGGTGTTCGTGCCGGTGTTCTTCGTGGTGGTGATGCGCCTGTCCGGGCGACGCCAGGCGCATGACAGCGACGGGCAACCCGTGCCCCGCGAAAGCTGACTGGACAAGCCCCGCCTGAAATGGGAGGGTTCTCGGCAATGACTGc < 8:24747/148‑1 (MQ=255) cGGTGGTGTTCGTGCCGGTGTTCTTCGTGGTGGTGATGCGCCTGTCCGGGCGACGCCAGGCGCATGACAGCGACGGGCAACCCGTGCCCCGCGAAAGCTGACTGGACAAGCCCCGCCTGAAATGGGAGGGTTCTCGGCAATGACTGccc < 6:142473/149‑1 (MQ=255) tgtTCGTGCCGGTGTTCTTCGTGGTGGTGATGCGCCTGTCCGGGCGACGCCAGGCGCATGACAGCGACGGGCAACCCGTGCCCCGCGAAAGCTGACTGGACAAGCCCCGCCTGAAATGGGAGGGTTCCGGGCAGTCATTGCCGAGAAcc < 6:193313/149‑1 (MQ=255) gtgtTCTTCGTGGTGGTGATGCGCCTGTCCGGGCGACGCCAGGCGCATGACAGCGACGGGCAACCCGTGCCCCGCGAAAGCTGACTGGACAAGCCCCGCCTGAAATGGGAGGGTTCTCGGCAATGACTGCCCGGAACCCTTTTTCATGc > 6:284083/1‑149 (MQ=255) gtgtTCTTCGTGGTGGTGATGCGCCTGTCCGGGCGACGCCAGGCGCATGACAGCGACGGGCAACCCGTGCCCCGCGAAAGCTGACTGGACAAGCCCCGCCTGAAATGGGAGGGTTCTCGGCAATGACTGCCCGGAACCCTTTTTCATGc > 4:223740/1‑149 (MQ=255) cttcGTGGTGGTGATGCGCCTGTCAGGGCGACGCCAGGCGCATGACAGCGACGGGCAACCCGGGCCCCGCGAAAGCTGACTGGACAAGCCCCGCCTGAAATGGGCGGGTTCTCGGCAATGACTGCCCGGACCCCTTTGTCATGCCTgcc > 8:231831/1‑147 (MQ=255) tcGTGGTGGTGATGCGCCTGTCCGGGCGACGCCAGGCGCATGACAGCGACGGGCAACCCGTGCCCCGCGAAAGCTGACTGGACAAGCCCCGCCTGAAATGGGAGGGTTCCGGGCAGTCATTGCCGAGAACCCt > 4:22205/1‑133 (MQ=255) tcGTGGTGGTGATGCGCCTGTCCGGGCGACGCCAGGCGCATGACAGCGACGGGCAACCCGTGCCCCGCGAAAGCTGACTGGACAAGCCCCGCCTGAAATGGGAGGGTTCCGGGCAGTCATTGCCGAGAACCCt < 3:22205/133‑1 (MQ=255) cGTGGTGGTGATGCGCCTGTCCGGGCGACGCCAGGCGCATGACAGCGACGGGCAACCCGTGCCCCGCGAAAGCTGACTGGACAAGCCCCGCCTGAAATGGGAGGGTTCTCGGCAATGACTGCCCGGAACCCTTTTTCATGCCTGAc < 3:181547/146‑1 (MQ=255) cGTGGTGGTGATGCGCCTGTCCGGGCGACGCCAGGCGCATGACAGCGACGGGCAACCCGTGCCCCGCGAAAGCTGACTGGACAAGCCCCGCCTGAAATGGGAGGGTTCTCGGCAATGACTGCCCGGAACCCTTTTTCATGCCTGACGc < 3:223740/148‑1 (MQ=255) tggtggtgATGCGCCTGTCCGGGCGACGCCAGGCGCATGACAGCGACGGGCAACCCGTGCCCCGCGAAAGCTGACTGGACAAGCCCCGCCTGAAATGGGAGGGTTCTCGGCAATGACTGCCCGGAACCCTTTTTCAt < 1:151849/137‑1 (MQ=255) tggtggtgATGCGCCTGTCCGGGCGACGCCAGGCGCATGACAGCGACGGGCAACCCGTGCCCCGCGAAAGCTGACTGGACAAGCCCCGCCTGAAATGGGAGGGTTCTCGGCAATGACTGCCCGGAACCCTTTTTCAt > 2:151849/1‑137 (MQ=255) ggtggtgATGCGCCTGTCCGGGCGACGCCAGGCGCATGACAGCGACGGGCAACCCGTGCCCCGCGAAAGCTGACTGGACAAGCCCCGCCTGAAATGGGAGGGTTCTCGGCAATGACTGCCCGGAACCCTTTTTCATGCCTGACGCACTc > 7:186328/1‑149 (MQ=255) gcgcCTGTCCGGGCGACGCCAGGCGCATGACAGCGACGGGCAACCCGTGCCCCGCGAAAGCTGACTGGACAAGCCCCGCCTGAAATGGGAGGGTTCTCGGCAATGACTGCCCGGAACCCTTTTTCATGCCTGACGCACTCCCCCCCTGc > 2:183605/1‑149 (MQ=255) gcCTGTCCGGGCGACGCCAGGCGCATGACAGCGACGGGCAACCCGTGCCCCGCGAAAGCTGACTGGACAAGCCCCGCCTGAAATGGGAGGGTTCTCGGCAATGACTGCCCGGAACCCTTTTTCATGCCTGACGCACTCCCCCCCTGCTc > 3:18319/1‑149 (MQ=255) ccTGTCCGGGCGACGCCAGGCGCATGACAGCGACGGGCAACCCGTGCCCCGCGAAAGCTGACTGGACAAGCCCCGCCTGAAATGGGAGGGTTCTCGGCAATGACTGCCCGGAACCCTTTTTCATGCCTGACGCACTCCCCCCCTGCTc > 2:191846/1‑148 (MQ=255) cGGGCGACGCCAGGCGCATGACAGCGACGGGCAACCCGTGCCCCGCGAAAGCTGACTGGACAAGCCCCGCCTGAAATGGGAGGGTTCTCGGCAATGACTGCCCGGAACCCTTTTTCATGCCTGACGCACTCCCCCCCTGCTCCATCCt > 3:217656/1‑148 (MQ=255) gggCGACGCCAGGCGCATGACAGCGACGGGCAACCCGTGCCCCGCGAAAGCTGACTGGACAAGCCCCGCCTGAAATGGGAGGGTTCTCGGCAATGACTGCCCGGAACCCTTTTTCATGCCTGACGCACTCCCCCCCTGCTCCATCCTCa > 7:88217/1‑149 (MQ=255) cGACGCCAGGCGCATGACAGCGACGGGCAACCCGTGCCCCGCGAAAGCTGACTGGACAAGCCCCGCCTGAAATGGGAGGGTTCTCGGCAATGACTGCCCGGAACCCTTTTTCATGCCTGACGCACTCCCCCCCTGCTCCATCCTCATtc < 8:186328/149‑1 (MQ=255) gCCAGGCGCATGACAGCGACGGGCAACCCGTGCCCCGCGAAAGCTGACTGGACAAGCCCCGCCTGAAATGGGAGGGTTCTCGGCAATGACTGCCCGGAACCCTTTTTCATGCCTGACGCACTCCCCCCCTGCTCCATCCTc < 6:4792/141‑1 (MQ=255) gCCAGGCGCATGACAGCGACGGGCAACCCGTGCCCCGCGAAAGCTGACTGGACAAGCCCCGCCTGAAATGGGAGGGTTCTCGGCAATGACTGCCCGGAACCCTTTTTCATGCCTGACGCACTCCCCCCCTGCTCCATCCTc > 5:4792/1‑141 (MQ=255) cGACGGGCAACCCGTGCCCCGCGAAAGCTGACTGGACAAGCCCCGCCTGAAATGGGAGGGTTCTCGGCAATGACTGCCCGGAACCCTTTTTCATGCCTGACGCACTCCCCCCCTGCTCCATCCTCATTCTCGCCGGTGGCCGCGGCCAg < 7:212040/149‑1 (MQ=255) gCACCCCGTGCCCCGCGAAAGCTGACTGGACAAGCCCCGCCTGAAATGGGAGGGTTCTCGGCAATGCCTGCCGGGCCCCCTTTTTCATGCCTGACGCACTCCCCCCCTGCACCATCCTCATTCTCCCCGGTGGCCGCCGACAGCGCATg > 5:95273/1‑149 (MQ=255) gCAACCCGTGCCCCGCGAAAGCTGACTGGACAAGCCCCGCCTGAAATGGGAGGGTTCTCGGCAATGACTGCCCGGAACCCTTTTTCATGCCTGACGCACTCCCCCCCTGCTCCATCCTCATTCTCGCCGGTGGCCGCGGCCAGCGCATg > 4:176997/1‑149 (MQ=255) aaCCCGTGCCCCGCGAAAGCTGACTGGACAAGCCCCGCCTGAAATGGGAGGGTTCTCGGCAATGACTGCCCGGAACCCTTTTTCATGCCTGACGCACTCCCCCCCTGCTCCATCCTCATTCTCGCCGGTGGCCGCGGCCAGCGCATggg < 4:18319/149‑1 (MQ=255) aaCCCGTGCCCCGCGAAAGCTGACTGGACAAGCCCCGCCTGAAATGGGAGGGTTCTCGGCAATGACTGCCCGGAACCCTTTTTCATGCCTGACGCACTCCCCCCCTGCTCCATCCTCATTCTCGCCGGTGGCCGCGGCCAGCGCATggg > 2:261049/1‑149 (MQ=255) ccccGCGAAAGCTGACTGGACAAGCCCCGCCTGAAATGGGAGGGTTCTCGGCAATGACTGCCCGGAACCCTTTTTCATGCCTGACGCACTCCCCCCCTGCTCCATCCTCATTCTCGCCGGTGGCCGCGGCCAGCGCATGGGCGGCCGt < 1:261049/148‑1 (MQ=255) gcgAAAGCTGACTGGACAAGCCCCGCCTGAAATGGGAGGGTTCTCGGCAATGACTGCCCGGAACCCTTTTTCATGCCTGACGCACTCCCCCCCTGCTCCATCCTCATTCTCGCCGGTGGCCGCGGCCAGCGCATGGGCGGCCGTGACaa < 5:284083/149‑1 (MQ=255) cgAAAGCTGACTGGACAAGCCCCGCCTGAAATGGGAGGGTTCTCGGCAATGACTGCCCGGAACCCTTTTTCATGCCTGACGCACTCCCCCCCTGCTCCATCCTCATTCTCGCCGGTGGCCGCGGCCAGCGCATGGGCGGCCGTGACAAg < 1:191846/149‑1 (MQ=255) aaGCTGACTGGACAAGCCCCGCCTGAAATGGGAGGGTTCTCGGCAATGACTGCCCGGAACCCTTTTTCATGCCTGACGCACTCCCCCCCTGCTCCATCCTCATTCTCGCCGGTGGCCGCGGCCAGCGCATGGGCGGCCGTGACAAGGGc < 3:176997/149‑1 (MQ=255) cTGACTGGAAAAGCCCCGCCTGAAATGGGAGGGTTCTCGGCAATGACTGCCCGGAACACTTTTTCATGACTGACGCACTCCCCCCCTGCTCCATCCACATACACGCCGGTGGCCGCGGCCAGCGCATGCGCGGCCGTGACAAGGGCCtg > 2:152530/1‑149 (MQ=255) aCTGGACAAGCCCCGCCTGAAATGGGAGGGTTCTCGGCAATGACTGCCCGGAACCCTTTTTCATGCCTGACGCACTCCCCCCCTGCTCCATCCTCATTCTCGCCGGTGGCCGCGGCCAGCGCATGGGCGGCCGTGACAAGGGCCtggtg < 8:88217/149‑1 (MQ=255) aCTCGACAAGCGCCGCCTGAAATGGGAGGGTTCTCGGCAATGACTGCCCGGAACCCTTTTTCATGCCTGACGCACTCCCCCCCTGCTCCATCCTGATTCTCGCCGGTGGCAGCGGCCAGCGCATGGGCGGCCGTGACAAGGGCCtggt < 1:152530/148‑1 (MQ=255) tGGACAAGCCCCGCCTGAAATGGGAGGGTTCTCGGCAATGACTGCCCGGAACCCTTTTTCATGCCTGACGCACTCCCCCCCTGCTCCATCCTCATTCTCGCCGGTGGCCGCGGCCAGCGCATGGGCGGCCGTg < 2:60944/133‑1 (MQ=255) tGGACAAGCCCCGCCTGAAATGGGAGGGTTCTCGGCAATGACTGCCCGGAACCCTTTTTCATGCCTGACGCACTCCCCCCCTGCTCCATCCTCATTCTCGCCGGTGGCCGCGGCCAGCGCATGGGCGGCCGTg > 1:60944/1‑133 (MQ=255) gCCCCGCCTGAAATGGGAGGGTTCTCGGCAATGACTGCCCGGAACCCTTTTTCATGCCTGACGCACTCCCACCCTGCTCCATCCTCATTATCGCCGGTGGCCGCGGCCAGCGCATGGGCGGCAGTGACAAGCGCCTGTTGGACTGGc > 6:205149/1‑147 (MQ=255) ccGCCTGAAATGGGAGGGTTCTCGGCAATGACTGCCCGGAACCCTTTTTCATGCCTGACGCACTCCCCCCCTGCTCCATCCTCATTCTCGCCGGTGGCCGCGGCCAGCGCATGGGCGGCCGTGACAAGGGCCTGGTGGACTGGCAAGGc > 6:171807/1‑149 (MQ=255) gaaaaaGGGTTCCGGGCAGTCATTGCCGAGAACCCTTTTTCATGCCTGACGCACTCCCCCCCTGCTCCATCCTCATTCTCGCCGGTGGCCGCGGCCAGCGCATGGGCGGCCGTGACAAGGGCCTGGTGGACTGGCAAGGCGAGCCGtt > 4:254742/6‑148 (MQ=38) ggTTCTCGGCAATGACTGCCCGGAACCCTTTTTCATGCCTGACGCACTCCCCCCCTGCTCCATCCTCATTCTCGCCGGTGGCCGCGGCCAGCGCATGGGCGGCCGTGACAAGGGCCTGGTGGACTGGCAAGGCGAGCCGTTGATTGcac < 2:105486/149‑1 (MQ=255) | TCAGCGCCACCCTGGCGGTGGTGTTCGTGCCGGTGTTCTTCGTGGTGGTGATGCGCCTGTCCGGGCGACGCCAGGCGCATGACAGCGACGGGCAACCCGTGCCCCGCGAAAGCTGACTGGACAAGCCCCGCCTGAAATGGGAGGGTTCTCGGCAATGACTGCCCGGAACCCTTTTTCATGCCTGACGCACTCCCCCCCTGCTCCATCCTCATTCTCGCCGGTGGCCGCGGCCAGCGCATGGGCGGCCGTGACAAGGGCCTGGTGGACTGGCAAGGCGAGCCGTTGATTGCAC > NC_002947/3919232‑3919523 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 21 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |