Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | NC_002947 | 409,081 | C→T | 70.6% | G529D (GGC→GAC) | aceE ← | pyruvate dehydrogenase E1 component |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_002947 | 409,081 | 0 | C | T | 70.6% | 65.7 / 33.4 | 51 | G529D (GGC→GAC) | aceE | pyruvate dehydrogenase E1 component |
Reads supporting (aligned to +/- strand): ref base C (10/5); new base T (20/16); total (30/21) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 5.43e-01 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.63e-01 |
TACGAGGTACCGGCAGCGATGAACGACGACATGGCGCCGGCTTCGTTGATGCCTTCCTCGAGGATCTGGCCCTTCTTGTCTTCGCGGTAGAACATCACCTGGTCTTTATCGACTGGCTCGTAGAGCTGGCCGACCGACGAGTAGATGCCCAACTGGCGGAACATGCCTTCCATACCGAAGGTACGGGCTTCGTCCGGGATGATCGGGACGATGCGCTGGCCAATGTCCTTGTCCTTGACCAGCTGCGCCAGAATACGCACGAAGGCCATGGTGGTGGAGATTTCGCGGTCGCCCGA > NC_002947/408934‑409229 | tACGAGGTACCGGCAGCGATGAACGACGACATGGCGCCGGCTTCGTTGATGCCTTCCTCGAGGATCTGGCCCTTCTTGTCTTCGCGGTAGAACATCACCTGGTCTTTATCGACTGGCTCGTAGAGCTGGCCGACCGACGAGTAGATGccc < 2:205104/150‑1 (MQ=255) aGGTACCGGCAGCGATGAACGACGACATGGCGCCGGCTTCGTTGATGCCTTCCTCGAGGATCTGGCCCTTCTTGTCTTCGCGGTAGAACATCACCTGGTCTTTATCGACTGGCTCGTAGAGCTGGCCGACCGACGAGTAGATGTCCAACt > 1:762681/1‑150 (MQ=255) gccGGCAGCGATGAACGACGACATGGCGCCGGCTTCGTTGATGCCTTCCTCGAGGATCTGGCCCTTCTTGTCTTCGCGGTAGAACATCACCTGGTCTTTATCGACTGGCTCGTAGAGCTGGCCGACCGACGAGTAGATGTCCAACTGGCg < 1:431081/149‑1 (MQ=255) gATGAACGACGACATGGCGCCGGCTTCGTTGATGCCTTCCTCGAGGATCTGGCCCTTCTTGTCTTCGCGGTAGAACATCACCTGGTCTTTATCGACTGGCTCGTAGAGCTGGCCGACCGACGAGTAGATGTCCAACTGGCGGAACATGcc > 1:274323/1‑150 (MQ=255) gAACGACGACATGGCGCCGGCTTCGTTGATGCCTTCCTCGAGGATCTGGCCCTTCTTGTCTTCGCGGTAGAACATCACCTGGTCTTTATCGACTGGCTCGTAGAGCTGGCCGACCGACGAGTAGATGTCCAACTGGCGGAACATGCCtt > 2:313198/1‑149 (MQ=255) gAACGACGACATGGCGCCGGCTTCGTTGATGCCTTCCTCGAGGATCTGGCCCTTCTTGTCTTCGCGGTAGAACATCACCTGGTCTTTATCGACTGGCTCGTAGAGCTGGCCGACCGACGAGTAGATGTCCAACTGGCGGAACATGCCtt < 1:313198/149‑1 (MQ=255) gAACGACGACATGGCGCCGGCTTCGTTGATGCCTTCCTCGAGGATCTGGCCCTTCTTGTCTTCGCGGTAGAACATCACCTGGTCTTTATCGACTGGCTCGTAGAGCTGGCCGACCGACGAGTAGATGCCCAACTGGCGGAACATGCCTTc > 2:1064158/1‑150 (MQ=255) cATGGCGCCGGCTTCGTTGATGCCTTCCTCGAGGATCTGGCCCTTCTTGTCTTCGCGGTAGAACATCACCTGGTCTTTATCGACTGGCTCGTAGAGCTGGCCGACCGACGAGTAGATGTCCAACTGGCGGAACATGCCTTCCATACCGaa < 1:963960/150‑1 (MQ=255) gcCGGCTTCGTTGATGCCTTCCTCGAGGATCTGGCCCTTCTTGTCTTCGCGGTAGAACATCACCTGGTCTTTATCGACTGGCTCGTAGAGCTGGCCGACCGACGAGTAGATGCCCAACTGGCGGAAc > 1:306734/1‑127 (MQ=255) gcCGGCTTCGTTGATGCCTTCCTCGAGGATCTGGCCCTTCTTGTCTTCGCGGTAGAACATCACCTGGTCTTTATCGACTGGCTCGTAGAGCTGGCCGACCGACGAGTAGATGCCCAACTGGCGGAAc < 2:306734/127‑1 (MQ=255) ggCTTCGTTGATGCCTTCCTCGAGGATCTGGCCCGCCTTGCCTTCGCGGTAGAACATCACCTGGGCGTTATCGACTGGCTCGTAGGGCTGGCAGACCGAGAAGTAGATGTCCCAATGGCGGAACATTCCTTCCATACCGAAGGTACGGGc > 2:75845/1‑150 (MQ=255) ttCGTTGATGCCTTCCTCGAGGATCTGGCCCTTCTTGTCTTCGCGGTAGAACATCACCTGGTCTTTATCGACTGGCTCGTAAAGCTGGCCGACCGACGAGTAGATGCCCAACTGGCGGAACATGCCTTCCATACCGAAGGTACGGGCTTc > 2:1001636/1‑150 (MQ=255) ttGATGCCTTCCTCGAGGATCTGGCCCTTCTTGTCTTCGCGGTAGAACATCACCTGGTCTTTATCGACTGGCTCGTAGAGCTGGCCGACCGACGAGTAGATGTCCAACTGGCGGAACATGCCTTCCATACCGAAGGTACGGGCTTCGTcc > 1:411356/1‑150 (MQ=255) gATGCCTTCCTCGAGGATCTGGCCCGTCTTGTCTTCGCGGTAGATCATCACCTGGTCTTTATCGACTGGCTCGTAGAGCTGGCCGACCGACGAGTAGATGTCCAACTGGCGGAACATGCCTTCCATACCGAAGGTACGGGCTTCGTCCgg < 1:75845/150‑1 (MQ=255) tGCCTTCCTCGAGGATCTGGCCCTTCTTGTCTTCGCGGTAGAACATCACCTGGTCTTTATCGACTGGCTCGTAGAGCTGGCCGACCGACGAGTAGATGTCCAACTGGCGGAACATGCCTTCCATACCGAAGGTACGGGCTTCGTCCGGga > 2:821376/1‑150 (MQ=255) tGCCTTCCTCGAGGATCTGGCCCTTCTTGTCTTCGCGGTAGAACATCACCTGGTCTTTATCGACTGGCTCGTAGAGCTGGCCGACCGACGAGTAGATGCCCAACTGGCGGAACATGCCTTCCATACCGAAGGTACGGGCTTCGTCCGGga > 1:756146/1‑150 (MQ=255) cTTCCTCGAGGATCTGGCCCTTCTTGTCTTCGCGGTAGAACATCACCTGGTCTTTATCGACTGGCTCGTAGAGCTGGCCGACCGACGAGTAGATGCCCAACTGGCGGAACATGCCTTCCATACCGAAGGTACGGGCTTCGTCCGGgatga > 2:1028294/1‑150 (MQ=255) ttCCTCGAGGATCTGGCCCTTCTTGTCTTCGCGGTAGAACATCACCTGGTCTTTATCGACTGGCTCGTAGAGCTGGCCGACCGACGAGTAGATGTCCAACTGGCGGAACATGCCTTCCATACCGAAGGTACGGGCTTCGTCCGGgatgat > 2:559606/1‑150 (MQ=255) ccTCGAGGATCTGGCCCTTCTTGTCTTCGCGGTAGAACATCACCTGGTCTTTATCGACTGGCTCGTAGAGCTGGCCGACCGACGAGTAGATGTCCAACTGGCGGAACATGCCTTCCATACCGAAGGTACGGGCTTCGTCCGGGATGATCg < 1:967161/150‑1 (MQ=255) tCGAGGATCTGGCCCTTCTTGTCTTCGCGGTAGAACATCACCTGGTCTTTATCGACTGGCTCGTAGAGCTGGCCGACCGACGAGTAGATGTCCAACTGGCGGAACATGCCTTCCATACCGAAGGTACGGGCTTCGTCCGGGATGATCggg > 2:318032/1‑150 (MQ=255) tCGAGGATCTGGCCCTTCTTGTCTTCGCGGTAGAACATCACCTGGTCTTTATCGACTGGCTCGTAGAGCTGGCCGACCGACGAGTAGATGCCCAACTGGCGGAACATGCCTTCCATACCGAAGGTACGGGCTTCGTCCGGGATGATCggg > 2:1010266/1‑150 (MQ=255) cGAGGATCTGGCCCTTCTTGTCTTCGCGGTAGAACATCACCTGGTCTTTATCGACTGGCTCGTAGAGCTGGCCGACCGACGAGTAGATGTCCAACTGGCGGAACATGCCTTCCATACCGAAGGTACGGGCTTCGTCCGGGATGATCGGGa < 2:762681/150‑1 (MQ=255) gAGGATCTGGCCCTTCTTGTCTTCGCGGTAGAACATCACCTGGTCTTTATCGACTGGCTCGTAGAGCTGGCCGACCGACGAGTAGATGTCCAACTGGCGGAACATGCCTTCCATACCGAAGGTACGGGCTTCGTCCGGGATGATCGGGAc > 1:794465/1‑150 (MQ=255) gAGGATCTGGCCCTTCTTGTCTTCGCGGTAGAACATCACCTGGTCTTTATCGACTGGCTCGTAGAGCTGGCCGACCGACGAGTAGATGCCCAACTGGCGGAACATGCCTTCCATACCGAAGGTACGGGCTTCGTCCGGGATGATCGGGa < 2:77943/149‑1 (MQ=255) gAGGATCTGGCCCTTCTTGTCTTCGCGGTAGAACATCACCTGGTCTTTATCGACTGGCTCGTAGAGCTGGCCGACCGACGAGTAGATGCCCAACTGGCGGAACATGCCTTCCATACCGAAGGTACGGGCTTCGTCCGGGATGATCGGGa > 1:77943/1‑149 (MQ=255) ggCCCTTCTTGTCTTCGCGGTAGAACATCACCTGGTCTTTATCGACTGGCTCGTAGAGCTGGCCGACCGACGAGTAGATGTCCAACTGGCGGAACATGCCTTCCATACCGAAGGTACGGGCTTCGTCCGGGATGATCGGGACGATGCGCt < 1:334751/150‑1 (MQ=255) ggCCCTTCTTGTCTTCGCGGTAGAACATCACCTGGTCTTTATCGACTGGCTCGTAGAGCTGGCCAACCGACGAGTAGATGTCCAACTGGCGGAACATGCCTTCCATACCGAAGGTACGGGCTTCGTCCGGGATGATCGGGACGATGCGCt > 2:66228/1‑150 (MQ=255) ccTTCTTGTCTTCGCGGTAGAACATCACCTGGTCTTTATCGACTGGCTCGTAGAGCTGGCCGACCGACGAGTAGATGTCCAACTGGCGGAACAt > 1:379297/1‑94 (MQ=255) ccTTCTTGTCTTCGCGGTAGAACATCACCTGGTCTTTATCGACTGGCTCGTAGAGCTGGCCGACCGACGAGTAGATGTCCAACTGGCGGAACAt < 2:379297/94‑1 (MQ=255) tCTTCGCGGTAGAACATCACCTGGTCTTTATCGACTGGCTCGTAGAGCTGGCCGACCGACGAGTAGATGCCCAACTGGCGGAACATGCCTTCCATACCGAAGGTACGGGCTTCGTCCGGGATGATCGGGACGATGCGCTGGCCAAtgtcc < 1:1028294/150‑1 (MQ=255) cATCACCTGGTCTTTATCGACTGGCTCGTAGAGCTGGCCGACCGACGAGTAGATGTCCAACTGGCGGAACATGCCTTCCATACCGAAGGTACGGGCTTCGTCCGGGATGATCGGGACGATGCGCTGGCCAAtgtccttgtc > 1:535001/1‑141 (MQ=255) cATCACCTGGTCTTTATCGACTGGCTCGTAGAGCTGGCCGACCGACGAGTAGATGTCCAACTGGCGGAACATGCCTTCCATACCGAAGGTACGGGCTTCGTCCGGGATGATCGGGACGATGCGCTGGCCAAtgtccttgtc < 2:535001/141‑1 (MQ=255) cATCACCTGGTCTTTATCGACTGGCTCGTAGAGCTGGCCGACCGACGAGTAGATGTCCAACTGGCGGAACATGCCTTCCATACCGAAGGTACGGGCTTCGTCCGGGATGATCGGGACGATGCGCTGGCCAATGTCCTTGTCCTTGACCAg < 1:559606/150‑1 (MQ=255) cACCTGGTCTTTATCGACTGGCTCGTAGAGCTGGCCGACCGACGAGTAGATGTCCAACTGGCGGAACATGCCTTCCATACCGAAGGTACGGGCTTCGTCCGGGATGATCGGGACGATGCGCTGGc < 2:752160/125‑1 (MQ=255) cACCTGGTCTTTATCGACTGGCTCGTAGAGCTGGCCGACCGACGAGTAGATGTCCAACTGGCGGAACATGCCTTCCATACCGAAGGTACGGGCTTCGTCCGGGATGATCGGGACGATGCGCTGGc > 1:752160/1‑125 (MQ=255) aCCTGGTCTTTATCGACTGGCTCGTAGAGCTGGCCGACCGACGAGTAGATGTCCAACTGGCGGAACATGCCTTCCATACCGAAGGTACGGGCTTCGTCCGGGATGATCGGGACGATGCGCTGGCCAATGTCCTTGTCCTTGACCAGCTgc > 1:10405/1‑150 (MQ=255) aTCGACTGGCTCGTAGAGCTGGCCGACCGACGAGTAGATGTCCAACTGGCGGAACATGCCTTCCATACCGAAGGTACGGGCTTCGTCCGGGATGATCGGGACGATGCGCTGGCCAATGTCCTTGTCCTTGACCAGCTg > 1:417639/1‑138 (MQ=255) aTCGACTGGCTCGTAGAGCTGGCCGACCGACGAGTAGATGTCCAACTGGCGGAACATGCCTTCCATACCGAAGGTACGGGCTTCGTCCGGGATGATCGGGACGATGCGCTGGCCAATGTCCTTGTCCTTGACCAGCTg < 2:417639/138‑1 (MQ=255) aTCGACTGGCTCGTAGAGCTGGCCGACCGACGAGTAGATGCCCAACTGGCGGAACATGCCTTCCATACCGAAGGTACGGGCTTCGTCCGGGATGATCGGGACGATGCGCTGGCCAAt > 1:229273/1‑117 (MQ=255) aTCGACTGGCTCGTAGAGCTGGCCGACCGACGAGTAGATGCCCAACTGGCGGAACATGCCTTCCATACCGAAGGTACGGGCTTCGTCCGGGATGATCGGGACGATGCGCTGGCCAAt < 2:229273/117‑1 (MQ=255) gACTGGCTCGTAGAGCTGGCCGACCGACGAGTAGATGTCCAACTGGCGGAACATGCCTTCCATACAGAAGGTACGGGCTTCGTCCGGGATGATCGGGACGATGCGCTGGCCAATGTCCTTGTCCTTGACCAGCTGCGCCAGAATACGCAc > 2:436404/1‑150 (MQ=255) aCTGGCTCGTAGAGCTGGCCGACCGACGAGTAGATGCCCAACTGGCGGAACATGCCTTCCATACCGAAGGTACGGGCTTCGTCCGGgatgat < 1:747985/92‑1 (MQ=255) aCTGGCTCGTAGAGCTGGCCGACCGACGAGTAGATGCCCAACTGGCGGAACATGCCTTCCATACCGAAGGTACGGGCTTCGTCCGGgatgat > 2:747985/1‑92 (MQ=255) tGGCTCGTAGAGCTGGCCGACCGACGAGTAGATGTCCAACTGGCGGAACATGCCTTCCATACCGAAGGTACGGGCTTCGTCCGGGATGATCGGGACGATGCGCTGGCCAATGTCCTTGTCCTTGACCAGCTGCGCCAGAATACGCACGaa > 2:777922/1‑150 (MQ=255) ggCTCGTAGAGCTGGCCGACCGACGAGTAGATGTCCAACTGGCGGAACATGCCTTCCATACCGAAGGTACGGGCTTCGTCCGGGATGATCGGGACGATGCGCTGGCCAATGTCCTTGTCCTTGACCAGCTGCGCCAGAATACGCACGAAg < 1:777922/150‑1 (MQ=255) tGGCCGACCGACGAGTAGATGTCCAACTGGCGGAACATGCCTTCCATACCGAAGGTACGGGCTTCGTCCGGGATGATCGGGACGATGCGCTGGCCAATGTCCTTGTCCTTGACCAGCTGCGCCAGAATACGCACGAAGGCCAtggtggtg < 1:821376/150‑1 (MQ=255) gCCGACCGACGAGTAGATGTCCAACTGGCGGAACATGCCTTCCATACCGAAGGTACGGGCTTCGTCCGGGATGATCGGGACGATGCGCTGGCCAATGTCCTTGTCCTTGACCAGCTGCGCCAGAATACGCACGAAGGCCATGGTGGTGga > 1:738660/1‑150 (MQ=255) gaccgacGAGTAGATGTCCAACTGGCGGAACATGCCTTCCATACCGAAGGTACGGGCTTCGTCCGGGATGATCGAGACGACGCGCTGGCCAATGTCCTTGTCCTTGACCAGCTGCGCCAGAATACGCACGAAGGCCATGGTGGTGGAGAt > 1:81033/1‑150 (MQ=255) ccgacGAGTAGATGTCCAACTGGCGGAACATGCCTTCCATACCGAAGGTACGGGCTTCGTCCGGGATGATCGGGACGATGCGCTGGCCAATGTCCTTGTCCTTGACCAGCTGCGCCAGAATACGCACGAAGGCCATGGTGGTGGAGAttt < 2:10405/150‑1 (MQ=255) cgacgaGTAGATGTCCAACTGGCGGAACATGCCTTCCATACCGAAGGTACGGGCTTCGTCCGGGATGATCGGGACGATGCGCTGGCCAATGTCCTTGTCCTAGACCAGCTGCGCCAGAATACGCACGAAGGCCATGGTGGTGGAGATTTc > 1:718246/1‑150 (MQ=255) cctaGTAGATGTCCAACTGGCGGAACATGCCTTCCATACCGAAGGTACGGGCTTCGTCCGGGATGATCGGGACGATGCGCTGGCCAATGTCCTTGTCCTTGACCAGCTGCGCCAGAATACGCACGAAGGCCATGGTGGTGGAGATTTcgc < 1:964932/147‑1 (MQ=255) gCCCAACTGGCGGAACATGCCTTCCATACCGAAGGTACGGGCTTCGTCCGGGATGATCGGGACGATGCGCTGGCCAATGTCCTTGTCCTTGACCAGCTGCGCCAGAATACGCACGAAGGCCATGGTGGTGGAGATTTCGCGGTCGCCCGa > 1:760178/1‑150 (MQ=255) | TACGAGGTACCGGCAGCGATGAACGACGACATGGCGCCGGCTTCGTTGATGCCTTCCTCGAGGATCTGGCCCTTCTTGTCTTCGCGGTAGAACATCACCTGGTCTTTATCGACTGGCTCGTAGAGCTGGCCGACCGACGAGTAGATGCCCAACTGGCGGAACATGCCTTCCATACCGAAGGTACGGGCTTCGTCCGGGATGATCGGGACGATGCGCTGGCCAATGTCCTTGTCCTTGACCAGCTGCGCCAGAATACGCACGAAGGCCATGGTGGTGGAGATTTCGCGGTCGCCCGA > NC_002947/408934‑409229 |
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: ‑ |