| 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 | 83.9% | 62.9 / 9.1 | 31 | I305V (ATT→GTT) | rpoS | RNA polymerase sigma factor RpoS |
| Reads supporting (aligned to +/- strand): ref base T (1/4); new base C (11/15); total (12/19) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 6.24e-01 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.51e-01 | |||||||||||
ACTTACTCGCGGAACAGCGCTTCGATATTCAGCCCCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCGAACACGTTCACGGGTGAGGCCAATTTCACGACCTACATCTTCCAGTGTTGCCGCTTCGTACCCCAGCAAACCGAATCGAC > NZ_CP009273/2859916‑2860084 | acttactCGCGGAACAGCGCTTCGATATTCAGCCCCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCGa < 2:595550/90‑1 (MQ=255)acttactCGCGGAACAGCGCTTCGATATTCAGCCCCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAACCTGGCGa < 2:36845/90‑1 (MQ=255)acttactCGCGGAACAGCGCTTCGATATTCAGCCCCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAACCTGGCGa < 2:670189/90‑1 (MQ=255) tactCGCGGAACAGCGCTTCGATATTCAGCCCCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAACCTGGCGAaca < 1:334007/90‑1 (MQ=255) ctCGCGGAACAGCGCTTCGATATTCAGCCCCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAACCTGGCGAACACg > 1:426847/1‑90 (MQ=255) ggAACAGCGCTTCGATATTCAGCCCCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAACCTGGCGAACACGTTCAc < 1:431328/90‑1 (MQ=255) aGCGCTTCGATATTCAGCCCCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAACCTGGCGAACACGTTCACGGGTg < 1:351129/90‑1 (MQ=255) tCGATATTCAGCCCCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAACCTGGCGAACACGTTCACGGGTGAGGCCa > 2:814279/1‑90 (MQ=255) tCAGCCCCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCGAACACGTTCACGGGTGAGGCCAATTTCAc < 1:678269/90‑1 (MQ=255) tCAGCCCCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAACCTGGCGAACACGTTCACGGGTGAGGCCAATTTCAc < 1:458469/90‑1 (MQ=255) cAGCCCCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAACCTGGCGAACACGTTCACGGGTGAGGCCAATTTCACg > 2:489886/1‑90 (MQ=255) gCCCCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAACCTGGCGAACACGTTCACGGGTGAGGCCAATTTCACGAc > 2:774462/1‑90 (MQ=255) gCCCCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAACCTGGCGAACACGTTCACGGGTGAGGCCAATTTCACGAc > 2:470530/1‑90 (MQ=255) ccccTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAACCTGGCGAACACGTTCACGGGTGAGGCCAATTTCACGAcc < 2:531873/90‑1 (MQ=255) tGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCGAACACGTTCACGGGTGAGGCCAATTTCAc < 1:812106/82‑1 (MQ=255) tGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCGAACACGTTCACGGGTGAGGCCAATTTCAc > 2:812106/1‑82 (MQ=255) tGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAACCTGGCGAACACGTTCACGGGTGAGGCCAATTTCACGACCTACa < 1:449951/90‑1 (MQ=255) tGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAACCTGGCGAACACGTTCACGGGTGAGGCCAATTTCACGACCTACa < 1:398953/90‑1 (MQ=255) ttGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAACCTGGCGAACACGTTCACGGGTGAGGCCAATTTCACGACCTACATCTTc > 1:244200/1‑90 (MQ=255) tCGCGCAAACGGCGCAGGCCTTCAACCTGAACCTGGCGAACACGTTCACGGGTGAGGCCAATTTCACGACCTACATCTTCCAGTGTTGcc > 1:702725/1‑90 (MQ=255) tCGCGCAAACGGCGCAGGCCTTCAACCTGAACCTGGCGAACACGTTCACGGGTGAGGCCAATTTCACGACCTACATCTTCCAGTGTTGcc < 2:470679/90‑1 (MQ=255) gcgcAAACGGCGCAGGCCTTCAACCTGAACCTGGCGAACACGTTCACGGGTGAGGCCAATTTCACGACCTACATCTTCCAGTGTTGCCGc < 2:89787/90‑1 (MQ=255) aaCGGCGCAGGCCTTCAACCTGAACCTGGCGAACACGTTCACGGGTGAGGCCAATTTCACGACCTACATCTTCCAGTGTTGCCGCTTCGt < 1:558268/90‑1 (MQ=255) gcgcAGGCCTTCAACCTGAACCTGGCGAACACGTTCACGGGTGAGGCCAATTTCACGACCTACATCTTCCAGTGTTGCCGCTTCGTAccc > 2:162262/1‑90 (MQ=255) gcgcAGGCCTTCAACCTGAACCTGGCGAACACGTTCACGGGTGAGGCCAATTTCACGACCTACATCTTCCAGTGTTGCCGCTTCGTAccc > 2:87098/1‑90 (MQ=255) cAGGCCTTCAACCTGAATCTGGCGAACACGTTCACGGGTGAGGCCAATTTCACGACCTACATCTTCCAGTGTTGCCGCTTCGTACCCCAg < 1:368052/90‑1 (MQ=255) tCAACCTGAACCTGGCGAACACGTTCACGGGTGAGGCCAATTTCACGACCTACATCTTCCAGTGTTGCCGCTTCGTACCCCAGCAAACCg > 1:828945/1‑90 (MQ=255) tCAACCTGAACCTGGCGAACACGTTCACGGGTGAGGCCAATTTCACGACCTACATCTTCCAGTGTTGCCGCTTCGTACCCCAGCAAACCg > 2:837317/1‑90 (MQ=255) gAACCTGGCGAACACGTTCACGGGTGAGGCCAATTTCACGACCTACATCTTCCAGTGTTGc < 2:684022/61‑1 (MQ=255) gAACCTGGCGAACACGTTCACGGGTGAGGCCAATTTCACGACCTACATCTTCCAGTGTTGCCGCTTCGTACCCCAGCAAACCGAATCGAc < 2:682271/90‑1 (MQ=255) gAACCTGGCGAACACGTTCACGGGTGAGGCCAATTTCACGACCTACATCTTCCAGTGTTGCCGCTTCGTACCCCAGCAAACCGAATCGAc < 1:246636/90‑1 (MQ=255) | ACTTACTCGCGGAACAGCGCTTCGATATTCAGCCCCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCGAACACGTTCACGGGTGAGGCCAATTTCACGACCTACATCTTCCAGTGTTGCCGCTTCGTACCCCAGCAAACCGAATCGAC > NZ_CP009273/2859916‑2860084 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 29 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |