Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | minE | 1,523,433 | G→T | 24.9% | A78S (GCC→TCC) | flk → | predicted flagella assembly protein |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | minE | 1,523,433 | 0 | G | T | 24.9% | 91.5 / 12.6 | 36 | A78S (GCC→TCC) | flk | predicted flagella assembly protein |
Reads supporting (aligned to +/- strand): ref base G (17/10); new base T (7/2); total (24/12) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 6.85e-01 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.29e-01 |
CTGTGGGCGGGTATCAGGCACGATATTGGTCTGTCGGGAGATTCACCGCTGCTTTCGCGTCACTTCCCTGCCGCTGAGCATAATCTGGCGCAACGTCTGCTGGCCGCGCAAAAAAGCCATTCTGCCCGCCAGCTTT > minE/1523364‑1523499 | cTGTGGGCGGGTATCAGGCACGATATTGGTCTGTCGGGAGATTCACCGCTGCTTTCGCGTCACTTCCCTGc > 1:101162/1‑71 (MQ=255) tATCAGGCACGATATTGGTCTGTCGGGAGATTCACCGCTGCTTTCGCGTCACTTCCCTGCCGCTGAGCAt > 1:2197560/1‑70 (MQ=255) tCAGGCACGATATTGGTCTGTCGGGAGATTCACCGCTGCTTTCGCGTCACTTCCCTGCCGCTGAGCATa < 1:878762/69‑1 (MQ=255) aGGCACGATATTGGTCTGTCGGGAGATTCACCGCTGCTTTCGCGTCACTTCCCTGCCGCTGAGCa < 1:1450519/65‑1 (MQ=255) aCGATATTGGTCTGTCGGGAGATTCACCGCTGCTTTCGCGTCACTTCCCTGCCGCTGAGCATAATCTGgcg > 1:792730/1‑71 (MQ=255) aCGATATTGGTCTGTCGGGAGATTCACCGCTGCTTTCGCGTCACTTCCCTGCCGCTGAGCATAATATGgcg > 1:1894131/1‑71 (MQ=255) gATATTGGTCTGTCGGGAGATTCACCGCTGCTTTCGCGTCACTTCCCTGCCGCTGAGCATAATCTGGCGCa < 1:1077721/71‑1 (MQ=255) tatTGGTCTGTCGGGAGATTCACCGCTGCTTTCGCGTCACTTCCCTGCCGCTGAGCATAATCTGGCGCAAc < 1:781086/71‑1 (MQ=255) tGGTCTGTCGGGAGATTCACCGGTGCTTTCGCGTCACTTCCCTGCCGCTGAGCATAATCTGGCGCAACGTc > 1:1331746/1‑71 (MQ=255) ggTCTGTCGGGAGATTCACCGCTGCTTTCGCGTCACTTCCCTGCCGCTGAGCATAATCTGGCGCAACGTct > 1:470285/1‑71 (MQ=255) gTCTGTCGGGAGATTCACCGCTGCTTTCGCGTCACTTCCCTTCCGCTGAGCATAATCTGGCGCAACGTctg < 1:855600/71‑1 (MQ=255) tGTCGGGAGATTCACCGCTGCTTTCGCGTCACTTCCCTTCCGCTGAGCATAATCTGGCGCAACGTctgctg > 1:1281654/1‑71 (MQ=255) tGTCGGGAGATTCACCGCTGCTTTCGCGTCACTTCCCTTCCGCTGAGCATAATCTGGCGCAACGTctgctg > 1:1141098/1‑71 (MQ=255) tGTCGGGAGATTCACCGCTGCTTTCGCGTCACTTCCCTGCCGCTGAGCATAAt > 1:1316640/1‑53 (MQ=255) gggAGATTCACCGCTGCTTTCGCGTCACTTCCCTGCCGCTGAGCATAATCTGGCGCAACGTCTGCTGGCcg > 1:987513/1‑71 (MQ=255) gaTTCACCGCTGCTTTCGCGTCACTTCCCTGCCGCTGAGCATa < 1:1525236/43‑1 (MQ=255) tCACCGCTGCTTTCGCGTCACTTCCCTGCCGCTGAGCATAATCTGGCGCAACGTCTGCTGGCCGCGCaaa > 1:1684358/1‑70 (MQ=255) aCCGCTGCTTTCGCGTCACTTCCCTGCCGCTGAGCATAATCTGGCGCAACGTCTGCTGGCCGCGCaaaaa < 1:1485817/70‑1 (MQ=255) cGCTGCTTTCGCGTCACTTCCCTGCCGCTGAGCATAATCTGGCGCAACGTCTGCTGGCCGCGCaaaaaa > 1:1401554/1‑69 (MQ=255) cGCTGCTTTCGCGTCACTTCCCTGCCGCTGAGCATAATCTGGCGCAACGTCTGCTGGCCGCGCaaaaaa > 1:279447/1‑69 (MQ=255) cGCTGCTTTCGCGTCACTTCCCTGCCGCTGAGCATAATCTGGCGCAACGTCTGCTGGCCGCGCAAAAAAGc > 1:1156866/1‑71 (MQ=255) gctgctTTCGCGTCACTTCCCTTCCGCTGAGCATAATCTGGCGCAACGTCTGCTGGCCGCGCAAAAAAg > 1:997944/1‑69 (MQ=255) gctgctTTCGCGTCACTTCCCTTCCGCTGAGCATAATCTGGCGCAACGTCTGCTGGCCGCGCAAAAAAGc > 1:1991758/1‑70 (MQ=255) ctTTCGCGTCACTTCCCTTCCGCTGAGCATAATCTGGCGAAACGTCTGCTGGCCGCGCAAAAAAGCCAt < 1:2084260/69‑1 (MQ=255) ttCGCGTCACTTCCCTGCCGCTGAGCATAATCTGGCGCAACGTCTGCTGGCCGCGCAAAAAAGCCATTCTg > 1:1202674/1‑71 (MQ=255) gTCACTTCCCTTCCGCTGAGCATAATCTGGCGCAACGTCTGCTGGCCGCGCAAAAAAGCCATTCTGCCCGc > 1:1474428/1‑71 (MQ=255) gTCACTTCCCTTCCGCTGAGCATAATCTGGCGCAACGTCTGCTGGCCGCGCAAAAAAGCCATTCTGCCCGc > 1:969321/1‑71 (MQ=255) tCACTTCCCTTCCGCTGAGCATAATCTGGCGCAACGTCTGCTGGCCGCGCAAAAAAGCCATTCTGCCCGc > 1:1472562/1‑70 (MQ=255) tCACTTCCCTGCCGCTGAGCATAATCTGGCGCAACGTCTGCTGGCCGCGCAAAAAAGCCATTCTGCCCGcc < 1:552271/71‑1 (MQ=255) aCTTCCCTGCCGCTGAGCATAATCTGGCGCAACGTCTGCTGGCCGCGCAAAAAAGCCATTCTGCCCGCCAg > 1:553459/1‑71 (MQ=255) aCTTCCCTGCCGCTGAGCATAATCTGGCGCAACGTCTGCTGGCCGCGCAAAAAAGCCATTCTGCCCGCCAg > 1:1971287/1‑71 (MQ=255) aCTTCCCTGCCGCTGAGCATAATCTGGCGCAACGTCTGCTGGCCGCGCAAAAAAGCCATTCTGCCCGCCAg > 1:1265149/1‑71 (MQ=255) ttCCCTGCCGCTGAGCATAATCTGGCGCAACGTCTGCTGGCCGCGCaaa < 1:850902/49‑1 (MQ=255) cccTGCCGCTGAGCATAATCTGGCGCAACGTCTGCTGGCCGCGCAAAAAAGCCATTCTGCCCGCCAGCttt > 1:7928/1‑71 (MQ=255) cTGCCGCTGAGCATAATCTGGCGCAACGTCTGCTGGCCGCGCAAAAAAGc < 1:538945/50‑1 (MQ=255) cTGCCGCTGAGCATAATCTGGCGCAACGTCTGCTGGCCGCGCAAAAAAGc < 1:1218433/50‑1 (MQ=255) | CTGTGGGCGGGTATCAGGCACGATATTGGTCTGTCGGGAGATTCACCGCTGCTTTCGCGTCACTTCCCTGCCGCTGAGCATAATCTGGCGCAACGTCTGCTGGCCGCGCAAAAAAGCCATTCTGCCCGCCAGCTTT > minE/1523364‑1523499 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |