Predicted mutation | ||||||
---|---|---|---|---|---|---|
evidence | seq id | position | mutation | annotation | gene | description |
RA | NC_002947 | 4,586,035 | +G | intergenic (+118/+97) | 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,033 | 1 | . | G | 95.7% | 68.2 / ‑2.9 | 23 | intergenic (+116/+99) | PP_4061/PP_4063 | hypothetical protein/long‑chain fatty acid‑‑CoA ligase |
Reads supporting (aligned to +/- strand): ref base . (0/1); new base G (7/15); total (7/16) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 1.00e+00 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.26e-01 |
GCCGAGCCAGGAGCCACCGGTGAAGGCTTAAGGGTAGACCACGGCACCGGCGTGAACAGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTCTG‑C‑GGTAAACCCGCTCCTA‑CCAGGG‑A‑CGTGCCGGCCTTGAATGCGGCGCTGCCCCTGGTAGGAGCGGGTTTACCCGCGAAGA > NC_002947/4585887‑4586112 | gccgAGCCAGGAGCCACCGGTGAAGGCTTAAGGGTAGACCACGGCACCGGCGTGAACAGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCG‑C‑gg < 1:276045/149‑1 (MQ=255) gccgAGCCAGGAGCCACCGGTGAAGGCTTAAGGGTAGACCACGGCACCGGCGTGAACAGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCG‑C‑gg < 5:21171/149‑1 (MQ=255) aGGAGCCACCGGTGAAGGCTTAAGGGTAGACCACGGCACCGGCGTGAACAGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCG‑CGGGTAAAccc > 8:290130/1‑149 (MQ=255) aGCCACCGGTGAAGGCTTAAGGGTAGACCACGGCACCGGCGTGAACAGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCG‑CGGGTAAACCCGCt < 2:68050/149‑1 (MQ=255) gCCACCGGTGAAGGCTTAAGGGTAGACCACGGCACCGGCGTGAACAGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCG‑CGGGTAAACCCGCTc < 1:59451/149‑1 (MQ=255) aCCGGTGAAGGCTTAAGGGTAGACCACGGCACCGGCGTGAACAGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCG‑CGGGTAAACCCGCTCCTa > 3:268703/1‑149 (MQ=255) aaGGCTTAAGGGTAGACCACGGCACCGGCGTGAACAGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCG‑CGGGTAAACCCGCTCCTA‑CCAGGG‑a < 2:314358/149‑1 (MQ=255) ggTAGACCACGTCACCGGCGTGAACAGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCG‑CGGGTAAACCCGCTCCTA‑CCAGGG‑ACCGTGCCGGc < 5:215398/149‑1 (MQ=38) ggCGTGAACAGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCG‑CGGGTAAACCCGCTCCTA‑CCAGGG‑ACCGTGCCGGCCTTGAATGCGGCGCTg < 4:410546/149‑1 (MQ=39) aaCAGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCG‑CGGGTAAACCCGCTCCTA‑CCAGGG‑ACCGTGCCGGCCTTGAATGCGGCGCTGCCCCt < 3:348912/148‑1 (MQ=39) aCAGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCG‑CGGGTAAACCCGCTCCTA‑CCAGGG‑ACCGTGCCGGCCTTGAATGCGGCGCTGCCCCTgg < 7:290130/149‑1 (MQ=39) aGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCG‑CGGGTAAACCCGCTCCTA‑CCAGGG‑ACCGTGCCGGCCTTGAATGCGGCGCTGCCCCTGGTa < 6:320928/149‑1 (MQ=39) cTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCG‑CGGGTAAACCCGc < 5:302807/100‑1 (MQ=39) cTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCG‑CGGGTAAACCCGc > 6:302807/1‑100 (MQ=39) aTCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTAG‑CGGGTAAAACCGCGCCGA‑CCAGGGAA‑CCTTGCGGCGGTGAAATCGGCGCTGCCCCAgggtggacg > 7:290354/1‑142 (MQ=255) cgcTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCG‑CGGGTAAACCCGCTCCTA‑CCAGGG‑ACCGTGCCGGCCTTGAATGCGGCGCTGCCCCTGGTAGGAGc < 4:147359/144‑1 (MQ=38) cgcTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCG‑CGGGTAAACCCGCTCCTA‑CCAGGG‑ACCGTGCCGGCCTTGAATGCGGCGCTGCCCCTGGTAGGAGc < 1:375179/144‑1 (MQ=38) cgcTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCG‑CGGGTAAACCCGCTCCTA‑CCAGGG‑ACCGTGCCGGCCTTGAATGCGGCGCTGCCCCTGGTAGGAGCGGGt < 6:138731/148‑1 (MQ=38) cgcTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCGCTTCG‑CGGGTAAACCCGCTCCTA‑CCAGGG‑ACCGTGCCGGCCTTGAATGCGGCGCTGCCCCTGGTAGGAGCggg < 1:391655/147‑1 (MQ=38) tgctCCTTTGATAGGTTGAGTATGAGAGGGTTCTTTGCTTGGTGAAGATCTTGCATTTTTTGGTCCTGGATCTTGGTG‑GGTTAATCTGCTCCTA‑CGAGTG‑ACCGTGGCGGTATTGTTTGAGGCGCTGCCCCTGGTAGGAGCGGGtttt < 4:108369/148‑2 (MQ=255) ctCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCG‑CGGGTAAACCCGATACTA‑CCAGGG‑ACCGTGCCGGCCTTGAATGCGGCTCTGCCCCTGGGAGGAGCCGGTTTAca > 3:108369/1‑147 (MQ=255) ggATGACCCGGTGCGTGGCAAGGGCAAGATCTTGCATTGATCGGGCCGGCCTCTTCG‑CGGGTAAACCCGCTCCTA‑CCCGGG‑ACCGGGCCGGCCTTGAATGCGGCGCTGCCCCTGGTAGGAGCGGGTTTACCCGCGAaga < 5:352953/139‑1 (MQ=25) ggATGACCAGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCG‑CGGGTAAAACCGCGCCGACCCGGGG‑C‑CGGGGCGGGCGTGGATGCGGGGCGGCGCCGGGTAGGGGGGGGTTTAACCGCGAagg > 6:352953/1‑138 (MQ=255) ccGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCG‑CGGGTAAACCCGCTCCTA‑CCAGGG‑ACCGTGCCGGCCTTGAATGCGGCGCt < 7:139012/101‑1 (MQ=25) ccGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCG‑CGGGTAAACCCGCTCCTA‑CCAGGG‑ACCGTGCCGGCCTTGAATGCGGCGCt > 8:139012/1‑101 (MQ=25) | GCCGAGCCAGGAGCCACCGGTGAAGGCTTAAGGGTAGACCACGGCACCGGCGTGAACAGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTCTG‑C‑GGTAAACCCGCTCCTA‑CCAGGG‑A‑CGTGCCGGCCTTGAATGCGGCGCTGCCCCTGGTAGGAGCGGGTTTACCCGCGAAGA > NC_002947/4585887‑4586112 |
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: ‑ |