Predicted mutation | ||||||
---|---|---|---|---|---|---|
evidence | seq id | position | mutation | annotation | gene | description |
RA | NZ_CP009273 | 2,859,998 | T→C | I305V (ATT→GTT) | rpoS ← | RNA polymerase sigma factor RpoS |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NZ_CP009273 | 2,859,998 | 0 | T | C | 92.6% | 71.6 / ‑0.6 | 27 | I305V (ATT→GTT) | rpoS | RNA polymerase sigma factor RpoS |
Reads supporting (aligned to +/- strand): ref base T (0/2); new base C (11/14); total (11/16) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 4.99e-01 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.64e-01 |
TACTTACTCGCGGAACAGCGCTTCGATATTCAGCCCCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCGAACACGTTCACGGGTGAGGCCAATTTCACGACCTACATCTTCCAGTGTTGCCGCTTCGTACCCCAGCAAACCGAATCGAC > NZ_CP009273/2859915‑2860084 | tacttactCGCGGAACAGCGCTTCGATATTCAGCCCCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAACCTGGCg > 1:212053/1‑90 (MQ=255) cttactCGCGGAACAGCGCTTCGATATTCAGCCCCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAACCTGGCGaa < 2:56808/90‑1 (MQ=255) tactCGCGGAACAGCGCTTCGATATTCAGCCCCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAACCTGGCGAaca < 1:409387/90‑1 (MQ=255) tactCGCGGAACAGCGCTTCGATATTCAGCCCCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAACCTGGCGAaca > 1:474565/1‑90 (MQ=255) tactCGCGGAACAGCGCTTCGATATTCAGCCCCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAACCTGGCGAaca < 2:575539/90‑1 (MQ=255) actCGCGGAACAGCGCTTCGATATTCAGCCCCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAACCTGGCGAacac > 1:33154/1‑90 (MQ=255) tCGATATTCAGCCCCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAACCTGGCGAaca > 2:524975/1‑72 (MQ=255) tCGATATTCAGCCCCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAACCTGGCGAaca < 1:524975/72‑1 (MQ=255) atTCAGCCCCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAACCTGGCGAACACGTTCACGGGTGAGGCCAATTTc < 1:199751/90‑1 (MQ=255) tCAGCCCCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAACCTGGCGAACACGTTCACGGGTGAGGCCAATTTCAc < 2:214323/90‑1 (MQ=255) ccccTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCGAACACGTTCACGGGTGAGGCCAATTTCACGAcc < 2:8104/90‑1 (MQ=255) ccTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAACCTGGCGAACACGTTCACGGGTGAGGCCAATTTCACGAcc < 1:566126/88‑1 (MQ=255) ccTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAACCTGGCGAACACGTTCACGGGTGAGGCCAATTTCACGAcc > 2:566126/1‑88 (MQ=255) ccTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAACCTGGCGAACACGTTCACGGGTGAGGCCAATTTCACGACCTa > 2:505805/1‑90 (MQ=255) ccTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAACCTGGCGAACACGTTCACGGGTGAGGCCAATTTCACGACCTa > 1:515520/1‑90 (MQ=255) tGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCGAACACGTTCACGGGTGAGGCCAATTTCACGACCTACa < 2:224400/90‑1 (MQ=255) gCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAACCTGGCGAACACGTTCACGGGTGAGGCCAATTTCACGACCTACAt < 2:234272/90‑1 (MQ=255) gCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAACCTGGCGAACACGTTCACGGGTGAGGCCAATTTCACGACCTACAt < 2:360558/90‑1 (MQ=255) ttGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAACCTGGCGAACACGTTCACGGGTGAGGCCAATTTCACGACCTACATCTTc > 2:90861/1‑90 (MQ=255) ttCGCGCAAACGGCGCAGGCCTTCAACCTGAACCTGGCGAACACGTTCACGGGTGAGGCCAATTTCACGACCTACATCTTCCAGTGTTGc < 1:195925/90‑1 (MQ=255) ttCGCGCAAACGGCGCAGGCCTTCAACCTGAACCTGGCGAACACGTTCACGGGTGAGGCCAATTTCACGACCTACATCTTCCAGTGTTGc < 1:426332/90‑1 (MQ=255) tCGCGCAAACGGCGCAGGCCTTCAACCTGAACCTGGCGAACACGTTCACGGGTGAGGCCAATTTCACGACCTACATCTTCCAGTGTTGcc > 2:471979/1‑90 (MQ=255) aaaCGGCGCAGGCCTTCAACCTGAACCTGGCGAACACGTTCACGGGTGAGGCCAATTTCACGACCTACATCTTCCAGTGTTGCCGc < 2:147312/86‑1 (MQ=255) aaaCGGCGCAGGCCTTCAACCTGAACCTGGCGAACACGTTCACGGGTGAGGCCAATTTCACGACCTACATCTTCCAGTGTTGCCGc > 1:147312/1‑86 (MQ=255) gAACCTGGCGAACACGTTCACGGGTGAGGCCAATTTCACGACCTACATCTTCCAGTGTTGCCGCTTCGTACCCCAGCAAACCGAATCGAc < 2:289705/90‑1 (MQ=255) gAACCTGGCGAACACGTTCACGGGTGAGGCCAATTTCACGACCTACATCTTCCAGTGTTGCCGCTTCGTACCCCAGCAAACCGAATCGAc < 2:212053/90‑1 (MQ=255) gAACCTGGCGAACACGTTCACGGGTGAGGCCAATTTCACGACCTACATCTTCCAGTGTTGCCGCTTCGTACCCCAGCAAACCGAATCGAc > 1:510471/1‑90 (MQ=255) | TACTTACTCGCGGAACAGCGCTTCGATATTCAGCCCCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCGAACACGTTCACGGGTGAGGCCAATTTCACGACCTACATCTTCCAGTGTTGCCGCTTCGTACCCCAGCAAACCGAATCGAC > NZ_CP009273/2859915‑2860084 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 31 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |