Predicted mutation | ||||||
---|---|---|---|---|---|---|
evidence | seq id | position | mutation | annotation | gene | description |
RA | NC_002947 | 4,586,057 | +C | intergenic (+140/+75) | 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 | 90.9% | 66.2 / 1.1 | 22 | intergenic (+139/+76) | PP_4061/PP_4063 | hypothetical protein/long‑chain fatty acid‑‑CoA ligase |
Reads supporting (aligned to +/- strand): ref base . (2/0); new base C (7/13); total (9/13) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 1.56e-01 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.21e-01 |
GCTTAAGGGTAGACCACGGCACCGGCGTGAACAGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTCTG‑C‑GG‑TAAACCCGC‑TCCTACCAGGG‑A‑CGTGCCGGCCTTGAATGCGGCGCTGCCCC‑TGGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCCGCGCCAGCACATCAACCCGCAGATGCAGCCGCAATCTCCGCCACACTGATCTCACGCATGCGGAACTTC > NC_002947/4585912‑4586189 | gCTTAAGGGTAGACCACGGCACCGGCGTGAACAGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCG‑CGGG‑TAAACCCGC‑TCCTACCAGGG‑A‑c < 6:204182/147‑1 (MQ=255) tAAGGGTAGACCACGGCACCGGCGTGAACAGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGAACGGGTGCTTGTCAAGGTGAACATCTTGCATTTAACGGTCCGGCCTCTTCG‑C‑GGCTAAACCCGC‑TCCTACCAGGG‑A‑ccgt > 8:351827/1‑144 (MQ=255) aGGGTAGACCACGGCACCGGCGTGAACAGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCG‑CGGG‑TAAACCCGC‑TCCTACCAGGG‑A‑ccgt < 5:128499/145‑4 (MQ=39) aCCACGGCACCGGCGTGAACAGCTGATCTGCGCTGCTCCCTTGAAAGGTTGATGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCG‑CGGG‑TAAACCCGC‑TCCTACCAGGG‑ACCGTGCCGGCCt < 8:196798/146‑1 (MQ=38) ccACGGCACCGGCGTGAACAGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCG‑CGGG‑TAAACCCGC‑TCCTACCAGGG‑ACCGTGCCGGCCt < 2:17457/145‑1 (MQ=38) ccACGGCACCGGCGTGAACAGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCG‑CGGG‑TAAACCCGC‑TCCTACCAGGG‑ACCGTGCCGGCCt > 1:17457/1‑145 (MQ=38) ccGGCGTGAACAGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCG‑CGGG‑TAAACCCGC‑TCCTACCAGGG‑ACCGTGCCGGCCTTGAATgcggc > 3:141041/1‑147 (MQ=39) aGCTGATCTGCGCTGCTCGATTGAATGGTTGATGATGAACGGGTGCTTGTCAATGTTAAGATCTTGCATTTATCGGGCCGGCCTCTTCG‑CGGG‑TAAACCCGC‑TCCTACCAGGG‑ACCGTGCCGGCCTTGAATGCGGCGCTGCCCC‑Tgg < 7:351827/147‑1 (MQ=25) aGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCG‑CGGG‑TAAACCCGC‑TCCTACCAGGG‑ACCGTGCCGGCCTTGAATGCGGCGCTGCCCC‑Tgg < 6:346927/147‑1 (MQ=39) gctgctCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCG‑CGGG‑TAAACCCGC‑TCCTACCAGGG‑ACCGTGCCGGCCTTGAATGCGGCGCTGCCCC‑TGGTAGGAGCggg < 2:182093/146‑1 (MQ=38) gctgctCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCG‑CGGG‑TAAACCCGC‑TCCTACCAGGG‑ACCGTGCCGGCCTTGAATGCGGCGCTGCCCC‑TGGTAGGAGCggg < 6:1938/146‑1 (MQ=38) gctgctCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCG‑CGGG‑TAAACCCGC‑TCCTACCAGGG‑ACCGTGCCGGCCTTGAATGCGGCGCTGCCCC‑TGGTAGGAGCgg < 1:210208/145‑1 (MQ=38) gctgctCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCG‑CGGG‑TAAACCCGC‑TCCTACCAGGG‑ACCGTGCCGGCCTTGAATGCGGCGCTGCCCC‑TGGTAGGAGCg < 4:141041/144‑1 (MQ=38) gctgctCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCG‑CGGG‑TAAACCCGC‑TCCTACCAGGG‑ACCGTGCCGGCCTTGAATGCGGCGCTGCCCC‑TGGTAGGAGCg < 3:27872/144‑1 (MQ=38) gctgctCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGAGCGGGCCGGCCTCTTCG‑CGGG‑TAAACCCGC‑TCCTACCAGGG‑ACCGTGCCGGCCTTGAATGCGGCGCTGCCCC‑TGGTAGGAg < 8:149560/142‑1 (MQ=38) tCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCG‑CGGG‑TAAACCCGC‑TCCTACCAGGG‑ACCGTGCCGGCCTTGAATGCGGCGCTGCCCC‑TGGTAGGAGCGGGTTTAcc < 4:250984/147‑1 (MQ=38) tCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCG‑CGGG‑TAAACCCGC‑TCCTACCAGGG‑ACCGTGCCGGCCTTGAATGCGGCGCTGCCCC‑TGGTAGGAGCGGGTTTAcc < 6:105345/147‑1 (MQ=38) tGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGGCAAGATCTTGCATTGATCGGGCCGGCCTCTTCG‑CGGG‑TAAACCCGC‑TCCTACCAGGG‑ACCGTGCCGGCCTTGAATGCGGCGCTGCCCC‑TGGTAGGAGCGGGTTTACCCGCGa < 7:75048/147‑1 (MQ=37) gTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCACTTAGCC‑GG‑TAAAACCGCAGCAAACCCGGC‑C‑CCTGGCGGGCGTGTATGCGGCGCGGCCCC‑GGGGAGGGGCGGGGTTAGCCGCCGAGGGGGCGGCCGCGCCCGc > 2:314094/1‑145 (MQ=255) gcggg‑TAAACCCGC‑TCCTACCAGGGCC‑CGTGCCGGCCTTGAATGCGGCGCTGCCCC‑TGGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCCGCGCCAGCACATCAACCCGCAGATGCAGCCGCAATCTCCGCCACACTGATCTCAc > 3:212663/4‑146 (MQ=255) g‑TAAACCCGC‑TCCTACCAGGG‑ACCGTGCCGGCCTTGAATGCGGCGCTGCCCCGTGGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCCGCGCCAGCACATCAACCCGCAGATGCAGCCGCAATCTCCGCCACACTGATCTCACGCAt > 1:101588/1‑147 (MQ=255) g‑TAAACCCGC‑TCCTACCAGGG‑ACCGTGCCGGCCTTGAATGCGGCGCTGCCCC‑TGGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCCGCGCCAGCACATCAACCCGCAGATGCAGCCGCAATCTCCGCCACACTGATCTCACGCATg > 8:173480/1‑147 (MQ=255) aaaCCCGC‑TCCTACCAGGG‑ACCGTGCCGGCCTTGAATGCGGCGCTGCCCC‑TGGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCCGCGCCAGCCCATCAACCCGCAGATGCAGCCGCAATCTCCGCCACACTGATCTCACGCATGc > 5:269061/1‑146 (MQ=255) gC‑TCCTACCAGGG‑ACCGTGCCGGCCTTGAATGCGGCGCTGCCCC‑TGGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCCGCGCCAGCACATCAACCCGCAGATGCAGCCGCAATCTCCGCCACACTGATCTCACGCATGCGGAACtt > 3:320086/1‑147 (MQ=255) c‑TCCTACCAGGG‑ACCGTGCCGGCCTTGAATGCGGCGCTGCCCC‑TGGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCCGCGCCAGCACATCAACCCGCAGATGCAGCCGCAATCTCCGCCACACTGATCTCACGCATGCGGAACTTc > 5:288169/1‑147 (MQ=255) | GCTTAAGGGTAGACCACGGCACCGGCGTGAACAGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTCTG‑C‑GG‑TAAACCCGC‑TCCTACCAGGG‑A‑CGTGCCGGCCTTGAATGCGGCGCTGCCCC‑TGGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCCGCGCCAGCACATCAACCCGCAGATGCAGCCGCAATCTCCGCCACACTGATCTCACGCATGCGGAACTTC > NC_002947/4585912‑4586189 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |