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 | 91.2% | 100.9 / 4.2 | 34 | intergenic (+139/+76) | PP_4061/PP_4063 | hypothetical protein/long‑chain fatty acid‑‑CoA ligase |
Reads supporting (aligned to +/- strand): ref base . (3/0); new base C (13/18); total (16/18) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 9.36e-02 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.13e-01 |
AAGGCTTAAGGGTAGACCACGGCACCGGCGTGAACAGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTCTGC‑GGTAAACCCGCTCCTACCA‑GGGA‑C‑GTGCCGGCCTTGAATG‑CGGCGCTGCCCCTGGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCCGCGCCAGCACATCAACCCGCAGATGCAGCCGCAATCTCCGCCACACTGATCTCACGCATGCGGAACTTCTGCACC > NC_002947/4585909‑4586195 | aaGGCTTAAGGGTAGACCACGGCACCGGCGTGAACAGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACCA‑GGGa < 3:142084/149‑1 (MQ=255) gCTTAAGGGTAGACCACGGCACCGGCGTGAACAGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACCA‑GGGA‑ccg > 8:87284/1‑147 (MQ=255) cTTAAGGGTAGACCACGGCACCGGCGTGAACAGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACCA‑GGGA‑ccgt > 3:152869/1‑146 (MQ=39) aGGGTAGACCACGGCACCGGCGTGAACAGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACCA‑GGGACC‑GTGCCg < 7:151337/149‑1 (MQ=39) gggTAGACCACGGCACCGGCGTGAACAGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCGCGGGTAAAACCGCTACTACCAGGGGC‑C‑GTGCCgg > 1:13786/1‑149 (MQ=38) gTAGACCACGGCACCGGCGTGAACAGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACCA‑GGGACC‑GTGCCGGcc < 4:220308/149‑1 (MQ=39) gACCACGGCACCGGCGTGAACAGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACCA‑GGGACC‑GTGCCGGCCtt < 7:201017/148‑1 (MQ=39) cGGCACCGGCGTGAACAGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCTGGCCTCTTCGCGGGTAAACCCGCTCCTACCA‑GGGACC‑GTGCCGGCCTTGAATg‑c < 6:25215/149‑1 (MQ=38) cGGCACCGGCGTGAACAGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACCA‑GGGACC‑GTGCCGGCCTTGAATg‑c < 4:152869/149‑1 (MQ=39) gAACAGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGTCTCTTCGCGGGTAAACCCGCTCCTACCA‑GGGACC‑GTGCCGGCCTTGAATG‑CGGCGCTGCCCCt > 2:129599/1‑149 (MQ=39) cAGCTGATCTGCGCTGCTCCATTGAAAGGGTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACCA‑GGGACC‑GTGCCGGCCTTGAATG‑CGGCGCTGCCCCTGGt < 8:365261/149‑1 (MQ=38) aGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACCA‑GGGACC‑GTGCCg < 3:346513/122‑1 (MQ=38) aGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACCA‑GGGACC‑GTGCCg > 4:346513/1‑122 (MQ=38) aTCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACCA‑GGGACC‑GTGCCGGCCTTGAATG‑CGGCGCTGCCCCTGGTAGGa < 5:147672/147‑1 (MQ=39) cgcTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACCA‑GGGACC‑GTGCCGGCCTTGAATG‑CGGCGCTGCCCCTGGTAGGTGCggg < 1:294304/147‑1 (MQ=38) cgcTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACCA‑GGGACC‑GTGCCGGCCTTGAATG‑CGGCGCTGCCCCTGGTAGGAGCggg < 7:87284/147‑1 (MQ=38) cgcTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACCA‑GGGACC‑GTGCCGGCCTTGAAAG‑CGGCGCTGCCCCTGGTAGGAGCggg < 7:67951/147‑1 (MQ=38) ctCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGGCTCTTCGCGGGTAAACCCGCTCCTACCA‑GGGACC‑GTGCCGGCCTTGAATG‑CGGCGCTGCCCCTGGTAGGAGCGGGTTTAccc < 1:129599/149‑1 (MQ=38) tCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACCA‑GGGACC‑GTGCCGGCCTTGAATG‑CGGCGCTGcccc < 6:17798/128‑1 (MQ=38) tCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACCA‑GGGACC‑GTGCCGGCCTTGAATG‑CGGCGCTGcccc > 5:17798/1‑128 (MQ=38) tCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACCA‑GGGACC‑GTGCCGGCCTTGAATG‑CGGCGCTGCCCCTGGTAGGAGCgg < 6:75495/140‑1 (MQ=38) tCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACCA‑GGGACC‑GTGCCGGCCTTGAATG‑CGGCGCTGCCCCTGGTAGGAGCgg > 5:75495/1‑140 (MQ=38) aTTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTAGCGGGTAAACCCCCAACTACCA‑GGGACC‑CAGCCCGCCGTGAATG‑CCGCGCGGCCCCAGGTAGGAGCGGGGTTACCcgcg > 6:70518/1‑148 (MQ=255) ttGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCGCGGGTAAAACCCCACCTAACA‑GGGACC‑GTGCCGGCCTTGAAAG‑CGGCGCTGCCCCTGGTAGGAGCGGGTTTTACCGCGa > 2:71262/1‑148 (MQ=255) gATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACCA‑GGGACC‑GTGCCGGCCTTGAATG‑CGGCGCTGCCCCTGGTAGGAGCggg > 5:386041/1‑124 (MQ=37) gATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACCA‑GGGACC‑GTGCCGGCCTTGAATG‑CGGCGCTGCCCCTGGTAGGAGCggg < 6:386041/124‑1 (MQ=37) tGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACCA‑GGGACC‑GTGCCGGCCTTGAATG‑CGGCGCTGCCCCTgg < 6:378323/112‑1 (MQ=37) tGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACCA‑GGGACC‑GTGCCGGCCTTGAATG‑CGGCGCTGCCCCTgg > 5:378323/1‑112 (MQ=37) ttGGCACGGTGAAGAGCTTGCATTGATCGGGCCGGCCGCTTCGCGGGTAAACCCGCTCCTACCA‑GGGACC‑GTGCCGGCCTTGAATG‑CGGCGCTGCCCCTGGTAGGAGCGGGTTTACCCGCGAAGAggc < 5:387788/128‑1 (MQ=25) ggCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACCA‑GGGACC‑GTGCCGGCCTTGAATG‑CGGCGCTGc < 2:214825/93‑1 (MQ=25) ggCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACCA‑GGGACC‑GTGCCGGCCTTGAATG‑CGGCGCTGc > 1:214825/1‑93 (MQ=25) tAAACCCGCTCCTACCA‑GGGG‑CAGCGCC‑GCATTCAAGGCCGGCACGGTCCCTGGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCCGCGCCAGCACATCAACCCGCAGATGCAGCCGCAATCTCCGCCACACTGATCTCACGCATGCg > 7:102289/1‑148 (MQ=18) aCCCGCTCCTACCA‑GGGACC‑GTGCCGGCCTTGAATG‑CGGCGCTGCCCCTGGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCCGCGCCAGCACATCAACCCGCAGATGCAGCCGCAATCTCCGCCACACTGATCTCACGCATGCGGAAc > 7:52049/1‑149 (MQ=255) tCCTACCA‑GGGG‑CAGTGCC‑GCATTCAAGGCCGGCACGGTCCCTGGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCCGCGCCAGCACATCAACCCGCAGATGCAGCCGCAATCTCCGCCACACTGATCTCACGCATGCGGAACTTCTg > 6:159845/1‑148 (MQ=21) tCCTACCA‑GGGACC‑GTGCCGGCCTTGAATG‑CGGCGCTGCCACTGGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGACGCGCCAGCACATCAACCCGCAGATGCCGACGCAATCTCCGCCACCATGCTCTCACGCATCCCGCACTTCtt > 6:337227/1‑147 (MQ=255) tACCA‑GGGACC‑GTGCCGGCCTTGAATG‑CGGCGCTGCCCCTGGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCCGCGCCAGCACATCAACCCGCAGATGCAGCCGCAATCTCCGCCACACTGATCTCACGCATGCGGAACTTCTGCAcc > 1:151870/1‑149 (MQ=255) tACCA‑GGGACC‑GTGCCGGCCTTGAATG‑CGGCGCTGCCCCTGGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCCGCGCCAGCACATCAACCCGCAGATGCAGCCGCAATCTCCGCCACACTGATCTCACGCATGCGGAACTTCTGCAcc > 8:193280/1‑149 (MQ=255) | AAGGCTTAAGGGTAGACCACGGCACCGGCGTGAACAGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTCTGC‑GGTAAACCCGCTCCTACCA‑GGGA‑C‑GTGCCGGCCTTGAATG‑CGGCGCTGCCCCTGGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCCGCGCCAGCACATCAACCCGCAGATGCAGCCGCAATCTCCGCCACACTGATCTCACGCATGCGGAACTTCTGCACC > NC_002947/4585909‑4586195 |
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: ‑ |