Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | NC_002947 | 85,731 | C→A | 8.7% | M329I (ATG→ATT) | PP_0075 ← | choline sulfate transporter |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_002947 | 85,731 | 0 | C | A | 8.7% | 136.9 / 3.7 | 46 | M329I (ATG→ATT) | PP_0075 | choline sulfate transporter |
Reads supporting (aligned to +/- strand): ref base C (21/21); new base A (2/2); total (23/23) | |||||||||||
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.15e-01 |
GAACAATGCACGGATCCCTCGGTGGTCCACCAACCCCCAGCAGATCAGCAGAATGCTGCCGGCCATGGCCGGTATCGGCAGGTGCGCGATCAGGCCGGCGCCGGTCAC‑GGCGAACAACGCCACCCACAGCGCTGAAAACACCCCGGCCATCGGTGAACGGGCGCCCGCCTCGTAGCTCAGGCCCGAGCGGGTGAAGGAACCGGCAGACAGGTAACCAGAGAAAAATGCACCGGCAATATTCGACAGGCCCTGTGCGCGTATTTCCTGGTCCGGGTCGATCAGTTGTTCTGAA > NC_002947/85584‑85875 | gAACAATGCACGGATCCCTCGGTGGTCCACCACCGCTCAGCAGATCAGCAGAATGCTGCCGGCCATGGCCGGTATCGGCAGGTGCGCGATCAGGCCGGCGCCGGTCAC‑GGCGAACAACGCCACCCACAGCGCTGAAAACACCCCGGCCAt > 1:640953/1‑150 (MQ=255) aaTGCACGGATCCCTCGGTGGTCCACCAACCCCCAGCAGATCAGCAGAATGCTGCCGGCCATGGCCGGTATCGGCAGGTGCGCGATCAGGCCGGCGCCGGTCAC‑GGCGAACAACGCCACCCACAGCGCTGAAAACACCCCGGCCATCGGt > 2:136282/1‑150 (MQ=255) tGCACGGATCCCTCGGTGGTCCACCAACCCCCAGCAGATCAGCAGAATGCTGCCGGCCATGGCCGGTATCGGCAGGTGCGCGATCAGGCCGGCGCCGGTCAC‑GGCGAACAACGCCACCCACAGCGCTGAAAACACCCCGGCCATCGGTGa < 1:269415/150‑1 (MQ=255) ggATCCCTCGGTGGTCCACCAACCCCCAGCAGATCAGCAGAATGCTGCCGGCCATGGCCGGTATCGGCAGGTGCGCGATCAGGCCGGCGCCGGTCAC‑GGCGAACAACGCCACCCACAGCGCTGAAAACACCCCGGCCATCGGTGAACggg > 1:360940/1‑150 (MQ=255) cTCGGTGGTCCACCAACCCCCAGCAGATCAGCAGAATGCTGCCGGCCATGGCCGGTATCGGCAGGTGCGCGATCAGGCCGGCGCCGGTCAC‑GGCGAACAACGCCACCCACAGCGCTGAAAACACCCCGGCCATCGGTGAACGGGcgcccg > 1:871523/1‑150 (MQ=255) accaACCCCCAGCAGATCAGCAGAATGCTGCCGGCCATGGCCGGTATCGGCAGGTGCGCGATCAGGCCGGCGCCGG‑CACAGGCGAACCCCGCGTCCTTCATCGCTGAAAACACCCCGGCCATCGGTGAACGGGCGCCCGCCTCGTAGCTc < 1:136282/150‑1 (MQ=255) aaCCCCCAGCAGATCAGCAGAATGCTGCCGGCCATGGCCGGTATCGGCAGGTGCGCGATCAGGCCGGCGCCGGTCAC‑GGCGAACAACGCCACCCACAGCGCTGAAAACACCCCGGCCATCGGTGAACGGGCGCCCGCCTCGTAGCTCAgg < 1:113959/150‑1 (MQ=255) aaCCCCCAGCAGATCAGCAGAATGCTGCCGGCCATGGCCGGTATCGGCAGGTGCGCGATCAGGCCGGCGCCGGTCAC‑GGCGAACAACGCCACCCACAGCGCTGAAAACACCCCGGCCATCGGTGAACGGGCGCCCGCCTCGTAGCTCAgg > 2:113959/1‑150 (MQ=255) ccAGCAGATCAGCAGAATGCTGCCGGCCATGGCCGGTATCGGCAGGTGCGCGATCAGGCCGGCGCCGGTCAC‑GGCGAACAACGCCACCCACAGAGCTGAAAACACCCCGGCCATCGGTGAACGGGCGCCCGCCTCGTAGCTCAGGCCCGa > 1:906683/1‑150 (MQ=255) ccAGCAGATCAGCAGAATGCTGCCGGCCATGGAATGTTTCGGCAGGTGCGCGATCAGGCCGGCGCCGGTCAC‑GGCGAACAACGCCACCCACAGCGCTGAAAACACCCCGGCCATCGGTGAACGGGCGCCCGCCTCGTAGCTCAGGCCCGa < 2:640506/150‑1 (MQ=255) gcggATCAGCAGACTGCTTGCAGCCATGGCCGGTAGCGGCGGGTGCGCGATCAGGCCGGCGCCGGTCAC‑GGCGAACAACGCCACCCACAGCGCTGAAAACACCCCGGCCATCGGTGAACGGGCGCCCGCCTCGTAGCTCAGGCCCGAGCg < 2:360940/147‑1 (MQ=255) cagATCAGCAGAATGCTGCCGGCCATGGCCGGTATCGGCAGGTGCGCGATCAGGCCGGCGCCGGTCAC‑GGCGAACAACGCCACCCACAGCGCTGAAAACACCCCGGCCATCGGTGAACGGGCGCCCGCCTCGTAGCTCAGGCCCGAGCgg < 2:871523/150‑1 (MQ=255) gATCAGCAGAATGCTGCCGGCCATGGCCGGTATCGGCAGGTGCGCGATCAGGCCGGCGCCGGTCAC‑GGCGAACAACGCCACCCACAGCGCTGAAAACACCCCGGCCATCGGTGAACGGGCGCCCGCCTCGTAGCTCAGGCCCGAGCGGGt < 2:35886/150‑1 (MQ=255) gATCAGCAGAATGCTGCCGGCCATGGCCGGTATCGGCAGGTGCGCGATCAGGCCGGCGCCGGTCAC‑GGCGAACAACGCCACCCACAGCGCTGAAAACACCCCGGCCATCGGTGAACGGGCGCCCGCCTCGTAGCTCAGGCCCGAGCGGGt > 2:683893/1‑150 (MQ=255) aTCAGCAGAATGCTGCCGGCCATGGCCGGTATCGGCAGGTGCGCGATCAGGCCGGCGCCGGTCAC‑GGCGAACAACGCCACCCACAGCGCTGAAAACACCCCGGCCATCGGTGAACGGGCGCCCGCCTCGTAGCTCAGGCCCGAGCGGGTg < 2:951689/150‑1 (MQ=255) agAATGCTGCCGGCCATGGCCGGTATCGGCAGGTGCGCGATCAGGCCGGCGCCGGTCAC‑GGCGAACAACGCCACCCACAGCGCTGAAAACACCCCGGCCATCGGTGAACGGGCGCCAGCCTCGTAGCTCAGGCCCGAGCGGGTGAAGGaa > 2:17382/1‑150 (MQ=255) gAATGCTGCCGGCCATGGCCGGTATCGGCAGGTGCGCGATCAGGCCGGCGCCGGTCAC‑GGCGAACAACGCCACCCACAGCGCTGAAAACACCCCGGCCATCGGTGAACGGGCGCCCGCCTCGTAGCTCAGGCCCGAGCGGGTGAAGGAAc < 2:766044/150‑1 (MQ=255) ccGGCCATGGCCGGTATCGGCAGGTGCGCGATCAGGCCGGCGCCGGTCAC‑GGCGAACAACGCCACCCACAGCGCTGAAAACACCCCGGCCATCGGTGAACGGGCGCCCGCCTCGTAGCTCAGGCCCGAGCGGGTGAAGGAACCGGCAGAc < 1:859487/150‑1 (MQ=255) cATGGCCGGTATCGGCAGGTGCGCGATCAGGCCGGCGCCGGTCAC‑GGCGAACAACGCCACCCACAGCGCTGAAAACACCCCGGCCATCGGTGAACGGGCGCCCGCCTCGTAGCTCAGGCCCGAGCGGGTGAAGGAACCGGCAGACAGGTa > 2:968962/1‑150 (MQ=255) tGGCCGGTATCGGCAGGTGCGCGATCAGGCCGGCGCCGGTCAC‑GGCGAACAACGCCACCCACAGCGCTGAAAACACCCCGGCCATCGGTGAACGGGCGCCCGCCTCGTAGCTCAGGCCCGAGCGGGTGAAGGAACCGGCAGCTAGGTAAc > 2:748505/1‑150 (MQ=255) ccGGTATCGGCAGGTGCGCGATCAGGCCGGCGCCGGTCAC‑GGCGAACAACGCCACCCACAGCGCTGAAAACACCCCGGCCATCGGTGAACGGGCGCCCGCCTCGTAGCTCAGGCCCGAGCGGGTGAAGGAACCGGCAGACAGGTAACCag < 1:328125/150‑1 (MQ=255) aTCGGCAGGTGCGCGATCAGGCCGGCGCCGGTCAC‑GGCGAACAACGCCACCCACAGCGCTGAAAACACCCCGGCCATCGGTGAACGGGCGCCCGCCTCGTAGCTCAGGCCCGAGCGGGTGAAGGAACCGGCAGACAGGTAACCAGAGaaa > 2:712280/1‑150 (MQ=255) ggCAGGTGCGCGATCAGGCCGGCGCCGGTCAC‑GGCGAACAACGCCACCCACAGCGCTGAAAACACCCCGGCCATCGGTGAACGGGCGCCCGCCTCGTAGCTCAGGCCCGAGCGGGTGAAGGAACCGGCAGACAGGTAACCAGAGAAAAAt < 2:241302/150‑1 (MQ=255) gCAGGTGCGCGATCAGGCCGGCGCCGGTCAC‑GGCGAACAACGCCACCCACAGCGCTGAAAACACCCCGGCCATCGGTGAACGGGCGCCCGCCTCGTAGCTCAGGCCCGGGCGGGTGAAGGAACCGGCAGACAGGTAACCAGCGAAAAATg > 1:1083895/1‑150 (MQ=255) gTGCGCGATCAGGCCGGCGCCGGTCAC‑GGCGAACAACGCCACCCACAGCGCTGAAAACACCCCGGCCATCGGTGAGCGGGCGCCCGCCTCGTAGCTCAGGCCCGAGCGGGTGAAGGAACCGGCAGACAGGTAACCAGAGAAAAATGCAcc > 1:105139/1‑150 (MQ=255) gcgcGATCAGGCCGGCGCCGGTCAC‑GGCGAACAACGCCACCCACAGCGCTGAAAACACCCCGGCCATCGGTGAACGGGCGCCCGCCTCGTAGCTCAGGCCCGAGCGGGTGAAgg > 1:799/1‑114 (MQ=255) gcgcGATCAGGCCGGCGCCGGTCAC‑GGCGAACAACGCCACCCACAGCGCTGAAAACACCCCGGCCATCGGTGAACGGGCGCCCGCCTCGTAGCTCAGGCCCGAGCGGGTGAAgg < 2:799/114‑1 (MQ=255) cAGGCCGGCGCCGGTCAC‑GGCGAACAACGCCACCCACAGCGCTGAAAACACCCCGGCCATCGGTGAACGGGCGCCAGCCTCGTAGCTCAGGCCCGAGCGGGTGAAGGAACCGGCAGACAGGTAACCAGAGAAAAATGCACCGGCAATAtt < 1:17382/150‑1 (MQ=255) aGGCCGGCGCCGGTCAC‑GGCGAACAACGCCACCCACAGCGCTGAAAACACCCCGGCCATCGGTGAACGGGCGCCCGCCTCGTAGCTCAGGCCCGAGCGGGTGAAGGAACCGGCAGACAGGTAACCATATAAAAATGCACCGGCAATATTc < 2:210625/150‑1 (MQ=255) aGGCCGGCGCCGGTCAC‑GGCGAACAACGCCACCCACAGCGCTGAAAACACCCCGGCAATCGGTGAACGGGCGCCCGCCTCGTAGCTCAGGCCCGAGCGGGTGAAGGAACCGGCAGACAGGTAACCAGAGAAAAATGCACCGGCAATATTc > 2:116047/1‑150 (MQ=255) aGGCCGGCGCCGGTCAC‑GGCGAACAACGCCACCCACAGCGCTGAAAACACCCCGGCAATCGGTGAACGGGCGCCCGCCTCGTAGCTCAGGCCCGAGCGGGTGAAGGAACCGGCAGACAGGTAACCAGAGAAAAATGCACCGGCAATATTc > 2:666816/1‑150 (MQ=255) gcCGGTCAC‑GGCGAACAACGCCACCCACAGCGCTGAAAACACCCCGGCCATCGGTGAACGGGcgcccgcc > 1:657884/1‑70 (MQ=255) gcCGGTCAC‑GGCGAACAACGCCACCCACAGCGCTGAAAACACCCCGGCCATCGGTGAACGGGcgcccgcc < 2:657884/70‑1 (MQ=255) gAACAACGCCACCCACAGCGCTGAAAACACCCCGGCCATCGGTGAACGGGCGCCCGCCTCGTAGCTCAGGCCCGAGCGGGTGAAGGAACCGGCAGACAGGTAACCAGAGAAAAACGCACCGGCAATATTGGACAGGCCCTGTGCGCGTaa > 2:9724/1‑149 (MQ=255) ccccaccctcAGCGCTGAAACCGCCTCGGCCATCGGTGAACGGGCGCCCGCCTCGTAGCTCAGGCCCGAGCGGGTGAAGGAACCGGCAGACAGGTAACCAGAGAAAAATGCACCGGCAATATTCGACAGGCCCTGTGCGCGTATTTCCTg < 2:906683/148‑1 (MQ=255) cccacAGCGCTGAAAACACCCCGGCCATCGGTGAACGGGCGCCCGCCTCGTAGCTCAGGCCCGAGCGGGTGAAGGAACCGGCAGACAGGTAACCAGAGAAAAATGCACCGGCAATATTCGACAGGCCCTGTGCGCGTATTTCCTGGTCCg < 2:713258/150‑1 (MQ=255) tGAAAACACCCCGGCCATCGGTGAACGGGCGCCCGCCTCGTAGCTCAGGCCCGAGCGGGTGaa < 1:63868/63‑1 (MQ=255) tGAAAACACCCCGGCCATCGGTGAACGGGCGCCCGCCTCGTAGCTCAGGCCCGAGCGGGTGaa > 2:63868/1‑63 (MQ=255) gAAAACACCCCGGCAATCGGTGAACGGGCGCCCGCCTCGTAGCTCAGGCCCGAGCGGGTGAAGGAACCGGCAGACAGGTAACCAGAGAAAAATGCACCGGCAATATTCGACAGGCCCTGTGCGCGTATTTCCTGGTCCGGGTCGATCAgt < 1:116047/150‑1 (MQ=255) gAAAACACCCCGGCAATCGGTGAACGGGCGCCCGCCTCGTAGCTCAGGCCCGAGCGGGTGAAGGAACCGGCAGACAGGTAACCAGAGAAAAATGCACCGGCAATATTCGACAGGCCCTGTGCGCGTATTTCCTGGTCCGGGTCGATCAgt < 1:666816/150‑1 (MQ=255) aaaCACCCCGGCCATCGGTGAACGGGCGCCCGCCTCGTAGCTCAGGCCCGAGCGGGTGAAGGAACCGGCAGACAGGTAACCAGAGAAAAATGCACCGGCAATATTCGACAGGCCCTGTGCGCGTATTTCCTGGTCCGGGTCGATCAgttg > 2:965455/1‑150 (MQ=255) ccccGGCCATCGGTGAACGGGCGCCCGGTGCGTAGCTCAGGCCCGAGCGGGTGAAGGAACCGGCAGACAGGTAACCAGAGAAAAATGCACCGGCAATATTCGACAGGCCCTGTGCGCGTATTTCCTGGTCCGGGTCGATCAGTTGTTCTg > 1:308950/1‑150 (MQ=255) cccGGCCATCGGTGAACGGGCGCCCGCCTCGTAGCTCAGGCCCGAGCGGGTGAAGGAACCGGCAGAc < 1:185516/67‑1 (MQ=255) cccGGCCATCGGTGAACGGGCGCCCGCCTCGTAGCTCAGGCCCGAGCGGGTGAAGGAACCGGCAGAc > 2:185516/1‑67 (MQ=255) ccGGCCATCGGTGAACGGGCGCCCGCCTCGTAGCTCAGGCCCGAGCGGGTGAAGGAACCGGCAGACAGGTAACCAGAGAAAAATGCACCGGCAATATTCGACAGGCCCTGTGCGCGTATTTCCTGGTCCGGGTCGATCAGTTGTTCTGaa > 1:1087719/1‑150 (MQ=255) ccGGCCATCGGTGAACGGGCGCCCGCCTCGTAGCTCAGGCCCGAGCGGGTGAAGGAACCGGCAGACAGGTAACCAGAGAAAAAGGCACCGGCAATATTCGACAGGCCCTGTGCGCGTATTTCCTGGTCCGGGTCGATCAGTTGTTCTGaa < 2:646206/150‑1 (MQ=255) | GAACAATGCACGGATCCCTCGGTGGTCCACCAACCCCCAGCAGATCAGCAGAATGCTGCCGGCCATGGCCGGTATCGGCAGGTGCGCGATCAGGCCGGCGCCGGTCAC‑GGCGAACAACGCCACCCACAGCGCTGAAAACACCCCGGCCATCGGTGAACGGGCGCCCGCCTCGTAGCTCAGGCCCGAGCGGGTGAAGGAACCGGCAGACAGGTAACCAGAGAAAAATGCACCGGCAATATTCGACAGGCCCTGTGCGCGTATTTCCTGGTCCGGGTCGATCAGTTGTTCTGAA > NC_002947/85584‑85875 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |