| Predicted mutation | ||||||
|---|---|---|---|---|---|---|
| evidence | seq id | position | mutation | annotation | gene | description |
| RA | NC_002947 | 4,586,030 | 2 bp→TC | intergenic (+113/+101) | 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,030 | 0 | C | T | 96.3% | 81.4 / ‑4.5 | 27 | intergenic (+113/+102) | PP_4061/PP_4063 | hypothetical protein/long‑chain fatty acid‑‑CoA ligase |
| Reads supporting (aligned to +/- strand): ref base C (1/0); new base T (12/14); total (13/14) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 4.81e-01 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.89e-01 | |||||||||||
| * | NC_002947 | 4,586,031 | 0 | T | C | 96.3% | 82.4 / ‑5.9 | 27 | intergenic (+114/+101) | PP_4061/PP_4063 | hypothetical protein/long‑chain fatty acid‑‑CoA ligase |
| Reads supporting (aligned to +/- strand): ref base T (0/0); major base C (12/14); minor base G (1/0); total (13/14) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 4.81e-01 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 | |||||||||||
CCACCGGTGAAGGCTTAAGGGTAGACCACGGCACCGGCGTGAACAGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCT‑CTCTGC‑GGTAAACCCGCTCCTACCAGGG‑A‑CGTGCCGGCCTTGAATGCGGCGCTGCCCCTGGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCCGCGCCAG > NC_002947/4585900‑4586127 || ccACCGGTGAAGGCTTAAGGGTAGACCACGGCACCGGCGTGAACAGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGCGAAGATCTTGCATTGATCGGGCCGGCCT‑CTTCGCGGGTAAACCCGCTc > 4:194365/1‑148 (MQ=255) aCCGGTGAAGGCTTAAGGGTAGACCACGGCACCGGCGTGAACAGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCT‑CTTCGCGGGTAAACCCGCTCCTa > 4:132241/1‑149 (MQ=255) ggTGAAGGCTTAAGGGTAGACCACGGCACCGGCGTGAACAGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCT‑CTTCGCGGGTAAACCCGCTCCTACCa < 8:129026/149‑1 (MQ=255) gTGAAGGCTTAAGGGTAGACCACGGCACCGGCGTGAACAGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCT‑CTTCGCGGGTAAACCCGCTCCTACCAg > 2:190104/1‑149 (MQ=255) gTGAAGGCTTAAGGGTAGACCACGGCACCGGCGTGAACAGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCT‑CTTCGCGGGTAAACCCGCTCCTACCAg > 4:5699/1‑149 (MQ=255) gAAGGCTTAAGGGTAGACCACGGCACCGGCGTGAACAGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCT‑CTTCGCGGGTAAACCCGCTCCTACCAgg > 4:392929/1‑148 (MQ=255) aaGGCTTAAGGGTAGACCACGGCACCGGCGTGAACAGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCT‑CTTCGCGGGTAAACCCGCTc > 8:27930/1‑139 (MQ=255) aaGGCTTAAGGGTAGACCACGGCACCGGCGTGAACAGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCT‑CTTCGCGGGTAAACCCGCTc < 7:27930/139‑1 (MQ=255) ttAAGGGTAGACCACGGCACCGGCGTGAACAGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTAGTCGAG‑GGTAAAACCGCTCCTACCAgggtccgtg > 7:170564/1‑143 (MQ=255) cGGCACCGGCGTGAACAGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCT‑CTTCGCGGGTAAACCCGCTCCTACCAGGG‑A‑cc > 8:91991/1‑131 (MQ=255) cGGCACCGGCGTGAACAGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCT‑CTTCGCGGGTAAACCCGCTCCTACCAGGG‑A‑cc < 7:91991/132‑2 (MQ=255) aGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCT‑CTTCGCGGGTAAACCCGCTCCTACCAGGG‑ACCGTGCCGGCCTTGAATGCGGCGCTGCCCCTgg > 4:35871/1‑147 (MQ=39) aGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCT‑CTTCGCGGGTAAACCCGCTCCTACCAGGG‑ACCGTGCCGGCCTTGAATGCGGCGCTGCCCCTGGTa < 2:138149/149‑1 (MQ=39) gCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCT‑CTTCGCGGGTAAACCCGCTCCTACCAGGG‑ACCGTGCCGGCCTTGAATGCGGCGCTGCCCCTGGTAg < 3:35871/149‑1 (MQ=39) cTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCT‑CTTCGCGGGTAAACCCGCTcc > 8:159410/1‑103 (MQ=38) cTGATCTGCGCTCCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCT‑CTTCGCGGGTAAACCCGCTcc < 7:159410/103‑1 (MQ=38) tGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGCGAAGATCTTGCATTGATCGGGCCGGCCT‑CTTCGCGGGTAAACCCGCTCCTACCAGGG‑ACCGTGCCGGCCTTGAATGCGGCGCTGCCCCTGGTAGGa < 3:194365/149‑1 (MQ=38) cgcTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCT‑CTTCGCGGGTAAACCCGCTCCTACCAGGG‑ACCGTGCCGGCCTTGAATGCGGCGCTGCCCCTGGTAGGAGCgg < 3:5699/146‑1 (MQ=38) cgcTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCT‑CTTCGCGGGTAAACCCGCTCCTACCAGGG‑ACCGTGCCGGCCTTGAATGCGGCGCTGCCCCTGGTAGGAGCgg < 2:93013/146‑1 (MQ=38) gctgctCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCT‑CTTCGCGGGTAAACCCGCTCCTACCAGGG‑ACCGTGCCGGCCTTGAATGCGGCGCTGCCCCTGGTAGGAGCGGGttt < 6:38891/149‑1 (MQ=38) ctgctCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCT‑CTTCGCGGGTAAACCCGCTCCTACCAGGG‑ACCGTGCCGGCCTTGAATGCGGCGCTGCCCCTGGTAGGAGCGGGTTTa < 1:190104/149‑1 (MQ=38) tCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCT‑CTTCGCGGGTAAACCCGCTCCTACCAGGG‑ACCGTGCCGGCCTTGAATgc > 7:354225/1‑117 (MQ=37) tCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCT‑CTTCGCGGGTAAACCCGCTCCTACCAGGG‑ACCGTGCCGGCCTTGAATgc < 8:354225/117‑1 (MQ=37) aTTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCT‑CTTCGCGGGTAAACCCGCTCCTACCAGGG‑ACCGTGCCGGCCTTGAATGCGGCGCTGCCCCTGGTAGGAGCGGGTTTACCCGCGa < 3:150104/149‑1 (MQ=38) aCCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCA‑ATTCGCGGGTAAAACCGCTCCAACCAGGGCC‑CGGGGCGGCCTGTAATGAGGCGGGGCCGCGGGGGGGGGCGGGTGTGCCCGCGCAGAGgccgggcggggccg > 6:393820/1‑143 (MQ=255) ttGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCT‑CTTCGCGGGTAAACCCGCTCCTACCAGGG‑ACCGTGCCGGCCTTGAATGCGGCGCTGCCCCTGGt > 4:350468/1‑102 (MQ=25) ttGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCT‑CTTCGCGGGTAAACCCGCTCCTACCAGGG‑ACCGTGCCGGCCTTGAATGCGGCGCTGCCCCTGGt < 3:350468/102‑1 (MQ=25) || CCACCGGTGAAGGCTTAAGGGTAGACCACGGCACCGGCGTGAACAGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCT‑CTCTGC‑GGTAAACCCGCTCCTACCAGGG‑A‑CGTGCCGGCCTTGAATGCGGCGCTGCCCCTGGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCCGCGCCAG > NC_002947/4585900‑4586127 |
| 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: ‑ |