Predicted mutation | ||||||
---|---|---|---|---|---|---|
evidence | seq id | position | mutation | annotation | gene | description |
RA | NC_002947 | 307,866 | +C | intergenic (+881/‑988) | hslO → / → PP_0254 | molecular chaperone HslO/hypothetical protein |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_002947 | 307,865 | 1 | . | C | 96.3% | 87.8 / ‑2.9 | 27 | intergenic (+880/‑989) | hslO/PP_0254 | molecular chaperone HslO/hypothetical protein |
Reads supporting (aligned to +/- strand): ref base . (0/1); new base C (15/11); total (15/12) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 4.44e-01 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.33e-01 |
GCCATGTTCTCGGTGCAGAACTTCCTGCTGCCGGCTGCCGATGTACTGCCTATGCACTGCGCGGCCAACATCGGCGAAGAAGGCGATGTGACCTTGTTCTTCGGCCTGTCCGGCACTGGCAAGACTACCCTTTCGGCTGACGAAAG‑CGTTACCTGATCGGTGACGACGAGCACGGTTGGGGCGAAGGCGTTGTCTTCAACATGGAAGGCGGCTGCTACGCCAAGTGCATCGACCTGTCCGAGAAGAACGAGCCGGTCATCTGGAAAGCCATCAAGCACGGCGCAGTAC > NC_002947/307720‑308007 | gCCATGTTCTCGGTGCAGAACTTCCTGCTGCCGGCTGCCGATGTACTGCCTATGCACTGCGCGGCCAACATCGGCGAAGAAGGCGATGTGACCTTGTTCTTCGGCCTGTCCGGCACTGGCAAGACTACCCTTTCGGCTGACGAAAG‑cc < 1:173354/148‑2 (MQ=255) ccATGTTCTCGGTGCAGAACTTCCTGCTGCCGGCTGCCGATGTACTGCCTATGCACTGCGCGGCCAACATCGGCGAAGAAGGCGATGTGACCTTGTTCTTCGGCCTGTCCGGCACTGGCAAGACTACCCTTTCGGCTGACGAAAG‑ccg > 5:149963/1‑146 (MQ=255) cATGTTCTCGGTGCAGAACTTCCTGCTGCCGGCTGCCGATGTACTGCCTATGCACTGCGCGGCCAACATCGGCGAAGAAGGCGATTTGACCTTGTTCTTCGGCCTGTCCGGCACTGGCAAGACTACCCTTTCGTCAGCCGAAAG‑GGta < 5:61869/148‑2 (MQ=255) cATGTTCTCGGTGCAGAACTTCCTGCTGCCGGCTGCCGATGTACTGCCTATGCACTGCGCGGCCAACATCGGCGAAGAAGGCGATGTGACCTTGTTCTTCGGCCTGTCCGGCACTGGCAAGACTACCCTTTCGGCTGACGAAAG‑ccgt < 5:77405/148‑4 (MQ=255) tctcGGTGCAGAACTTCCTGCTGCCGGCTGCCGATGTACTGCCTATGCACTGCGCGGCCAACATCGGCGAAGAAGGCGATGTGACCTTGTTCTTCGGCCTGTCCGGCACTGGCAAGACTACCCTTTCGGCTGACGAAAGCCGTTAcc < 5:49365/147‑1 (MQ=255) aGAACTTCCTGCTGCCGGCTGCCGATGTACTGCCTATGCACTGCGCGGCCAACATCGGCGAAGAAGGCGATGTGACCTTGTTCTTCGGCCTGTCCGGCACTGGCAAGACTACCCTTTCGGCTGACGAAAGCCGTTACCTGATCGGTg < 4:143007/147‑1 (MQ=255) tgccggctgccgATGTACTGCCTATGCACTGCGCGGCCAACATCGGCGAAGAAGGCGATGTGACCTTGTTCTTCGGCCTGTCCGGCACTGGCAAGACTACCCTTTCGGCTGACGAAAGCCGTTACCTGATCGGTGACGACGAGCACgg < 5:148278/148‑1 (MQ=255) ccgATGTACTGCCTATGCACTGCGCGGCCAACATCGGCGAAGAAGGCGATGTGACCTTGTTCTTCGGCCTGTCCGGCACTGGCAAGACTACCCTTTCGGCTGACGAAAGCCGTTACCTGATCGGTGACGACGAGCACGGTTGGGGc > 2:119829/1‑146 (MQ=255) cgATGTACTGCCTATGCACTGCGCGGCCAACATCGGCGAAGAAGGCGATGTGACCTTGTTCTTCGGCCTGTCCGGCACTGGCAAGACTACCCTTTCGGCTGACGAAAGCCGTTACCTGATCGGTGCCGACGAGCACGGTTGGGTCGaa > 6:123888/1‑148 (MQ=255) aCTGCCTATGCACTGCGCGGCCAACATCGGCGAAGAAGGCGATGTGACCTTGTTCTTCGGCCTGTCCGGCACTGGCAAGACTACCCTTTCGGCTGACGAAAGCCGTTACCTGATCGGTGACGACGAGCACGGTTGGGGCGAAGGCGtt < 8:79372/148‑1 (MQ=255) aCTGCCTATGCACTGCGCGGCCAACATCGGCGAAGAAGGCGATGTGACCTTGTTCTTCGGCCTGTCCGGCACTGGCAAGACTACCCTTTCGGCTGACGAAAGCCGTTACCTGATCGGTGACGACGAGCACGGTTGGGGCGAAGGCGtt < 1:152543/148‑1 (MQ=255) cTATGCACTGCGCGGCCAACATCGGCGAAGAAGGCGATGTGACCTTGTTCTTCGGCCTGTCCGGCACTGGCAAGACTACCCTTTCGGCTGACGAAAGCCGTTACCTGATCGGTGACGACGAGCACGGTTGGGGCGAAGGCGTTGTCt > 6:117898/1‑147 (MQ=255) ggCCAACATCGGCGAAGAAGGCGATGTGACCTTGTTCTTCGGCCTGTCCGGCACTGGCAAGACTACCCTTTCGGCTGACGAAAGCCGTTACCTGATCGGTGACGACGAGCACGGTTGGGGCGAAGGCGTTGTCTTCAACATGGAAgga > 3:106815/1‑147 (MQ=255) ccAACATCGGCGAAGAAGGCGATGTGACCTTGTTCTTCGGCCTGTCCGGCACTGGCAAGACTACCCTTTCGGCTGACGAAAGCCGTTACCTGATCGGTGCCGACGAGCACGGTTGGGGCGAAGGCGTTGTCTTCAACATGGAAggcgg > 1:47091/1‑148 (MQ=255) aaCATCGGCGAAGAAGGCGATGTGACCTTGTTCTTCGGCCTGTCCGGCACTGGCAAGACTACCCTTTCGGCTGACGAAAGCCGTTACCTGATCGGTGACGACGAGCACGGTTGGGGCGAAGGCGTTGTCTTCAACATGGAAGGCGgct > 8:129571/1‑148 (MQ=255) cTTGTTCTTCGGCCTGTCCGGCACTGGCAAGACTACCCTTTCGGCTGACGAAAGCCGTTACCTGATCGGTGACGACGAGCACGGTTGGGGCGAAGGCGTTGTCTTCAACATGGAAGGCGGCTGCTACGCCAAGTGCATCGACCTGTcc > 3:90490/1‑148 (MQ=255) ttcttcGGCCTGTCCGGCACTGGCAAGACTACCCTTTCGGCTGACGAAAGCCGTTACCTGATCGGTGACGACGAGCACGGTTGGGGCGAAGGCGTTGTCTTCAACATGGAAGGCGGCTGCTACGCCAAGTGCATCGACCTGTCCgag < 5:6638/147‑1 (MQ=255) cttcGGCCTGTCCGGCACTGGCAAGACTACCCTTTCGGCTGACGAAAGCCGTTACCTGATCGGTGACGACGAGCACGGTTGGGGCGAAGGCGTTGTCTTCAACATGGAAGGCGGCTGCTACGCCAAGTGCATCGACCTGTCCGagaa > 3:102450/1‑147 (MQ=255) ggCCTGTCCGGCACTGGCAAGACTACCCTTTCGGCTGACGAAAGCCGTTACCTGATCGGTGACGACGAGCACGGTTGGGGCGAAGGCGTTGTCTTCAACATGGAAGGCGGCTGCTACGCCAAGTGCATCGACCTGTCCGAGAAGAAc < 3:80716/147‑1 (MQ=255) gCCTGTCCGGCACTGGCAAGACTACCCTTTCGGGTGACGAAAGCCGTTACCTGATCGGTGACGACGAGCACGGTTAGGGCGAAGGCGTTGTCTTCAACATGGAAGGCGGCTGCTACGCCAAGTGCATCGACCTGTCCGAGAAGAACGa > 4:150578/1‑148 (MQ=255) gCCTGTCCGGCACTGGCAAGACTACCCTTTCGGCTGACGAAAGCCGTTACCTGATCGGTGACGACGAGCACGGTTGGGGCGAAGGCGTTGTCTTCAACATGGAAGGCGGCTGCTACGCCAAGTGCATCGACCTGTCCGAGAAGAACGa > 5:86997/1‑148 (MQ=255) gCCTGTCCGGCACTGGCAAGACTACCCTTTCGGCTGACGAAAGCCGTTACCTGATCGGTGACGACGAGCACGGTTGGGGCGAAGGCGTTGTCTTCAACATGGAAGGCGGCTGCTACGCCAAGTGCATCGACCTGTCCGAGAAGAACGa > 3:153631/1‑148 (MQ=255) catGTCCGGCACTGGCAAGACTACCCTTTAGGCTGACGAAAGCCGTTAAATGATAGGTGAAGACGAGCACGGGTGGGGCGAAGGCGTTGTCTTCAACATGGAAGGCGGCTGCTACGCCAAGTGCATAGACCTGTCCGAGAAGAACGAg > 6:150935/3‑148 (MQ=255) cTGTCCGGCACTGGCAAGACTACCCTTTCGGCTGACGAAAGCCGTTACCTGATCGGTGACGACGAGCACGGTTGGGGCGAAGGCGTTGTCTTCAACATGGAAGGCGGCTGCTACGCCAAGTGCATCGACCTGTCCGAGAAGAACGAGc < 4:102450/148‑1 (MQ=255) cTGTCCGGCACTGGCAAGACTACCCTTTCGGCTGACGAAAGCCGTTACCTGATCGGTGACGACGAGCACGGTTGGGGCGAAGGCGTTGTCTTCAACATGGAAGGCGGCTGCTACGCCAAGTGCATCGACCTGTCCGAGAAGAACGAGc < 6:83481/148‑1 (MQ=255) cAAGACTCCCCTTTCGGCCGACGAAAGCCGTTACCTGATCGGTGACGCCGAGCACGGTTGGCGCGAAGCCGTTGTCTTCAACATGGAAGGCGGCGGCTACGCCAACGGCATCGACCTGTCCGCGAAGCACGAGCCGGTCATCTGGaaa < 8:71438/148‑1 (MQ=255) tACCCTTTCGGCTGACGAAAGCCGTTACCTGATCGGTGACGACGAGCACGGTTGGGGCGAAGGCGTTGTCTTCAACATGGAAGGCGGCTGCTACGCCAAGTGCATCGACCTGTCCGAGAAGAACGAGCCGGTCATCTGGAAAGCCAt > 3:1441/1‑147 (MQ=255) tGACGAAAGCCGTTACCTGATCGGTGACGACGAGCACGGTTGGGGCGAAGGCGTTGTCTTCAACATGGAAGGCGGCTGCTACGCCAAGTGCATCGACCTGTCCGAGAAGAACGAGCCGGTCATCTGGAAAGCCATCAAGCACGGCGCa < 1:38382/148‑1 (MQ=255) tGACGAAAGCCGTTACCTGATCGGTGACGACGAGCACGGTTGGGGCGAAGGCGTTGTCTTCAACATGGAAGGCGGCTGCTACGCCAAGTGCATCGACCTGTCCGAGAAGAACGAGCCGGTCATCTGGAAAGCCATCAAGCACGGCGCa > 5:35279/1‑148 (MQ=255) gAAAGCCGTTACCTGATCGGTGACGACGAGCACGGTTGGGGCGAAGGCGTTGTCTTCAACATGGAAGGCGGCTGCTACGCCAAGTGCATCGACCTGTCCGAGAAGAACGAGCCGGTCATCTGGAAAGCCATCAAGCACGGCGCAGTAc > 3:16126/1‑148 (MQ=255) | GCCATGTTCTCGGTGCAGAACTTCCTGCTGCCGGCTGCCGATGTACTGCCTATGCACTGCGCGGCCAACATCGGCGAAGAAGGCGATGTGACCTTGTTCTTCGGCCTGTCCGGCACTGGCAAGACTACCCTTTCGGCTGACGAAAG‑CGTTACCTGATCGGTGACGACGAGCACGGTTGGGGCGAAGGCGTTGTCTTCAACATGGAAGGCGGCTGCTACGCCAAGTGCATCGACCTGTCCGAGAAGAACGAGCCGGTCATCTGGAAAGCCATCAAGCACGGCGCAGTAC > NC_002947/307720‑308007 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 21 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |