Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | NC_002947 | 5,478,830 | G→C | 17.4% | I33M (ATC→ATG) | PP_4816 ← | hypothetical protein |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_002947 | 5,478,830 | 0 | G | C | 17.4% | 54.0 / 6.6 | 23 | I33M (ATC→ATG) | PP_4816 | hypothetical protein |
Reads supporting (aligned to +/- strand): ref base G (8/11); new base C (2/2); total (10/13) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 1.00e+00 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.29e-01 |
TTGGCATATACCATCACCCCGACTGCCGGCGGCACCAGCAGCAGCGGGAGATAGCGCAGCAAGCTGCTTGCGGCCTCGTTCAGCGGCTTGCCGACCTCACCGCGGGCCATGAGGAAGAGCAGCAGCAACAACAGGCCGATGATCGGCCCGGGTAAAAACGGCACGAACAGGTGATTGATCGCCGTCCCCAGCAGCTGGAACAGTACCAGCCAGGTCAAACCACGCAACAGCATAAGCATCTCCCTCGGGCATGGCCGCCATTATAAGCACGCTCGG > NC_002947/5478690‑5478965 | ttGGCATATACCATCACCCCGACTGCCGGCGGCACCAGCAGCAGCGGGAGATAGCGCAGCAAGCTGCTTGCGGCCTCGTTCAGCGGCTTGCCGACCTCACCGCGGGCCATGAGGAAGAGCAGCAGCAACAACAGGCCGATGATCGGCCCg > 1:118476/1‑150 (MQ=255) gACTGCCGGCGGCACCAGCAGCAGCGGGAGATAGCGCAGCAAGCTGCTTGCGGCCTCGTTCAGCGGCTTGCCGACCTCACCGCGGGCCATGAGGAAGAGCAGCAGCAACAACAGGCCGATGATCGGCCCGGGTAAAAACGGCACGAACAg < 1:443630/150‑1 (MQ=255) ggcACCAGCAGCAGCGGGAGATAGCGCAGCAAGCTGCTTGCGGCCTCGTTCAGCGGCTTGCCGACCTCACCGCGGGCCATGAGGAAGAGCAGCAGCAACAACAGGCCGATGATCGGCCCGGGTAAAAACGGCACGAACAGGTGATTGATc < 1:279018/150‑1 (MQ=255) ggcACCAGCAGCAGCGGGAGATAGCGCAGCAAGCTGCTTGCGGCCTCGTTCAGCGGCTTGCCGACCTCACCGCGGGCCATGAGGAAGAGCAGCAGCAACAACAGGCCGATGATCGGCCCGGGTAAAAACGGCACGAACAGGTGATTGATc < 1:102152/150‑1 (MQ=255) aCCAGCAGCAGCGGGAGATAGCGCAGCAAGCTGCTTGCGGCCTCGTTCAGCGGCTTGCCGACCTCACCGCGGGCCATGAGGAAGAGCAGCAGCAACAACAGGCCGATCATCGGcctgtt > 1:602673/1‑115 (MQ=255) aCCAGCAGCAGCGGGAGATAGCGCAGCAAGCTGCTTGCGGCCTCGTTCAGCGGCTTGCCGACCTCACCGCGGGCCATGAGGAAGAGCAGCAGCAACAACAGGCCGATCATCGGcctgtt < 2:602673/119‑5 (MQ=255) cagcagcagCGGGAGATAGCGCAGCAAGCTGCTTGCGGCCTCGTTCAGCGGCTTGCCGACCTCACCGCGGGCCATGAGGAAGAGCAGCAGCAACAACAGGCCGATGATCGGCCCGGGTAAAAACGGCACGAACAg > 1:128684/1‑135 (MQ=255) cagcagcagCGGGAGATAGCGCAGCAAGCTGCTTGCGGCCTCGTTCAGCGGCTTGCCGACCTCACCGCGGGCCATGAGGAAGAGCAGCAGCAACAACAGGCCGATGATCGGCCCGGGTAAAAACGGCACGAACAg < 2:128684/135‑1 (MQ=255) ccTCGTTCAGCGGCTTGCCGACCTCACCGCGGGCCATGAGGAAGAGCAGCAGCAACAACAGGCCGATGATCGGCCCGGGTAAAAACGGCACGAACAGGTGATTGATCGCCGTCCCCAGCAGCTGGAACAGTACCAGCCAGGTCAAACCAc < 2:515254/150‑1 (MQ=255) aGCGGCTTGCCGACCTAACCGCGGGCGATGAGGAAGAGCGGCAGCAACAACAGGCCGATGATCGGCCCGGGTAAAAACGGCACGAACAGGTGATTGATCGCCGTCCCCAGCAGCTGGAACAGTACCAGCCAGGTCAAACCACGCAACAGc > 2:441858/1‑150 (MQ=255) gCTTGCCGACCTCACCGCGGGCCATGAGGAAGAGCAGCAGCAACAACAGGCCGATGATCGGcgc < 2:352431/64‑3 (MQ=255) gCTTGCCGACCTCACCGCGGGCCATGAGGAAGAGCAGCAGCAACAACAGGACGATGATCGGcgc > 1:352431/1‑62 (MQ=255) cGACCTCACCGCGGGCCATGAGGAAGAGCAGCAGCAACAACAGGCCGATGATCGGCCCGGGTAAAAACGGCACGAACAGGTGATTGATCGCCGTCCCCAGCAGCTGGAACAGTACCAGCCAGGTCAAACCACGCAACAGCATaa > 1:496671/1‑144 (MQ=255) cGACCTCACCGCGGGCCATGAGGAAGAGCAGCAGCAACAACAGGCCGATGATCGGCCCGGGTAAAAACGGCACGAACAGGTGATTGATCGCCGTCCCCAGCAGCTGGAACAGTACCAGCCAGGTCAAACCACGCAACAGCATaa < 2:496671/144‑1 (MQ=255) aCCGCGGGCCATGAGGAAGAGCAGCAGCAACAACAGGCCGATGATCGGCCCGGGTAAAAACGGCACGAACAGGTGATTGATCGCCGTCCCCAGCAGCTGGAACAGTACCAGCCAGGTCAAACCAc < 2:410314/125‑1 (MQ=255) aCCGCGGGCCATGAGGAAGAGCAGCAGCAACAACAGGCCGATGATCGGCCCGGGTAAAAACGGCACGAACAGGTGATTGATCGCCGTCCCCAGCAGCTGGAACAGTACCAGCCAGGTCAAACCAc > 1:410314/1‑125 (MQ=255) ccGCGGGCCATGAGGAAGAGCAGCAGCAACAACAGGCCGATGATCGGCCCGGGTAAAACCGGCACGAACAGGTGATTGATCGCCGTCCCCAGCAGCTGGAACAGTACCAGCCAGGTCAAACCACGCAACAGCATAAGCATCTCCCTCggg < 1:339106/150‑1 (MQ=255) gagCAGCAGCAACAACAGGCCGATGATCGGCCCGGGTAAAAACGGCACGAACAGGTGATTGATCGCCGTCCCCAGCAGCTGGAACAGTACCAGCCAGGTCAAACCACGCAACAGCATAAGCATCTCCCTCGGGCATGGCCGCCATTATaa > 1:49028/1‑150 (MQ=255) agcagcaACAACAGGCCGATGATCGGCCCGGGTAAAAACGGCACGAACAGGTGATTGATCGCCGTCCCCAGCAGCTGGAACAGTACCAGCCAGGTCAAACCACGCAACAGCATAAGCATCTCCCTCGGGCATGGCCGCCATTATAAGCAc < 2:118476/150‑1 (MQ=255) aacaacAGGCCGATGATCGGCCCGGGTAAAAACGGCACGAACAGGTGATTGATCGCCGTCCCCAGCAGCTGGAACAGTACCAGCCAGGTCAAACCACGCAACAGCATAAGCATCTCCCTCGGGCATGGCCGCCATTATAAGCACGCTCgg < 1:111549/150‑1 (MQ=255) aacaacAGGCCGATGATCGGCCCGGGTAAAAACGGCACGAACAGGTGATTGATCGCCGTCCCCAGCAGCTGGAACAGTACCAGCCAGGTCAAACCACGCAACAGCATAAGCATCTCCCTCGGGCATGGCCGCCATTATAAGCAAGCTCgg > 1:588932/1‑150 (MQ=255) ccgggCCGATCATCGGCCCGGGTAAAAACGGCACGAACAGGTGATTGATCGCCGTCCCCAGCAGCTGGAACAGTACCAGCCAGGTCAAACCACGCAACAGCATAAGCATCTCCCTCGGGCAt > 1:75007/4‑122 (MQ=255) ccgggCCGATCATCGGCCCGGGTAAAAACGGCACGAACAGGTGATTGATCGCCGTCCCCAGCAGCTGGAACAGTACCAGCCAGGTCAAACCACGCAACAGCATAAGCATCTCCCTCGGGCAt < 2:75007/119‑1 (MQ=255) | TTGGCATATACCATCACCCCGACTGCCGGCGGCACCAGCAGCAGCGGGAGATAGCGCAGCAAGCTGCTTGCGGCCTCGTTCAGCGGCTTGCCGACCTCACCGCGGGCCATGAGGAAGAGCAGCAGCAACAACAGGCCGATGATCGGCCCGGGTAAAAACGGCACGAACAGGTGATTGATCGCCGTCCCCAGCAGCTGGAACAGTACCAGCCAGGTCAAACCACGCAACAGCATAAGCATCTCCCTCGGGCATGGCCGCCATTATAAGCACGCTCGG > NC_002947/5478690‑5478965 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |