Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | CP009273 | 3,303,385 | G→T | 8.9% | N381K (AAC→AAA) | pnp ← | polynucleotide phosphorylase/polyadenylase |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | CP009273 | 3,303,385 | 0 | G | T | 8.9% | 185.5 / 9.9 | 56 | N381K (AAC→AAA) | pnp | polynucleotide phosphorylase/polyadenylase |
Reads supporting (aligned to +/- strand): ref base G (35/16); new base T (3/2); total (38/18) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 6.52e-01 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.61e-01 |
TCAGACACAACACGTACGGTGTACGGGAATTTGTCCATATCCGGCATGACTGCCAGCACGCCGCGCTTCGCCAGACGACCGTGACCAATTTCACGACGCTTCGGAGAACCGACCATGCCGGTTTCGCCTACGGAGTACGGAGGGAAGTTGTAGTGGAACAGGAAGGTATCGGTACGTTCGCCCATCAGTTCATCAAGAACCTGCGCGTCACGAGCAGTACCCAGCGTTGCGGTAACCAGCGCCTGCGTTTCACCACGGGTGAACAGCGCAGAACCGTGAGTACGCGGCAGCACG > CP009273/3303239‑3303532 | tCAGACACAACACGTACGGTGTACGGGAATTTGTCCATATCCGGCATGACTGCCAGCACGCCGCGCTTCGCCAGACGACCGTGACCAATTTCACGACGCTTCGGAGAACCGACCGTGCCGGTTTCGCCTACGGAGTACGGAGGGAAGtt < 1:394012/149‑1 (MQ=255) aGACACAACACGTACGGTGTACGGGAATTTGTCCATATCCGGCATGACTGCCAGCACGCCGCGCTTCGCCAGACGACCGTGACCAATTTCACGACGCTTCGGAGAACCGACCATGCCGGTTTCGCCTACGGAGTACGGAGGGAAGTTg > 2:895478/1‑148 (MQ=255) cacaACACGTACGGTGTACGGGAATTTGTCCATATCCGGCATGACTGCCAGCACGCCGCGCTTCGCCAGACGACCGTGACCAATTTCACGACGCTTCGGAGAACCGACCATGCCGGTTTCGCCTACGGAGTACGGAGGGAAGTTGTAGTg > 1:854909/1‑150 (MQ=255) gTACGGTGTACGGGAATTTGTCCATATCCGGCATGACTGCCAGCACGCCGCGCTTCGCCAGACGACCGTGACCAATTTCACGACGCTTCGGAGAACCGACCATGCCGGTTTCGCCTACGGAGTACGGAGGGAAGTTGTAGTGGAACAgg > 2:297994/1‑149 (MQ=255) cGGTGTACGGGAATTTGTCCATATCCGGCATGACTGCCAGCACGCCGCGCTTCGCCAGACGACCGTGACCAATTTCACGACGCTTCGGAGAACCGACCATGCCGGTTTCGCCTACGGAGTACGGAGGGAAGTTGTAGTGGAACAGGAAgg > 1:973040/1‑150 (MQ=255) gtgtACGGGAATTTGTCCATATCCGGCATGACTGCCAGCACGCCGCGCTTCGCCAGACGACCGTGACCAATTTCACGACGCTTCGGAGAACCGACCATGCCGGTTTCGCCTACGGAGTACGGAGGGAAGTTGTAGTGGAACAGGAAGGTa < 2:757236/150‑1 (MQ=255) tgtACGGGAATTTGTCCATATCCGGCATGACTGCCAGCACGCCGCGCTTCGCCAGACGACCGTGACCAATTTCACGACGCTTCGGAGAACCGACCATGCCGGTTTCGCCTACGGAGTACGGAGGGAAGTTGTAGTGGAACAGGAAGGTAt > 2:930359/1‑150 (MQ=255) gtACGGGAATTTGTCCATATCCGGCATGACTGCCAGCACGCCGCGCTTCGCCAGACGACCGTGACCAATTTCACGACGCTTCGGAGAACCGACCATGCCGGTTTCGCCTACGGAGTACGGAGGGAAGTTGTAGTGGAACAGGAAGGTATc > 1:668282/1‑150 (MQ=255) cGGGAATTTGTCCATATCCGGCATGACTGCCAGCACGCCGCGCTTCGCCAGACGACCGTGACCAATTTCACGACGCTTCGGAGAACCGACCATGCCGGTTTCGCCTACGGAGTACGGAGGGAAGTTGTAGTGGAACAGGAAGGTATCgg > 1:696513/1‑149 (MQ=255) aaTTTGTCCATATCCGGCATGACTGCCAGCACGCCGCGCTTCGCCAGACGACCGTGACCAATTTCACGACGCTTCGGAGAACCGACCATGCCGGTTTCGCCTACGGAGTACGGAGGGAAGTTGTAGTGGAACAGGAAGGTATCGGTACGt > 2:85065/1‑150 (MQ=255) tttGTCCATATCCGGCATGACTGCCAGCACGCCGCGCTTCGCCAGACGACCGTGACCAATTTCACGACGCTTCGGAGAACCGACCATGCCGGTTTCGCCTACGGAGTACGGAGGGAAGTTGTAGTGGAACAGGAAGGTATCGGTACGTTc < 1:329549/150‑1 (MQ=255) tCCATATCCGGCATGACTGCCAGCACGCCGCGCTTCGCCAGACGACCGTGACCAATTTCACGACGCTTCGGAGAACCGACCATGCCGGTTTCGCCTACGGAGTACGGAGGGAAGTTGTAGTGGAACAGGAAGGTATCGGTACGTTCGcc < 1:803834/149‑1 (MQ=255) ccATATCCGGCATGACTGCCAGCACGCCGCGCTTCGCCAGACGACCGTGACCAATTTCACGACGCTTCGGAGAACCGACCATGCCGGTTTCGCCTACGGAGTACGGAGGGAAGTTGTAGTGGAACAGGAAGGTATCGGTACGTTCGCCCa > 2:950284/1‑150 (MQ=255) ccATATCCGGCATGACTGCCAGCACGCCGCGCTTCGCCAGACGACCGTGACCAATTTCACGACGCTTCGGAGAACCGACCATGCCGGTTTCGCCTACGGAGTACGGAGGGAAGTTGTAGTGGAACAGGAAGGTATCGGTACGTTCGCCCa > 2:567303/1‑150 (MQ=255) tCCGGCATGACTGCCAGCACGCCGCGCTTCGCCAGACGACCGTGACCAATTTCACGACGCTTCGGAGAACCGACCATGCCGGTTTCGCCTACGGAGTACGGAGGGAAGTTGTAGTGGAACAGGAAGGTATCGGTACGTTCGCCCATCAGt < 1:995276/150‑1 (MQ=255) ccGGCATGACTGCCAGCACGCCGCGCTTCGCCAGACGACCGTGACCAATTTCACGACGCTTCGGAGAACCGACCATGCCGGTTTCGCCTACGGAGTACGGAGGGAAGTTGTAGTGGAACAGGAAGGTATCGGTACGTTCGCCCATCAGtt > 1:342471/1‑150 (MQ=255) gCATGACTGCCAGCACGCCGCGCTTCGCCAGACGACCGTGACCAATTTCACGACGCTTCGGAGAACCGACCATGCCGGTTTCGCCTACGGAGTACGGAGGGAAGTTGTAGTGGAACAGGAAGGTATCGGTACGTTCGCCCATCAGTtcat > 1:586107/1‑150 (MQ=255) gCATGACTGCCAGCACGCCGCGCTTCGCCAGACGACCGTGACCAATTTCACGACGCTTCGGAGAACCGACCATGCCGGTTTCGCCTACGGAGTACGGAGGGAAGTTGTAGTGGAACAGGAAGGTATCGGTACGTTCGCCCATCAGTtcat > 1:464662/1‑150 (MQ=255) aTGACTGCCAGCACGCCGCGCTTCGCCAGACGACCGTGACCAATTTCACGACGCTTCGGAGAACCGACCATGCCGGTTTCGCCTACGGAGTACGGAGGGAATTTGTAGTGGAACAGGAAGGTATCGGTACGTTCGCCCATCAGTtcatc > 2:695917/1‑149 (MQ=255) gCCAGCACGCCGCGCTTCGCCAGACGACCGTGACCAATTTCACGACGCTTCGGAGAACCGACCATGCCGGTTTCGCCTACGGAGTACGGAGGGAATTTGTAGTGGAACAGGAAGGTATCGGTACGTTCGCCCATCAGTTCATCAAGAAcc > 2:785089/1‑150 (MQ=255) ttCGCCAGACGACCGTGACCAATTTCACGACGCTTCGGAGAACCGACCATGCCGGTTTCGCCTACGGAGTACGGAGGGAAGTTGTAGTGGAACAGGAAGGTATCGGTACGTTCGCCCATCAGTTCATCAAGAACCTGCGCGTCACGAg > 1:856643/1‑148 (MQ=255) cgacCGTGACCAATTTCACGACGCTTCGGAGAACCGACCATGCCGGTTTCGCCTACGGAGTACGGAGGGAAGTTGTAGTGGAACAGGAAGGTATCGGTACGTTCGCCCATCAGTTCATCAAGAACCTGCGCGTCACGAGCAGTACCCAGc < 2:757423/150‑1 (MQ=255) ccGTGACCAATTTCACGACGCTTCGGAGAACCGACCATGCCGGTTTCGCCTACGGAGTACGGAGGGAAGTTGTAGTGGAACAGGAAGGTATCGGTACGTTCGCCCATCAGTTCATCAAGAACCTGCGCGTCACGAGCAGTACCCAGCGtt < 2:737040/150‑1 (MQ=255) aCCAATTTCACGACGCTTCGGAGAACCGACCATGCCGGTTTCGCCTACGGAGTACGGAGGGAAGTTGTAGTGGAACAGGAAGGTATCGGTACGTTCGCCCATCAGTTCATCAAGAACCTGCGCGTCACGAGCAGTACCCAGCGTTGCgg < 1:884718/149‑1 (MQ=255) cAATTTCACGACGCTTCGGAGAACCGACCATGCCGGTTTCGCCTACGGAGTACGGAGGGAAGTTGTAGTGGAACAGGAAGGTATCGGTACGTTCGCCCATCAGTTCATCAAGAACCTGCGCGTCACGAGCAGTACCCAGCGTTGCGGTaa > 2:793935/1‑150 (MQ=255) aaTTTCACGACGCTTCGGAGAACCGACCATGCCGGTTTCGCCTACGGAGTACGGAGGGAATTTGTAGTGGAACAGGAAGGTATCGGTACGTTCGCCCATCAGTTCATCAAGAACCTGCGCGTCACGAGCAGTACCCAGCGTTGCGGTAAc > 2:829734/1‑150 (MQ=255) aTTTCACGACGCTTCGGAGAACCGACCATGCCGGTTTCGCCTACGGAGTACGGAGGGAAGTTGTAGTGGAACAGGAAGGTATCGGTACGTTCGCCCATCAGTTCATCAAGAACCTGCGCGTCACGAGCAGTACCCAGCGTTGCGGTAAcc < 1:851488/150‑1 (MQ=255) tttCACGACGCTTCGGAGAACCGACCATGCCGGTTTCGCCTACGGAGTACGGAGGGAAGTTGTAGTGGAACAGGAAGGTATCGGTACGTTCGCCCATCAGTTCATCAAGAACCTGCGCGTCACGAGCAGTACCCAGCGTTGCGGTAAc > 1:299297/1‑148 (MQ=255) tCACGACGCTTCGGAGAACCGACCATGCCGGTTTCGCCTACGGAGTACGGAGGGAAGTTGTAGTGGAACAGGAAGGTATCGGTACGTTCGCCCATCAGTTCATCAAGAACCTGCGCGTCACGAGCAGTACCCAGCGTTGCGGTAACCAgc > 2:877998/1‑150 (MQ=255) acgacgCTTCGGAGAACCGACCATGCCGGTTTCGCCTACGGAGTACGGAGGGAAGTTGTAGTGGAACAGGAAGGTATCGGTACGTTCGCCCATCAGTTCATCAAGAACCTGCGCGTCACGAGCAGTACCCAGCGTTGCGGTAACCAgcgc < 1:589052/150‑1 (MQ=255) gacgCTTCGGAGAACCGACCATGCCGGTTTCGCCTACGGAGTACGGAGGGAAGTTGTAGTGGAACAGGAAGGTATCGGTACGTTCGCCCATCAGTTCATCAAGAACCTGCGCGTCACGAGCAGTACCCAGCGTTGCGGTAACCAGCGCCt > 1:822470/1‑150 (MQ=255) cgCTTCGGAGAACCGACCATGCCGGTTTCGCCTACGGAGTACGGAGGGAAGTTGTAGTGGAACAGGAAGGTATCGGTACGTTCGCCCATCAGTTCATCAAGAACCTGCGCGTCACGAGCAGTACCCAGCGTTGCGGTAACCAGCGCCTGc > 2:578356/1‑150 (MQ=255) cgCTTCGGAGAACCGACCATGCCGGTTTCGCCTACGGAGTACGGAGGGAAGTTGTAGTGGAACAGGAAGGTATCGGTACGTTCGCCCATCAGTTCATCAAGAACCTGCGCGTCACGAGCAGTACCCAGCGTTGCGGTAACCAGCGCCTGc > 1:941232/1‑150 (MQ=255) gagaACCGACCATGCCGGTTTCGCCTACGGAGTACGGAGGGAAGTTGTAGTGGAACAGGAAGGTATCGGTACGTTCGCCCATCAGTTCATCAAGAACCTGCGCGTCACGAGCAGTACCCAGCGTTGCGGTAACCAGCGCCTGCGTTTcac < 1:84389/150‑1 (MQ=255) agaACCGACCATGCCGGTTTCGCCTACGGAGTACGGAGGGAAGTTGTAGTGGAACAGGAAGGTATCGGTACGTTCGCCCATCAGTTCATCAAGAACCTGCGCGTCACGAGCAGTACCCAGCGTTGCGGTAACCAGCGCCTGCGTTTcac > 1:917384/1‑149 (MQ=255) gaACCGACCATGCCGGTTTCGCCTACGGAGTACGGAGGGAAGTTGTAGTGGAACAGGAAGGTATCGGTACGTTCGCCCATCAGTTCATCAAGAACCTGCGCGTCACGAGCAGTACCCAGCGTTGCGGTAACCAGCGCCTGCGTTTcacca > 2:628750/1‑150 (MQ=255) aaCCGACCATGCCGGTTTCGCCTACGGAGTACGGAGGGAATTTGTAGTGGAACAGGAAGGTATCGGTACGTTCGCCCATCAGTTCATCAAGAACCTGCGCGTCACGAGCAGTACCCAGCGTTGCGGTAACCAGCGCCTGCGTTTcaccac < 1:695917/150‑1 (MQ=255) aaCCGACCATGCCGGTTTCGCCTACGGAGTACGGAGGGAAGTTGTAGTGGAACAGGAAGGTATCGGTACGTTCGCCCATCAGTTCATCAAGAACCTGCGCGTCACGAGCAGTACCCAGCGTTGCGGTAACCAGCGCCTGCGTTTcaccac < 1:950284/150‑1 (MQ=255) cGACCATGCCGGTTTCGCCTACGGAGTACGGAGGGAAGTTGTAGTGGAACAGGAAGGTATCGGTACGTTCGCCCATCAGTTCATCAAGAACCTGCGCGTCACGAGCAGTACCCAGCGTTGCGGTAACCAGCGCCTGCGTTTCACCACggg > 1:291099/1‑150 (MQ=255) gACCATGCCGGTTTCGCCTACGGAGTACGGAGGGAAGTTGTAGTGGAACAGGAAGGTATCGGTACGTTCGCCCATCAGTTCATCAAGAACCTGCGCGTCACGAGCAGTACCCAGCGTTGCGGTAACCAGCGCCTGCGTTTCACCACGGGt > 2:535458/1‑150 (MQ=255) gACCATGCCGGTTTCGCCTACGGAGTACGGAGGGAAGTTGTAGTGGAACAGGAAGGTATCGGTACGTTCGCCCATCAGTTCATCAAGAACCTGCGCGTCACGAGCAGTACCCAGCGTTGCGGTAACCAGCGCCTGCGTTTCACCACGGGt < 1:793935/150‑1 (MQ=255) aCCATGCCGGTTTCGCCTACGGAGTACGGAGGGAATTTGTAGTGGAACAGGAAGGTATCGGTACGTTCGCCCATCAGTTCATCAAGAACCTGCGCGTCACGAGCAGTACCCAGCGTTGCGGTAACCAGCGCCTGCGTTTCACCACGGGTg < 1:783375/150‑1 (MQ=255) aCCATGCCGGTTTCGCCTACGGAGTACGGAGGGAAGTTGTAGTGGAACAGGAAGGTATCGGTACGTTCGCCCATCAGTTCATCAAGAACCTGCGCGTCACGAGCAGTACCCAGCGTTGCGGTAACCAGCGCCTGCGTTTCACCACGGGt > 2:54452/1‑149 (MQ=255) aCCATGCCGGTTTCGCCTACGGAGTACGGAGGGAAGTTGTAGTGGAACAGGAAGGTATCGGTACGTTCGCCCATCAGTTCATCAAGAACCTGCGCGTCACGAGCAGTACCCAGCGTTGCGGTAACCAGCGCCTGCGTTTCACCACGGGt > 2:870965/1‑149 (MQ=255) cATGCCGGTTTCGCCTACGGAGTACGGAGGGAAGTTGTAGTGGAACAGGAAGGTATCGGTACGTTCGCCCATCAGTTCATCAAGAACCTGCGCGTCACGAGCAGTACCCAGCGTTGCGGTAACCAgcgc > 1:563992/1‑129 (MQ=255) aTGCCGGTTTCGCCTACGGAGTACGGAGGGAAGTTGTAGTGGAACAGGAAGGTATCGGTACGTTCGCCCATCAGTTCATCAAGAACCTGCGCGTCACGAGCAGTACCCAGCGTTGCGGTAACCAGCGcc < 2:563992/129‑1 (MQ=255) aTGCCGGTTTCGCCTACGGAGTACGGAGGGAAGTTGTAGTGGAACAGGAAGGTATCGGTACGTTCGCCCATCAGTTCATCAAGAACCTGCGCGTCACGAGCAGTACCCAGCGTTGCGGTAACCAGCGCCTGCGTTTCACCACGGGTGaa > 2:498281/1‑149 (MQ=255) gCCGGTTTCGCCTACGGAGTACGGAGGGAAGTTGTAGTGGAACAGGAAGGTATCGGTACGTTCGCCCATCAGTTCATCAAGAACCTGCGCGTCACGAGCAGTACCCAGCGTTGCGGTAACCAGCGCCTGCGTTTCACCACGGGTGAACAg > 1:329531/1‑150 (MQ=255) cGGTTTCGCCTACGGAGTACGGAGGGAAGTTGTAGTGGAACAGGAAGGTATCGGTACGTTCGCCCATCAGTTCATCAAGAACCTGCGCGTCACGAGCAGTACCCAGCGTTGCGGTAACCAGCGCCTGCGTTTCACCACGGGTGAACAgcg < 1:886457/150‑1 (MQ=255) gTTTCGCCTACGGAGTACGGAGGGAAGTTGTAGTGGAACAGGAAGGTATCGGTACGTTCGCCCATCAGTTCATCAAGAACCTGCGCGTCACGAGCAGTACCCAGCGTTGCGGTAACCAGCGCCTGCGTTTCACCACGGGTGAACAgcg > 2:15419/1‑148 (MQ=255) cGCCTACGGAGTACGGAGGGAAGTTGTAGTGGAACAGGAAGGTATCGGTACGTTCGCCCATCAGTTCATCAAGAACCTGCGCGTCACGAGCAGTACCCAGCGTTGCGGTAACCAGCGCCTGCGTTTCACCACGGGTGAACAGCGCAGAAc > 1:1003543/1‑150 (MQ=255) ccTACGGAGTACGGAGGGAAGTTGTAGTGGAACAGGAAGGTATCGGTACGTTCGCCCATCAGTTCATCAAGAACCTGCGCGTCACGAGCAGTACCCAGCGTTGCGGTAACCAGCGCCTGCGTTTCACCACGGGTGAACAGCGCAGAACCg > 2:578933/1‑150 (MQ=255) cTACGGAGTACGGAGGGAAGTTGTAGTGGAACAGGAAGGTATCGGTACGTTCGCCCATCAGTTCATCAAGAACCTGCGCGTCACGAGCAGTACCCAGCGTTGCGGTAACCAGCGCCTGCGTTTCACCACGGGTGAACAGCGCAGAACCGt > 2:238710/1‑150 (MQ=255) cggagtacggagGGAAGTTGTAGTGGAACAGGAAGGTATCGGTACGTTCGCCCATCAGTTCATCAAGAACCTGCGCGTCACGAGCAGTACCCAGCGTTGCGGTAACCAGCGCCTGCGTTTCACCACGGGTGAACAGCGCAGAAACGTg < 1:281608/148‑1 (MQ=255) ggacGGAAGTTGTAGTGGAACAGGAAGGTATCGGTACGTTCGCCCATCAGTTCATCAAGAACCTGCGCGTCACGAGCAGTACCCAGCGTTGCGGTAACCAGCGCCTGCGTTTCACCACGGGTGAACAGCGCAGAACCGTGAGTACGCgg > 1:701711/1‑149 (MQ=255) aaGTTGTAGTGGAACAGGAAGGTATCGGTACGTTCGCCCATCAGTTCATCAAGAACCTGCGCGTCACGAGCAGTACCCAGCGTTGCGGTAACCAGCGCCTGCGTTTCACCACGGGTGAACAGCGCAGAACCGTGAGTACGCGGCAGCACg > 2:363479/1‑150 (MQ=255) | TCAGACACAACACGTACGGTGTACGGGAATTTGTCCATATCCGGCATGACTGCCAGCACGCCGCGCTTCGCCAGACGACCGTGACCAATTTCACGACGCTTCGGAGAACCGACCATGCCGGTTTCGCCTACGGAGTACGGAGGGAAGTTGTAGTGGAACAGGAAGGTATCGGTACGTTCGCCCATCAGTTCATCAAGAACCTGCGCGTCACGAGCAGTACCCAGCGTTGCGGTAACCAGCGCCTGCGTTTCACCACGGGTGAACAGCGCAGAACCGTGAGTACGCGGCAGCACG > CP009273/3303239‑3303532 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 29 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |