Predicted mutation | ||||||
---|---|---|---|---|---|---|
evidence | seq id | position | mutation | annotation | gene | description |
RA | NZ_CP009273 | 2,122,528 | T→G | T202P (ACC→CCC) | wcaE ← | colanic acid biosynthesis glycosyltransferase WcaE |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NZ_CP009273 | 2,122,528 | 0 | T | G | 86.4% | 63.5 / 3.1 | 22 | T202P (ACC→CCC) | wcaE | colanic acid biosynthesis glycosyltransferase WcaE |
Reads supporting (aligned to +/- strand): ref base T (1/2); new base G (11/8); total (12/10) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 5.71e-01 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.44e-01 |
ATTCAGCCCAAAAGCCAGGCACATGTAATATTTGTCGTTGGACTTTTTTCGCGTCAGCACACAATTCCATATTATTGGTGGTAGATACCCCACCCATGGAAAATTCAGACACCAGGCCATTGAGTTTTTTAAATGCATAACCGGCTTTATACATTTTGGCTGCCAG > NZ_CP009273/2122447‑2122612 | aTTCAGCCCAAAAGCCAGGCACATGTAATATTTGTCGTTGGACTTTTTTCGCGTCAGCACACAATTCCATATTATTGGTGGTAGATAccc < 2:196455/90‑1 (MQ=255) cccAAAAGCCAGGCACATGTAATATTTGTCGTTGGACTTTTTTCGCGTCAGCACACAATTCCATATTATTGGTGGTAGATACcccaccca > 1:166501/1‑90 (MQ=255) aaGCCAGGCACATGTAATATTTGTCGTTGGACTTTTTTCGCGTCAGCACACAATTCCATATTATTGGTGGGAGATACCCCACCCATGGaa < 2:33594/90‑1 (MQ=255) aaGCCAGGCACATGTAATATTTGTCGTTGGACTTTTTTCGCGTCAGCACACAATTCCATATTATTGGTGGGAGATACCCCACCCATGGaa < 2:151423/90‑1 (MQ=255) gCACATGTAATATTTGTCGTTGGACTTTTTTCGCGTCAGCACACAATTCCATATTATTGGTGGGAGATACCCCACCCATGGAAAATTCAg < 1:781186/90‑1 (MQ=255) gCACATGTAATATTTGTCGTTGGACTTTTTTCGCGTCAGCACACAATTCCATATTATTGGTGGGAGATACCCCACCCATGGAAAATTCAg < 2:274036/90‑1 (MQ=255) aTGTAATATTTGTCGTTGGACTTTTTTCGCGTCAGCACACAATTCCATATTATTGGTGGGAGATACCCCACCCATGGAAAATTCAGacac > 1:321504/1‑90 (MQ=255) ttGTCGTTGGACTTTTTTCGCGTCAGCACACAATTCCATATTATTGTTGGGAGAAACCCCACCCATGGAAAATTCAGACACCAGGCCAtt > 1:281757/1‑90 (MQ=255) ttGTCGTTGGACTTTTTTCGCGTCAGCACACAATTCCATATTATTGGTGGGAGATACCCCACCCATGGAAAATGCAGACACCAGGCCAtt > 2:639562/1‑90 (MQ=255) gTCGTTGGACTTTTTTCGCGTCAGCACACAATTCCATATTATTGGTGGGAGATACCCCACCCATGGAAAATTCAGACACCAGGCCATTGa < 1:35491/90‑1 (MQ=255) ggACTTTTTTCGCGTCAGCACACAATTCCATATTATTGGTGGGAg < 2:734771/45‑1 (MQ=255) ttttCGCGTCAGCACACAATTCCATATTATTGGTGGTAGATACCCCACCCATGGAAAATTCAGACACCAGGCCATTGAGTTTTTTAAATg < 2:166501/90‑1 (MQ=255) gTCAGCACACAATTCCATATTATTGGTGGGAGATACCCCACCCATGGAAAATTCAGACACCAGGCCATTGAGTTTTTTAAATGCATAAcc > 1:685428/1‑90 (MQ=255) tCAGCACACAATTCCATATTATTGGTGGGAGATACCCCACCCATGGAAAATTCAGACACCAGGCCATTGAGTTTTTTAAATGCATAACCg < 1:647694/90‑1 (MQ=255) caATTCCATATTATTGGTGGGAGATACCCCACCCATGGAAAATTCAGACACCAGGCCATTGAGTTTTTTAAATGCATAACCGGCTTtata > 2:69420/1‑90 (MQ=255) ccATATTATTGGTGGGAGATACCCCACCCATGGAAAATTCAGACACCAGGCCATTGAGTTTTTTAAATGCATAACCGGCTTTATACAttt > 1:120933/1‑90 (MQ=255) ttattGGTGGGAGATACCCCACCCATGGAAAATTCAGACACCAGGCCATTGAGTTTTTTAAATGCATAACCGGCTTTATACATTTTGGCt > 2:882998/1‑90 (MQ=255) ttattGGTGGGAGATACCCCACCCATGGAAAATTCAGACACCAGGCCATTGAGTTTTTTAAATGCATAACCGGCTTTATACATTTTGGCt > 1:576710/1‑90 (MQ=255) attGGTGGGAGATACCCCACCCATGGAAAATTCAGACACCAGGCCATTGAGTTTTTTa < 2:448933/58‑1 (MQ=255) attGGTGGGAGATACCCCACCCATGGAAAATTCAGACACCAGGCCATTGAGTTTTTTa > 1:448933/1‑58 (MQ=255) attGGTGGGAGATACCCCACCCATGGAAAATTCAGACACCAGGCCATTGAGTTTTTTAAATGCATAACCGGCTTTATACATTTTGGCTGc > 1:304755/1‑90 (MQ=255) attGGTGGGAGATACCCCACCCATGGAAAATTCAGACACCAGGCCATTGAGTTTTTTAAATGCATAACCGGCTTTATACATTTTGGCTGc > 1:892388/1‑90 (MQ=255) ggtgggAGATACCCCACCCATGGAAAATTCAGACACCAGGCCATTGAGTTTTTTAAATGCATAACCGGCTTTATACATTTTGGCTGCCAg > 1:392863/1‑90 (MQ=255) | ATTCAGCCCAAAAGCCAGGCACATGTAATATTTGTCGTTGGACTTTTTTCGCGTCAGCACACAATTCCATATTATTGGTGGTAGATACCCCACCCATGGAAAATTCAGACACCAGGCCATTGAGTTTTTTAAATGCATAACCGGCTTTATACATTTTGGCTGCCAG > NZ_CP009273/2122447‑2122612 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |