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 | 93.5% | 94.1 / 0.7 | 31 | 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 (17/12); total (19/12) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 5.10e-01 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.82e-01 |
TAGACCACGGCACCGGCGTGAACAGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTCTGC‑GGTAAACCCGCTCCTACCA‑GGG‑A‑CGTGCCGGCCTTGAATGCGGCGCTGCCCCTGGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCCGCGCCAGCACATCAACCCGCAGATGCAGCCGCAATCTCCGCCACACTGATCTCACGCATGCGGAACTTCTGCACCTTG > NC_002947/4585921‑4586198 | tAGACCACGGCACCGGCGTGAACAGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACCA‑GGG‑ACCGTGCCGGCCt > 5:37486/1‑149 (MQ=39) tAGACCACGGCACCGGCGTGAACAGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACCA‑GGG‑ACCGTGCCGGCCt < 6:197528/149‑1 (MQ=39) gACCACGGCACCGGCGTGAACAGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACCA‑GGG‑ACCGTGCCGGCCtt > 8:205315/1‑148 (MQ=39) cACGGCACCGGCGTGAACAGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACCA‑GGG‑ACCGTGCCGGCCTTGAAt < 7:89227/149‑1 (MQ=39) cACGGCACCGGCGTGAACAGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTAGCGGGTAAAACCCCTACTAACA‑TGG‑ACCCTGCCGGCCTTGAAt > 5:196996/1‑149 (MQ=255) ggCACCGGCGTGAACAGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACCA‑GGG‑ACCGTGCCGGCCTTGAATgcg < 7:350124/149‑1 (MQ=39) ggCACCGGCGTGAACAGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACCA‑GGG‑ACCGTGCCGGCCTTGAATgcg < 6:37486/149‑1 (MQ=39) aGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCGCGGGTTAAACCCCTCCAACCA‑GGGCC‑CGTGCCGGCCGTGAATGCGGCGCGGCCCCAggga > 1:77217/1‑147 (MQ=255) aGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACCA‑GGG‑ACCGTGCCGGCCTTGAATGCGGCGCTGCCCCTGGTa < 6:247868/149‑1 (MQ=39) aTCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCGCGGGTTAAACCGCGGCTACCAGGGG‑C‑CGGGGCGGGCATGAATGCGGCGCGGGCCCgggggggggc > 5:115311/1‑142 (MQ=255) cgcTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACCA‑GGG‑ACCGTGCCGGCCTTGAATGCGGCGCTGCCCCTGGTAGGAGCgg < 7:205315/146‑1 (MQ=38) ctCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACCA‑GGG‑ACCGTGCCGGCCTTGAATGCGGCGCTGCCCCTGGTAGGAGCGGGTTTAccc < 3:402952/149‑1 (MQ=38) cATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACCA‑GGG‑ACCGTGCCGGCCTTGAATGCGGCGCTGCCCCTGGTAGGAGCgg > 5:25073/1‑138 (MQ=38) cATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACCA‑GGG‑ACCGTGCCGGCCTTGAATGCGGCGCTGCCCCTGGTAGGAGCgg < 6:25073/138‑1 (MQ=38) gAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACCA‑GGG‑ACCGTGCCGGCCTTGAATGCGGCGCTGCCCCTGGTAGGAGCGGGTTTACCCGCGAag < 6:189103/148‑1 (MQ=38) gAAGATCTTGCATTGATCGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACCA‑GGG‑ACCGTGCCGGCCGTGAATGCGGCGCTGCCCCGGGTAGGAGCGGGGTTTCCCGCGAAGCGGCCGGCCGCGCCCGGACCTCCACCCGCAGCTg > 5:403880/1‑148 (MQ=25) gATCTTGCATTGATCGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACCA‑GGG‑ACCGTGCCGGCCTTGAATGCGGCGCTGCCCCTGGTAGGAGCGGGt < 1:359472/99‑1 (MQ=25) gATCTTGCATTGATCGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACCA‑GGG‑ACCGTGCCGGCCTTGAATGCGGCGCTGCCCCTGGTAGGAGCGGGt > 2:359472/1‑99 (MQ=25) ttcgcgggTAAACCCGCTCCTACCA‑GGG‑ACCGTGCCCGCCTTGAATGCGGCGCGGCCCCGGGTAGGAGCCGGTTTACCCGCGAACAGGCCGGCCCCGCCAGCACAACAACCCCCAGATGCCGCCGCAAACTCCGCCCCACGGATCTCAc > 2:359738/7‑149 (MQ=25) tcgcgggTAAACCCGCACCTACCA‑GGG‑ACCGTGCCGGCCTTGAAAGCGGCGCTGACACTGGTAGGAGCGGGTTTACCCGCGAAGAGGCCCGCCGCGCCAGCACAACAACCCGCAGAAGCAGACGCAATATCCGCCACACTGATCTAACg > 8:159760/6‑149 (MQ=37) gTAAACCCGCTCCTACCA‑GGG‑ACCGTGCCGGCCTTGAATGCGGCGCTGCCCCTGGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCCGCGCCAGCACATCAACCCGCAGATGCAGCCGCAATCTCCGCCACACTGATCTCACGCATGCg > 6:141717/1‑149 (MQ=255) tCCTACCA‑GGG‑ACCGTGCCGGCCTTGAATGCGGCGCTGcccc > 7:229977/1‑42 (MQ=255) tCCTACCA‑GGG‑ACCGTGCCGGCCTTGAATGCGGCGCTGcccc < 8:229977/42‑1 (MQ=255) tACCA‑GGG‑ACCGTGCCGGCCTTGAATGCGGCGCTGCCCCTGGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCCGCGCCAGCACATCAACCCGCAGATGCAGCCGCAATCTCCGCCACACTGATCTCACGCATGCGGAACTTCTGCAcc > 7:276709/1‑149 (MQ=255) tACCA‑GGG‑ACCGTGCCGGCCTTGAATGCGGCGCTGCCCCTGGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCCGCGCCAGCACATCAACCCGCAGATGCAGCCGCAATCTCCGCCACACTGATCTCACGCATGCGGAACTTCTGCAcc > 2:103364/1‑149 (MQ=255) tACCA‑GGG‑ACCGTGCCGGCCTTGAATGCGGCGCTGCCCCTGGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCCGCGCCAGCACATCAACCCGCAGATGCAGCCGCAATCTCCGCCACACTGATCTCACGCATGCGGAACTTCTGCAcc > 1:344353/1‑149 (MQ=255) aCCA‑GGG‑ACCGTGCCGGCCTTGAATGCGGCGCTGCCCCTGGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCCGCGCCAGCACATCAACCCGCAGATGCAGCCGCAATCTCCGCCACACTGATCTCACGCATGCGGAACTTCTGCACCt > 4:32302/1‑149 (MQ=255) aCCA‑GGG‑ACCGTGCCGGCCTTGAATGCGGCGCTGCCCCTGGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCCGCGCCAGCACATCAACCCGCAGATGCAGCCGCAATCTCCGCCACACTGATCTCACGCATGCGGAACTTCTGCACCt > 3:296630/1‑149 (MQ=255) cA‑GGG‑ACCGTGCCGGCCTTGAATGCGGCGCTGCCCCTGGTAGGAGCGGGTTTATCCGCGAAGAGGCCGGCCGCGCCAGCACATCAACCCGCAGATGCAGCCGCAATCTCCGCCACACTGATCTCACGCATGCGGAACTTCTGCACCTTg > 2:113465/1‑149 (MQ=255) cA‑GGG‑ACCGTGCCGGCCTTGAATGCGGCGCTGCCCCTGGTAGGAGCGGGTTTACCCGCGaaagg > 6:123327/1‑61 (MQ=25) cA‑GGG‑ACCGTGCCGGCCTTGAATGCGGCGCTGCCCCTGGTAGGAGCGGGTTTACCCGCGaaagg < 5:123327/64‑4 (MQ=25) | TAGACCACGGCACCGGCGTGAACAGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTCTGC‑GGTAAACCCGCTCCTACCA‑GGG‑A‑CGTGCCGGCCTTGAATGCGGCGCTGCCCCTGGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCCGCGCCAGCACATCAACCCGCAGATGCAGCCGCAATCTCCGCCACACTGATCTCACGCATGCGGAACTTCTGCACCTTG > NC_002947/4585921‑4586198 |
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: ‑ |