Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | CP009273 | 1,123,224 | A→T | 32.2% | intergenic (+39/‑71) | yceM → / → murJ | putative oxidoreductase/putative peptidoglycan lipid II flippase |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | CP009273 | 1,123,224 | 0 | A | T | 32.2% | 30.5 / 20.2 | 31 | intergenic (+39/‑71) | yceM/murJ | putative oxidoreductase/putative peptidoglycan lipid II flippase |
Reads supporting (aligned to +/- strand): ref base A (12/9); new base T (4/6); total (16/15) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 4.58e-01 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.95e-01 |
TTGAATGTGTGCAAAACCAGACAGTTCCGCAAACCGCCGGCGAACAGGCCGTGCTGGCGCAACGTATCGTTGACAAGATCTGGCGCGATGCGATGAGTGAATAAACCCTGTAACATCTGGCGGTAGCAATTCACCATAATCCAGGTAGACTATTCGCCTCTTTCAGCGCCTGCCTTGCAGGCGTTTTGCCCGTGGGTCTGGATTAGAACACCGATGAATTTATTAAAATCGCTGGCCGCCGTCAGCTCGATGACCATGTTTTCGCGTGTGCTTGGCTT > CP009273/1123082‑1123359 | ttGAATGTGTGCAAAACCAGACAGTTCCGCAAACCGCCGGCGAACAGGCCGTGCTGGCGCAACGTATCGTTGACAAGATCTGGCGCGATGCGATGAGTGAATAAACCCTGTAACATCTGGCGGTAGCAATTCACCATAATCCTGGTAGa < 2:295127/149‑1 (MQ=255) ttGAATGTGTGCAAAACCAGACAGTTCCGCAAACCGCCGGCGAACAGGCCGTGCTGGCGCAACGTATCGTTGACAAGATCTGGCGCGATGCGATGAGTGAATAAACCCTGTAACATCTGGCGGTAGCAATTCACCATAATCCAGGTAGa < 1:217204/149‑1 (MQ=255) tgCAAAACCAGACAGTTCCGCAAACCGCCGGCGAACAGGCCGTGCTGGCGCAACGTATCGTTGCCAAGATCTGGCGCGATGCGATCAGTGCCCCATCGCGGTATCATCGGGCGGTAGCCATTCACCATAATCCTGGTAGACTATTCGcc > 1:76325/1‑149 (MQ=255) gCAAAACCAGACAGTTCCGCAAACCGCCGGCGAACAGGCCGTGCTGGCGCAACGTATCGTTGACAAGACCTGGCGCGATGCGATGAGTGAATAAACCATGTAACATCTGGCGGGAGCAATTCACCATAATCCTGGTAGACTAGTCGcca > 2:156425/1‑148 (MQ=255) gTTCCGCAAACCGCCGGCGAACAGGCCGTGCTGGCGCAACGTATCGTTGACAAGATCTGGCGCGATGCGATGAGTGAATAAACCCTGTAACATCTGGCGGTAGCAATTCACCATAATCCTGGTAGACTATTCGCCCCTTTCAGCgcctg > 2:101216/1‑149 (MQ=255) gCAAGCCGCCGGCGAACAGGCCGTGCGGGCGCAACGTATCGTTGACAAGATCTGGCGCGATGCGATGAGTGAATAAACCCTGTACCATCTGGCGGTAGCAATTCACCATAATCCAGGTAGACTATTCGCCTCTTTCAGCGCCTGCCTTg < 2:140621/149‑1 (MQ=255) cAAACCGCCGGCGAACAGGCCGTGCTGGCGCAACGTATCGTTGACAAGATCTGGCGCGATGCGATGAGTGAATAAACCCTGTAACATCTGGCGGTAGCAATTCACCATAATCCAGGTAGACTATTCGCCTCTTTCAGCGCCTGCCTTGc > 1:249033/1‑149 (MQ=255) cAAACCGCCGGCGAACAGGCCGTGCTGGCGCAACGTATAGCTAAGCAGATCTGGCGCGATGCGCTGAGTGAATCAACCCTGTAACATCTGGCGGAAGCAATTCACCATAATCCAGGTAGCCTATTCGGCTCTTTCAGCGCCTGCCTTGc > 1:93621/1‑149 (MQ=255) gAACAGGCCGTGCTGGCGCAACGTATCGTTGACAAGATCTGGTGCGATGCGATGAGTGAATAAACCCTGTAACATCTGGCGGTAGCAATTCACCATAATCCAGGTAGACTATTCGCCTCTTTCAGCGCCTGCCTTGCAGGCGTTTTGcc < 1:290182/149‑1 (MQ=255) cAGGCCGTGCTGGCGGAACGTATCGTTGACAAGATCTGGCGCGATGCGATGAGTGAATAAACCCGGTAACATCTGGCGGTAGCAATTCACCATAATCAAGGTAGACTATTCGCCTCGTTCAGCGCCTGCCTTGCAGGCGTTTTGCCCGt > 1:225701/1‑149 (MQ=255) gCCGTGCTGGCGCAACGTATCGTTGACAAGATCTGGCGCGATGCGATGAGTGAATAAACCCTGTAACATCTGGCGGTAGCAATTCACCATAATCCTGGTAGACTATTCGCCTCTTTCAGCGCCTGCCTTGCAGGCGTTTTGCCCGTggg > 1:351057/1‑149 (MQ=255) cgcAACGTACCGTTGACAAGATCTGGCGCGATGCGATGAGTGAATAAACCCTGTAACATCTGGCGGTAGCAATTCACCATAATCCTGGTAGACTATTCGCCTCTTTCAGCGCCTGCCTTGCAGGCGTTTTGCCCGTGGGTCTGGATTAg < 2:73368/149‑1 (MQ=255) aaCGTATCGTTGACAAGATCTGGCGCGATGCGATGAGTGAATAAACCCTGTAACATCTGGCGGTAGCAATTCACCATAATCCAGGTAGACTATTCGCCTCTTTCAGCGCCTGCCTTGCAGGCGTTTTGCCCGTGGGTCTGGATTAGAac > 2:193005/1‑149 (MQ=255) atccgtgactacagcTGGCGCAATGTGAGGAGTGAATAAACCCTGCAACATCTGGCGGTAGCATTTCACCATAAGCCTGGTAGACTGTGCGCCTCTTTCAGCGCCTGCCTTGCTGGCGGTTTGCCCGTGGGTCTGGATTAGAACACCGa < 1:71946/135‑1 (MQ=255) cAAGATCTGGCGCGATGCGATGAGTGAATAAACCCTGTAACATCTGGCGGTAGCACTTCACCATAATCCAGGTCGACTATTCGCCTCTTTCAGCGCCTGCCTGGCAGGCGTTTTGCCCGTGGGTCGGGATTAGAACACCGATGCGTtta > 2:89973/1‑149 (MQ=255) cAAGATCTGGCGCGATGCGATGAGTGAATAAACCCTGTAACATCTGGCGGTAGCAATTCACCATAATCCAGGTAGACTATTCGCCTCTTTCAGCGCCTGCCTTGCAGGCGTTTTGCCCGTGGGTCTGGATTAGAACACCGATGAATtta > 1:36978/1‑149 (MQ=255) cAAGATCTGGCGCGATGCGATGAGGGAATAAACCCTGTAACATCTGGCGGTAGCAATTCACCATAATCCAGGTAGACTATTCGCCTCTTTCAGCGCCTGCCTTGCAGGCGTTTTGCCCGTGGGTCTGGATTAGAACACCGATGAATtta < 1:152383/149‑1 (MQ=255) cAAGAGCTGGAGCGATGCGATGAGTGAATAAACCCTGTAACATCTGGCGGTAGCAATTCACCATAATCCAGGTAGACTATTCGCCTCTTTCAGCGCCTGCCTTGCAGGCGTTTTGCCCGTGGGTCTGGATTAGAACCCCGATTAATtta < 1:151151/149‑1 (MQ=255) gtCTGGCGCGATGCGATGAGTGTATCAACCCCGTAACATCTGGCGGTAGCAATTCACCATAATCCTGGTAGACTATTCGCCTCTTTCAGCGCCTGCCTTGCAGGCGTTTTGCCCGTGGGTCTGGATTAGAACACCGATGAATTTATTaa < 1:234087/148‑1 (MQ=255) aTCTGGCGCGATGCGATGAGTGAATAAACCCTGTAACATCTGGCGGTAGCAATTCACCATAATCCTGGTAGACTATTCGCCTCTTTAAGCGCCTGCCTTGCAGGCGTTTTGCCCGTGGGTCTGGATTAGAACACCGATGAATTTATTaa < 2:140638/149‑1 (MQ=255) ccgcgATGCGATGCGTCAATCGGCCCGAGAACATCTGGCGGTAGCAATTCACCATAATCCAGGTAGACTATTCGCCTCTTTCAGCGCCTGCCTTGCAGGCGTTTTGCCCGTGGGTCTGGATTAGAACACCGATGAATTTATTAAAATCg < 1:250435/148‑1 (MQ=255) gatgcgatgAGTGAATAAACCCTGTAACATCTGGCGGTAGCAATTCACCATAATCCAGGTAGACTATTCGCCTCTTTCAGCGCCTGCCTTGCAGGCGTTTTGCCCGTGGGTCTGGATTAGAACCCCGATGAATTTATTAAAATCGCTgg > 2:316111/1‑149 (MQ=255) atgcgatgAGTGAATAAACCCTGTAACATCTGGCGGTAGCAATTCACCATAATACAGGTACAGTATTGGCCTCTGTCGGGGCCTGCCTTGGAGGCGTTTTGCCCGGTGGTCTGGATTAGAACCCCGCTGAAGTTATTAAGTTCGCTGgc > 1:234500/1‑149 (MQ=255) tgtgatgAGTGCATAAACCCTTTAACCTCTGGCGGTAGCAATTCACCATAATCCAGGTAGACTATTCGCCTCTTTCAGCGCCTGCCTTGCAGGCGTTTTGCCCGTGGGTCTGGATTAGAACACCGATGAATTTATTAAAATCGCTGgcc < 1:179520/146‑1 (MQ=255) tGAATAAACCCTGTAACATCTGGCGGTAGCAATTCACCATAATCCAGGTGGACTATTCGCCTCTTTCAGCGCCTGCCTTGCAGGCGTGTGGCCCGTGGGTCTGGATTAGAACACCGATGAGTTTATTAAAATCGCTGGCCGCCGTCAGc > 2:305351/1‑149 (MQ=255) tGAATAAACCCTGTAACATCTGGCGGTAGCAATGCACCATAATCCAGGTGGCCTATTAGCCTCTTTCAGCGCCTGCCTTGCAGGCGTTTTGCCCGTGGGTCTGGACTAGACCAGCGATGAATTTATTAAAATCTCCGGGCGCCCTCAGc > 1:315359/1‑149 (MQ=255) aaTAAACCCTGTAACATCTGGCGGTAGCAATTCACCATAATCCAGGTAGACTATTCGCCTCTTTCAGCGCCTGCCTTGCAGGCGTTTTGCCCGTGGGTCTGGATTAGAACACCGATGAATTTATTAAAATCGCTGGCCGCCGTCAGCTc > 1:37232/1‑149 (MQ=255) ccTGTATCATCTGTCGGTAGCAATTCACCATAATCCAGGTAGACTATTTGCCGCTTTCAGCGCCTGCCTTGCAGGCGTTTTGCCCGTGGGTCTGGATTAGAACTCCGATGAATTTATTAAAATCGCTGGCCGCCGTCAGCTCGATGAcc < 1:252043/149‑1 (MQ=255) cTGTAACATCTGGCGGTAGCAATTCACCATAATCCAGGTAGACTATTCGCCTCTTTCAGCGCCTGCCTTGCAGGCGTTTTGCCCGTGGGTCTGGATTAGAACACCGATGAATTTATTAAAATCGCTGGCCGCCGTCAGCTCGATGCCCa > 2:65155/1‑149 (MQ=255) cAATTCACCATAATCCAGGTAGACTATTCGCCTCTTTCAGCGCCTGCCTTGCAGGCGTTTTGCCCGTGGGTCTGGATTAGAACACCGATGAATTTATTAAAATCGCTGGCCGCCGTCAGCTCGATGACCATGTTTTCGCGTGTGCTTgg < 2:37232/149‑1 (MQ=255) ttCACCATAATCCTGGTAGACTATTCGCCTCTTTCAGCGCCTGCCTTGCAGGCGTGTATCCGTTGGGTCTGGATTAGAACACCGATGAATTTATTAAAATCGCTGGCCGCCGTCAGCTCGATGACCATGTTTTCGCGTGTGCTTGGCtt < 2:351057/149‑1 (MQ=255) | TTGAATGTGTGCAAAACCAGACAGTTCCGCAAACCGCCGGCGAACAGGCCGTGCTGGCGCAACGTATCGTTGACAAGATCTGGCGCGATGCGATGAGTGAATAAACCCTGTAACATCTGGCGGTAGCAATTCACCATAATCCAGGTAGACTATTCGCCTCTTTCAGCGCCTGCCTTGCAGGCGTTTTGCCCGTGGGTCTGGATTAGAACACCGATGAATTTATTAAAATCGCTGGCCGCCGTCAGCTCGATGACCATGTTTTCGCGTGTGCTTGGCTT > CP009273/1123082‑1123359 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 8 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |