Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | NC_002947 | 4,586,056:1 | +C | 55.6% | intergenic (+139/+76) | PP_4061 → / ← PP_4063 | hypothetical protein/long‑chain fatty acid‑‑CoA ligase |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_002947 | 4,586,056 | 1 | . | C | 55.6% | 45.9 / 42.8 | 27 | intergenic (+139/+76) | PP_4061/PP_4063 | hypothetical protein/long‑chain fatty acid‑‑CoA ligase |
Reads supporting (aligned to +/- strand): ref base . (6/6); new base C (3/12); total (9/18) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 1.27e-01 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 |
GGCTTAAGGGTAGACCACGGCACCGGCGTGAACAGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTCTGC‑GGTAAACCCGCTCCTACCAGGG‑A‑C‑GTGCCGGCCTTGAATGC‑GGCGCTGCCCCTGGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCCGCGCCAGCACATCAACCCGCAGATGCAGCCGCAATCTCCGCCACACTGATCTCACGCATGCGGAACTTCTGCACCTTG > NC_002947/4585911‑4586198 | ggCTTAAGGGTAGACCACGGCACCGGCGTGAACAGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACCAGGG‑A‑cc > 1:55212/1‑148 (MQ=255) ggCTTAAGGGTAGACCACGGCACCGGCGTGAACAGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACCAGGG‑A‑cc > 6:201213/1‑148 (MQ=255) ggCTTAAGGGTAGACCACGGCACCGGCGTGAACAGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACCAGGG‑A‑cc < 4:75361/149‑2 (MQ=255) ggTAGACCACGGCACCGGCGTGAACAGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACCAGGG‑ACC‑GTGCCGGc < 7:87539/149‑1 (MQ=39) ggaGACCACGGCACCGGCGTGAACAGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACCAGGG‑ACC‑GTGCCGGcc < 4:193374/147‑1 (MQ=38) gTAGACCACGGCACCGGCGTGAACAGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACCAGGG‑ACC‑GTGCCGGcc < 6:214984/149‑1 (MQ=39) gTAGACCACGGCACCGGCGTGAACAGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACCAGGG‑ACC‑GTGCCGGc < 8:261577/148‑1 (MQ=39) gACCACGGCACCGGCGTGAACAGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACCAGGG‑ACC‑GTGCCGGCCtt < 2:55212/148‑1 (MQ=39) gACCACGGCACCGGCGTGAACAGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACCAGGG‑ACC‑GTGCCGGCCtt < 3:86293/148‑1 (MQ=39) cGGCGTGAACAGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACCAGGG‑ACC‑GTGCCGGCCTTGAATgc‑g < 5:270477/144‑1 (MQ=38) cGGCGTGAACAGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACCAGGG‑ACC‑GTGCCGGCCTTGAATgc‑g > 6:270477/1‑144 (MQ=38) ccggcaaggtgcagaagttccgcatgcgtgagatcagtgtggcggagattgcggctgcatctgcgggttgatgtgctggcgcggccggcctcttcgcgGGTAAACCCGCTCCTACCAGGG‑A‑CCGTGCCGGCCTTGAATGC‑GGCGCTGc > 3:202028‑M1/99‑148 (MQ=255) aaCAGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACCAGGG‑ACC‑GTGCCGGCCTTGAATGC‑GGCGCTGCCCCt < 1:127593/148‑1 (MQ=39) aaCAGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACCAGGG‑A‑ccgtgc < 1:221756/123‑6 (MQ=38) aaCAGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACCAGGG‑A‑ccgtgc > 2:221756/1‑118 (MQ=38) cAGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACCAGGG‑ACC‑GTGCCGGCCTTGAATGC‑GGCGGTGCCCCTGGt < 6:143017/149‑1 (MQ=38) cAGCTGATCTGCGCTGCTCCATTGAAAGGGTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACCAGGG‑ACC‑GTGCCGGCCTTGAATGC‑GGCGCTGCCCCTGGt < 5:201213/149‑1 (MQ=38) gtgcagaagttccgcatgcgtgagatcagtgtggcggagattgcggctgcatctgcgggttgatgtgctggcgcggccggcctcttcgcgGGTAAACCCGCTCCTACCAGGG‑A‑CCGTGCCGGCCTTGAATGC‑GGCGCTGCCCCTGGTAg > 1:51126‑M1/91‑149 (MQ=255) tgcagaagttccgcatgcgtgagatcagtgtggcggagattgcggctgcatctgcggggtgatgtgctggcgcggccggcctcttcgcgGGTAAACCCGCTCCTACCAGGG‑A‑CCGTGCCGGCCTTGAATGC‑GGCGCTGCCCCTGGTAgg > 8:66416‑M1/90‑149 (MQ=255) cTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCGCGGGTAAAACCGCTCCTACCAGGGAA‑C‑GTGCCGGCCTTGAATGC‑GGCGCTGCCCCGGGGAg > 7:147917/1‑148 (MQ=38) gctgctCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACCAGGG‑ACC‑GTGCCGGCCTTGAATGC‑GGCGCTGccc < 7:84050/132‑1 (MQ=38) gctgctCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACCAGGG‑ACC‑GTGCCGGCCTTGAATGC‑GGCGCTGccc > 8:84050/1‑132 (MQ=38) tgcgtgagatcagtgtggcggagattgcggctgcatctgcgggttgatgtgctggcgcggccggcctcttcgcgGGTAAACCCGCTCCTACCAGGG‑A‑CCGTGCCGGCCTTGAATGC‑GGCGCTGCCCCTg < 2:112221‑M1/55‑1 (MQ=255) tgcgtgagatcagtgtggcggagattgcggctgcatctgcgggttgatgtgctggcgcggccggcctcttcgcgGGTAAACCCGCTCCTACCAGGG‑A‑CCGTGCCGGCCTTGAATGC‑GGCGCTGCCCCTg > 1:112221‑M1/75‑129 (MQ=255) ggTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACCAGGG‑ACC‑GTGCCGGCCTTGAATGC‑GGCGCt < 8:215116/98‑1 (MQ=25) cggccTCTTCGCGGGTAAAACCGCTGCTACCAGGG‑ACC‑GTGCCGGCCTTGAATGC‑GGCGCTGCCCCGGGTAGGAGCGGGTTTACCCGCGAAGCGGCCGGCCGCGCCAGCACATCACCCCGCAGCTGCAGCCGCACTCTCCGCCacacg > 4:97129/1‑147 (MQ=37) cgcgGGTAAACCCGCTCCTACCAGGG‑G‑CAGCGCCG‑CATTCAAGGCCGGCACGGTCCCAGGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCCGCGCCAGCACATCAACCCGCAGATGCAGCCGCAATCTCCGCCACACTGATCTCACg < 4:63374‑M1/144‑1 (MQ=255) cgcgGGTAAACCCGCTCCTACCAGGG‑A‑CAGTGCCGGCCTTGAATGC‑GGCGCTGCCCCGGGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCCGCGCCAGCACATCAACCCGCACATGCAGCCGCAATCTCCGCCACACTGAACACACg < 2:175487‑M1/144‑1 (MQ=255) GGTAAACCCGCTCCTACCAGGG‑G‑CAGCGCCG‑CATTCAAGGCCGGCACGGTCCCTGGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCCGCGCCAGCACATCAACCCGCAGATGCAGCCGCAATCTCCGCCACACTGATCTCACGCATGc > 3:63374‑M1/1‑149 (MQ=255) CTACCAGGG‑A‑CCGTGCCGGCCTTGAATGC‑GGCGCTGCCCCTGGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCCGCGCCAGCACATCAACCCGCAGATGCAGCCGCAATCTCCGCCACACTGATCTCACGCATGCGGAACTTCTGCAc < 4:169577‑M1/149‑1 (MQ=255) ACCAGGG‑A‑CCGTGCCGGCCTTGAATGC‑GGCGCTGCCCCTGGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCCGCGCCAGCACATCAACCCGCAGATGCAGCCGCAATCTCCGCCACACTGATCTCACGCATGCGGAACTTCTGCACCt < 1:76664‑M1/149‑1 (MQ=255) CAGGG‑A‑CCGTGCCGGCCTTGAATGC‑GGCGCTGCCCCTGGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCCGCGCCAGCACATCAACCCGCAGATGCAGCCGCAATCTCCGCCACACTGATCTCACGCATGCGGAACTTCTGCACCTTg < 4:234680‑M1/149‑1 (MQ=255) | GGCTTAAGGGTAGACCACGGCACCGGCGTGAACAGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTCTGC‑GGTAAACCCGCTCCTACCAGGG‑A‑C‑GTGCCGGCCTTGAATGC‑GGCGCTGCCCCTGGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCCGCGCCAGCACATCAACCCGCAGATGCAGCCGCAATCTCCGCCACACTGATCTCACGCATGCGGAACTTCTGCACCTTG > NC_002947/4585911‑4586198 |
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: ‑ |