Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | NC_002947 | 5,059,508:1 | +G | 11.9% | coding (1326/1581 nt) | PP_4458 ← | opine ABC transporter substrate‑binding protein |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_002947 | 5,059,508 | 1 | . | G | 11.9% | 206.5 / 25.6 | 59 | coding (1326/1581 nt) | PP_4458 | opine ABC transporter substrate‑binding protein |
Reads supporting (aligned to +/- strand): ref base . (25/27); new base G (2/5); total (27/32) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 4.37e-01 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.32e-01 |
TGTTTCAGACCAAACAAAGGGATGGCGCATACGTCGTCGTGAATCTGAACCTGCGCCTTTTTCCACAAATCGAGTTGAGTCTGCGGATCAGCCTCGACCCTGGCTTTCCGGATGGCATCATCAGCTACTGAGCAATGTCCGAAATTC‑GACATTGCAGCGGGTGCGCCGATGGCAGAAGCCGAATCGTAAAACTCTGTCAGCCAGTAGTCGGCGCCCGGATATCGAGCAGCACCGTAAAACACGATCGCGCTCTGGTCTTGACGGCTTTTGGCCTGGTAAGTCGCGTG > NC_002947/5059362‑5059648 | tGTTTCAGACCAAACAAAGGGATGGCGCATACGTCGTCGTGAATCTGAACCTGCGCCTTTTTCCACAAATCGAGTTGAGTCTGCGGATCAGCCTCGACCCTGGCTTTCCGGATGGCATCATCAGCTACTGAGCAATGTCCGAAATTC‑Ga > 5:34384/1‑149 (MQ=255) tttCAGACCAAACAAAGGGATGGCGCATACGTCGTCGTGAATCTGAACCTGCGCCTTTTTCCACAAATCGAGTTGAGTCTGCGGATCAGCCTCGACCCTGGCTTTCCGGATGGCATCATCAGCTACTGAGCAATGTCCGAAATTC‑Ga < 7:279494/147‑1 (MQ=255) tttCAGACCAAACAAAGGGATGGCGCATACGTCGTCGTGAATCTGAACCTGCGCCTTTTTCCACAAATCGAGTTGAGTCTGCGGATCAGCCTCGACCCTGGCTTTCCGGATGGCATCATCAGCTACTGAGCAATGTCCGAAATTC‑GACa > 3:87758/1‑149 (MQ=255) gACCAAACAAAGGGATGGCGCATACGTCGTCGTGAATCTGAACCTGCGCCTTTTTCCACAAATCGAGTTGAGTCTGCGGATCAGCCTCGACCCTGGCTTTCCGGATGGCATCATCAGCTACTGAGCAATGT‑CGAATTTCGGACATTGct < 5:138912/149‑2 (MQ=255) gACCAAACAAAGGGATGGCGCATACGTCGTCGTGAATCTGAACCTGCGCCTTTTTCCACAAATCGAGTTGAGTCTGCGGATCAGCCTCGACCCTGGCTTTCCGGATGGCATCATCAGCTACTGAGCAATGT‑CGAATTTCGGACATTGct < 3:321043/149‑2 (MQ=255) acaaaGGGATGGCGCATACGTCGTCGTGAATCTGAACCTGCGCCTTTTTCCACAAATCGAGTTGAGTCTGCGGATCAGCCTCGACCCTGGCTTTCCGGATGGCATCATCAGCTACTGAGCAATGT‑CGAATTTCGGACATTGctcagtag < 2:32114/149‑8 (MQ=255) aaaGGGATGGCGCATATGTCGTCGTGAATCTGAACCTGCGCCTTTTTCCGCAAATCGAGTTGAGTCTGCGGATCAGCCTCGACCCTGGCTTTCCGGATGGCATCATCAGCTACTGAGCAATGTCCGAAATTC‑GACATTGAAGCGGGt > 5:261563/1‑147 (MQ=255) aaaGGGATGGCGCATACGTCGTCGTGAATCTGAACCTGCGCCTTTTTCCACAAATCGAGTTGAGTCTGCGGATAAGCCTCGACCCTGGATTTCCGGGTGGCATCATCAGATACTGAGCAATGTGCGAATTTG‑GACATTGCAGCGGGTgc > 2:327681/1‑149 (MQ=255) aGGGATGGCGCATACGTCGTCGTGAATCTGAACCTGCGCCTTTTTCCACAAATCGAGTTGAGTCTGCGGATCAGCCTCGACCCTGGCTTTCCGGATGGCATCATCAGCTACTGAGCAATGTCCGAAATTC‑GACATTGCAGCGGGTgcgc < 4:87758/149‑1 (MQ=255) gggATGGCGCATACGTCGTCGTGAATCTGAACCTGCGCCTTTTTCCACAAATCGAGTTGAGTCTGCGGATCAGCCTCGACCCTGGCTTTCCGGATGGCATCATCAGCTACTGAGCAATGTCCGAAATTC‑GACATTGCAGCgg > 1:276962/1‑142 (MQ=255) gggATGGCGCATACGTCGTCGTGAATCTGAACCTGCGCCTTTTTCCACAAATCGAGTTGAGTCTGCGGATCAGCCTCGACCCTGGCTTTCCGGATGGCATCATCAGCTACTGAGCAATGTCCGAAATTC‑GACATTGCAGCgg < 2:276962/142‑1 (MQ=255) taTGGCGCATACGTCGTCGTGAATCTGAACCTGCGCCTTTTTCCACAAATCGAGTTGAGTCTGCGGATCAGCCTCGACCCTGGCTTTCCGGATGGCATCATCAGCTACTGAGCAATGTCCGAAATTC‑GACATTGCAGCGGGTGCGCCGa < 7:302584/148‑1 (MQ=255) tGGCGCATACGTCGTCGTGAATCTGAACCTGCGCCTTTTTCCACAAATCGAGTTGAGTCTGCGGATCAGCCTCGACCCTGGCTTTCCGGATGGCATCATCAGCTACTGAGCAATGTCCGAAATTC‑GACATTGCAGCGGGTGCGCCGATg < 1:28360/149‑1 (MQ=255) cgcATACGTCGTCGTGAATCTGAACCTGCGCCTTTTTCCACAAATCGAGTTGAGTCTGCGGATCAGCCTCGACCCTGGCTTTCCGGATGGCATCATCAGCTACTGAGCAATGTCCGAAATTC‑GACATTGCAGCGGGTGCGCCGATGGCa < 7:251710/149‑1 (MQ=255) tcgtcgtGAATCTGAACCTGCGCCTTTTTCCACAAATCGAGTTGAGTCTGCGGATCAGCCTCGACCCTGGCTTTCCGGATGGCATCATCAGCTACTGAGCAATGTCCGAAATTC‑GACATTGCAGCGGGTGCGCCGATGGCAGAAGCCg > 7:54443/1‑148 (MQ=255) tGAATCTGAACCTGCGCCTTTTTCCACAAATCGAGTTGAGTCTGCGGATCAGCCTCGACCCTGGCTTTCCGGATGGCATCATCAGCTACTGAGCAATGTCCGAAATTC‑GACATTGCAGCGGGTGCGCCGATGGCAGAAGCCGAATCGTa < 8:54443/149‑1 (MQ=255) gcgcCTTTTTCCACAAATCGAGTTGAGTCTGCGGATCAGCCTCGACCCTGGCTTTCCGGATGGCATCATCAGCTACTGAGCAATGTCCGAAATTC‑GACATTGCAGCGGGTGCGCCGa < 5:344361/117‑1 (MQ=255) gcgcCTTTTTCCACAAATCGAGTTGAGTCTGCGGATCAGCCTCGACCCTGGCTTTCCGGATGGCATCATCAGCTACTGAGCAATGTCCGAAATTC‑GACATTGCAGCGGGTGCGCCGa > 6:344361/1‑117 (MQ=255) tttCCACAAATCGAGTTGAGTCTGCGGATCAGCCTCGACCCTGGCTTTCCGGATGGCATCATCAGCTACTGAGCAATGTCCGAAATTC‑GACATTGCAGCGGGTGCGCCGATGGCAGAAGCCGAATCGTAAAACTCTGTCAGCCagtag < 6:34384/148‑1 (MQ=255) aaaTCGAGTTGAGTCTGCGGATCAGCCTCGACCCTGGCTTTCCGGATGGCATCATCAGCTACTGAGCAATGTCCGAAATTC‑GACATTGCAGCGGGTGCGCCGATGGCAGAAGCCGAATCGTAAAACTCTGTCAGCCAGTAGTCGGCGcc > 6:89450/1‑149 (MQ=255) aGTTGAGTCTGCGGATCAGCCTCGACCCTGGCTTTCCGGATGGCATCATCAGCTACTGAGCAATGTCCGAAATTC‑GACATTGCAGCGGGTGCGCCGATGGCAGAAGCCGAATCGTAAAACTCTGTCAGCCAGTAGTCGGCGCCCGGata < 5:203039/149‑1 (MQ=255) gTTGAGTCTGCGGATCAGCCTCGACCCTGGCTTTCCGGATGGCATCATCAGCTACTGAGCAATGTCCGAAATTC‑GACATTGCAGCGGGTGCGCCGATGGCAGAAGCCGAATCGTAAAACTCTGTCAGCCAGTAGTCg > 6:346287/1‑137 (MQ=255) gTTGAGTCTGCGGATCAGCCTCGACCCTGGCTTTCCGGATGGCATCATCAGCTACTGAGCAATGTCCGAAATTC‑GACATTGCAGCGGGTGCGCCGATGGCAGAAGCCGAATCGTAAAACTCTGTCAGCCAGTAGTCg < 5:346287/137‑1 (MQ=255) gTCTGCGGATCAGACTCGACCCTGGCTTTACGGATGTCATCATCAGCTACTGAGCAATGTCAGAAATTC‑GACATTGCAGCGGGTGCGCCGATGGCAGAATCCGAATCGTAAAACTCTGTCAGCCAGTAGTCGGCTCCCGGATATCGatc < 1:327681/149‑3 (MQ=255) aGCCTCGACCCTGGCTTTCCGGATGGCATCATCAGCTACTGAGCAATGTCCGAAATTC‑GACATTGCAGCGGGTGCGCCGATGGCAGAAGCCGAATCGTAAAACTCTGTCAGCCAGTAGTCg > 6:361083/1‑121 (MQ=255) aGCCTCGACCCTGGCTTTCCGGATGGCATCATCAGCTACTGAGCAATGTCCGAAATTC‑GACATTGCAGCGGGTGCGCCGATGGCAGAAGCCGAATCGTAAAACTCTGTCAGCCAGTAGTCg < 5:361083/121‑1 (MQ=255) gCCTCGACCCTGGCTTTCCGGATGGCATCATCAGCTACTGAGCAATGTCCGAAATTC‑GACATTGCAGCGGGTGCGCCGATGGCAGAAGCCGAATCGTAAAACTCTGTCAGCCAGTAGTCGGCGCCCGGATATCGAGCAGCACCGTaaaa > 1:320404/1‑149 (MQ=255) gCCTCGACCCTGGCTTTCCGGATGGCATCATCAGCTACTGAGCAATGTCCGAAATTC‑GACATTGCAGCGGGTGCGCCGATGGCAGAAGCCGAATCGTAAAACTCTGTCAGCCAGTAGTCGGCGCCCGGATATCGAGCAGCACCGTaaaa > 5:70637/1‑149 (MQ=255) ccTCGACCCTGGCTTTCCGGATGGCATCATCAGCTACTGAGCAATGTCCGAAATTC‑GACATTGCAGCGGGTGCGCCGATGGCAGACGCCGAATCGTAAATCTCTGTCAGCCAGTAGTCGGCGCCCGGATATCGAGCAGCACCGTaaaa > 3:296095/1‑148 (MQ=255) ccTCGACCCTGGCTTTCCGGATGGCATCATCAGCTACTGAGCAATGTCCGAAATTC‑GACATTGCAGCGGGTGCGCCGATGGCAGAAGCCGAATCGTAAAACTCTGTCAGCCAGTAGTCGGCGCCCGGATATCGAGCAGCACCGTAAAac < 5:89450/149‑1 (MQ=255) cTCGACCCTGGCTTTCCGGATGGCATCATCAGCTACTGAGCAATGTCCGAAATTC‑GACATTGCAGCGGGTGCGCCGATGGCAGAAGCCGAATCGTAAAACTCTGTCAGCCAGTAGTCg < 6:166667/118‑1 (MQ=255) cTCGACCCTGGCTTTCCGGATGGCATCATCAGCTACTGAGCAATGTCCGAAATTC‑GACATTGCAGCGGGTGCGCCGATGGCAGAAGCCGAATCGTAAAACTCTGTCAGCCAGTAGTCGGCGCCCGGATATCGAGCAGCACCGTAAAaca > 6:70021/1‑149 (MQ=255) tCGACCCTGGCTTTCCGGATGGCATCATCAGCTACTGAGCAATGTCCGAAATTC‑GACATTGCAGCGGGTGCGCCGATGGCAGAAGCCGAATCGTAAAACTCTGTCAGCCAGTAGTCGGCGCCCGGATATCGAGCAGCACCGTAAAacac > 7:281478/1‑149 (MQ=255) cTGGCTTTCCGGATGGCATCATCAGCTACTGAGCAATGTCCGAAATTC‑GACATTGCAGCGGGTGCGCCGATGGCAGAAGCCGAATCGTAAAACTCTGTCAGCCAGTAGTCGGCGCCCGGATATCGAGCAGCACCGTAAAACACGATcgc < 8:281478/149‑1 (MQ=255) cTTTCCGGATGGCATCATCAGCTACTGAGCAATGTCCGAAATTC‑GACATTGCAGCGGGTGCGCCGATGGCAGAAGCCGAATCGTAAAACTCTGTCAGCCAGTAGTCGGCGCCCGGATATCGAGCAGCACCGTAAAACACGATCGCGctc > 7:105046/1‑149 (MQ=255) tCCGGATGGAATCATCAGATACTTACCAATGTCCTAAATTC‑GACATTGCAGCGGGTGCGCCTATGGCATAAGCCGAATCGTAAAACTCTGTCAGCCATTAGTCGGCGCCCGTATATCGAGCAGCACCGTAAAACACGATCGCGCTCTg < 6:70637/148‑1 (MQ=255) ccGGATGGCATCATCAGCTACTGAGCAATGTCCGAAATTC‑GACATTGCAGCGGGTGCGCCGATGGCAGAAGCTGAATCGTAAAACTCTGTCAGCCAGTAGTCGGCGCCCGGATATc < 3:299727/116‑1 (MQ=255) ccGGATGGCATCATCAGCTACTGAGCAATGTCCGAAATTC‑GACATTGCAGCGGGTGCGCCGATGGCAGAAGCTGAATCGTAAAACTCTGTCAGCCAGTAGTCGGCGCCCGGATATc > 4:299727/1‑116 (MQ=255) tGGCATCATCAGCTACTGAGCAATGTCCGAAATTC‑GACATTGCAGCGGGTGCGCCGATGGCAGAAGCCGAATCGTAAAACTCTGTCAGCCAGTAGTCGGCGCCCGGATATCGAGCAGCACCGTa < 7:292692/124‑1 (MQ=255) tGGCATCATCAGCTACTGAGCAATGTCCGAAATTC‑GACATTGCAGCGGGTGCGCCGATGGCAGAAGCCGAATCGTAAAACTCTGTCAGCCAGTAGTCGGCGCCCGGATATCGAGCAGCACCGTa > 8:292692/1‑124 (MQ=255) catcatCAGCTACTGAGCAATGTCCGAAATTC‑GACATTGCAGCGGGTGCGCCGATGGCAGAAGCCGAATCGTAAAACTCTGTCAGCCAGTAGTCGGCGCCCGGATATCGAGCAGCACCGTAAAACACGATCGCGCTCTGGTc < 7:275595/142‑1 (MQ=255) aatcatcAGCTACTGAGCAATGTCAGAAATTC‑GACATTGCAGCGGGTGCGCCGATGGCAGAAGACGAATCGTAAAACTCTGTCAGCCAGTAGTCGGCGCCCGGATATCTAGCAGAACCGTAAAACACGATCGCGCTCTggac > 8:275595/2‑140 (MQ=255) atcatcAGCTACTGAGCAATGTCCGAAATTC‑GACATTGCAGCGGGTGCGCCGATGGCAGAAGCCGAATCGTAAAACTCTGTCAGCCAGTAGTCGGCGCCCGGATATCGAGCAGCACCGTAAAACACGATCGCGCTCTGGTCTTGACGGc > 7:10052/1‑149 (MQ=255) cTGAGCAATGTCCGAAATTC‑GACATTGCAGCGGGTGCGCCGATGGCAGAAGCCGAATCGTAAAACTCTGTCAGCCAGTAGTCGGCGCCCGGATATCGAGCAGCACCGTAAAACACGATCGCGCTCTGGTCTTGACGGCTTTTGGCCTg < 8:105046/148‑1 (MQ=255) gAGCAATGTCCGAAATTC‑GACATTGCAGCGGGTGCGCCGATGGCAGAAGCCGAATCGTAAAACTCTGTCAGCCAGTAGTCGGCGCCCGGATATCGAGCAGCACCGTAAAACACGATCGCGCTCTGGTCTTGACGGCTTTTGGCCTGGTa > 5:220070/1‑149 (MQ=255) tgCAATGT‑CGAATTTCGGACATTGCAGCGGGTGCGCCGATGGCAGAAGCCGAAt > 3:128650/2‑54 (MQ=18) tgCAATGT‑CGAATTTCGGACATTGCAGCGGGTGCGCCGATGGCAGAAGCCGAAt < 4:128650/53‑1 (MQ=18) tgCAATGT‑CGAATTTCGGACATTGCAGCGGGTGCGCCGATGGCAGAAGCCGAATCGTAAAACTCTGTCAGCCAGTAGTCGGCGCCCGGATATCGAGCAGCACCGTAAAACACGATCGCGCTCTGGTCTTGACGGCTTTTg > 6:196205/2‑140 (MQ=255) tgCAATGT‑CGAATTTCGGACATTGCAGCGGGTGCGCCGATGGCAGAAGCCGAATCGTAAAACTCTGTCAGCCAGTAGTCGGCGCCCGGATATCGAGCAGCACCGTAAAACACGATCGCGCTCTGGTCTTGACGGCTTTTg < 5:196205/139‑1 (MQ=255) aGCAATGTCCGAAATTC‑GACATTGCAGCGGGTGCGCCGATGGCAGAAGCCGAATCGTAAAACTCTGTCAGCCAGTAGTCGGCGCCCGGATATCGAGCAGCACCGTAAAACACGATCGCGCTCTGGTCTTGACGGCTTTTGGCCTGGTaa < 6:220070/149‑1 (MQ=255) gCAATGTCCGAAATTC‑GACATTGCAGCGGGTGCGCCGATGGCAGAAGCCGAATCGTAAAACTCTGTCAGCCAGTAGTCGGCGCCCGGATATCGag < 6:339342/95‑1 (MQ=255) gCAATGTCCGAAATTC‑GACATTGCAGCGGGTGCGCCGATGGCAGAAGCCGAATCGTAAAACTCTGTCAGCCAGTAGTCGGCGCCCGGATATCGag > 5:339342/1‑95 (MQ=255) aaTGTCCGAAATTC‑GACATTGCAGCGGGTGCGCCGATGGCAGAAGCCGAATCGTAAAACTCTGTCAGCCAGTAGTCGGCGCCCGGATATCGAGCAGCACCGTAAAACACGATCGCGCTCTGGTCTTGACGGCTTTTGGCCTGGTaa > 6:16545/1‑146 (MQ=255) aaTGTCCGAAATTC‑GACATTGCAGCGGGTGCGCCGATGGCAGAAGCCGAATCGTAAAACTCTGTCAGCCAGTAGTCGGCGCCCGGATATCGAGCAGCACCGTAAAACACGATCGCGCTCTGGTCTTGACGGCTTTTGGCCTGGTaa < 5:16545/146‑1 (MQ=255) cGAAATTC‑GACATTGCAGCGGGTGCGCCGATGGCAGAAGCCGAATCGTAAAACTCTGTCAGCCAGTAGTCGGCGCCCGGATATCGAGCAGCACCGTAAAACACGATCGCGCTCTGGTCTTGACGGCTTTTGGCCTGGTAAGTCGCGTg < 3:169499/148‑1 (MQ=255) gAAATTC‑GACATTGCAGCGGGTGCGCCGATGGCAGAAGCCGAATCGTAAAACTCTGTCAGCCAGTAGTCGGCGCCCGGATATCGAGCAGCACCGTAAAACACGATCGCGCTCTGGTCTTGACGGCTTTTGGCCTg > 1:271029/1‑135 (MQ=255) gAAATTC‑GACATTGCAGCGGGTGCGCCGATGGCAGAAGCCGAATCGTAAAACTCTGTCAGCCAGTAGTCGGCGCCCGGATATCGAGCAGCACCGTAAAACACGATCGCGCTCTGGTCTTGACGGCTTTTGGCCTg < 2:271029/135‑1 (MQ=255) ttC‑GACATTGCAGCGGGTGCGCCGATGGCAGAAGCCGAATCGTAAAACTCTGTCAGCCAGTAGTCGGCGCCCGGATATCGAGCAGCACCGTAAAACACGATCGCGCTCTGGTCTTGACGGCTTTTGGCCTGGt > 3:268626/1‑133 (MQ=255) ttC‑GACATTGCAGCGGGTGCGCCGATGGCAGAAGCCGAATCGTAAAACTCTGTCAGCCAGTAGTCGGCGCCCGGATATCGAGCAGCACCGTAAAACACGATCGCGCTCTGGTCTTGACGGCTTTTGGCCTGGt < 4:268626/133‑1 (MQ=255) | TGTTTCAGACCAAACAAAGGGATGGCGCATACGTCGTCGTGAATCTGAACCTGCGCCTTTTTCCACAAATCGAGTTGAGTCTGCGGATCAGCCTCGACCCTGGCTTTCCGGATGGCATCATCAGCTACTGAGCAATGTCCGAAATTC‑GACATTGCAGCGGGTGCGCCGATGGCAGAAGCCGAATCGTAAAACTCTGTCAGCCAGTAGTCGGCGCCCGGATATCGAGCAGCACCGTAAAACACGATCGCGCTCTGGTCTTGACGGCTTTTGGCCTGGTAAGTCGCGTG > NC_002947/5059362‑5059648 |
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: ‑ |