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 | 97.1% | 109.8 / ‑3.1 | 34 | intergenic (+139/+76) | PP_4061/PP_4063 | hypothetical protein/long‑chain fatty acid‑‑CoA ligase |
Reads supporting (aligned to +/- strand): ref base . (1/0); new base C (14/19); total (15/19) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 4.41e-01 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.56e-01 |
ACCACGGCACCGGCGTGAACAGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTCTGC‑GGTAAACCCGCTC‑CTACCAGGGA‑C‑GTGCCGGCCTTGAATG‑CGGCGCTGCCCCTGGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCCGCGCCAGCACATCAACCCGCAGATGCAGCCGCAATCTCCGCCACACTGATCTCACGCATGCGGAACTT > NC_002947/4585924‑4586188 | aCCACGGCACCGGCGTGAACAGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCGCGGGTAAACCCGCTC‑CTACCAGGGACC‑GTGCCGGCCTTg > 1:185933/1‑148 (MQ=39) aCCACGGCACCGGCGTGAACAGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCGCGGGTAAACCCGCTC‑CTACCAGGGACC‑GTGCCGGCCTTg < 2:185933/148‑1 (MQ=39) cACGGCACCGGCGTGAACAGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCGCGGGTAAACCCGCTC‑CTACCAGGGACC‑GTGCCGGCCTTGCAt > 2:315151/1‑149 (MQ=38) gCACCGGCGTGAACAGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCGCGGGTAAACCCGCTC‑CTACCAGGGACC‑GTGCCGGCCTTGAATg‑cgg < 2:225651/149‑1 (MQ=39) gTGAACAGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCGCGGGTAAACCCGCTC‑CTACCAGGGACC‑GTGCCGGCCtt < 8:90847/133‑1 (MQ=38) gTGAACAGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCGCGGGTAAACCCGCTC‑CTACCAGGGACC‑GTGCCGGCCtt > 7:90847/1‑133 (MQ=38) tGAACAGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGTTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCGCGGGTAAACCCGCTC‑CTACCAGGGACC‑GTGCCGGCCTTGAATG‑CGGCGCTGcccc < 6:244956/149‑1 (MQ=38) tGAACAGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCGCGGGTAAACCCGCTC‑CTACCAGGGACC‑GTGCCGGCCtt < 8:146248/132‑1 (MQ=38) tGAACAGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCGCGGGTAAACCCGCTC‑CTACCAGGGACC‑GTGCCGGCCtt > 7:146248/1‑132 (MQ=38) tGAACAGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCGTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCGCGGGTAAACCCGCTC‑CTACCAGGGACC‑GTGCCGGCCTTGAATG‑CGGCGCTGcccc < 6:250657/149‑1 (MQ=38) cgcTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCGCGGGTAAACCCGCTC‑CTACCAGGGACC‑GTGCCGGCCTTGAATG‑CGGCGCTGCCCCTGGTAGGAGCgg < 1:182215/146‑1 (MQ=38) ctCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCGCGGGTAAACCCGCTC‑CTACCAGGGACC‑GTGCCGGCCTTGAATG‑CGGCGCTGCCCCTg < 3:288151/131‑1 (MQ=38) ctCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCGCGGGTAAACCCGCTC‑CTACCAGGGACC‑GTGCCGGCCTTGAATG‑CGGCGCTGCCCCTg > 4:288151/1‑131 (MQ=38) cATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCTGCCTCTTCGCGGGTAAACCCGCTC‑CTACCAGGGA‑ccgtgc < 7:140223/102‑6 (MQ=37) cATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCGCGGGTAAACCCGCTC‑CTACCAGGGA‑ccgtgc > 8:140223/1‑97 (MQ=37) cATTGAAACGTTGAGGAAGACCGGGGGCCTGGCAAGGTGCATATCATGCATTGCTCGGGCCGGCCGCTTCGCGGGTAAACCCGCTCATTACCAGGGACC‑GTGCCGGCCTTGAATG‑CGGCGCTGCCCCTGGTAGGAGCGGGTTTAcc < 1:315151/146‑1 (MQ=25) aTTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCGCGGGTAAACCCGCTC‑CTACCAGGGACC‑GTGCCGGCCTTGAATG‑CGGCGCTGCCCCTGGTAGGAGCgg < 8:65266/137‑1 (MQ=38) aTTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCGCGGGTAAACCCGCTC‑CTACCAGGGACC‑GTGCCGGCCTTGAATG‑CGGCGCTGCCCCTGGTAGGAGCgg > 7:65266/1‑137 (MQ=38) aTTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCGCGGGTAAACCCGCTC‑CTACCAGGGACC‑GTGCCGGCCTTGAATG‑CGGCGCTGCCCCTGGTAGGAGCGGGTTTACCCGCGa < 5:157743/149‑1 (MQ=38) aaaGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCGCGGGTAAACCCGCTC‑CTACCAGGGACC‑GTGCCGGCCTTGAATG‑CGGCGCTGCCCCTg < 7:81116/123‑1 (MQ=38) aaaGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCGCGGGTAAACCCGCTC‑CTACCAGGGACC‑GTGCCGGCCTTGAATG‑CGGCGCTGCCCCTg > 8:81116/1‑123 (MQ=38) ggTTCAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCGCGGGTAAACCCGCTC‑CTACCAGGGACC‑GTGCCGGCCTTGAATG‑CGGCGCTGCCCCTGGTAGGAGCGGGTTTACCCGCGAAGAggcc < 1:312896/149‑1 (MQ=38) ttGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCGCGGGTAAACCCGCTC‑CTACCAGGGACC‑GTGCCGGCCTTGAATG‑CGGCGCTGCCCCTGGTAGGAGCgg < 4:90552/128‑1 (MQ=38) ttGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCGCGGGTAAACCCGCTC‑CTACCAGGGACC‑GTGCCGGCCTTGAATG‑CGGCGCTGCCCCTGGTAGGAGCgg > 3:90552/1‑128 (MQ=38) gATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCGCGGGTAAACCCGCTC‑CTACCAGGGACC‑GTGCCGGCCTTGAATG‑AGGCGCTGCCCCTGGTAGGAGc < 6:143581/121‑1 (MQ=37) gATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCGCGGGTAAACCCGCTC‑CTACCAGGGACC‑GTGCCGGCCTTGAATG‑AGGCGCTGCCCCTGGTAGGAGc > 5:143581/1‑121 (MQ=37) cGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCGCGGGTAAACCCGCTC‑CTACCAGGGACC‑GTGCCGGCCTTGAATG‑CGGCGCTGCCCCTGGTa < 2:125018/110‑1 (MQ=37) cGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCGCGGGTAAACCCGCTC‑CTACCAGGGACC‑GTGCCGGCCTTGAATG‑CGGCGCTGCCCCTGGTa > 1:125018/1‑110 (MQ=37) ggTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCGCGGGTAAACCCGCTC‑CTACCAGGGACC‑GTGCCGGCCTTGAATG‑CGGCGCTGcccc > 4:140153/1‑103 (MQ=25) ggTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCGCGGGTAAACCCGCTC‑CTACCAGGGACC‑GTGCCGGCCTTGAATG‑CGGCGCTGcccc < 3:140153/103‑1 (MQ=25) tcgcgggTAAACCCGCTC‑CTACCAGGGACC‑GTGCCGGCCTTGAATG‑CGGCGCTGCCCCTGGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCCGCGCCAGCACATCAACCCGCAGATGCAGCCGCAATCTCCGCCACACTGATCTCAc < 2:112240/143‑1 (MQ=255) tcgcgggTAAACCCGCTC‑CTACCAGGGACC‑GTGCCGGCCTTGAATG‑CGGCGCTGCCCCTGGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCCGCGCCAGCACATCAACCCGCAGATGCAGCCGCAATCTCCGCCACACTGATCTCACg > 1:112240/6‑149 (MQ=255) gggTAAACCCGCTC‑CTACCAGGGACC‑GTGCCGGCCTTGAATG‑CGGCGCTGCCCCTGGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCCGCGCCAGCACATCAACCCGCAGATGCAGCCGCAATCTcc < 5:293647/127‑1 (MQ=255) gggTAAACCCGCTC‑CTACCAGGGACC‑GTGCCGGCCTTGAATG‑CGGCGCTGCCCCTGGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCCGCGCCAGCACATCAACCCGCAGATGCAGCCGCAATCTcc > 6:293647/2‑128 (MQ=255) ggTAAACCCGCTC‑CTACCAGGGG‑CAGCGCC‑GCATTCAAGGCCGGCACGGTCCCTGGTAGGAGCGGGTTTACCCGCGAAGCGGCCGGCCGCGCCAGCCCATCAACCCGCAGATGCAGCCGCACTCTCCGCCCCACTGATCTCACGCATGc > 3:97952/1‑149 (MQ=17) ccGCTC‑CTACCAGGGACC‑GTGCCGGCCTTGAATG‑CGGCGCTGCCCCTGGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCCGCGCCAGCACATCAACCCGCAGATGCAGCCGCAATCTCCGCCACACTGATCTCACGCATGCGGAACtt > 2:116381/1‑149 (MQ=255) | ACCACGGCACCGGCGTGAACAGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTCTGC‑GGTAAACCCGCTC‑CTACCAGGGA‑C‑GTGCCGGCCTTGAATG‑CGGCGCTGCCCCTGGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCCGCGCCAGCACATCAACCCGCAGATGCAGCCGCAATCTCCGCCACACTGATCTCACGCATGCGGAACTT > NC_002947/4585924‑4586188 |
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: ‑ |