Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | NC_002947 | 5,883,651 | C→G | 7.0% | A60G (GCC→GGC) | PP_5157 → | hypothetical protein |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_002947 | 5,883,651 | 0 | C | G | 7.0% | 181.8 / 5.3 | 57 | A60G (GCC→GGC) | PP_5157 | hypothetical protein |
Reads supporting (aligned to +/- strand): ref base C (29/24); new base G (2/2); total (31/26) | |||||||||||
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 |
GTACTGATCAGTCCATTGGCTCAGGCGGAGCTGATCGACGAAATCGCCGACCGGGGCGAACTGCGCATTGCCGTACAGGCCGAGAACTCGCCATACGCCTTCAAGCAAGACGACCACCTGACAGGTTTCGACATCGAATTTGGCCAGGCCCTGGCCAAGGAACTGGATTTGCGCGCCGAGTTCGTCGAAGCCCCGGCTACCGAAGTGCTGGCGGGTGTCGAGAGTGGCAAGTACGACATCGCATTCACCCCGTCGAGTGAATCGCTCACAGGCGACGGCCCGCTCGATGTAAGCCTGC > NC_002947/5883503‑5883800 | gTACTGATCAGTCCATTGGCTCAGGCGGAGCTGATCGACGAAATCGCCGACCGGGGCGAACTGCGCATTGCCGTACAGGCCGAGAACTCGCCATACGCCTTCAAGCAAGACGACCACCTGACAGGTTTCGACATCGAATTTGGCCAGGcc > 2:249178/1‑150 (MQ=255) cTGGTCAGTCGAGTGGCTCAGGCGGAGCGGAGCGACGAAATCGCCGACCGGGGCGAACTGCGCATTGCCGTACAGGCCGAGAACTCGCCATACGCCTTCAAGCAAGACGACCACCTGACAGGTTTCGACAGCGAATTTGGCCAGGCCCTg < 2:657584/150‑1 (MQ=255) tGATCAGTCCATTGGCTCAGGCGGAGCTGATCGACGAAATCGCCGACCGGGGCGAACTGCGCATTGCCGTACAGGCCGAGAACTCGCCATACGCCTTCAAGCAAGACGACCACCTGACAGGTTTCGACATCGAATTTGGCCAGGCCCTgg > 2:168860/1‑150 (MQ=255) tCCATTGGCTCAGGCGGAGCTGATCGACGAAATCGCCGACCGGGGCGAACTGCGCATTGCCGTACAGGCCGAGAACTCGCCATACTCCTTCAAGCAAGACGACCACCTGACAGGTTTCGACATCGAATTTGGCCAGGCCCTGGCCAAGGa > 2:1051455/1‑150 (MQ=255) ccATTGGCTCAGGCGGAGCTGATCGACGAAATCGCCGACCGGGGCGAACTGCGCATTGCCGTACAGGCCGAGAACTCGCCATACGCCTTCAAGCAAGACGACCACCTGACAGGTTTCGACATCGAATTTGGCCAGGCCCTGGCCAAGGaa < 1:17987/150‑1 (MQ=255) cATTGGCTCAGGCGGAGCTGATCGACGAAATCGCCGACCGGGGCGAACTGCGCATTGCCGTACAGGCCGAGAACTCGCCATACTCCTTCAAGCAAGACGACCACCTGACAGGTTTCGACATCGAATTTGGCCAGGCCCTGGCCAAGGAAc < 1:1051455/150‑1 (MQ=255) gCGGAGCTGATCGACGAAATCGCCGACCGGGGCGAACTGCGCATTGCCGTACAGGCCGAGAACTCGCCATACGCCTTCAAGCAAGACGACCACCTGACAGGTTTCGACATCGAATTTGGCCAGGCCCTGGCCAAGGAACTGGATTTgcgc < 1:1021929/150‑1 (MQ=255) gCTGATCGACGAAATCGCCGACCGGGGCGAACTGCGCATTGCCGTACAGGCCGAGAACTCGCCATACGCCTTCAAGCAAGACGACCACCTGACAGGTTTCGACATCGAATTTGGCCAGGCCCTGGCCAAGGAACTGGATTTGCGCGCCGa < 1:168860/150‑1 (MQ=255) gATCGACGAAATCGCCGACCGGGGCGAACTGCGCATTGCCGTACAGGCCGAGAACTCGCCATACGCCTTCAAGCAAGACGACCACCTGACAGGTTTCGACATCGAATTTGGCCAGGCCCTGGCCAAGGAACTGGATTTGCGCGCCGAGtt > 2:390628/1‑150 (MQ=255) cgacgaAATCGCCGACCGGGGCGAACTGCGCATTGCCGTACAGGCCGAGAACTCGCCATACGCCTTCAAGCAAGACGACCACCTGACAGGTTTCGACATCGAATTTGGCCAGGCCCTGGCCAAGGAACTGGATTTGCGCGCCGAGTtcgt < 2:902469/150‑1 (MQ=255) acgaAATCGCCGACCGGGGCGAACTGCGCATTGCCGTACAGGCCGAGAACTCGCCATACGCCTTCAAGCAAGACGACCACCTGACAGGTTTCGACATCGAATTTGGCCAGGCCCTGGCCAAGGAACTGGATTTGCGCGCCGAGTtcgtcg > 1:472178/1‑150 (MQ=255) gCCGACCGGGGCGAACTGCGCATTGCCGTACAGGCCGAGAACTCGCCATACGCCTTCAAGCAAGACGACCACCTGACAGGTTTCGACATCGAATTTGGCCAGGCCCTGGCCAAGGAACTGGATTTGCGCGCCGAGTTCGTCGAAGCCCCg > 1:1098175/1‑150 (MQ=255) ccGACCGGGGCGAACTGCGCATTGCCGTACAGGCCGAGAACTCGCCATACGCATAACAGCAACTCGTCCACCTGACAGGTTTCGACATCGAATTTGGCCAGGCCCTGGCCAAGGAACTGGATTTGCGCGCCGAGTTCGTCGAAGCCCCgg > 2:221968/1‑150 (MQ=255) gggCGAACTGCGCATTGCCGTACAGGCCGAGAACTCGCCATACGCCTTCAAGCAAGACGACCACCTGAAAGGTTTGGACATCGAATTTGGCCAGGCACTGGACAAGGAACTGGATTTGCGCGCCGAGTTCGTCGAAGTCCCGGCTagata > 1:758807/1‑146 (MQ=255) ggCGAACTGCGCATTGCCGTACAGGCCGAGAACTCGCCATACGCCTTCAAGCAAGACGACCACCTGACAGGTTTCGACATCGAATTTGGCCAGGCCCTGGCCAAGGAACTGGATTTGCGCGCCGAGTTCGTCGAAGCCCCGGCTACCGaa < 2:771016/150‑1 (MQ=255) tGCGCATTGCGGTACAGGGCTAGAACGCGCCATACGCCTTCAAGCAAGACGACCACCTGACAGGTTTCGACATCGAATTGGGCCAGGCCCTGGCCAAGGAACTGGATTTGCGCGCCGAGTTCGTCGAAGCCCCGGCTACCGAAGTGCTgg < 2:707043/150‑1 (MQ=255) tGCGCATTGCCGTACAGGCCGAGAACTCGCCATACGCCTTCAAGCAAGACGACCACCTGACAGGTTTCGACATCGAATTTGGCCAGGCCCTGGCCaa < 2:1255500/97‑1 (MQ=255) tGCGCATTGCCGTACAGGCCGAGAACTCGCCATACGCCTTCAAGCAAGACGACCACCTGACAGGTTTCGACATCGAATTTGGCCAGGCCCTGGCCaa > 1:1255500/1‑97 (MQ=255) gcATTGCCGTACAGGCCGAGAACTCGCCATACGCCTTCAAGCAAGACGACCACCTGACAGGTTTCGACATCGAATTTGGCCAGGCCCTGGCCAAGGAACTGGATTTGCGCGCCGAGTTCGTCGAAGCCCCGGCTACCGAAGTGCTGGCgg < 1:390628/150‑1 (MQ=255) aTTGCCGTACAGGCCGAGAACTCGCCATACGCCTTCAAGCAAGACGACCACCTGACAGGTTTCGACATCGAATTTGGCCAGGCCCTGGCCAAGGAACTGGATTTGCGCGCCGAGTTCGTCGAAGCCCCGGCTACCGAAGTGCTGGCGGgt > 1:693804/1‑150 (MQ=255) cGTACAGGCCGAGAACTCGCCATACGCCTTCAAGCAAGACGACCACCTGACAGGTTTCGACATCGAATTTGGCCAGGCCCTGGCCAAGGAACTGGATTTGCGCGCCGAGTTCGTCGAAGCCCCGGCTACCGAAGTGCTGGCGGGTGTCga > 1:954650/1‑150 (MQ=255) cGTACAGGCCGAGAACTCGCCATACGCCTTCAAGCAAGACGACCACCTGACAGGTTTCGACATCGAATTTGGCCAGGCCCTGGCCAAGGAACTGGATTTGCGCGCCGAGTTCGTCGAAGCCCCGGCTACCGAAGTGCTGGCGGGTGTCga > 1:174325/1‑150 (MQ=255) tACAGGCCGAGAACTCGCCATACGCCTTCAAGCAAGACGACCACCTGACAGGTTTCGACATCGAATTTGGCCAGGCCCTGGCCAAGGAACTGGATTTGCGCGCCGAGTTCGTCGAAGCCCCGGCTACCGAAGTGCTGGCGGGTGTCgaga < 2:693804/150‑1 (MQ=255) aGGCCGAGAACTCGCCATACGCCTTCAAGCAAGACGACCACCTGACAGGTTTCGACATCGAATTTGGCCAGGCCCTGGCCAAGGAACTGGATTTGCGCGCCGAGTTCGTCGAAGCCCCGGCTACCGAAGTGCTGGCGGGTGTCGAGAGTg < 1:1538494/150‑1 (MQ=255) aGGCCGAGAACTCGCCATACGCCTTCAAGCAAGACGACCACCTGACAGGTTTCGACATCGAATTTGGCCAGGCCCTGGCCAAGGAACTGGATTTGCGCGCCGAGTTCGTCGAAGCCCCGGCTACCGAAGTGCTGGCGGGTGTCGAGAGTg > 1:1578260/1‑150 (MQ=255) gCCGAGAACTCGCCATACGCCTTCAAGCAAGACGACCACCTGACAGGTTTCGACATCGAATTTGGCCAGGCCCTGGCCAAGGAACTGGATTTGCGCGCCGAGTTCGTCGAAGCCCCGGCTACCGAAGTGCTGGCGGGTGTCGAGAGTGGc < 2:954650/150‑1 (MQ=255) gCCGAGAACTCGCCATACGCCTTCAAGCAAGACGACCACCTGACAGGTTTCGACATCGAATTTGGCCAGGCCCTGGCCAAGGAACTGGATTTGCGCGCCGAGTTCGTCGAAGCCCCGGCTACCGAAGTGCTGGCGGGTGTCGAGAGTGGc > 1:857141/1‑150 (MQ=255) aTACGCCTTCAAGCAAGACGACCACCTGACAGGTTTCGACATCGAATTTGGCCAGGCCCTGGCCAAGGAACTGGATTTGCGCGCCGAGTTCGTCGAAGCCCCGGCTACCGAAGTGCTGGCGGGTGTCGAGAGTGGCAAGTACGACATCGc > 2:1103643/1‑150 (MQ=255) cGCCTTCAAGCAAGACGACCACCTGACAGGTTTCGACATCGAATTTGGCCAGGCCCTGGCCAAGGAACTGGATTTGCGCGCCGAGTTCGTCGAAGCCCCGGCTACCGAAGTGCTGGCGGGTGTCGAGAGTGGCAAGTACGACATCGCAtt < 2:857141/150‑1 (MQ=255) cTTCAAGCAAGACGACCACCTGACAGGTTTCGACATCGAATTTGGCCAGGCCCTGGCCAAGGAACTGGATTTGCGCGCCGAGTTCGTCGaa > 1:1484658/1‑91 (MQ=255) cTTCAAGCAAGACGACCACCTGACAGGTTTCGACATCGAATTTGGCCAGGCCCTGGCCAAGGAACTGGATTTGCGCGCCGAGTTCGTCGaa < 2:1484658/91‑1 (MQ=255) cTTCAAGCAAGACGACCACCTGACAGGTTTCGACATCGAATTTGGCCAGGCCCTGGCCAAGGAACTGGATTTGCGCGCCGAGTTCGTCGAAGCCCCGGCTACCGAAGTGCTGGCGGGTGTCGAGAGTGGCAAGTACGACATCGCATTCAc < 2:472178/150‑1 (MQ=255) aagcaagACGACCACCTGACAGGTTTCGACATCGAATTTGGCCAGGCCCTGGCCAAGGAACTGGATTTGCGCGCCGAGTTCGTCGAAGCCCCGGCTACCGAAGTGCTGGCGGGTGTCGAGAGTGGCAAGTACGACATCGCATTCACCCCg > 2:284597/1‑150 (MQ=255) gcaagACGACCACCTGACAGGTTTCGACATCGAATTTGGCCAGGCCCTGGCCAAGGAACTGGATTTGCGCGCCGAGTTCGTCGAAGCCCCGGCTACCGAAGTGCTGGCGGGTGTCga > 1:451033/1‑117 (MQ=255) gcaagACGACCACCTGACAGGTTTCGACATCGAATTTGGCCAGGCCCTGGCCAAGGAACTGGATTTGCGCGCCGAGTTCGTCGAAGCCCCGGCTACCGAAGTGCTGGCGGGTGTCga < 2:451033/117‑1 (MQ=255) gcaagACGACCACCTGACAGGTTTCGACATCGAATTTGGCCAGGCCCTGGCCAAGGAACTGGATTTGCGCGCCGAGTTCGTCGAAGCCCCGGCTACCGAAGTGCTGGCGGGTGTCGagag > 1:175818/1‑120 (MQ=255) gcaagACGACCACCTGACAGGTTTCGACATCGAATTTGGCCAGGCCCTGGCCAAGGAACTGGATTTGCGCGCCGAGTTCGTCGAAGCCCCGGCTACCGAAGTGCTGGCGGGTGTCGagag < 2:175818/120‑1 (MQ=255) aagACGACCACCTGACAGGTTTCGACATCGAATTTGGCCAGGCCCTGGCCAAGGAACTGGATTTGCGCGCCGAGTTCGTCGAAGCCCCGGCTACCGAAGTGCTTGCGGGTGTTGAGAGTGGCAAGTACGACATCGCATTCACCCCGTCGa > 1:513485/1‑150 (MQ=255) cTGACAGGTTTCGACATCGAATTTGGCCAGGCCCTGGCCAAGGAACTGGATTTGCGCGCCGAGTTCGTCGAAGCCCCGGCTACCGAAGTGCTGGCGGGTGTCGAGAGTGGCAAGTACGACATCGCATTCACCCCGTCGAGTGaa < 2:947591/144‑1 (MQ=255) cTGACAGGTTTCGACATCGAATTTGGCCAGGCCCTGGCCAAGGAACTGGATTTGCGCGCCGAGTTCGTCGAAGCCCCGGCTACCGAAGTGCTGGCGGGTGTCGAGAGTGGCAAGTACGACATCGCATTCACCCCGTCGAGTGaa > 1:947591/1‑144 (MQ=255) aCAGGTTTCGACATCGAATTTGGCCAGGCCCTGGCCAAGGAACTGGATTTGCGCGCCGAGTTCGTCGAAGCCCCGGCTACCGAAGTGCTGGCGGGTGTCGAGAGTGGCAAGTACGACATCGCATTCACCCCGTCGAGTGAATCGCTcaca > 1:1358464/1‑150 (MQ=255) tttCGACATCGAATTTGGCCAGGCCCTGGCCAAGGAACTGGATTTGCGCGCCGAGTTCGTCGAAGCCCCGGCTACCGaa < 2:1427962/79‑1 (MQ=255) tttCGACATCGAATTTGGCCAGGCCCTGGCCAAGGAACTGGATTTGCGCGCCGAGTTCGTCGAAGCCCCGGCTACCGaa > 1:1427962/1‑79 (MQ=255) tttCGACATCGAATTTGGCCAGGCCCTGGCCAAGGAACTGGATTTGCGCGCCGAGTTCGTCGAAGCCCCGGCTACCGAAGTGCTGGCGGGTGTCGAGAGTGGCAAGTACGACATCGCATTCACCCCGTCGAGTGAATCGCTCACAGGCGa > 2:1566502/1‑150 (MQ=255) aCATCGAATTTGGCCAGGCCCTGGCCAAGGAACTGGATTTGCGCGCCGAGTTCGTCGAAGCCCCGGCTACCGAAGTGCTGGCGGGTGTCGAGAGTgg > 2:1086073/1‑97 (MQ=255) aCATCGAATTTGGCCAGGCCCTGGCCAAGGAACTGGATTTGCGCGCCGAGTTCGTCGAAGCCCCGGCTACCGAAGTGCTGGCGGGTGTCGAGAGTgg < 1:1086073/97‑1 (MQ=255) aTCGAATTTGGCCAGGCCCTGGCCAAGGAACTGGATTTGCGCGCCGAGTTCGTCGAAGCCCCGGCTACCGAAGTGCTGGCGGGTGTCGAGAGTGGCAAGTACGACATCGCATTCACCCCGTCGAGTGAATCGCTCACAGGCGACGGCCCg > 2:216379/1‑150 (MQ=255) tCGAATTTGGCCAGGCCCTGGCCAAGGAACTGGATTTGCGCGCCGAGTTCGTCGAAGCCCCGGCTACCGAAGTGCTGGCGGGTGTCGAGAGTGGCAAGTACGACATCGCATTCACCCCGTCGAGTGAATCGCTcaca > 1:346963/1‑137 (MQ=255) tCGAATTTGGCCAGGCCCTGGCCAAGGAACTGGATTTGCGCGCCGAGTTCGTCGAAGCCCCGGCTACCGAAGTGCTGGCGGGTGTCGAGAGTGGCAAGTACGACATCGCATTCACCCCGTCGAGTGAATCGCTcaca < 2:346963/137‑1 (MQ=255) gAATTTGGCCAGGCCCTGGCCAAGGAACTGGATTTGCGCGCCGAGTTCGTCGAAGCCCCGGCTACCGAAGTGCTGGCGGGTGTCGAGAGTGGCAAGTACGACATCGCATTCACCCCGTCGAGTGAATCGCTCACAGGCGACGGCCCGCTc > 1:656759/1‑150 (MQ=255) ccttGGCCAGGGCCTGGCCAAGGAACTGGATTTGCGCGCCGAGTTCGTCGAAGCCCCGGCTACCGAAGTGCTGGCGGGTGTCGAGAGTGGCAAGTACGACATCGCAtt < 1:242982/106‑1 (MQ=255) ccttGGCCAGGGCCTGGCCAAGGAACTGGATTTGCGCGCCGAGTTCGTCGAAGCCCCGGCTACCGAAGTGCTGGCGGGTGTCGAGAGTGGCAAGTACGACATCGCAtt > 2:242982/3‑108 (MQ=255) ttGGCCAGGGCCTGGCCAAGGAACTGGATTTGCGCGCCGAGTTCGTCGAAGCCCCGGCTACCGAAGTGCTGGCGGgtg > 1:18165/1‑78 (MQ=255) ttGGCCAGGGCCTGGCCAAGGAACTGGATTTGCGCGCCGAGTTCGTCGAAGCCCCGGCTACCGAAGTGCTGGCGGgtg < 2:18165/78‑1 (MQ=255) tGGCCAGGCCCTGGCCAAGGAACTGGATTTGCGCGCCGAGTTCGTCGAAGCCCCGGCTACCGAAGTGCTGGCGGGTGTCGAGAGTGGCAAGTACGACATCGCATTCACCCCGTCGAGTGAATCGCTCACAGGCGACGGCCCGCTCGATGt < 2:1525715/150‑1 (MQ=255) cAGGCCCTGGCCAAGGAACTGGATTTGCGCGCCGAGTTCGTCGAAGCCCCGGCTACCGAAGTGCTGGCGGGTGTc < 2:302839/75‑1 (MQ=255) cAGGCCCTGGCCAAGGAACTGGATTTGCGCGCCGAGTTCGTCGAAGCCCCGGCTACCGAAGTGCTGGCGGGTGTc > 1:302839/1‑75 (MQ=255) cAGGCCCTGGCCAAGGAACTGGATTTGCGCGCCGAGTTCGTCGAAGCCCCGGCTACCGAAGTGCTGGCGGGTGTCGAGAGTGGCAAGTACGACATCGCATTCACCCCGTCGAGTGAATCGCTCACAGGCGACGGCCCGCTCGATGTAAGc > 1:1215658/1‑150 (MQ=255) cccTGGCCAAGGAACTGGATTTGCGCGCCGGGTTCGTCGAAGCCCCGGCTACCGAAGTGCTGGCGGGTCAGGAGAGTGGCAAGTACGACATCGCATTCACCCCGTCGAGTGAATCGCTCACAGGCGACGGCCCGCTCGATGTAAGCCTGc < 1:1103643/150‑1 (MQ=255) cccTGGCCAAGGAACTGGATTTGCGCGCCGAGTTCGTCGAAGCCCCGGCTACCGAAGTGCTGGCGGGTGTCGAGAGTGGCAAGTACGACATCGCATTCCCCCCGCCGATTGAATCGCTCACAGGCGACGGCCCGCTCGATGTAAGCCTGc < 1:284597/150‑1 (MQ=255) | GTACTGATCAGTCCATTGGCTCAGGCGGAGCTGATCGACGAAATCGCCGACCGGGGCGAACTGCGCATTGCCGTACAGGCCGAGAACTCGCCATACGCCTTCAAGCAAGACGACCACCTGACAGGTTTCGACATCGAATTTGGCCAGGCCCTGGCCAAGGAACTGGATTTGCGCGCCGAGTTCGTCGAAGCCCCGGCTACCGAAGTGCTGGCGGGTGTCGAGAGTGGCAAGTACGACATCGCATTCACCCCGTCGAGTGAATCGCTCACAGGCGACGGCCCGCTCGATGTAAGCCTGC > NC_002947/5883503‑5883800 |
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: ‑ |