Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | NC_000913 | 759,237 | A→G | 15.8% | T178A (ACC→GCC) | sucA → | 2‑oxoglutarate decarboxylase, thiamine triphosphate‑binding |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_000913 | 759,237 | 0 | A | G | 15.8% | 161.7 / 25.8 | 76 | T178A (ACC→GCC) | sucA | 2‑oxoglutarate decarboxylase, thiamine triphosphate‑binding |
Reads supporting (aligned to +/- strand): ref base A (35/29); new base G (6/6); total (41/35) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 1.00e+00 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.97e-02 |
GCCAGCGGCAAAGAAACCATGAAACTCGGCGAGCTGCTGGAAGCCCTCAAGCAAACCTACTGCGGCCCGATTGGTGCCGAGTATATGCACATTACCAGCACCGAAGAAAAACGCTGGATCCAACAGCGTATCGAGTCTGGTCGCGCGACTTTCAATAGCGAAGAGAAAAAACGCTTCTTAAGCGAACTGACCGCCGCT > NC_000913/759138‑759335 | gCCAGCGGCAAAGAAACCATGAAACTCGGCGAGCTGCCGGAAGCCCTCAAGCAAACCTACTGCGGCCCGATTGGTGCCGAGTATATGCACATTACCAGCa < 2:1614452/100‑1 (MQ=255) cAGCGGCAAAGAAACCATGAATCTCGGCGAGCTGCTGGAAGCCCTCAAGCAAACCTACTGCGGCCCGATTGGTGCCGAGTATATGCACATTACCAGCAcc < 2:2067121/100‑1 (MQ=255) cAGCGGCAAAGAAACCATGAAACTCGGCGAGCTGCTGGAAGCCCTCAAGCAAACCTACTGCGGCCCGATTGGTGCCGAGTATATGCACATTACCAGCAcc > 2:778114/1‑100 (MQ=255) cAGCGGCAAAGAAACCATGAAACTCGGCGAGCTGCTGGAAGCCCTCAAGCAAACCTACTGCGGCCCGATTGGTGCCGAGTATATGCACATTACCAGCAcc < 2:713313/100‑1 (MQ=255) aGCGGCAAAGAAACCATGAAACTCGGCGAGCTGCTGGAAGCCCTCAAGCAAACCTACTGCGGCCCGATTGGTGCCGAGTATATGCACATTACCAGCGCCg > 2:273682/1‑100 (MQ=255) gCGGCAAAGAAACCATGAAACTCGGCGAGCTGCTGGAAGCCCTCAAGCAAACCTACTGCGGCCCGATTGGTGCCGAGTATATGCACATTACCAGCACCga < 1:2379028/100‑1 (MQ=255) cGGCAAAGAAACCATGAAACTCGGCGAGCTGCTGGAAGCCCTCAAGCAAACCTACTGCGGCCCGATTGGTGCCGAGTATATGCACATTACCAGCGCCgaa > 1:1760314/1‑100 (MQ=255) cGGCAAAGAAACCATGAAACTCGGCGAGCTGCTGGAAGCCCTCAAGCAAACCTACTGCGGCCCGATTGGTGCCGAGTATATGCACATTACCAGCACCgaa < 2:891425/100‑1 (MQ=255) ggCAAAGAAACCATGAAACTCGGCGAGCTGCTGGAAGCCCTCAAGCAAACCTACTGCGGCCCGATTGGTGCCGAGTATATGCACATTACCAGCACCgaag > 1:1027973/1‑100 (MQ=255) cAAAGAAACCATGAAACTCGGCGAGTTGCTGGAAGCCCTCAAGCAAACATACTGCGGCCCGATTGGTGCCGAGTATATGCACATTACCAGCACCgaagaa > 2:898189/1‑100 (MQ=255) aaaGAAACCATGAAACTCGGCGAGCTGCTGGAAGCCCACAAGCAAACCTACTGCGGCCCGATTGGTGCCGAGTATATGCACATTACCAGCACCGAAGaaa > 2:367333/1‑100 (MQ=255) aaaCCATGAAACTCGGCGAGCTGCTGGAAGCCCTCAAGCAAACCTACTGCGGCCCGATTGGTGCCGAGTATATGCACATTACCAGCACCGAAGAAAAACg < 2:390927/100‑1 (MQ=255) ccATGAAACTCGGCGAGCTGCTGGAAGCCCTCAAGCAAACCTACTGCGGCCCGATTGGTGCCGAGTATATGCACATTACCAGCGCCGAAGAAAAACGCTg < 2:767629/100‑1 (MQ=255) cATGAAACTCGGCGAGCTGCTGGAAGCCCTCAAGCAAACCTACTGCGGCCCGATTGGTGCCGAGTATATGCACATTACCAGCACCGAAGAAAAACGCTgg > 1:416026/1‑100 (MQ=255) cATGAAACTCGGCGAGCTGCTGGAAGCCCTCAAGCAAACCTACTGCGGCCCGATTGGTGCCGAGTATATGCACATTACCAGCACCGAAGAAAAACGCTgg > 2:880844/1‑100 (MQ=255) tGAAACTCGGCGGGCTGCTGGAAGCCCTCAAGCAAACCTACTGCGGCCCGATTGGTGCCGAGTATATGCACATTACCAGCACCGAAGAAAAACGCTGGAt < 2:158323/100‑1 (MQ=255) gAAACTCGGCGAGCTGCTGGAAGCCCTCAAGCAAACCTACTGCGGCCCGATTGGTGCCGAGTATATGCACATTACCAGCACCGAAGAAAAACGCTGGATc < 1:523579/100‑1 (MQ=255) cGGCGAGCTGCTGGAAGCCCTCAAGCAAACCTACTGCGGCCCGCTTGGTGCCGAGTATATGCACATTACCAGCGCCGAAGAAAAACGCTGGATCCAACAg > 2:2176938/1‑100 (MQ=255) cGAGCTGCTGGAAGCACTCAAGCAACCCCACAGCGGCCCGATTGATGCCGAGTATATGCACATTACCAGCACCGAAGAAACACGCTAGATCCAACAGCGt < 2:705312/100‑1 (MQ=255) gctgctGGAAGCCCTCAAGCAAACCTACTGCGGCCCGATTGGTGCCGAGTATATGCACATTACCAGCACCGAAGAAAAACGCTGGATCCAACAGCGTATc < 2:404186/100‑1 (MQ=255) ctgctgGAAGCCCTCAAGCAAACCTACTGCGGCCCGATTGGTGCCGAGTATATGCACATTACCAGCGCCGAAGAAAAACGCTGGATCCAACAGCGTATCg < 1:1892713/100‑1 (MQ=255) ctgctgGAAGCCCTCAAGCAAACCTACTGCGGCCCGATTAGTGCCGAGTATATGCACATTACCAGCACCGAAGAAAAACGCTGGATCCAACAGCGTACCg > 2:160347/1‑100 (MQ=255) tgctgGAAGCCCTCAAGCAAACCTACTGCGGCCCGATTGGTGCCGAGTATATGCACATTACCAGCGCCGAAGAAAAACGCTGGATCCAACAGCGTATCGa < 1:1386151/100‑1 (MQ=255) tgctgGAAGCCCTCAAGCAAACCTACTGCGGCCCGATTGGTGCCGAGTATATGCACATTACCAGCACCGAAGAAAAACGCTGGATCCAACAGCGTATCGa > 1:2155788/1‑100 (MQ=255) gctgGAAGCCCTCAAGCAAACCTACTGCGGCCCGATTGGTGCCGAGTATATGCACATTACCAGCACCGAAGAAAAACGCTGGATCCAACAGCGTATCGAg > 2:2337607/1‑100 (MQ=255) gctgGAAGCCCTCAAGCAAACCTACTGCGGCCCGATTGGTGCCGAGTATATGCACATTACCAGCACCGAAGAAAAACGCTGGATCCAACAGCGTATCGAg > 2:1025942/1‑100 (MQ=255) gctgGAAGCCCTCAAGCAAACCTACTGCGGCCCGATTGGTGCCGAGTATATGCACATTACCAGCACCGAAGAAAAACGCTGGATCCAACAGCGTATCGAg > 1:865717/1‑100 (MQ=255) gctgGAAGCCCTCAAGCAAACCTACTGCGGCCCGATAGGTGCCGAGTATATGCACATTACCAGCACCGAAGAAAAACGCTGGATCCAACAGCGTATCGAg > 1:1689908/1‑100 (MQ=255) ctgGAAGCCCTCAAGCAAACCTACTGCGGCCCGATTGGTGCCGAGTATATGCACATTACCAGCACCGAAGAAAAACGCTGGATCCAACAGCGTATCGAGt < 1:1185276/100‑1 (MQ=255) ggAAGCCCTCAAGCAAACCTACTGCGGCCCGATTGGTGCCGAGTATATGCACATTACCAGCACCGAAGAAAAACGCTGGATCCAACAGCGTATCGAGTCt < 2:1013812/100‑1 (MQ=255) aGCCCTCAAGCAAACCTACTGCGGCCCGATTGGTGCCGAGTATATGCACATTACCAGCACCGAAGAAAAACGCTGGATCCAACAGCGTATCGAGTCTGGt < 1:2223780/100‑1 (MQ=255) gCCCTCAAGCAAACCTACTGCGGCCCGATTGGTGCCGAGTATATGCACATTACCAGCACCGAAGAAAAACGCTGGATCCAACAGCGTATCGAGTCTGGTc > 1:910284/1‑100 (MQ=255) ccTCAAGCAAACCTACTGCGGCCCGATTGGTGCCGAGTATATGCACATTACCAGCACCGCAGAAAAACGCTGGATCTAACAGCGCATCAAGTCTAGTagc > 1:760909/1‑100 (MQ=255) ccTCAAGCAAACCTACTGCGGCCCGATTGGTGCCGAGTATATGCACATTACCAGCACCGAAGAAAAACGCTGGATCCAACAGCGTATCGAGTCTGGTcgc > 1:617737/1‑100 (MQ=255) ccTCAAGCAAACCTACTGCGGCCCGATAACTGCCGCGTATAAGCAAACTCCCTGCACCGCAGAAAAACGCTGGAACCAACAGCTTAGCGCGTATGGt > 2:1790493/1‑97 (MQ=255) cTCAAGCAAACCTACTGCGGCCCGATTGGTGGCGAGTATATGCACATTACCAGCACCGAAGAAAAACGCTGGATCCAACAGCGTATCGAGTCTGGTcgcg < 1:1548777/100‑1 (MQ=255) cTCAAGCAAACCTACTGCGGCCCGATTGGTGCCGAGTATATGCACATTACCAGCGCCGAAGAAAAACGCTGGATCCAACAGCGTATCGAGTCTGGTcgcg > 1:1769307/1‑100 (MQ=255) cTCAAGCAAACCTACTGCGGCCCGATTGGTGCCGAGTATATGCACATTACCAGCACCGAAGAAAAACGCTGGATCCAACAGCGTATCGAGTCTGGTcgcg < 1:1548952/100‑1 (MQ=255) cAAGCAAACCTACTGCGGCCCGATTGGTGCCGAGTATATGCACATTACCAGCGCCGAAGAAAAACGCTGGATCCAACAGCGTATCGAGTCTGGTcgcgcg < 2:1760293/100‑1 (MQ=255) cAAGCAAACCTACTGCGGCCCGATTGGTGCCGAGTATATGCACATTACCAGCACCGAAGACAAACGCTGGCTCCAACAGCGTATCGAGTCTGGTagcgcg > 2:446926/1‑100 (MQ=255) cAAGCAAACCTACTGCGGCCCGATTGGTGCCGAGTATATGCACATTACCAGCACCGAAGAAAAACGATGGATCCAACAGCGTATCGAGTCTGGTcgcgcg > 1:1196336/1‑100 (MQ=255) aGCAAACCTACTGCGGCCCGATTGGTGCCGAGTATATGCACATTACCAGCACCGAAGAAAAACGCTGGATCCAACAGCGTATCGAGTCTGGTCGCGCGAc > 1:1057170/1‑100 (MQ=255) cAAACCTACTGCGGCCCGATTGGTGCCGAGTATATGCACATTACCAGCGCCGAAGAAAAACGCTGGATCCAACAGCGTATCGAGTCTGGTCGCGCGACtt < 1:1499008/100‑1 (MQ=255) aCCTACTGCGGCCCGCTTGGTGCCGAGTATATGCACATTACCAGCACCGAAGAAAAACGCTGGATCCGACAGCGTATCGAGTCTGGTCGCGCGACTTTCa < 1:2281151/100‑1 (MQ=255) cTACTGCGGCCCGATTGGTGCCGAGTATATGCACATTACCAGCACCGAAGAAAAACGCTGGATCCAACAGCGTATCGAGTCTGGTCGCGCGACTTTCAAt < 1:1463137/100‑1 (MQ=255) cTGCGGCCCGATTGGTGCCGAGTATATGCACATTACCAGCACCGAAGAAAAACGCTGGATCCAACAGCGTATCGAGTCTGGTCGCGCGACTTTCAATAGc < 1:642680/100‑1 (MQ=255) cTGCGGACCGATTGGTGCCGAGTATATGCACATTACCAGCACCGAAGAAAAACGCTGGATCCAACAGCGTATCGAGTATGGTCGAGCGACTTTCAATa > 2:1937381/1‑98 (MQ=255) tGCGGCCCGATTGGTGCCGAGTATATGCACATTACCAGCACCGAAGAAAAACGCTGGATCCAACAGCGTATCGAGTCTGGTCGCGCGACTTTCAATAGCg > 2:738153/1‑100 (MQ=255) tGCGGCCCGATTGGTGCCGAGTATATGCACATTACCAGCACCGAAGAAAAACGCTGGATCCAAAAGCGTATCGAGTCTGGTCGCGCGACTTTCAATAGCg > 2:984577/1‑100 (MQ=255) gCGGCCCGATTGGTGCCGAGTATATGCACATTACCAGCACCGAAGAAAAACGCTGGATCCAACAGCGTATCGAGTCTGGTCGCGCGACTTTCAATAGCGa < 1:2336110/100‑1 (MQ=255) ggCCCGATTGGTGCCGAGTATATGCACATTACCAGCACCGAAGAAAAACGCTGGATCCAACAGCGTATCGAGTCTGGTCGCGCGACTTTCAATAGCGAag > 2:155592/1‑100 (MQ=255) gCCCGATTGGTGCCGAGTATATGCACATTACCAGCACCGAAGAAAAACGCTGGATCCAACAGCGTATCGAGTCTGGTCGCGCGACTTTCAATAGCGAaga > 2:43838/1‑100 (MQ=255) ccGATTGGTGCCGAGTATATGCACATTACCATCACCGAAGAAAAACGCTGGATCCAACAGCGTATCGAGTCTGGTCGCGCGACTTTCAATAGCGAAgaga < 1:432490/100‑1 (MQ=255) cGATTGGTGCCGAGTATATGCACATTACCAGCACCGAAGAAAAACGCTGGATCCAACAGCGTATCGAGTCTGGTCGCGCGACTTTCAATAGCGAAGAGaa > 1:1245263/1‑100 (MQ=255) aTTGGTGCCGAGTATATGCACATTACCAGCACCGAAGAAAAACGCTGGATCCAACAGCGTATCGAGTCTGGTCGCGCGACTTTCAATAGCGAAGAGaaaa < 1:778123/100‑1 (MQ=255) tGGTGCCGAGTATATGCACATTACCAGCACCGAAGAAAAACGCTGGATCCAACAGCGTATCGAGTCTGGTCGCGCGACTTTCAATAGCGAAGAGaaaaaa < 1:1473782/100‑1 (MQ=255) tGCCGAGTATATGCACATTACCAGCGCCGAAGAAAAACGCTGGATCCAACAGCGTATCGAGTCTGGTCGCGCGACTTTCAATAGCGAAGAGAAAAAACg > 1:596254/1‑99 (MQ=255) ccGAGTATCTGCACATTACCAGCACCGAAGAAAACCGCTGGCTCCACCAGCGTATCGAGTCTGGTCGCGCGACTTCTAATTGCGAAGTGAAAAAACGc < 2:1394366/98‑1 (MQ=255) gTATATGCACATTACCAGCACCGAAGAAAAACGCTGGATCCAACAGCGTATCGAGTCTGGTCGCGCGACTTTCAATAGCGAAGAGAAAAAACGCTTCTTa > 1:1827567/1‑100 (MQ=255) tataTGCACATTACCAGCACCGAAGAAAAACGCTGGATCCAACAGCGTATCGAGTCTGGTCGCGCGACTTTCAATAGCGAAGAGAAAAAACGCTTCTTa < 1:1102816/99‑1 (MQ=255) atatGCACATTACCAGCACCGAAGAAAAACGCTGGATCCAACAGCGTATCGAGTCTGGTCGCGCGACTTTCAATAGCGAAGAGAAAAAACGCTTCTTAAg > 1:1254200/1‑100 (MQ=255) tatGCACATTACCAGCACCGAAGAAAAACGCTGGATCCAACAGCGTATCGAGTCTGGTCGCGCGACTTTCAATAGCGAAGAGAAAAAACGCTTCTTAAGc < 1:288222/100‑1 (MQ=255) tatGCACATTACCAGCACCGAAGAAAAACGCTGGATCCAACAGCGTATCGAGTCTGGTCGCGCGACTTTCAATAGCGAAGAGAAAAAACGCTTCTTAAGc > 1:81951/1‑100 (MQ=255) tatGCACATTACCAGCACCGAAGAAAAACGCTGGATCCAACAGCGTATCGAGTCTGGTCGCGCGACTTTCAATAGCGAAGAGAAAAAACGCTTCTTAAGc < 1:1633800/100‑1 (MQ=255) tGCACAAGACCAGAACAGAAGAAAAACGCTGGACACAACAGTGTATCGAGTCTGGTCGCGCGACTTTTAAAAAAAAAAAAAAAAAACGCTTCTTAAGCg < 2:45217/99‑1 (MQ=255) gCACATTACCAGCACCGAAGAAAAACGCTGGATCCAACAGCGTATCGAGTCTGGTCGCGCGACTTTCAATAGCGAAGAGAAAAAACGCTTCTTAAGCgga > 1:1235448/1‑98 (MQ=255) cacaTTACCAGCGCCGAAGAAAAACGCTGGATCCAACAGCGTATCGAGTCTGGTCGCGCGACTTTCAATAGCGAAGAGAAAAAACGCTTCTTAAGCGAAc < 2:1207084/100‑1 (MQ=255) cacaTTACCAGCACCGAAGAAAAACGCTGGATCCAACAGCGTATCGAGTCTGGTCGCGCGACTTTCAATAGCGAAGAGAAAAAACGCTTCTTAAGCGAAc > 1:2163277/1‑100 (MQ=255) caTTACCAGCGCCGAAGAAAAACGCTGGATCCAACAGCGTATCGAGTCTGGTCGCGCGACTTTCAATAGCGAAGAGAAAAAACGCTTCTTAAGCGAACTg > 2:414595/1‑100 (MQ=255) ttACCAGCACCGAAGAAAAACGCTGGATCCAACAGCGTATCGAGTCTGGTCGCGCGACTTTCAATAGCGAAGAGAAAAAACGCTTCTTAAGCGAACTGAc > 1:2249879/1‑100 (MQ=255) ccAGCACCGAAGAAAAACGCTGGATCCAACAGCGTATCGAGTCTGGTCGCGCGACTTTCAATAGCGAAGAGAAAAAACGCTTCTTAAGCGAACTGAccgc > 2:1345083/1‑100 (MQ=255) gCACCGAAGAAAAACGCTGGATCCAACAGCGTATCGAGTCTGGTCGCGCGACTTTCAATAGCGAAGAGAAAAAACGCTTCTTAAGCGAACTGACcgccgc < 2:416023/100‑1 (MQ=255) gCACCGAAGAAAAACGCTGGATCCAACAGCGTATCGAGTCTGGTCGCGCGACTTTCAATAGCGAAGAGAAAAAACGCTTCTTAAGCGAACTGACcgccgc < 1:271923/100‑1 (MQ=255) gCACCGAAGAAAAACGCTGGATCCAACAGCGTATCGAGTCTGGTCGCGCGACTTTCAATAGCGAAGAGAAAAAACGCTTCTTAAGCGAACTGACcgccgc > 2:1933917/1‑100 (MQ=255) cACCGAAGAAAAACGCTGGATCCAACAGCGTATCGAGTCTGGTCGCGCGACTTTCAATAGCGAAGAGAAAAAACGCTTCTTAAGCGAACTGAccgccg > 2:1892263/1‑98 (MQ=255) cACCGAAGAAAAACGCTGGATCCAACAGCGTATCGAGTCTGGTCGCGCGACTTTCAATAGCGAAGAGAAAAAACGCTTCTTAAGCGAACTGACCGCCGCt < 1:1933939/100‑1 (MQ=255) cACCGAAGAAAAACGCTGGATCCAACAGCGTATCGAGTCTGGTCGCGCGACTTTCAATAGCGAAGAGAAAAAACGCTTCTTAAGCGAACTGACCGCCGCt > 1:1025812/1‑100 (MQ=255) | GCCAGCGGCAAAGAAACCATGAAACTCGGCGAGCTGCTGGAAGCCCTCAAGCAAACCTACTGCGGCCCGATTGGTGCCGAGTATATGCACATTACCAGCACCGAAGAAAAACGCTGGATCCAACAGCGTATCGAGTCTGGTCGCGCGACTTTCAATAGCGAAGAGAAAAAACGCTTCTTAAGCGAACTGACCGCCGCT > NC_000913/759138‑759335 |
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: ‑ |