Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | NC_002947 | 975,290 | A→T | 8.2% | intergenic (+136/+33) | ycfJ → / ← suhB | flagella biosyntheis regulator/bifunctional inositol monophosphatase/glycerol‑2‑phosphatase |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_002947 | 975,290 | 0 | A | T | 8.2% | 141.5 / 4.6 | 49 | intergenic (+136/+33) | ycfJ/suhB | flagella biosyntheis regulator/bifunctional inositol monophosphatase/glycerol‑2‑phosphatase |
Reads supporting (aligned to +/- strand): ref base A (20/25); new base T (2/2); total (22/27) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 1.00e+00 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.50e-01 |
AGTGATTTGCGGTGTCTGTGCCGGCCTCTTCGCGGGTAAACCCGATCCTACAGGGGCAATGCGATCCATGTAGGAGCGGGTTTACCCGCGAAGAGGTCGGTACAGGCCACGCAAAACCACCAGGCACAAAAAAGCACCTCAAGAGGTGC‑TTTTTGTTTGCCGGCAATCTGCGCTTAGCGCTTCATGCTTTCTGGCAGGTGCGGCTGGATGGCAGTCAGTACCGCCTTGAAGCACTTGATGTTGCCTGCAACGATGTGGCCTTTGTCGAGGAAGTCGTGGCCACCGTTGAA > NC_002947/975150‑975439 | aGTGATTTGCGGTGTCTGTGCCGGCCTCTTCGCGGGTAAACCCGATCCTACAGGGGCAATGCGATCCATGTAGGAGCGGGTTTACCCGCGAAGAGGTCGGTACAGGCCACGCAAAACCACCAGGCACAAAAAAGCACCTCAAGAGGTGC‑t > 2:535442/1‑150 (MQ=255) gTGATTTGCGGTGTCTGTGCCGGCCTCTTCGCGGGTAAACCCGATCCTACAGGGGCAATGCGATCCATGTAGGAGCGGGTTTACCCGCGAAGAGGTCGGTACAGGCCACGCAAAACCACCAGGCACAAAAAAGCACCTCAAGAGGTGC‑tt > 2:75287/1‑150 (MQ=255) tttGCGGTGTCTGTGCCGGCCTCTTCGCGGGTAAACCCGATCCTACAGGGGCAATGCGATCCATGTAGGAGCGGGTTTACCCGCGAAGAGGTCGGTACAGGCCACGCAAAACCACCAGGCACAAAAAAGCACCTCAAGAGGTGC‑TTtttg > 1:775045/1‑150 (MQ=255) ttGCGGTGTCTGTGCCGGCCTCTTCGCGGGTAAACCCGATCCTACAGGGGCAATGCGATCCATGTAGGAGCGGGTTTACCCGCGAAGAGGTCGGTACAGGCCACGCAAAACCACCAGGCAC‑AAAAAGCACCTCTTGAGGTGC‑ttttttgt < 2:81131/150‑4 (MQ=38) tgtgCCGGCCTCTTCGCGGGTAAACCCGATCCTACAGGGGCAATGCGATCCATGTAGGAGCGGGTTTACCCGCGAAGAGGTCGGTACAGGCCACGCAAAACCACCAGGCACAAAAAAGCACCTCAAGAGGTGC‑TTtttgtt > 2:443970/1‑141 (MQ=255) tgtgCCGGCCTCTTCGCGGGTAAACCCGATCCTACAGGGGCAATGCGATCCATGTAGGAGCGGGTTTACCCGCGAAGAGGTCGGTACAGGCCACGCAAAACCACCAGGCACAAAAAAGCACCTCAAGAGGTGC‑TTtttgtt < 1:443970/141‑1 (MQ=255) cGGCCTCTTCGCGGGTAAACCCGATCCTACAGGGGCAATGCGATCCATGTAGGAGCGGGTTTACCCGCGAAGAGGTCGGTACAGGCCACGCAAAACCACCAGGCACAAAAAAGCACCTCAAGAGGTGC‑TTTTTGTTTGCCGGCAATCTgc < 1:78635/150‑1 (MQ=255) gCCTCTTCGCGGGTAAACCCGATCCTACAGGGGCAATGCGATCCATGTAGGAGCGGGTTTACCCGCGAAGAGGTCGGTACAGGCCACGCAAAACCACCAGGCACAAAAAAGCACCTCAAGAGGTGC‑TTTTTGTTTGCCGGCAATCTgcgc < 2:101056/150‑1 (MQ=255) ttCGCGGGTAAACCCGATCCTACAGGGGCAATGCGATCCATGTAGGAGCGGGTTTACCCGCGAAGAGGTCGGTACAGGCCACGCAAAACCACCAGGCACAAAAAAGCCAAGTAGGAGGTGC‑TTTTTGTTTGCCGGCAATCTGCGCTTAgg > 1:8392/1‑149 (MQ=39) gcgGGTAAACCCGATCCTACAGGGGCAATGCGATCCATGTAGGAGCGGGTTTACCCGCGAAGAGGTCGGTACAGGCCACGCAAAACCACCAGGCACAAAAAAGCACCTCAAGAGGTGC‑TTTTTGTTTGCCGGCAATCTGCGCTTAGCGCt < 2:122112/150‑1 (MQ=255) gTAAACCCGATCCTACAGGGGCAATGCGATCCGTTTTCGAACGGGTTTACCCGCGAAGAGGTCGGTACAGGCCACGCAAAACCACCAGGCACAAAAAAGCACCTCAAGAGGTGC‑TTTTTGTTTGCCGGCAATCTGCGCTTAGCGCTTCAt < 2:496164/150‑1 (MQ=255) aaaCCCGATCCTACAGGGGCAATGCGATCCATGTAGGAGCGGGTTTACCCGCGAAGAGGTCGGTACAGGCCACGCAAAACCACCAGGCACAAAAAAGCACCTCAAGAGGTGC‑TTTTTGTTTGCCGGCAATCTGCGCTTAGCGCTTCATGc < 2:433965/150‑1 (MQ=255) aCCCGATCCTACAGGGGCAATGCGATCCATGTAGGAGCGGGTTTACCCGCGAAGAGGTCGGTACAGGCCACGCAAAACCACCAGGCACAAAAAAGCACCTCAAGAGGTGC‑TTTTTGTTTGCCGGCAATCTGCGCTTAGCGCTTCa < 1:718972/145‑1 (MQ=255) aCCCGATCCTACAGGGGCAATGCGATCCATGTAGGAGCGGGTTTACCCGCGAAGAGGTCGGTACAGGCCACGCAAAACCACCAGGCACAAAAAAGCACCTCAAGAGGTGC‑TTTTTGTTTGCCGGCAATCTGCGCTTAGCGCTTCa > 2:718972/1‑145 (MQ=255) cccGATCCTACAGGGGCAATGCGATCCATGTAGGAGCGGGTTTACCCGCGAAGAGGTCGGTACAGGCCACGCAAAACCACCAGGCACAAAAAAGCACCTCAAGAGGTGC‑TTTTTGTTTGCCGGCAATCTGCGCTTAGCGCTTCATGCttt > 1:128001/1‑150 (MQ=255) cGATCCTATAGGGGCAATGCGATCCATGTAGGAGCGGGTTTACCCGCGAAGAGGTCGGTACAGGCCACGCAAAACCACCAGGCACAAAAAAGCACCTCAAGAGGTGC‑TTTTTGTTTGCCGGCAATCTGCGCTTAGCGCTTCATGCTTTCt < 2:168030/150‑1 (MQ=255) aaTGCGATCCATGTAGGAGCGGGTTTACCCGCGAAGAGGTCGGTACAGGCCACGCAAAACCACCAGGCACAAAAAAGCACCTCAAGAGGTGC‑TTTTTGTTTGCCGGCAATCTGCGCTTAGCGCTTCATGCTTTCTGGCAGGTGCGGCTgg < 2:775045/150‑1 (MQ=255) gTAGGAGCGGGTTTACCCGCGAAGAGGTCGGTACAGGCCACGCAAAACCACCAGGCACAAAAAAGCACCTCAAGAGGTGC‑TTTTTGTTTGCCGGCAATCTGCGCTTAGCGCTTCATGCTTTCTGGCAGGTGCGGCTGGATGGcagtcagt < 1:428837/150‑1 (MQ=255) gAGCGGGTTTACCCGCGAAGAGGTCGGTACAGGCCACGCAAAACCACCAGGCACAAAAAAGCACCTCAAGAGGTGC‑TTTTTGTTTGCCGGCAATCTGCGCTTAGCGCTTCATGCTTTCTGGCAGGTGCGGCTGGATGGCAGTCAGTACCg > 2:271577/1‑150 (MQ=255) gCGGGTTTACCCGCGAAGAGGTCGGTACAGGCCACGCAAAACCACCAGGCACAAAAAAGCACCTCAAGAGGTGC‑TTTTTGTTTGCCGGCAATCTgc < 1:788112/96‑1 (MQ=255) gCGGGTTTACCCGCGAAGAGGTCGGTACAGGCCACGCAAAACCACCAGGCACAAAAAAGCACCTCAAGAGGTGC‑TTTTTGTTTGCCGGCAATCTgc > 2:788112/1‑96 (MQ=255) gTTTACCCGCGAAGAGGTCGGTACAGGCCACGCAAAACCACCAGGCACAAAAAAGCACCTCAAGAGGTGC‑TTTTTGTTTGCCGGCAATCTGCGCTTAGCGCTTCATGCTTTCTGGCAGGTGCGGCTGGATGGCAGTCAGTACCGCCTTGa > 1:852923/1‑150 (MQ=255) ccGCGAAGAGGTCGGTACAGGCCACGCAAAACCACCAGGCACAAAAAAGCACCTCAAGAGGTGC‑ttttt < 1:790857/69‑1 (MQ=255) ccGCGAAGAGGTCGGTACAGGCCACGCAAAACCACCAGGCACAAAAAAGCACCTCAAGAGGTGC‑ttttt > 2:790857/1‑69 (MQ=255) gcgAAGAGGTCGGTACAGGCCACGCAAAACCACCAGGCACAAAAAAGCACCTCAAGAGGTGC‑TTTTTGTTTGCCGGCAATCTGCGCTTAGCGCTTCATGCTTTCTGGCAGGTGCGGCTGGATGGCAGTCAGTACCGCCTTGAAGCACTTg < 2:128001/150‑1 (MQ=255) tCGGTACAGGCCACGCAAAACCACCAGGCACAAAAAAGCACCTCAAGAGGTGC‑TTTTTGTTTGCCGGCAATCTGCGCTTAGCGCTTCATGCTTTCTGGCAGGTGCGGCTGGATGGCAGTCAGTACCGCCTTGAAGCACTTGATGTTGcc < 2:786201/149‑1 (MQ=255) tCGGTACAGGCCACGCAAAACCACCAGGCACAAAAAAGCACCTCAAGAGGTGC‑TTTTTGTTTGCCGGCAATCTGCGCTTAGCGCTTCATGCTTTCTGGCAGGTGCGGCTGGATGGCAGTCAGTACCGCCTTGAAGCACTTGATGTTGcc > 1:786201/1‑149 (MQ=255) aGGCCACGCAAAACCACCAGGCACAAAAAAGCACCTCAAGAGGTGC‑TTTTTGTTTGCCGGCAATCTGCGCTTAGCGCTTCATGCTTTCt < 2:681519/89‑1 (MQ=255) aGGCCACGCAAAACCACCAGGCACAAAAAAGCACCTCAAGAGGTGC‑TTTTTGTTTGCCGGCAATCTGCGCTTAGCGCTTCATGCTTTCt > 1:681519/1‑89 (MQ=255) ggCCACGCAAAACCACCAGGCACAAAAAAGCACCTCAAGAGGTGC‑TTTTTGTTTGCCGGCAATCTGCGCTTAGCGCTTCATGCTTTCTGGCAGGTGCGGCTGGATGGCAGTCAGTACCGCCTTGAAGCACTTGATGTTGCCTGCAACGAt < 2:653892/150‑1 (MQ=255) ggCCACGCAAAACCACCAGGCACAAAAAAGCACCTCAAGAGGTGC‑TTTTTGTTTGCCGGCAATCTGCGCTTAGCGCTTCATGCTTTCTGGCAGGTGCGGCTGGATGGCAGTCAGTACCGCCTTGAAGCACTTGATGTTGCCTGCAACGAt < 1:535442/150‑1 (MQ=255) aaCCACCAGGCACAAAAAAGCACCTCAAGAGGTGC‑TTTTTGTTTGCCGGCAATCTGCGCTTAGCGCTTCATGCTTTCTGGCAGGTGCGGCTGGATGGCAGTCAGTACCGCCTTGAAGCACTTGATGTTGCCTGCAACGATGTGGCCTTTg > 1:105993/1‑150 (MQ=255) accaGGCACAAAAAAGCACCTCAAGAGGTGC‑TTTTTGTTTGCCGGCAATCTGCGCTTAGCGCTTCATGCTTTCTGGCAGGTGCGGCTGGATGGCAGTCAGTACCGCCTTGAAGCACTTGATGTTGCCTGCAACGATGTGGCCTTTGTCGa < 2:804696/150‑1 (MQ=255) accaGGCACAAAAAAGCACCTCAAGAGGTGC‑TTTTTGTTTGCCGGCAATCTGCGCTTAGCGCTTCATGCTTTCTGGCAGGTGCGGCTGGATGGCAGTCAGTACCGCCTTGAAGCACTTGATGTTGCCTGCAACGATGTGGCCTTTGTCGa < 2:852923/150‑1 (MQ=255) ccaGGCACAAAAAAGAGCTTCAAGAGGTGC‑TTTTTGTTTGCCGGCAATCTGCGCTTAGCGCTTCATGCTTTCTGGCAGGTGCGGCTGGATGGCAGTCAGTACCGCCTTGAAGCACTTGATGTTGCCTGCAACGATGTGGCCTTTGTCGAg < 2:326471/150‑1 (MQ=255) ggCACAAAAAAGCACCTCAAGAGGTGC‑TTTTTGTTTGCCGGCAATCTGCGCTTAGCGCTTCATGCTTTCTGGCAGGTGCGGCTGGATGGCAGTCAGTACCGCCTTGAAGCACTTGATGTTGCCTGCAACGATGTGGCCTTTGTCGAGGaa < 1:75287/150‑1 (MQ=255) gCACAAAAAAGCACCTCAAGAGGTGC‑TTTTTGTTTGCCGGCAATCTGCGCTTAGCGCTTCATGCTTTCTGGCAGGTGCGGCTgg > 2:571403/1‑84 (MQ=255) gCACAAAAAAGCACCTCAAGAGGTGC‑TTTTTGTTTGCCGGCAATCTGCGCTTAGCGCTTCATGCTTTCTGGCAGGTGCGGCTgg < 1:571403/84‑1 (MQ=255) gCACAAAAAAGCACCTCAAGAGGTGC‑TTTTTGTTTGCCGGCAATCTGCGCTTAGCGCTTCATGCTTTCTGGCAGGTGCGGCTGGATGGCAGTCAGTACCGCCTTGAAGCACTTGATGTTGCCTGCAACGATGTGGCCTTTGTCGAGGAAg > 2:371257/1‑150 (MQ=255) acaAAAAAGCACCTCAAGAGGTGC‑TTTTTGTTTGCCGGCAATCTGCGCTTAGCGCTTCATGCTTTCTGGCAGGTGCGGCTGGATGGCAGTCAGTACCGCCTTGAAGCACTTGATGTTGCCTGCAACGATGTGGCCTTTGTCGAGGAAGTc > 2:433956/1‑150 (MQ=255) aacaaaaaGCACCTCTTGAGGTGCTTTTTTGTTTGCCGGCAATCTGCGCTTAGCGCTTCATGCTTTCTGGCAGGTGCGGCTGGATGGCAGTCAGTACCGCCTTGAAGCACTTGATGTTGCCTGCAACGATGTGGCCTTTGTCGAGGAAGt > 1:6697/4‑150 (MQ=255) caaaaaGCACCTCTTGAGGTGC‑TTTTTGTTTGCCGGCAATCTGCGCTTAGCGCTTCATGCTTTCTg > 2:111612/2‑66 (MQ=37) caaaaaGCACCTCTTGAGGTGC‑TTTTTGTTTGCCGGCAATCTGCGCTTAGCGCTTCATGCTTTCTg < 1:111612/65‑1 (MQ=37) aaaaaaGCACCTCAAGAGGTGC‑TTTTTGTTTGCCGGCAATCTGCGCTTAGCGCTTCATGCTTTCTGGCAGGTGCGGCTGGATGGCAGTCAGTACCGCCTTGAAGCACTTGATGTTGCCTGCAACGAt > 2:176542/1‑127 (MQ=255) aaaaaaGCACCTCAAGAGGTGC‑TTTTTGTTTGCCGGCAATCTGCGCTTAGCGCTTCATGCTTTCTGGCAGGTGCGGCTGGATGGCAGTCAGTACCGCCTTGAAGCACTTGATGTTGCCTGCAACGAt < 1:176542/127‑1 (MQ=255) aaaaaGCACCTCAAGAGGTGC‑TTTTTGTTTGCCGGCAATCTGCGCTTAGCGCTTCATGCTTTCTGGCAGGTGCGGCTGGATGGCAGTCAGTACCGCCTTGAAGCACTTGATGTTGCCTGCAACGATGTGGCCTTTGTCGAGGAAGTCGTg < 1:433956/150‑1 (MQ=255) aaaaGCACCTCAAGAGGTGC‑TTTTTGTTTGCCGGCAATCTGCGCTTAGCGCTTCATGCTTTCTGGCAGGTGCGGCTGGATGGCAGTCAGTACCGCCTTGAAGCACTTGATGTTGCCTGCAACGATGTGGCCTTTGTCGAGGAAGTCGTgg < 2:105993/150‑1 (MQ=255) cTCAAGAGGTGC‑TTTTTGTTTGCCGGCAATCTGCGCTTAGCGCTTCATGCTTTCTGGCAGGTGCGGCTGGATGGCAGTCAGTACCGCCTTGAAGCACTTGATGTTGCCTGCAACGATGTGGCCTTTGTCGAGGAAGTCGTGGCCACCGtt > 1:453373/1‑150 (MQ=255) cAAGAGATGC‑TTTTTGTTTGCCGGCAATCTGCGCTTAGCGCTTAATGATTTCTGGAACGAGAGGAGAAATGGCTGTCAGTACCGGCTTGGAGAAGGTGAGTTTGTTTAGATGGATGTGGCCTTTGTCGAGGAAGTCGTGGCCGCCGTTGa > 2:735753/1‑150 (MQ=255) aaGAGGTGC‑TTTTTGTTTGCCGGCAATCTGCGCTTAGCGCTTCATGCTTTCTGGCAGGTGCGGCTGGATGGCAGTCAGTACCGCCTTGAAGCACTTGATGTTGCCTGCAACGATGTGGCCTTTGTCGAGGAAGTCGTGGCCACCGTTGaa < 1:501307/150‑1 (MQ=255) | AGTGATTTGCGGTGTCTGTGCCGGCCTCTTCGCGGGTAAACCCGATCCTACAGGGGCAATGCGATCCATGTAGGAGCGGGTTTACCCGCGAAGAGGTCGGTACAGGCCACGCAAAACCACCAGGCACAAAAAAGCACCTCAAGAGGTGC‑TTTTTGTTTGCCGGCAATCTGCGCTTAGCGCTTCATGCTTTCTGGCAGGTGCGGCTGGATGGCAGTCAGTACCGCCTTGAAGCACTTGATGTTGCCTGCAACGATGTGGCCTTTGTCGAGGAAGTCGTGGCCACCGTTGAA > NC_002947/975150‑975439 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |