Predicted mutation | ||||||
---|---|---|---|---|---|---|
evidence | seq id | position | mutation | annotation | gene | description |
RA | NZ_CP009273 | 1,423,602 | T→C | R99R (CGT→CGC) | stfR → | prophage tail fiber N‑terminal domain‑containing protein |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NZ_CP009273 | 1,423,602 | 0 | T | C | 85.7% | 40.6 / 2.8 | 21 | R99R (CGT→CGC) | stfR | prophage tail fiber N‑terminal domain‑containing protein |
Reads supporting (aligned to +/- strand): ref base T (1/2); new base C (10/8); total (11/10) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 5.86e-01 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.38e-01 |
GACCATTACCGTGTATGAAGATTCTCAACCCGGTACGCTGAATGATTTTCTCGGTGCCATGACGGAGGATGATGCCCGTCCGGAGGCACTGCGCCGTTTTGAACTGATGGTGGAAGAGGTGGCGCGTAACGCGTCC > NZ_CP009273/1423524‑1423659 | gACCATTACCGTGTATGAAGATTCTCAACCCGGTACGCTGAATGATTTTCTCGGTGCCATGACGGAGGATGATGCCCGTCCGGAGGCACt > 1:724396/1‑90 (MQ=255) aTTACCGTGTATGAAGATTCTCAACCCGGTACGCTGAATGATTTTCTCGGTGCCATGACGGAGGATGATGCCCGTCCGGAGGCACTgcgc < 2:521310/90‑1 (MQ=255) aTTACCGTGTATGAAGATTCTCAACCCGGTACGCTGAATGATTTTCTCGGTGCCATGACGGAGGATGATGCCCGCCCGGAGGCACTgcgc < 2:796538/90‑1 (MQ=255) aTTACCGTGTATGAAGATTCTCAACCCGGTACGCTGAATGATTTTCTCGGTGCCATGACGGAGGATGATGCCCGCCCGGAGGCACTgcgc > 1:796538/1‑90 (MQ=255) aTTACCGTGTATGAAGATTCTCAACCCGGTACGCTGAATGATTTTCTCGGTGCCATGACGGAGGATGATGCCCGCCCGGAGGCACTgcgc < 1:828459/90‑1 (MQ=255) aTTACCGTGTATGAAGATTCTCAACCCGGTACGCTGAATGATTTTCTCGGTGCCATGACGGAGGATGATGCCCGCCCGGAGGCACTgcgc < 2:381773/90‑1 (MQ=255) gtgtATGAGGATTCTCAACCCGGTACGCTGAATGATTTTCTCGGTGCCATGACGGAGGATGATGCCCGCCCGGAGGCACTGCGCCGtttt > 1:117495/1‑90 (MQ=255) gtgtATGAAGATTCTCAACCCGGTACGCTGAATGATTTTCTCGGTGCCATGACGGAGGATGATGCCCGCCCGGAGGCACTGCGCCGtttt > 2:328608/1‑90 (MQ=255) gAAGATTCTCAACCCGGTACGCTGAATGATTTTCTCGGTGCCATGACGGAGGATGATGCCCGCCCGGAGGCACTGCGCCGTTTTGAACTg < 1:382738/90‑1 (MQ=255) aGATTCTCAACCCGGTACGCTGAATGATTTTCTCGGTGCCATGACGGAGGATGATGCCCGCCCGGAGGCACTGCGCCGTTTTGAACTGAt < 1:316618/90‑1 (MQ=255) aGATTCTCAACCCGGTACGCTGAATGATTTTCTCGGTGCCATGACGGAGGATGATGCCCGCCCGGAGGCACTGCGCCGTTTTGAACTGAt < 1:581368/90‑1 (MQ=255) tctcAACCCGGTACGCTGAATGATTTTCTCGGTGCCATGCCGGAGGTTAAAACACGCCCGGAGGCACTGCGCCGTTTGGAACTGAtggtg > 2:826288/1‑90 (MQ=255) tctcAACCCGGTACGCTGAATGATTTTCTCGGTGCCATGACGGAGGATGATGCCCGCCCGGAGGCACTGCGCCGTTTTGAAc < 1:614386/82‑1 (MQ=255) tctcAACCCGGTACGCTGAATGATTTTCTCGGTGCCATGACGGAGGATGATGCCCGCCCGGAGGCACTGCGCCGTTTTGAAc > 2:614386/1‑82 (MQ=255) tctcAACCCGGTACGCTGAATGATTTTCTCGGTGCCATGACGGAGGATGATGCCCGCCCGGAGGCACTGCGCCGTTTTGAACTGAtggtg > 1:267063/1‑90 (MQ=255) tctcAACCCGGTACGCTGAATGATTTTCTCGGTGCCATGACGGAGGATGATGCCCGCCCGGAGGCACTGCGCCGTTTTGAACTGAtggtg > 1:832529/1‑90 (MQ=255) tctcAACCCGGTACGCTGAATGATTTTCTCGGTGCCATGACGGAGGATGATGCCCGCCCGGAGGCACTGCGCCGTTTTGAACTGAtggtg > 1:535319/1‑90 (MQ=255) ggTACGCTGAATGATTTTCTCGGTGCCATGACGGAGGATGATGCCCGCCCGGAGGCACTGCGCCGTTTTGAACTGATGGTGGAAGAGGTg > 1:724815/1‑90 (MQ=255) tttCTCGGTGCCATGACGGAGGATGATGCCCGTCCGGAGGCACTGCGCCGTTTTGAACTGATGGTGGAAGAGGTGGCGCGTAACGCGTcc < 1:476215/90‑1 (MQ=255) atgCCCGCCCGGAGGCACTGCGCCGTTTTGAACTGATGGTGGAAGAGGTGGCGCGTAACGCGTcc > 1:7563/1‑65 (MQ=255) atgCCCGCCCGGAGGCACTGCGCCGTTTTGAACTGATGGTGGAAGAGGTGGCGCGTAACGCGTcc < 2:7563/65‑1 (MQ=255) | GACCATTACCGTGTATGAAGATTCTCAACCCGGTACGCTGAATGATTTTCTCGGTGCCATGACGGAGGATGATGCCCGTCCGGAGGCACTGCGCCGTTTTGAACTGATGGTGGAAGAGGTGGCGCGTAACGCGTCC > NZ_CP009273/1423524‑1423659 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |