| Predicted mutation | ||||||
|---|---|---|---|---|---|---|
| evidence | seq id | position | mutation | annotation | gene | description |
| RA | NC_000913 | 2,793,867 | C→A | A378D (GCT→GAT) | gabT → | 4‑aminobutyrate aminotransferase GabT |
| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | NC_000913 | 2,793,867 | 0 | C | A | 93.8% | 336.0 / 13.0 | 96 | A378D (GCT→GAT) | gabT | 4‑aminobutyrate aminotransferase GabT |
| Reads supporting (aligned to +/- strand): ref base C (3/3); new base A (57/33); total (60/36) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 6.69e-01 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.91e-01 | |||||||||||
CGGCGATCACAACAAGCCGGACGCCAAACTCACCGCCGAGATCGTGGCTCGCGCCCGCGATAAAGGCCTGATTCTTCTCTCCTGCGGCCCGTATTA > NC_000913/2793820‑2793915 | cGGCGATCACAACAAGCCGGACGCCAAACTCACCGCCGAGATCGTGGATc < 1:1125058/50‑1 (MQ=255)cGGCGATCACAACAAGCCGGACGCCAAACTCACCGCCGAGATCGTGGATc < 2:5442589/50‑1 (MQ=255)cGGCGATCACAACAAGCCGGACGCCAAACTCACCGCCGAGATCGTGGATc > 1:5736849/1‑50 (MQ=255) ggCGATCACAACAAGCCGGACGCCAAACTCACCGCCGAGATCGTGGATcg > 1:4033841/1‑50 (MQ=255) ggCGATCACAACAAGCCGGACGCCAAACTCACCGCCGAGATCGTGGATcg > 1:2929010/1‑50 (MQ=255) ggCGATCACAACAAGCCGGACGCCAAACTCACCGCCGAGATCGTGGATcg > 2:1579884/1‑50 (MQ=255) ggCGATCACAACAAGCCGGACGCCAAACTCACCGCCGAGATCGTGGATcg > 1:1930806/1‑50 (MQ=255) ggCGATCACAACAAGCCGGACGCCAAACTCACCGCCGAGATCGTGGATcg > 2:3281181/1‑50 (MQ=255) cGATCACAACAAGCCGGACGCCAAACTCACCGCCGAGATCGTGGCTcgcg > 1:2446732/1‑50 (MQ=255) gATCACAACAAGCCGGACGCCAAACTCACCGCCGAGATCGTGGATCgcgc > 1:3118322/1‑50 (MQ=255) gATCACAACAAGCCGGACGCCAAACTCACCGCCGAGATCGTGGATCgcgc > 2:2800821/1‑50 (MQ=255) tCACAACAAGCCGGACGCCAAACTCACCGCCGAGATCGTGGATCGCGccc > 2:3694818/1‑50 (MQ=255) cacaACAAGCCGGACGCCAAACTCACCGCCGAGATCGTGGATCGCGCCcg > 2:359018/1‑50 (MQ=255) cacaACAAGCCGGACGCCAAACTCACCGCCGAGATCGTGGATCGCGCCcg > 1:2743741/1‑50 (MQ=255) acaacaAGCCGGACGCCAAACTCACCGCCGAGATCGTGGATCGCGCCcgc > 2:2420830/1‑50 (MQ=255) acaacaAGCCGGACGCCAAACTCACCGCCGAGATCGTGGATCGCGCCcgc > 2:2283763/1‑50 (MQ=255) acaacaAGCCGGACGCCAAACTCACCGCCGAGATCGTGGATCGCGCCcgc > 2:2035860/1‑50 (MQ=255) acaacaAGCCGGACGCCAAACTCACCGCCGAGATCGTGGATCGCGCCcgc > 2:134692/1‑50 (MQ=255) caacaaGCCGGACGCCAAACTCCCCGCCGAGATCGTGGATCGCGCCcgcg < 2:3712642/50‑1 (MQ=255) caacaaGCCGGACGCCAAACTCACCGCCGAGATCGTGGATCGCGCCcgc > 1:3519282/1‑49 (MQ=255) aaacatGCCGGACGCCAAACTGTCCGCCGAGATCGTGGATCGCGCCcgcg < 2:391147/49‑1 (MQ=255) aacaaGCCGGACGCCAAACTCACCGCCGAGATCGTGGATCGCGCCCGCGa > 2:1288380/1‑50 (MQ=255) aacaaGCCGGACGCCAAACTCACCGCCGAGATCGTGGATCGCGCCCGCGa > 1:1376937/1‑50 (MQ=255) caaGCCGGACGCCAAACTCACCGCCGAGATCGTGGATCGCGCCCGCGAt > 1:4381002/1‑49 (MQ=255) aGCCGGACGCCAAACTCACCGCCGAGATCGTGGATCGCGCCCGCGATa > 1:2212364/1‑48 (MQ=255) aGCCGGACGCCAAACTCACCGCCGAGATCGTGGATCGCGCCCGCGATa > 2:5265158/1‑48 (MQ=255) aGCCGGACGCAAAACTAACCGCCGAGATCGTGGATCGCGCCCGCGATaaa < 2:399197/50‑1 (MQ=255) gCCGGACGCAAAACTAACCGCCGAGATCGTGGATCGCGCCCGCGATAAAg < 1:1780887/50‑1 (MQ=255) ccGGACGCCAAACTCACCGCCGAGATCGTGGATCGCGCCCGCGATAAAgg > 1:5710738/1‑50 (MQ=255) ccGGACGCCAAACTCACCGCCGAGATCGTGGATCGCGCCCGCGATAAAgg > 1:3791567/1‑50 (MQ=255) cGGACGCAAAACTCACCGCCGAGATCGTGGATCGCGCCCGCGATAAAGGc < 2:858139/50‑1 (MQ=255) ggACGCCAAACTCACCGCCGAGATCGTGGATCGCGCCCGCGATAAAGGcc < 2:5525286/50‑1 (MQ=255) ggACGCCAAACTCACCGCCGAGATCGTGGATCGCGCCCGCGATAAAGGcc > 2:5317776/1‑50 (MQ=255) aCGCCAACCTCACCGCCGAGATCGTGGATCGCGCCCGCGATAAAGGCCTg < 2:1328472/50‑1 (MQ=255) aCGCCAAACTCACCGCCGAGATCGTGGATCGCGCCCGCGATAAAGGCCt > 1:70207/1‑49 (MQ=255) aCGCCAAACTCACCGCCGAGATCGTGGATCGCGCCCGCGATAAAGGCCt > 1:2133104/1‑49 (MQ=255) aCGCCAAACTCACCGCCGAGATCGTGGATCGCGCCCGCGATAAAGGCCTg < 1:101996/50‑1 (MQ=255) aCGCCAAACTCACCGCCGAGATCGTGGATCGCGCCCGCGATAAAGGCCTg < 1:2496848/50‑1 (MQ=255) aCGCCAAACTCACCGCCGAGATCGTGGATCGCGCCCGCGATAAAGGACt > 2:3254143/1‑49 (MQ=39) cGCCAAACTCACCGCCGAGATCGTGGATCGCGCCCGCGATAAAGGCCTGa > 2:1964152/1‑50 (MQ=255) cGCAAAACTCACCGCCGAGATCGTGGATCGCGCCCGCGATAAAGGCCTGa < 2:4573147/50‑1 (MQ=255) gCCAAACTCACCGCCGAGATCGTGGCTCGCGCCCGCGATAAAGGCCTGAt > 1:3198037/1‑50 (MQ=255) gCCAAACTCACCGCCGAGATCGTGGATCGCGCCCGCGATAAAGGCCTGAt > 2:3483897/1‑50 (MQ=255) gCCAAACTCACCGCCGAGATCGTGGATCGCGCCCGCGATAAAGGCCTGAt > 1:2334859/1‑50 (MQ=255) ccAAACTCACCGCCGAGATCGTGGATCGCGCCCGCGATAAAGGCCTGAtt > 1:3206661/1‑50 (MQ=255) ccAAACTCACCGCCGAGATCGTGGATCGCGCCCGCGATAAAGGCCTGAtt > 2:3558930/1‑50 (MQ=255) ccAAACTCACCGCCGAGATCGTGGATCGCGCCCGCGATAAAGGCCTGAtt > 2:4428147/1‑50 (MQ=255) ccAAACTCACCGCCGAGATCGTGGATCGCGCCCGCGATAAAGGCCTGAtt > 1:5562522/1‑50 (MQ=255) cAAACTCACCGCCGAGATCGTGGCTCGCGCCCGCGATAAAGGCCTGAttc < 2:5246227/50‑1 (MQ=255) cAAACTCACCGCCGAGATCGTGGATCGCGCCCGCGATAAAGGCCTGAttc < 2:1429609/50‑1 (MQ=255) cAAACTCACCGCCGAGATCGTGGATCGCGCCCGCGATAAAGGCCTGAttc < 1:2498075/50‑1 (MQ=255) cAAACTCACCGCCGAGATCGTGGATCGCGCCCGCGATAAAGGCCTGAttc < 1:4651750/50‑1 (MQ=255) cAAACTCACCGCCGAGATCGTGGATCGCGCCCGCGATAAAGGCCTGAttc < 1:4056056/50‑1 (MQ=255) cAAACTCACCGCCGAGATCGTGGATCGCGCCCGCGATAAAGGCCTGAttc > 2:1108016/1‑50 (MQ=255) aCTCACCGCCGAGATCGTGGCTCGCGCCCGCGATAAAGGCCTGAttcttc < 1:4115404/50‑1 (MQ=255) aCTCACCGCCGAGATCGTGGATCGCGCCCGCGATAAAGGCCTGAttcttc < 2:167172/50‑1 (MQ=255) aCTCACCGCCGAGATCGTGGATCGCGCCCGCGATAAAGGCCTGAttcttc > 1:2165956/1‑50 (MQ=255) cTCACCGCCGAGATCGTGGATCGCGCCCGCGATAAAGGCCTGATtcttct < 1:18165/50‑1 (MQ=255) cTCACCGCCGAGATCGTGGATCGCGCCCGCGATAAAGGCCTGATtcttct < 1:2993766/50‑1 (MQ=255) tCACCGCCGAGATCGTGGATCGCGCCCGCGATAAAGGCCTGATTCTtctc > 1:3989897/1‑50 (MQ=255) tCACCGCCGAGATCGTGGATCGCGCCCGCGATAAAGGCCTGATTCTtctc > 1:5635732/1‑50 (MQ=255) cACCGCCGAGATCGTGGATCGCGCCCGCGATAAAGGCCTGATTCTtctct > 1:2846644/1‑50 (MQ=255) aCCGCCGAGATCGTGGCTCGCGCCCGCGATAAAGGCCTGATTCTtctc < 2:4453138/48‑1 (MQ=255) aCCGCCGAGATCGTGGATCGCGCCCGCGATAAAGGCCTGATTCTtctctc < 1:210703/50‑1 (MQ=255) aCCGCCGAGATCGTGGATCGCGCCCGCGATAAAGGCCTGATTCTtctctc > 1:5091988/1‑50 (MQ=255) cgccgAGATCGTGGATCGCGCCCGCGATAAAGGCCTGATTCTTCTCTCCt > 1:396331/1‑50 (MQ=255) cgccgAGATCGTGGATCGCGCCCGCGATAAAGGCCTGATTCTTCTCTCCt < 2:1346079/50‑1 (MQ=255) cgccgAGATCGTGGATCGCGCCCGCGATAAAGGCCTGATTCTTCTCTCCt < 1:4075262/50‑1 (MQ=255) cgccgAGATCGTGGATCGCGCCCGCGATAAAGGCCTGATTCTTCTCTCCt < 1:1490270/50‑1 (MQ=255) cgccgAGATCGTGGATCGCGCCCGCGATAAAGGCCTGATTCTTCTCTCCt < 1:1203699/50‑1 (MQ=255) cgccgAGATCGTGGATCGCGCCCGCGATAAAGGCCTGATTCTTCTCTCCt > 2:1190670/1‑50 (MQ=255) gccgAGATCGTGGATCGCGCCCGCGATAAAGGCCTGATTCTTCTCTCCTg > 1:2426446/1‑50 (MQ=255) ccgAGATCGTGGATCGCGCCCGCGATAAAGGCCTGATTCTTCTCTCCTGc > 1:874649/1‑50 (MQ=255) ccgAGATCGTGGATCGCGCCCGCGATAAAGGCCTGATTCTTCTCTCCTGc < 1:2275372/50‑1 (MQ=255) cgAGATCGTGGATCGCGCCCGCGATAAAGGCCTGATTCTTCTCTCCTGCg > 1:3997506/1‑50 (MQ=255) gagaTCGTGGATCGCGCCCGCGATAAAGGCCTGATTCTTCTCTCCTGCgg < 2:3691157/50‑1 (MQ=255) gagaTCGTGGATCGCGCCCGCGATAAAGGCCTGATTCTTCTCTCCTGCgg > 2:2124734/1‑50 (MQ=255) agaTCGTGGATCGCGCCCGCGATAAAGGCCTGATTCTTCTCTCCTGCGGc > 1:701498/1‑50 (MQ=255) agaTCGTGGATCGCGCCCGCGATAAAGGCCTGATTCTTCTCTCCTGCGGc > 2:3444969/1‑50 (MQ=255) agaTCGTGGATCGCGCCCGCGATAAAGGCCTGATTCTTCTCTCCTGCGGc > 1:4768694/1‑50 (MQ=255) agaTCGTGGATCGCGCCCGCGATAAAGGCCTGATTCTTCTCTCCTGCGGc > 2:41172/1‑50 (MQ=255) aTCGTGGCTCGCGCCCGCGATAAAGGCCTGATTCTTCTCTCCTGCGGccc > 1:3128379/1‑50 (MQ=255) aTCGTGGATCGCGCCCGCGATAAAGGCCTGATTCTTCTCTCCTGCGGccc > 1:2370421/1‑50 (MQ=255) aTCGTGGATCGCGCCCGCGATAAAGGCCTGATTCTTCTCTCCTGCGGccc > 1:456636/1‑50 (MQ=255) aTCGTGGATCGCGCCCGCGATAAAGGCCTGATTCTTCTCTCCTGCGGccc < 1:1241850/50‑1 (MQ=255) aTCGTGGATCGCGCCCGCGATAAAGGCCTGATTCTTCTCTCCTGCGGccc < 1:3113715/50‑1 (MQ=255) aTCGTGGATCGCGCCCGCGATAAAGGCCTGATTCTTCTCTCCTGCGGccc > 1:2786778/1‑50 (MQ=255) aTCGTGGATCGCGCCCGCGATAAAGGCCTGATTCTTCTCTCCTGCGGccc < 1:5257778/50‑1 (MQ=255) tCGTGGATCGCGCCCGCGATAAAGGCCTGATTCTTCTCTCCTGCGGCCCg > 1:143759/1‑50 (MQ=255) tCGTGGATCGCGCCCGCGATAAAGGCCTGATTCTTCTCTCCTGCGGCCCg < 1:161917/50‑1 (MQ=255) tCGTGGATCGCGCCCGCGATAAAGGCCTGATTCTTCTCTCCTGCGGCCCg > 2:150838/1‑50 (MQ=255) cGTGGATCGCGCCCGCGATAAAGGCCTGATTCTTCTCTCCTGCGGCCCg > 1:160362/1‑49 (MQ=255) gTGGATCGCGCCCGCGATAAAGGCCTGATTCTTCTCTCCTGCGGCCCGTa < 1:361956/50‑1 (MQ=255) gTGGATCGAGCCCGCGATAAAGGCCTGATTCTTCTCTCCTGCGGCCCGt < 1:4047356/49‑1 (MQ=39) tGGATCGCGCCCGCGATAAAGGCCTGATTCTTCTCTCCTGCGGCCCGTAt > 2:3925027/1‑50 (MQ=255) tGGATCGCGCCCGCGATAAAGGCCTGATTCTTCTCTCCTGCGGCCCGTAt < 2:4194377/50‑1 (MQ=255) cTCGCGCCCGCGATAAAGGCCTGATTCTTCTCTCCTGCGGCCCGTATTa > 2:5146050/1‑49 (MQ=255) | CGGCGATCACAACAAGCCGGACGCCAAACTCACCGCCGAGATCGTGGCTCGCGCCCGCGATAAAGGCCTGATTCTTCTCTCCTGCGGCCCGTATTA > NC_000913/2793820‑2793915 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 17 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |