Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | NC_002947 | 538,400 | A→G | 14.6% | I259V (ATC→GTC) | rpoB → | DNA‑directed RNA polymerase subunit beta |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_002947 | 538,400 | 0 | A | G | 14.6% | 96.3 / 8.8 | 41 | I259V (ATC→GTC) | rpoB | DNA‑directed RNA polymerase subunit beta |
Reads supporting (aligned to +/- strand): ref base A (19/16); new base G (3/3); total (22/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 = 9.20e-01 |
GCCTCGGTGCTGCTGCGCGCCCTGGGTTACAGCACTGAAGAAGTGCTGAACACCTTCTACACCACCAACGTGTTCCACATTTCCGGCGAAAAACTCAGCCTGGAACTGGTACCGCAGCGTCTGCGTGGTGAAGTTGCAGTCATGGACATCCATGACGAAACCGGCAAAGTCATCGTCGAGCAAGGCCGCCGTATTACCGCGCGCCACATCAACCAGCTCGAGAAGGCCGGCGTCAAGCAACTGGACGTTCCAATGGAATACGTCCTGGGCCGCACTACCGCCAAGGCTATC > NC_002947/538253‑538543 | gCCTCGGTGCTGCTGCGCGCCCTGGGTTACAGCACTGAAGAAGTGCTGAACACCTTCTACACCACCAACGTGTTCCACATTTCCGGCGAAAAACTCAGCCTGGAACTGGTACCGCAGCGTCTGCGTGGTGAAGTTGCAGTCATGGACATc > 2:509733/1‑150 (MQ=255) cTCGGTGCTGCTGCGCGCCCTGGGTTACAGCACTGAAGAAGTGCTGAACACCTTCTACACCACCAACGTGTTCCACATTTCCGGCGAAAAACTCAGCCTGGAACTGGTACCGCAGCGTCTGCGTGGTGAAGTTGCAGTCATGGACATCCa > 2:468933/1‑150 (MQ=255) ggTGCTGCTGCGCGCCCTGGGTTACAGCACTGAAGAAGTGCTGAACACCTTCTACACCACCAACGTGTTCCACATTTCCGGCGAAAAACTCAGCCTGGAACTGGTACCGCAGCGTCTGCGTGGTGAAGTTGCAGTCATGGATGTCCATGa > 2:485752/1‑150 (MQ=255) ggTGCTGCTGCGCGCCCTGGGTTACAGCACTGAAGAAGTGCTGAACACCTTCTACACCACCAACGTGTTCCACATTTCCGGCGAAAAACTCAGCCTGGAACTGGTACCGCAGCGTCTGCGTGGTGAAGTTGCAGTCATGGACATCCATGa < 1:468933/150‑1 (MQ=255) tgctgctgcGCGCCCTGGGTTACAGCACTGAAGAAGTGCTGAACACCTTCTACACCACCAACGTGTTCCACATTTCCGGCGAAAAACTCAGCCTGGAACTGGTACCGCAGCGTCTGCGTGGTGAAGTTGCAGTCATGGATGTCCATGAct < 1:485752/150‑2 (MQ=255) ccTGGGTTACAGCACTGAAGAAGTGCTGAACACCTTCTACACCACCAACGTGTTCCACATTTCCGGGGAAAAACTCAGCCTGGAACTGGTACCGCAGCGTCTGCGTGGGGAAGTTGGAGGAATGGACATACATAACGAAACCGGCAAGGt > 2:815192/1‑150 (MQ=255) aGCACTGAAGAAGTGCTGAACACCTTCTACACCACCAACGTGTTCCACATTTCCGGCGAAAAACTCAGCCTGGAACTGGTACCGCAGCGTCTGCGTGGTGAAGTTGCAGTCATGGACATCCATGACGAAACCGGCAAAGTCATCGTCGAg > 1:277405/1‑150 (MQ=255) gaagaaGTGCTGAACACCTTCTACACCACCAACGTGTTCCACATTTCCGGCGAAAAACTCAGCCTGGAACTGGTACCGCAGCGTCTGCGTGGTGAAGTTGCAGTCATGGATGTCCATGActgc < 1:507062/123‑4 (MQ=255) gaagaaGTGCTGAACACCTTCTACACCACCAACGTGTTCCACATTTCCGGCGAAAAACTCAGCCTGGAACTGGTACCGCAGCGTCTGCGTGGTGAAGTTGCAGTCATGGATGTCCATGActgc > 2:507062/1‑120 (MQ=255) gTGCTGAACACCTTCTACACCACCAACGTGTTCCACATTTCCGGCGAAAAACTCAGCCTGGAACTGGTACCGCAGCGTCTGCGTGGTGAAGTTGCAGTCATGGACATCCa > 2:7400/1‑110 (MQ=255) gTGCTGAACACCTTCTACACCACCAACGTGTTCCACATTTCCGGCGAAAAACTCAGCCTGGAACTGGTACCGCAGCGTCTGCGTGGTGAAGTTGCAGTCATGGACATCCa < 1:7400/110‑1 (MQ=255) acCTTCTACACCACCAACGTGTTCCACATTTCCGGCGAAAAACTCAGCCTGGAACTGGTACCGCAGCGTCTGCGTGGTGAAGTTGCAGTCATGGTCATCCATGACGAAACCGGCAAAGTCAtcgtc < 1:582249/126‑1 (MQ=255) acCTTCTACACCACCAACGTGTTCCACATTTCCGGCGAAAAACTCAGCCTGGAACTGGTACCGCAGCGTCTGCGTGGTGAAGTTGCAGTCATGGACATCCATGACGAAACCGGCAAAGTCAtcgtc > 2:582249/1‑126 (MQ=255) acCTTCTACACCACCAACGTGTTCCACATTTCCGGCGAAAAACTCAGCCTGGAACTGGTACCGCAGCGTCTGCGTGGTGAAGTTGCAGTCATGGACATCCATGACGAAACCGGCAAAGTCATCGTCGAGCAAGGCCGCCGTATTACcgcg > 2:583931/1‑150 (MQ=255) ccTTCTCCACCCCCCACTTGTTCCACATCTCCGGCGAAAAACTCAGCCTGGAACTGGTACCGCAGCGTCTGCGTGGTGAAGTTGCAGTCATGGACATCCATGACGAAACCGGCAAAGTCATCGTCGAGCAAGGCCGCCGTATTACcgcgc < 1:815192/150‑1 (MQ=255) ttCTACACCACCAACGTGTTCCACATTTCCGGCGAAAAACTCAGCCTGGAACTGGTACCGCAGCGTCTGCGTGGTGAAGTTGCAGTCATGGACATCCATGACGAAACCGGCAAAGTCATCGTCGAGCAAGGCCGCCGTATTACCgcgcgc > 2:697428/1‑150 (MQ=255) caccacCAACGTGTTCCACATTTCCGGCGAAAAACTCAGCCTGGAACTGGTACCGCAGCGTCTGCGTGGTGAAGTTGCAGTCATGGACATCCATGACGAAACCGGCAAAGTCATCGTCGAGCAAGgcc > 2:702251/1‑128 (MQ=255) caccacCAACGTGTTCCACATTTCCGGCGAAAAACTCAGCCTGGAACTGGTACCGCAGCGTCTGCGTGGTGAAGTTGCAGTCATGGACATCCATGACGAAACCGGCAAAGTCATCGTCGAGCAAGgcc < 1:702251/128‑1 (MQ=255) ccaccaACGTGTTCCACATTTCCGGCGAAAAACTCAGCCTGGAACTGGTACCGCAGCGTCTGCGTGGTGAAGTTGCAGTCATGGACATCCATGACGAAACCGGCAAAGTCATCGTCGAGCAAGGCCGCCGTATTACCGCGCGCCACATCa > 1:766329/1‑150 (MQ=255) caccaACGTGTTCCACATTTCCGGCGAAAAACTCAGCCTGGAACTGGTACCGCAGCGTCTGCGTGGTGAAGTTGCAGTCATGGACATCCATGACGAAACCGGCAAAGTCATCGTCGAGCAAGGCCGCCGTATTACcgcgc > 2:24562/1‑140 (MQ=255) caccaACGTGTTCCACATTTCCGGCGAAAAACTCAGCCTGGAACTGGTACCGCAGCGTCTGCGTGGTGAAGTTGCAGTCATGGACATCCATGACGAAACCGGCAAAGTCATCGTCGAGCAAGGCCGCCGTATTACcgcgc < 1:24562/140‑1 (MQ=255) accaACGTGTTCCACATTTCCGGCGAAAAACTCAGCCTGGAACTGGTACCGCAGCGTCTGCGTGGTGAAGTTGCAGTCATGGACATCCAt > 1:102217/1‑90 (MQ=255) accaACGTGTTCCACATTTCCGGCGAAAAACTCAGCCTGGAACTGGTACCGCAGCGTCTGCGTGGTGAAGTTGCAGTCATGGACATCCAt < 2:102217/90‑1 (MQ=255) caACGTGTTCCACATTTCCGGCGAAAAACTCAGCCTGGAACTGGTACCGCAGCGTCTGCGTGGTGAAGTTGCAGTCATGGACATCCATGACGAAACCGGCAAAGTCATCGTCGAGCAAGGCCGCCGTATTACCGCGCGCCACATCAACCa > 2:507070/1‑150 (MQ=255) gtgtTCCACATTTCCGGCGAAAAACTCAGCCTGGAACTGGTACCGCAGCGTCTGCGTGGTGAAGTTGCAGTCATGGACATCCATGACGAAACCGGCAAAGTCATCGTCGAGCAAGGCCGCCGTATTACCGCGCGCCACATCAACCAGCTc < 1:583931/150‑1 (MQ=255) ccACATTTCCGGCGAAAAACTCAGCCTGGAACTGGTACCGCAGCGTCTGCGTGGTGAAGTTGCAGTCATGGACATCCATGACGAAACCGGCAAAGTCATCGTCGAGCAAGGCCGCCGTATTACCGCGCGCCACATCAACCAGCTCGAGaa < 1:369963/150‑1 (MQ=255) cacaTTTCCGGCGAAAAACTCAGCCTGGAACTGGTACCGCAGCGTCTGCGTGGTGAAGTTGCAGTCATGGACATCCATGACGAAACCGGCAAAGTCATCGTCGAGCAAGGCCGCCGTATTACCGCGCGCCACATCAACCAGCTCGAGAAg > 2:137164/1‑150 (MQ=255) cGAAAAACTCAGCCTGGAACTGGTACCGCAGCGTCTGCGTGGTGAAGTTGCAGTCATGGACATCCATGACGAAACCGGCAAAGTCATCGTCGAGCAAGGCCGCCGTATTACCGCGCGCCACATCAACCAGCTCGAGAAGGCCGGCGTCaa < 2:277405/150‑1 (MQ=255) aCTCAGCCTGGAACTGGTACCGCAGCGTCTGCGTGGTGAAGTTGCAGTCATGGACATCCATGACGAAACCGGCAAAGTCATCGTCGAGCAAGGCCGCCGTATTACCGCGCGCCACATCAACCAGCTCGAGAAGGCCGGCGTCAAGCAACt < 1:507070/150‑1 (MQ=255) gCCTGGAACTGGTACCGCAGCGTCTGCGTGGTGAAGTTGCAGTCATGGACATCCATGACGAAACCGGCAAAGTCATCGTCGAGCAAGGCCGCCGTa > 1:762087/1‑96 (MQ=255) gCCTGGAACTGGTACCGCAGCGTCTGCGTGGTGAAGTTGCAGTCATGGACATCCATGACGAAACCGGCAAAGTCATCGTCGAGCAAGGCCGCCGTa < 2:762087/96‑1 (MQ=255) tGGAACTGGTACCGCAGCGTCTGCGTGGTGAAGTTGCAGTCATGGACATCCATGACGAAACCGGCAAAGTCATCGTCGAGCAAGGCCGCCGTATTACCGCGCGCCACATCAACCAGCTCGAGAAGGCCGGCGTCAAGCAACTGGACGTTc < 1:137164/150‑1 (MQ=255) aCTGGTACCGCAGCGTCTGCGTGGTGAAGTTGCAGTCATGGACATCCATGACGAAACCGGCAAAGTCATCGTCGAGCAAGGCCGCCGTATTACCGCGCGCCACATCAACCAGCTCGAGAAGGCCGGCGTCAAGCAACTGGACGTTCCAAt > 1:757296/1‑150 (MQ=255) aCCGCAGCGTCTGCGTGGTGAAGTTGCAGTCATGGACATCCATGACGAAACCGGCAAAGTCATCGTCGAGCAAGGCCGCCGTATTACCGCGCGCCACATCAACCAGCTCGAGAAGGCCGGCGTCAAGCAACTGGACGTTCCAATGGAATa > 1:106922/1‑150 (MQ=255) gCAGCGTCTGCGTGGTGAAGTTGCAGTCATGGACATCCATGACGAAACCGGCAAAGTCATCGTCGAGCAAGGCCGCCGTATTACCGCGCGCCACATCAACCAGCTCGAGAAGGCCGGCGTCAAGCAACTGGACGTTCCAATGGAATACGt > 1:18280/1‑150 (MQ=255) aGCGTCTGCGTGGTGAAGTTGCAGTCATGGACATCCATGACGAAACCGGCAAAGTCATCGTCGAGCAAGGCCGCCGTATTACCGCGCGCCACATCAACCAGCTCGAGAAGGCCGGCGTCAAGCAACTGGACGTTCCAATGGAATACGTcc < 1:509733/150‑1 (MQ=255) tGCAGTCATGGACATCCATGACGAAACCGGCAAAGTCATCGTCGAGCAAGGCCGCCGTATTACCGCGCGCCACAACAACCAGCTCGAGAAGGCCGGCGTCAAGCAACTGGACGTTCCAATGGAATACGTCCTGGGCCGCACTACCGCCaa < 1:697428/150‑1 (MQ=255) gCAGTCATGGACATCCATGACGAAACCGGCAAAGTCATCGTCGAGCAAGGCCGCCGTATTACCGCGCGCCACATCAACCAGCTCGAGAAGGCCGGCGTCAAGCAACTGGACGTTCCAATGGAATACGTCCTGGGCCGCACTACCGCCAAg < 2:18280/150‑1 (MQ=255) tCATGGATGTCCATGACGAAACCGGCAAAGTCATCGTCGAGCAAGGCCGCCGTATTACCGCGCGCCACATCAACCAGCTCGAGAAgg > 1:390304/1‑87 (MQ=255) tCATGGATGTCCATGACGAAACCGGCAAAGTCATCGTCGAGCAAGGCCGCCGTATTACCGCGCGCCACATCAACCAGCTCGAGAAgg < 2:390304/87‑1 (MQ=255) aTGGACATCCATGACGAAACCGGCAAAGTCATCGTCGAGCAAGGCCGCCGTATTACCGCGCGCCACATCAACCAGCTCGAGAAGGCCGGCGTCAAGCAACTGGACGTTCCAATGGAATACGTCCTGAGCAGCACTACGGCCAAGGCTaag > 2:688614/1‑148 (MQ=255) | GCCTCGGTGCTGCTGCGCGCCCTGGGTTACAGCACTGAAGAAGTGCTGAACACCTTCTACACCACCAACGTGTTCCACATTTCCGGCGAAAAACTCAGCCTGGAACTGGTACCGCAGCGTCTGCGTGGTGAAGTTGCAGTCATGGACATCCATGACGAAACCGGCAAAGTCATCGTCGAGCAAGGCCGCCGTATTACCGCGCGCCACATCAACCAGCTCGAGAAGGCCGGCGTCAAGCAACTGGACGTTCCAATGGAATACGTCCTGGGCCGCACTACCGCCAAGGCTATC > NC_002947/538253‑538543 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 22 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |