Predicted mutation | ||||||
---|---|---|---|---|---|---|
evidence | seq id | position | mutation | annotation | gene | description |
RA | NZ_CP009273 | 3,103,871 | (A)7→6 | intergenic (+71/+78) | BW25113_RS15410 → / ← yghD | tRNA‑Phe/GspM family type II secretion system protein YghD |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NZ_CP009273 | 3,103,865 | 0 | A | . | 94.1% | 58.1 / ‑1.0 | 17 | intergenic (+65/+84) | BW25113_RS15410/yghD | tRNA‑Phe/GspM family type II secretion system protein YghD |
Reads supporting (aligned to +/- strand): ref base A (1/0); new base . (2/14); total (3/14) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 1.76e-01 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.47e-01 |
AGATTCCGAGTCCGGGCACCACTAATTCTTAAGAACCCGCCCACAAGGCGGGTTTTTGCTTTTGGATCTGACAATAACCTTCACGAAAAAAATTAGCTTATAAAGTCTGGGGGAATTACTCTCGCCACGTTAACGAGAGTAATTTTATTGATATTAATCTCCTGATACTT > NZ_CP009273/3103780‑3103949 | aGATTCCGAGTCCGGGCACCACTAATTCTTAAGAACCCGCCCACAAGGCGGGTTTTTGCTTTTGGATCTGACAATAACCTTCACGaaaaa > 1:344762/1‑90 (MQ=255) aaGAACCCGCCCACAAGGCGGGTTTTTGCTTTTGGATCTGACAATAACCTTCACG‑AAAAAATTAGCTTATAAAGTCTGGGGGAATTActc < 2:199101/90‑1 (MQ=255) cacaAGGCGGGTTTTTGCTTTTGGATCTGACAATAACCTTCACG‑AAAAAATTAGCTTATAAAGTCTGGGGGAATTACTCTCGCCACGTTa < 1:197653/90‑1 (MQ=255) aaGGCGGGTTTTTGCTTTTGGATCTGACAATAACCTTCACG‑AAAAAATTAGCTTATAAAGTCTGGGGGAATTACTCTCGCCACGTTAACg < 2:291280/90‑1 (MQ=255) gggTTTTTGCTTTTGGATCTGACAATAACCTTCACG‑AAAAAATTAGCTTATAAAGTCTGGGGGAATTACTCTCGCCACGTTAACGAGAGt < 2:225172/90‑1 (MQ=255) cTTTTGGATCTGACAATAACCTTCACG‑AAAAAATTAGCTTATAAAGTCTGGGGGAACTACTCTCGCCACGTTAACGAGAGTAATTTTAtt < 1:515312/90‑1 (MQ=255) ttGGATCTGACAATAACCTTCACG‑AAAAAATTAGCTTATAAAGTCTGGGGGAATTACTCTCGCCACGTTAACGAGAGTAATTTTATTGat < 2:123427/90‑1 (MQ=255) ggATCTGACAATAACCTTCACG‑AAAAAATTAGCTTATAAAGTCTGGGGGAATTACTCTCGCCACGTTAACGAGAGTAATTTTATTGatat < 1:327050/90‑1 (MQ=255) cTGACAATAACCTTCACG‑AAAAAATTAGCTTATAAAGTCTGGGGGAATTACTCTCGCCACGTTAACGAGAGTAATTTTATTGATATTAAt < 2:681555/90‑1 (MQ=255) tGACAATAACCTTCACG‑AAAAAATTAGCTTATAAAGTCTGGGGGAATTACTCTCGCCACGTTAACGAGAGTAATTTTATTGATATTAAtc < 1:460148/90‑1 (MQ=255) aTAACCTTCACG‑AAAAAATTAGCTTATAAAGTCTGGGGGAATTACTCTCGCCACGTTAACGAGAGTAATTTTATTGATATTAATCTCCTg < 1:113704/90‑1 (MQ=255) tAACCTTCACG‑AAAAAATTAGCTTATAAAGTCTGGGGGAATTACTCTCGCCACGTTAACGAGAGTAATTTTATTGATATTAATCTCCTGa < 2:344762/90‑1 (MQ=255) aaCCTTCACG‑AAAAAATTAGCTTATAAAGTCTGGGGGAATTACTCTCGCCACGTTAACGAGAGTAATTTTATTGATATTAATCTCCTGAt < 1:583250/90‑1 (MQ=255) aaCCTTCACG‑AAAAAATTAGCTTATAAAGTCTGGGGGAATTACTCTCGCCACGTTAACGAGAGTAATTTTATTGATATTAATCTCCTGAt < 2:340367/90‑1 (MQ=255) ttCACG‑AAAAAATTAGCTTATAAAGTCTGGGGGAATTACTCTCGCCACGTTAACGAGAGTAATTTTATTGATATTAATCTCCTGATACtt > 1:823716/1‑90 (MQ=255) tCACGAAAAAAATTAGCTTATAAAGTCTGGGGGAATTACTCTCGCCACGTTAACGAGAGTAATTTTATTGATATTAATCTCCTGATACtt > 2:364588/1‑90 (MQ=255) tCACG‑AAAAAATTAGCTTATAAAGTCTGGGGGAATTACTCTCGCCa < 1:310870/46‑1 (MQ=255) tCACG‑AAAAAATTAGCTTATAAAGTCTGGGGGAATTACTCTCGCCa > 2:310870/1‑46 (MQ=255) | AGATTCCGAGTCCGGGCACCACTAATTCTTAAGAACCCGCCCACAAGGCGGGTTTTTGCTTTTGGATCTGACAATAACCTTCACGAAAAAAATTAGCTTATAAAGTCTGGGGGAATTACTCTCGCCACGTTAACGAGAGTAATTTTATTGATATTAATCTCCTGATACTT > NZ_CP009273/3103780‑3103949 |
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: ‑ |