Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | NC_002947 | 4,586,030 | 2 bp→TC | 100% | intergenic (+113/+101) | 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,030 | 0 | C | T | 96.6% | 83.8 / ‑4.7 | 29 | intergenic (+113/+102) | PP_4061/PP_4063 | hypothetical protein/long‑chain fatty acid‑‑CoA ligase |
Reads supporting (aligned to +/- strand): ref base C (1/0); new base T (12/16); total (13/16) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 4.48e-01 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 | |||||||||||
Rejected as polymorphism: E-value score below prediction cutoff. | |||||||||||
Rejected as polymorphism: Frequency below/above cutoff threshold. | |||||||||||
Rejected as polymorphism: Variant not supported by required number of reads on each strand. | |||||||||||
* | NC_002947 | 4,586,031 | 0 | T | C | 100.0% | 90.8 / NA | 29 | intergenic (+114/+101) | PP_4061/PP_4063 | hypothetical protein/long‑chain fatty acid‑‑CoA ligase |
Reads supporting (aligned to +/- strand): ref base T (0/0); new base C (13/16); total (13/16) | |||||||||||
Rejected as polymorphism: Frequency below/above cutoff threshold. | |||||||||||
Rejected as polymorphism: Variant not supported by required number of reads on each strand. |
AGGAGCCACCGGTGAAGGCTTAAGGGTAGACCACGGCACCGGCGTGAACAGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTC‑TCTGC‑GGTAAACCCGCTCC‑TACCAGGGA‑C‑GTGCCGGCC‑TTGAATG‑CGGCGCTGCCCCTGGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCCGCGCCA > NC_002947/4585895‑4586126 || aGGAGCCACCGGTGAAGGCTTAAGGGTAGACCACGGCACCGGCGTGAACAGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTC‑TTCGC‑gg < 4:172984/141‑1 (MQ=255) aGGAGCCACCGGTGAAGGCTTAAGGGTAGACCACGGCACCGGCGTGAACAGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTC‑TTCGC‑gg > 3:172984/1‑141 (MQ=255) ccACCGGTGAAGGCTTAAGGGTAGACCACGGCACCGGCGTGAACAGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTC‑TTCGCGGGTAAACCCGCTcc > 5:140589/1‑149 (MQ=255) aCCGGTGAAGGCTTAAGGGTAGACCACGGCACCGGCGTGAACAGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTCCATTGATCGGGCCGGCCAC‑TTCGCGGGTAAACCCGCTccaa > 5:280481/1‑147 (MQ=255) tGAAGGCTTAAGGGTAGACCACGGCACCGGCGTGAACAGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTC‑TTCGCGGGTAAACCCGCTCC‑TACCAgg < 5:156592/149‑1 (MQ=255) tggCGTAAGGGTAGACCACGGCACCGGCGTGAACAGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTC‑TTCGCGGGTAAACCCGCTCC‑TACCAGGGA‑c < 8:66013/148‑1 (MQ=39) ttAAGGGTAGACCACGGCACCGGCGTGAACAGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTC‑TTCGCGGGTAAACCCGCTCC‑TACCAGGGA‑ccgtg > 4:117978/1‑145 (MQ=39) ttAAGGGTAGACCACGGCACCGGCGTGAACAGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTC‑TTCGCGGGTAAACCCGCTCC‑TACCAGGGA‑ccgtg > 2:67947/1‑145 (MQ=39) gggTAGACCACGGCACCGGCGTGAACAGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTC‑TTCGCGGGTAAACCCGCTCC‑TACCAGGGACC‑GTGCCgg > 1:245817/1‑149 (MQ=39) aCCACGGCACCGGCGTGAACAGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTC‑TTCGCGGGTAAACCCGCTCC‑TACCAGGGACC‑GTGCCGGCC‑TTg > 5:152136/1‑148 (MQ=39) aCCAAGGCACCGGCGTGAACAGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTC‑TTCGCGGGTAAACCCGCTCC‑TACCAGGGACC‑GTGCCGGCC‑TTg > 7:237844/1‑148 (MQ=38) cGGCACCGGCGTGAACAGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTC‑TTCGCGGGTAAACCCGCTCC‑TACCAGGGACC‑GTGCCGGCC‑TTGAATg‑c > 1:121069/1‑149 (MQ=39) gCACCGGCGTGAACAGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTC‑TTCGCGGGTAAACCCGCTCC‑TACCAGGGACC‑GTGCCGGCC‑TTGAATg‑cgg < 2:121069/149‑1 (MQ=39) gCACCGGCGTGAACAGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTC‑TTCGCGGGTAAACCCGCTCC‑TACCAGGGACC‑GTGCCGGCC‑TTGAATg‑cgg < 8:237844/149‑1 (MQ=39) cgcTGCTCCATTGAAAGGTTGAGGATGTCCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCAGGGCTC‑TTCGCGGGTAAACCCGCTCC‑TACCAGGGACC‑GTGCCGGCC‑TTGAATG‑CGGCGCTGCCCCTGGTAGGAGCgg < 2:247765/146‑1 (MQ=255) cgcTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTC‑TTCGCGGGTAAACCCGCTCC‑TACCAGGGG‑CAGCGCC‑GCA‑TTCAAGGCCGGCACGGTCCCTGGTAGGAGc < 6:140589/144‑1 (MQ=14) cgcTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTC‑TTCGCGGGTAAACCCGCTCC‑TACCAGGGACC‑GTGCCGGCC‑TTGAATG‑CGGCGCTGCCCCTGGTAGGAGc < 2:245817/144‑1 (MQ=38) cgcTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTC‑TTCGCGGGTAAACCCGCTCC‑TACCAGGGACC‑GTGCCGGCC‑TTGAATG‑CGGCGCTGCCCCTGGTAGGAGCgg < 1:67947/146‑1 (MQ=38) cgcTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTC‑TTCGCGGGTAAACCCGCTCC‑TACCAGGGACC‑GTGCCGGCC‑TTGAATG‑CGGCGCTGCCCCTGGTAGGAGCg < 3:117978/145‑1 (MQ=38) gctgctCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTC‑TTCGCGGGTAAACCCGCTCC‑TACCAGGGACC‑GTGCCGGCC‑TTGAATG‑CGGCGCTGCCCCTGGTAGGAGCGGGttt < 6:152136/149‑1 (MQ=38) cATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTC‑TTCGCGGGTAAACCCGCTCC‑TACCAGGGACC‑GTGCCGGCC‑TTGAATG‑CGGCGCTGCCCCTGGTAg > 8:193665/1‑132 (MQ=38) cATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTC‑TTCGCGGGTAAACCCGCTCC‑TACCAGGGACC‑GTGCCGGCC‑TTGAATG‑CGGCGCTGCCCCTGGTAg < 7:193665/132‑1 (MQ=38) aacggTTGAGGATGAGCGGGTGGTTGGCATGGTGAAGATCTTGCATTGCTCGGGCGGGCCTC‑TTCGCGGGTAAACCCGCTCC‑TACCAGGGACC‑GTGCCGGCC‑TTGAATG‑CGGCGGTGCCGCTGGTAGGAGCGGGTTTACCCGCGAaga < 6:280481/145‑1 (MQ=25) aaaGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTC‑TTCGCGGGTAAACCCGCTCC‑TACCAGGGACC‑GTGCCGGCC‑TTGAATG‑CGGCGCTGCCCCTGGt > 1:239291/1‑125 (MQ=38) aaaGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTC‑TTCGCGGGTAAACCCGCTCC‑TACCAGGGACC‑GTGCCGGCC‑TTGAATG‑CGGCGCTGCCCCTGGt < 2:239291/125‑1 (MQ=38) gACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCCCC‑GGTAAAACCCCACCAAACAACGAA‑C‑CAGCCGCCCATTAAATA‑C‑CCGCAGCCCCGGGTAGGAGCGGGTTTTACCGCCAAGAGTCCGGCCCCCCCa > 6:150769/1‑149 (MQ=255) gTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTC‑TTCGCGGGTAAACCCGCTCC‑TACCAGGGACC‑GTGCCGGCC‑TTGAATG‑CGGCGCTGcccc > 7:280160/1‑102 (MQ=25) gTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTC‑TTCGCGGGTAAACCCGCTCC‑TACCAGGGACC‑GTGCCGGCC‑TTGAATG‑CGGCGCTGcccc < 8:280160/102‑1 (MQ=25) gTGAAGATCTTGCATTGATCGGGCCGGCCTC‑TTCGCGGGTAAACCCGCTCC‑TACCAGGGACC‑GTGCCGGCC‑TTGAATg‑cgg < 6:104310/81‑1 (MQ=21) || AGGAGCCACCGGTGAAGGCTTAAGGGTAGACCACGGCACCGGCGTGAACAGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTC‑TCTGC‑GGTAAACCCGCTCC‑TACCAGGGA‑C‑GTGCCGGCC‑TTGAATG‑CGGCGCTGCCCCTGGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCCGCGCCA > NC_002947/4585895‑4586126 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |