| Predicted mutation | |||||||
|---|---|---|---|---|---|---|---|
| evidence | seq id | position | mutation | freq | annotation | gene | description |
| RA | NC_000913 | 2,291,284 | C→A | 14.3% | intergenic (+137/+74) | narP → / ← ccmH | response regulator in two‑component regulatory system with NarQ/heme lyase, CcmH subunit |
| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | NC_000913 | 2,291,284 | 0 | C | A | 14.3% | 108.7 / 9.5 | 35 | intergenic (+137/+74) | narP/ccmH | response regulator in two‑component regulatory system with NarQ/heme lyase, CcmH subunit |
| Reads supporting (aligned to +/- strand): ref base C (19/11); new base A (2/3); total (21/14) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 3.69e-01 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.59e-01 | |||||||||||
AAATAGCCCGATGGATTTATCATCGGGCTGAGATTTATGACAAACGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGATCGATCAATTTATTGAATTTACACAATTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATAAACTAAGCGCACTTTGTCAACAGTCTAGCCCGATGGCATCACCATCGGGCCTCTTTTTATTTACTCTCCTGCGGCGACAAATGTTGCATCGCCTGCGCGATACTACGTTCAATCAC > NC_000913/2291148‑2291417 | aaaTAGCCCGATAGATGTATCATCGAGTTGAGAATTATGACAAACGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGATCGGTCAATTTATTGAATTTACACAATTTTCTAGGCCGGATGCGGCGTGAACGCCTTATCCg < 4:15148/149‑1 (MQ=18) gATGGATTTATCATCGGGCTGAGATTTATGACAAACGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGATCGATCAATTTATTGAATTTACACAATTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATAAACt > 4:703258/1‑149 (MQ=255) agaTTTATGACAAACGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGATCGATCAATTTATTGAATTTACACAATTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATAAACTAAGCGCACTTTGTCAACAGt < 1:68249/148‑1 (MQ=255) aCAAACGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGATCGATCAATTTATTGAATTTACACAATTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATAAACTAAGCGCACTTTGTCAACAGTCTAGCCCGAt < 4:841677/149‑1 (MQ=255) gCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGATCGATCAATTTATTGAATTTACACAATTTTGTAGGCCGGATAAGGCGTTAACGCCGCATCCGGCATAAACTAAGCGCACTTTGTCAACAGTCTAGCCCGAt > 4:419390/1‑143 (MQ=255) gCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGATCGATCAATTTATTGAATTTACACAATTTTGTAGGCCGGATAAGGCGTTAACGCCGCATCCGGCATAAACTAAGCGCACTTTGTCAACAGTCTAGCCCGAt < 3:419390/143‑1 (MQ=255) gcTACGCTTATCAGGCCTACGTGGATCGATCAATTTATTGAATTTACACAATTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATAAACTAAGCGCACTTTGTCAACAGTCTAGCCCGATGGCATCACCATCGGGCCTCtttt > 4:735674/1‑149 (MQ=255) cTACGCTTATCAGGCCTACGTGGATCGATCAATTTATTGAATTTACACAATTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATAAACTAAGCGCACTTTGTCAACAGTCTAGCCCGATGGCATCACCATCGGGCCTCtttt > 1:3343/1‑148 (MQ=255) cTACGCTTATCAGGCCTACGTGGATCGATCAATTTATTGAATTTACACAATTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATAAACTAAGCGCACTTTGTCAACAGTCTAGCCCGATGGCATCACCATCGGGCCTCtttt > 2:3150/1‑148 (MQ=255) tACGCTTATCAGGCCTACGTGGATCGATCAATTTATTGAATTTACACAATTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATAAACTAAGCGCACTTTGTCAACAGTCTAGCCCGATGGCATCACCATCGGGCCTCttttt > 2:524667/1‑148 (MQ=255) aCGCTTATCAGGCCTACGTGGATCGATCAATTTATTGAATTTACACAATTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATAAACTAAGCGCACTTTGTCAACAGTCTAGCCCGATGGCATCACCATCGGGCCTCtttt > 4:332926/1‑146 (MQ=255) aCGCTTATCAGGCCTACGTGGATCGATCAATTTATTGAATTTACACAATTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATAAACTAAGCGCACTTTGTCAACAGTCTAGCCCGATGGCATCACCATCGGGCCTCtttt < 3:332926/146‑1 (MQ=255) aCGCTTATCAGGCCTACGTGGATCGATCAATTTATTGAATTTACACAATTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATAAACTAAGCGCACTTTGTCAACAGTCTAGCCCGATGGCATCACCATCGGGCCTCTTtttat < 3:114529/149‑1 (MQ=255) gCTTATCAGGCCTACGTGGATCGATCAATTTATTGAATTTACACAATTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATAAACTAAGCGCACTTTGTCAACAGTCTAGCCCGATGGCATCACCATCGGGCCTCTTtttattt > 4:517166/1‑149 (MQ=255) cTTATCAGGCCTACGTGGATCGATCAATTTATTGAATTTACACAATTTTGTAGGCCGGATGCGGCGTGAACGCCTTATCCgg < 2:741271/82‑1 (MQ=25) cTTATCAGGCCTACGTGGATCGATCAATTTATTGAATTTACACAATTTTGTAGGCCGGATGCGGCGTGAACGCCTTATCCgg > 1:741271/1‑82 (MQ=21) gcaTCAGGCCTACGTGGATCGATCAATTTATTGAATTTACACAATTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATAAACTAAGCGCACTTTGTCAACAGTCTAGCCCGATGGCATCACCATCGGGCCTCTTtttattta < 1:257625/146‑1 (MQ=255) gTGGATCGATCAATTTATTGAATTTACACAATTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATAAACTAAGCGCACTTTGTCAACAGTCTAGCCCGATGGCATCACCATCGGGCCTCTTTTTATTTACTCTCCTgcggcg < 1:832778/148‑1 (MQ=255) atcgatcAATTTATTGAATTTACACAATTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATAAACTAAGCGCACTTTGTCAACAGTCTAGCCCGATGGc < 3:524090/105‑1 (MQ=255) atcgatcAATTTATTGAATTTACACAATTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATAAACTAAGCGCACTTTGTCAACAGTCTAGCCCGATGGc > 4:524090/1‑105 (MQ=255) gatcAATTTATTGAATTTACACAATTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATAAACTAAGCGCACTTTGTCAACAGTCTAGCCCGATGGCATCACCATCGGGCCTCTTTTTATTTACTCTCCTGCGGCGACAAATg < 1:524667/148‑1 (MQ=255) tcAATTTATTGAATTTACACAATTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATAAACTAAGCGCACTTTGTCAACAGTCTAGCCCGATGGCATCACCATCGGGCCTCTTTTTATTTACTCTCCTGCGGCGACAAATGTTg > 3:454855/1‑149 (MQ=255) aTTTATTGAATTTACACAATTTTGTAGGCCGGATAAGGCGTTCACGCCGCATTCGGCATAAACTAAGCGCACTTTGTCAACAGTCTAGCCCGATGGCATCACCATCGGGCCTCTTTTTATTTACTCTCCTGCGGCGACAAATGTTGCAt > 3:77144/1‑149 (MQ=255) aTTTATTGAATTTACACAATTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATAAACTAAGCGCACTTTGTCAACAGTCTAGCCCGATGGCATCACCATCGGGCCTCTTTTTATTTACTCTCCTGCGGCGACAAATGTTGCAt > 4:713228/1‑149 (MQ=255) aTTTACACAATTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATAAACTAAGCGCACTTTGTCAACAGTCTAGCCCGATGGCATCACCATCGGGCCTCTTTTTATTTACTCTCCTGCGGCGACAAATGTTGCATCGCCTgcgc > 4:55345/1‑149 (MQ=255) aTTTACACAATTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATAAACTAAGCGCACTTTGTCAACAGTCTAGCCCGATGGCATCACCATCGGGCCTCTTTTTATTTACTCTCCTGCGGCGACAAATGTTGCATCGCCTgcg > 1:464837/1‑148 (MQ=255) ttACACAATTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATAAACTAAGCGCACTTTGTCAACAGTCTAGCCCGATGGCATCACCATCGGGCCTCTTTTTATTTACTCTCCTGCGGCGACAAATGTTGCATCGCCTGcgcg < 1:870377/148‑1 (MQ=255) ttACACAATTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATAAACTAAGCGCACTTTGTCAACAGTCTAGCCCGATGGCATCACCATCGGGCCTCTTTTTATTTACTCTCCTGCGGCGACAAATGTTGCATCGCCTGcgcg > 1:817286/1‑148 (MQ=255) cacaATTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATAAACTAAGCGCACTTTGTCAACAGTCTAGCCCGATGGCATCACCATCGGGCCTCTTTTTATTTACTCTCCTGCGGCGACAAATGTTGCATCGCCTGCGCGAta > 3:647764/1‑148 (MQ=255) cacaATTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATAAACTAAGCGCACTTTGTCAACAGTCTAGCCCGATGGCATCACCATCGGGCCTCTTTTTATTTACTCTCCTGCGGCGACAAATGTTGCATCGCCTGCGCGAta < 4:647764/148‑1 (MQ=255) acaATTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATAAACTAAGCGCACTTTGTCAACAGTCTAGCCCGATGGCATCACCATCGGGCCTCTTTTTATTTACTCTCCTGCGGCGACAAATGTTGCATCGCCTGCGCGAtac > 2:663524/1‑148 (MQ=255) aaTTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATAAACTAAGCGCACTTTGTCAACAGTCTAGCCCGATGGCATCACCATCGGGCCTCTTTTTATTTACTCTCCTGCGGCGACAAATGTTGCATCGCCTGCGCGAtacta < 2:817286/148‑1 (MQ=255) aTTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATAAACTAAGCGCACTTTGTCAACAGTCTAGCCCGATGGCATCACCATCGGGCCTCTTTTTATTTACTCTCCTGCGGCGACAAATGTTGCATCGCCTGCGCGAtactac > 2:492506/1‑148 (MQ=255) tttGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATAAACTAAGCGCACTTTGTCAACAGTCTAGCCCGATGGCATCACCATCGGGCCTCTTTTTATTGACTCTCCTGCGGCGACAAATGTTGCATCGCCTGCGCGATACTACGt > 1:267391/1‑148 (MQ=255) ccGGATAAGGCGTTCACGCCGCATCCGGCATAAACTAAGCGCACTTTGTCAACAGTCTAGCCCGATGGCATCACCATCGGGCCTCTTTTTATTTACTCTCCTGCGGCGACAAATGTTGCATCGCCTGCGCGATACTACGTTCAATCAc > 2:878450/1‑148 (MQ=255) | AAATAGCCCGATGGATTTATCATCGGGCTGAGATTTATGACAAACGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGATCGATCAATTTATTGAATTTACACAATTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATAAACTAAGCGCACTTTGTCAACAGTCTAGCCCGATGGCATCACCATCGGGCCTCTTTTTATTTACTCTCCTGCGGCGACAAATGTTGCATCGCCTGCGCGATACTACGTTCAATCAC > NC_000913/2291148‑2291417 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |