Predicted mutation | ||||||
---|---|---|---|---|---|---|
evidence | seq id | position | mutation | annotation | gene | description |
RA | NZ_CP009273 | 3,736,824 | A→G | A244A (GCA→GCG) | yiaK → | 3‑dehydro‑L‑gulonate 2‑dehydrogenase |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NZ_CP009273 | 3,736,824 | 0 | A | G | 97.9% | 166.4 / ‑4.2 | 48 | A244A (GCA→GCG) | yiaK | 3‑dehydro‑L‑gulonate 2‑dehydrogenase |
Reads supporting (aligned to +/- strand): ref base A (0/1); new base G (27/20); total (27/21) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 4.38e-01 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.75e-01 |
CCGCATTTTGCCGATGGGCTACTGGAAAGGTTCTGGCATGTCGATTGTGCTGGATATGATCGCTACTCTCCTTTCCGACGGCGCATCCGTTGCCGAAGTCACCCAGGACAACAGCGACGAATACGGCATTTCACAAATTTTTATTGCCATTGAAGTGGACAAGCTTA > NZ_CP009273/3736740‑3736906 | ccgcATTTTGCCGATGGGCTACTGGAAAGGTTCTGGCATGTCGATTGTGCTGGATATGATCGCTACTCTCCTTTCCGACGGCGCGTCCGt < 2:969417/90‑1 (MQ=255) aTTTTGCCGATGGGCTACTGGAAAGGTTCTGGCATGTCGATTGTGCTGGATATGATCGCTACTCTCCTTTCCGACGGCGCGTCCGTTGcc < 2:448881/90‑1 (MQ=255) aTTTTGCCGATGGGCTACTGGAAAGGTTCTGGCATGTCGATTGTGCTGGATATGATCGCTACTCTCCTTTCCGACGGCGCGTCCGTTGcc < 1:986329/90‑1 (MQ=255) ttGCCGATGGGCTACTGGAAAGGTTCTGGCATGTCGATTGTGCTGGATATGATCGCTACTCTCCTTTCCGACGGCGCGTCCGTTGCCGaa > 1:624299/1‑90 (MQ=255) gATGGGCTACTGGAAAGGTTCTGGCATGTCGATTGTGCTGGATATGATCGCTACTCTCCTTTCCGACGGCGCGTCCGTTGCCGAAGTCAc < 1:595895/90‑1 (MQ=255) tACTGGAAAGGTTCTGGCATGTCGATTGTGCTGGATATGATCGCTACTCTCCTTTCCGACGGCGCGTCCGTTGCCGAAGTCACCCAGGac > 1:369309/1‑90 (MQ=255) tGGAACGGTTCTGGCATGTCGATTGTGCTGGATATGATCGCTACTCTCCTTTCCGACGGCGCGTCCGTTGCCGAAGTCACCCAGGacaac < 2:498258/90‑1 (MQ=255) aGGTTCTGGCATGTCGATTGTGCTGGATATGATCGCTACTCTCCTTTCCGACGGCGCATCCGTTGCCGAAGTCACCCAGGACAACAGcga < 1:1060360/90‑1 (MQ=255) ttCTGGCATGTCGATTGTGCTGGATATGATCGCTACTCTCCTTTCCGACGGCGCGTCCGTTGCCGAAGTCACCCAGGACAACAGcgacga > 1:347637/1‑90 (MQ=255) tGTCGATTGTGCTGGATATGATCGCTACTCTCCTTTCCGACGGCGCGTCCGTTGCCGAAGTCACCCAgg < 2:981727/69‑1 (MQ=255) tGTCGATTGTGCTGGATATGATCGCTACTCTCCTTTCCGACGGCGCGTCCGTTGCCGAAGTCACCCAgg > 1:981727/1‑69 (MQ=255) tGTCGATTGTGCTGGATATGATCGCTACTCTCCTTTCCGACGGCGCGTCCGTTGCCGAAGTCACCCAGGACAACAGCGACGAATACGGCa > 2:905260/1‑90 (MQ=255) tGTCGATTGTGCTGGATATGATCGCTACTCTCCTTTCCGACGGCGCGTCCGTTGCCGAAGTCACCCAGGACAACAGCGACGAATACGGCa > 2:222441/1‑90 (MQ=255) gTCGATTGTGCTGGATATGATCGCTACTCTCCTTTCCGACGGCGCGTCCGTTGCCGAAGTCACCCAGGACAACAGCGACGAATACGGCAt > 1:504174/1‑90 (MQ=255) gTCGATTGTGCTGGATATGATCGCTACTCTCCTTTCCGACGGCGCGTCCGTTGCCGAAGTCACCCAGGACAACAGCGACGAATACGGCAt > 2:959165/1‑90 (MQ=255) gTCGATTGTGCTGGATATGATCGCTACTCTCCTTTCCGACGGCGCGTCCGTTGCCGAAGTCACCCAGGACAACAGCGACGAATACGGCAt > 1:809620/1‑90 (MQ=255) gTCGATTGTGCTGGATATGATCGCTACTCTCCTTTCCGACGGCGCGTCCGTTGCCGAAGTCACCCAGGACAACAGCGACGAATACGGCAt > 1:92220/1‑90 (MQ=255) gTCGATTGTGCTGGATATGATCGCTACACTCCTTTCCGACGGCGCGTCCGTTGCCGAAGTCAccc < 1:210301/65‑1 (MQ=255) gTCGATTGTGCTGGATATGATCGCTACACTCCTTTCCGACGGCGCGTCCGTTGCCGAAGTCAccc > 2:210301/1‑65 (MQ=255) atGATCGCTACTCTCCTTTCCGACGGCGCGTCCGTTGCCGAAGTCACCCAGGACAACAGCGACGAATACGGCATTTCACAAATTTTTAtt < 2:809620/90‑1 (MQ=255) tGATCGCTACTCTCCTTTCCGACGGCGCGTCCGTTGCCGAAGTCAccc < 2:217222/48‑1 (MQ=255) tGATCGCTACTCTCCTTTCCGACGGCGCGTCCGTTGCCGAAGTCAccc > 1:217222/1‑48 (MQ=255) tGATCGCTACTCTCCTTTCCGACGGCGCGTCCGTTGCCGAAGTCACCCAGGACAACAGCGACGAATACGTCATTTCACAAATTTTTATTg < 2:369309/90‑1 (MQ=255) tGATCGCTACTCTCCTTTCCGACGGCGCGTCCGTTGCCGAAGTCACCCAGGACAACAGCGACGAATACGGCATTTCACAAATTTTTATTg > 1:364686/1‑90 (MQ=255) tGATCGCTACTCTCCTTTCCGACGGCGCGTCCGTTGCCGAAGTCACCCAGGACAACAGCGACGAATACGGCATTTCACAAATTTTTATTg > 1:948906/1‑90 (MQ=255) tGATCGCTACTCTCCTTTCCGACGGCGCGTCCGTTGCCGAAGTCACCCAGGACAACAGCGACGAATACGGCATTTCACAAATTTTTATTg > 2:457554/1‑90 (MQ=255) tGATCGCTACTCTCCTTTCCGACGGCGCGTCCGTTGCCGAAGTCACCCAGGACAACAGCGACGAATACGGCATTTCACAAATTTTTATTg > 2:813972/1‑90 (MQ=255) gCTACTCTCCTTTCCGACGGCGCGTCCGTTGCCGAAGt < 1:973244/38‑1 (MQ=38) gCTACTCTCCTTTCCGACGGCGCGTCCGTTGCCGAAGt > 2:973244/1‑38 (MQ=38) gCTACTCTCCTTTCCGACGGCGCGTCCGTTGCCGAAGt < 1:868396/38‑1 (MQ=38) gCTACTCTCCTTTCCGACGGCGCGTCCGTTGCCGAAGt > 2:868396/1‑38 (MQ=38) gCTACTCTCCTTTCCGACGGCGCGTCCGTTGCCGAAGTCACCCAGGACAACAGCGACGAATACGGCATTTCACAAATTTTTATTGCCAtt > 1:601772/1‑90 (MQ=255) gCTACTCTCCTTTCCGACGGCGCGTCCGTTGCCGAAGTCACCCAGGACAACAGCGACGAATACGGCATTTCACAAATTTTTATTGCCAtt > 1:251778/1‑90 (MQ=255) aCTCTCCTTTCCGACGGCGCGTCCGTTGCCGAAGTCACCCAGGACAACAGCGACGAATACGGCATTTCACAAATTTTTATTGCCATTGaa < 1:959165/90‑1 (MQ=255) ctcCTTTCCGACGGCGCGTCCGTTGCCGAAGTCACCCAGGACAACAGCGACGAATACGGCATTTCACAAATTTTTATTGCCATTGAAGTg < 1:376736/90‑1 (MQ=255) cTTTCCGACGGCGCGTCCGTTGCCGAAGTCACCCAGGACAACAGCGACGAATACGGCATTTCACAAATTTTTATTGCCATTGAAGTGGAc < 2:9540/90‑1 (MQ=255) ttCCGACGGCGCGTCCGTTGCCGAAGTCACCCAGGACAACAGCGACGAATAc > 1:698524/1‑52 (MQ=255) ttCCGACGGCGCGTCCGTTGCCGAAGTCACCCAGGACAACAGCGACGAATAc < 2:698524/52‑1 (MQ=255) ttCCGACGGCGCGTCCGTTGCCGAAGTCACCCAGGACAACAGCGACGAATACGGCATTTCACAAATTTTTATTGCCATTGAAGTGGACaa > 1:396116/1‑90 (MQ=255) ttCCGACGGCGCGTCCGTTGCCGAAGTCACCCAGGACAACAGCGACGAATACGGCATTTCACAAATTTTTATTGCCATTGAAGTGGACaa > 1:621311/1‑90 (MQ=255) ttCCGACGGCGCGTCCGTTGCCGAAGTCACCCAGGACAACAGCGACGAATACGGCATTTCACAAATTTTTATTGCCATTGAAGTGGACaa < 1:610821/90‑1 (MQ=255) ttCCGACGGCGCGTCCGTTGCCGAAGTCACCCAGGACAACAGCGACGAATACGGCATTTCACAAATTTTTATTGCCATTGAAGTGGACaa > 1:973746/1‑90 (MQ=255) ttCCGACGGCGCGTCCGTTGCCGAAGTCACCCAGGACAACAGCGACGAATACGGCATTTCACAAATTTTTATTGCCATTGAAGTGGACaa < 1:1005324/90‑1 (MQ=255) ttCCGAAGGCGCGTCCGTTGCCGAAGTCACCCAGGACAACAGCGACGAATACGGCATTTCACAAATTTTTATTGCCATTGAAGTGGACaa > 1:575201/1‑90 (MQ=255) aCGGCGCGTCCGTTGCCGAAGTCACCCAGGACAACAGCGACGAATACGGCGTTTCACAAATTTTTATTGCCATTGAAGTGGACAAGCTTa < 2:93992/90‑1 (MQ=255) aCGGCGCGTACGTTGCCGAAGTCACCCAGGACAACAGCGACGAATACGGCATTTCACAAATTTTTATTGCCATTGAAGTGGACAAGCTTa > 2:36804/1‑90 (MQ=255) ggCGCGTCCGTTGCCGAAGTCACCCAGGACAACAGCGACGAATACGGCATTTCACaa > 1:976344/1‑57 (MQ=255) ggCGCGTCCGTTGCCGAAGTCACCCAGGACAACAGCGACGAATACGGCATTTCACaa < 2:976344/57‑1 (MQ=255) | CCGCATTTTGCCGATGGGCTACTGGAAAGGTTCTGGCATGTCGATTGTGCTGGATATGATCGCTACTCTCCTTTCCGACGGCGCATCCGTTGCCGAAGTCACCCAGGACAACAGCGACGAATACGGCATTTCACAAATTTTTATTGCCATTGAAGTGGACAAGCTTA > NZ_CP009273/3736740‑3736906 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 26 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |