Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | AP009048 | 2,951,877 | G→A | 6.8% | A928V (GCT→GTT) | recB ← | exonuclease V (RecBCD complex), beta subunit |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | AP009048 | 2,951,877 | 0 | G | A | 6.8% | 233.8 / 4.9 | 73 | A928V (GCT→GTT) | recB | exonuclease V (RecBCD complex), beta subunit |
Reads supporting (aligned to +/- strand): ref base G (36/32); new base A (2/3); total (38/35) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 6.66e-01 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.92e-01 |
GGTGACGCACCGCGCGGAAACTGATGTGGTGTTAACGTCGGTTCTTCAACGACGCTGGCAACGCCTGCAGCATCGACATCCAGCCGAGGCATCAAATCCTGGGCGATACCGTGACCACGCTGTTGCAAACCAGAGTA > AP009048/2951808‑2951944 | ggTGACGCACCGCGCGGAAACTGATGTGGTGTTAACGTCGGTTCTTCAAAGACGCTGGCAACGCCTgcagc < 1:3226302/71‑1 (MQ=255) gTGACGCACCGCGCGGAAACTGATGTGGTGTTAACGTCGGTTCTTCAACGACGCTGGCAACGCCTgcagca < 1:1288309/71‑1 (MQ=255) gTGACGCACCGCGCGGAAACTGATGTGGTGTTAACGTCGGTTCTTCAACGACGCTGGCAACGCCTgcaaca < 1:3232997/71‑1 (MQ=255) gACGCACCGCGCGGAAACTGATGTGGTGTTAACGTCGGTTCTTCAACGACGCTGGCAACGCCTGCAGCATc > 1:1513721/1‑71 (MQ=255) aCCGCGCGGAAACTGATGTGGTGTTAACGTCGGTTCTTCAACGACGCTGGCAACGCCTGCAGCATCGACAt > 1:1673003/1‑71 (MQ=255) gcgcgGAAACTGATGTGGTGTTAACGTCGGTTCTTCAACGACGCTGGCAACGCCTGCAGCATCGACATc > 1:2512268/1‑69 (MQ=255) ggAAACTGATGTGGTGTTAACGTCGGTTCTTCAACGACGCTGGCAACGCCTGCAGCATCGACATCCAGcc < 1:2305506/70‑1 (MQ=255) aaaCTGATGTGGTGTTAACGTCGGTTCTTCAACGACGCTGGCAACGCCTGCAGCATCGACATCCAGCCGAg > 1:3045005/1‑71 (MQ=255) tGATGTGGTGTTAACGTCGGTTCTTCAACGACGCTGGCAACGCCTGCAGCATCGACAt > 1:1837598/1‑58 (MQ=255) aTGTGGTGTTAACGTCGGTTCTTCAACGACGCTGGCAACGCCTGCAGCATCGACATCCAGCCGAGGCATc > 1:1940694/1‑70 (MQ=255) aTGTGGTGTTAACGTCGGTTCTTCAACGACGCTGGCAACGCCTGCAGCATCGACATCCAGCCGAGGCATc < 1:4128131/70‑1 (MQ=255) tgtgGTGTTAACGTCGGTTCTTCAACGACGCTGGCAACGCCTGCAGCATCGACATCCAGCCGAGGCATCa > 1:1236620/1‑70 (MQ=255) tgtTAACGTCGGTTCTTCAACGACGCTGGCAACGCCTGCAGCATCGACATCCAGCCGAGGCATCAAATCCt > 1:1100537/1‑71 (MQ=255) tgtTAACGTCGGTTCTTCAACGACGCTGGCAACGCCTGCAGCATCGACATCCAGCCGAGGCATCAAATCCt > 1:3063817/1‑71 (MQ=255) tgtTAACGTCGGTTCTTCAACGACGCTGGCAACGCCTGCAGCATCGACATCCAGCCGAGGCATCAAATCCt < 1:921802/71‑1 (MQ=255) tgtTAACGTCGGTTCTTCAACGACGCTGGCAACGCCTGCAACATCGACATCCAGCCGAGGCATCAAATCCt < 1:597560/71‑1 (MQ=255) gtTAACGTCGGTTCTTCAACGACGCTGGCAACGCCTGCAGCATCGACATCCAGCCGAg < 1:949905/58‑1 (MQ=255) ttAACGTCGGTTCTTCAACGACGCTGGCAACGCCTGCAGCATCGACATCCAGCCGAGGCATCa < 1:2104561/63‑1 (MQ=255) aaCGTCGGTTCTTCAACGACGCTGGCAACGCCTGCAGCATCGACATCCAGCCGa > 1:1998254/1‑54 (MQ=255) aaCGTCGGTTCTTCAACGACGCTGGCAACGCCTGCAGCATCGACATCCAGCCGAGGCATCAAATCCTGGGc > 1:862568/1‑71 (MQ=255) aaCGTCGGTTCTTCAACGACGCTGGCAACGCCTGCAGCATCGACATCCACCCGAGGCATCAAATCCTGGGc > 1:410971/1‑71 (MQ=255) cGTCGGTTCTTCAACGACGCTGGCAACGCCTGCAGCATCGACATCCAGCCGAGGCATCAAATCCTGGGCGa > 1:4310832/1‑71 (MQ=255) gTCGGTTCTTCAACGACGCTGGCAACGCCTGCAGCATCGACATCCAGCCGAGGCATCAAATCCTGGGCGAt < 1:3744279/71‑1 (MQ=255) cGGTTCTTCAACGACGCTGGCAACGCCTGCAGCATCGACATCCAGCCGAGGCATCAAATCCTGGGCGATAc > 1:3837762/1‑71 (MQ=255) ggTTCTTCAACGACGCTGGCAACGCCTGCAGCATCGACATCCAGCCGAGGCATCaa > 1:2732000/1‑56 (MQ=255) ggTTCTTCAACGACGCTGGCAACGCCTGCAGCATCGACATCCAGCCGAGGCATCAAATCCTGGGCGATAcc > 1:1485063/1‑71 (MQ=255) ggTTCTTCAACGACGCTGGCAACGCCTGCAACATCGACATCCAGCCGAGGCATCAAATCCTg < 1:3268711/62‑1 (MQ=255) gTTCTTCAACGACGCTGGCAACGCCTGCAGCATCGACATCCAGCCGAGGCATCAAATCCTGGGCGATACCg < 1:3844001/71‑1 (MQ=255) gTTCTTCAACGACGCTGGCAACGCCTGCAACATCGACATCCAGCCGAGGCATCAAATc > 1:3065334/1‑58 (MQ=255) ttcttcAACGACGCTGGCAACGCCTGCAGCATCGACATCCAGCCGAGGCATCAAATc < 1:2008777/57‑1 (MQ=255) ttcttcAACGACGCTGGCAACGCCTGCAGCATCGACATCCAGCCGAGGCATCAAATCCTGGGCGATACCGt < 1:4419994/71‑1 (MQ=255) ttcttcAACGACGCTGGCAACGCCTGCAGCATCGACATCCAGCCGAGGCATCAAATCCTGGGCGATACCGt < 1:933213/71‑1 (MQ=255) ttcttcAACGACGCTGGCAACGCCTGCAGCATCGACATCCAGCCGAGGCATCAAATCCTGGGCGATACCGt < 1:4434778/71‑1 (MQ=255) ttcttcAACGACGCTGGCAACGCCTGCAGCATCGACATCCAGCCGAGGCATCAAATCCTGGGCGATACCGt < 1:3484423/71‑1 (MQ=255) tcttcAACGACGCTGGCAACGCCTGCAGCATCGACATCCAGCCGAGGCATCAAATCCTGGGCGATACCGt > 1:2783076/1‑70 (MQ=255) tcttcAACGACGCTGGCAACGCCTGCAGCATCGACATCCAGCCGAGGCATCAAATCCTGGGCGATACCGt > 1:4468824/1‑70 (MQ=255) cttcAACGACGCTGGCAACGCCTGCAGCATCGACATCCAGCCGa > 1:448918/1‑44 (MQ=255) tcAACGACGCTGGCAACGCCTGCAGCATCGACATCCAGCCGAGGCATCAAATCCTGGGCGATACCGTGAc > 1:1719323/1‑70 (MQ=255) cAACGACGCTGGCAACGCCTGCAGCATCGACATCCAGCCGAGGCATCAAATCCTg < 1:3695344/55‑1 (MQ=255) cAACGACGCTGGCAACGCCTGCAGCATCGACATCCAGCCGAGGCATCAAATCCTGGGCGATACCGTGAcc < 1:1336076/70‑1 (MQ=255) cAACGACGCTGGCAACGCCTGCAGCATCGACATCCAGCCGAGGCATCAAATCCTGGGCGATACCGTGAc < 1:159247/69‑1 (MQ=255) cAACGACGCTGGCAACGCCTGCAGCATCGACATACAGCCGAGGCATCAAATCCTg < 1:1686209/55‑1 (MQ=255) aaCGACGCTGGCAACGCCTGCAGCATCGACATCCAGCCGAGGCATCAAATCCt < 1:2845138/53‑1 (MQ=255) acgacgCTGGCAACGCCTGCAGCATCGACATCAAGCCGAGGCATCAAATCCTGGGCGATACCGTGACCACg > 1:1390546/1‑71 (MQ=255) cgacgCTGGCAACGCCTGCAGCATCGACATCCAGCCGAGGCATCAAATCCTGGGCGa > 1:842361/1‑57 (MQ=255) cgacgCTGGCAACGCCTGCAGCATCGACATCCAGCCGAGGCATCAAATCCTGGGCGATACCGTGACCACGc > 1:942023/1‑71 (MQ=255) cgacgCTGGCAACGCCTGCAGCATCGACATCCAGCCGAGGCATCAAATCCTGGGCGATACCGTGACCACGc < 1:238887/71‑1 (MQ=255) gacgCTGGCAACGCCTGCAGCATCGACATCCAGCCGAgga > 1:3227710/1‑39 (MQ=255) gacgCTGGCAACGCCTGCAGCATCGACATCCAGCCGAGGCATCAAATCCTGGGCGATACCGTGACCACGCt < 1:2813839/71‑1 (MQ=255) gacgCTGGCAACGCCTGCAGCATCGACATCCAGCCGAGGCATCAAATCCTGGGCGATACCGTGACCACGCt < 1:996566/71‑1 (MQ=255) acgCTGGCAACGCCTGCAGCATCGACATCCAGCCGAGGc > 1:4468753/1‑39 (MQ=255) acgCTGGCAACGCCTGCAGCATCGACATCCAGCCGAGGCATCAAATCCTGGGCGATACCGTGACCACGCTg < 1:868808/71‑1 (MQ=255) acgCTGGCAACGCCTGCAGCATCGACATCCAGCCGAGGCATCAAATCCTGGGCGATACCGTGACCACGCTg < 1:1143069/71‑1 (MQ=255) acgCTGGCAACGCCTGCAACATCGACATCCAGCCGAGGCATCAAATCCTGGGCGATACCGTGACCACGCTg < 1:3541470/71‑1 (MQ=255) cgCTGGCAACGCCTGCAGCATCGACATCCAGCCGAGGCATCAAATCCTGGGCGATACCGTGACCACGCTGt < 1:3578991/71‑1 (MQ=255) cgCTGGCAACGCCTGCAGCATCGACATCCAGCCGAGGCATCAAATCCTGGGCGATACCGTGACCACGCTGt < 1:3195975/71‑1 (MQ=255) gCTGGCAACGCCTGCAGCATCGACATCCAGCCGAGGCATCAAATcc > 1:3954324/1‑46 (MQ=255) cTGGCAACGCCTGCAGCATCGACATCCAGCCGAGGCATCAAATCCTGGGCGATACCGTGACCACGCTGTTg < 1:1542637/71‑1 (MQ=255) cTGGCAACGCCTGCAGCATCGACATCCAGCCGAGGCATCAAATCCTGGGCGATACCGTGACCACGCTGTTg > 1:1864046/1‑71 (MQ=255) tGGCAACGCCTGCAGCATCGACATCCAGCCGAGGCATCAAATCCTGGGCGATACCg > 1:1888447/1‑56 (MQ=255) gCAACGCCTGCAGCATCGACATCCAGCCGAGGCATCAAATCCTGGGCGATACCg < 1:4530258/54‑1 (MQ=255) gCAACGCCTGCAGCATCGACATCCAGCCGAGGCATCAAATCCTGGGCGATACCGTGACCACGCTGTTGCa > 1:3955084/1‑70 (MQ=255) gCAACGCCTGCAGCATCGACATCCAGCCGAGGCATCAAATCCTGGGCGATACCGTGACCACGCTGTTGCa > 1:4075723/1‑70 (MQ=255) aaCGCCTGCAGCATCGACATCCAGCCGAGGCATCAAATCCTGGGCGATACCGTg > 1:1020327/1‑54 (MQ=255) cGCCTGCAGCATCGACATCCAGCCGAGGCATCAAATCCTGGGCGATACCGTGACCACGCTGTTGCaaa > 1:220708/1‑68 (MQ=255) cGCCTGCAGCATCGACATCCAGCCGAGGCATCAAATCCTGGGCGATACCGTGACCACGCTGTTGCAAACCa < 1:4314943/71‑1 (MQ=255) cGCCTGCAGCATCGACATCCAGCCGAGGCATCAAATCCTGGGCGATACCGTGACCACGCTGTTGCAAACCa < 1:2621859/71‑1 (MQ=255) gCCTGCAGCATCGACATCCAGCCGAGGCATCAAATCCTGGGCGATACCGTGACCACGCTGTTGCAAACCag < 1:1941353/71‑1 (MQ=255) ccTGCAGCATCGACATCCAGCCGAGGCATCAAATCCTGGGCGATACCGTGACCACGCTGt > 1:1719946/1‑60 (MQ=255) ccTGCAGCATCGACATCCAGCCGAGGCATCAAATCCTGGGCGATACCGTGACCACGCTGTTGCAAACCaga < 1:4300510/71‑1 (MQ=255) cTGCAGCATCGACATCCAGCCGAGGCATCAAATCCTGGGCGATAcc < 1:2821657/46‑1 (MQ=255) cTGCAGCATCGACATCCAGCCGAGGCATCAAATCCTGGGCGATACCGTGACCACGCTGTTGCa > 1:634707/1‑63 (MQ=255) cTGCAGCATCGACATCCAGCCGAGGCATCAAATCCTGGGCGATACCGTGACCACGCTGTTGCAAACCagag > 1:3638036/1‑71 (MQ=255) cTGCAACATCGACATCCAGCCGAGGCATCAAATCCTGGGCGATACCGTGACCACGCTGTTGCAAACCaga > 1:1336670/1‑70 (MQ=255) tGCAGCATCGACATCCAGCCGAGGCATCAAATCCTGGGCGATACCGTGACCACGCTGtt > 1:2102890/1‑59 (MQ=255) tGCAGCATCGACATCCAGCCGAGGCATCAAATCCTGGGCGATACCGTGACCACGCTGTTGCa < 1:3730080/62‑1 (MQ=255) gcagcaTCGACATCCAGCCAAGGCATCAAATCCTGGGCGATACCGTGACCACGCTGTTGCaaa < 1:1521639/63‑1 (MQ=255) gcaacaTCGACATCCAGCCGAGGCATCAAATCCTGGGCGATACCGTGACCACGCTGtt > 1:2908658/1‑58 (MQ=255) cagcaTCGACATCCAGCCGAGGCATCAAATCCTGGGCGATACCGTGACCACGCTGTTGCAAACCAGAGTa < 1:688587/70‑1 (MQ=255) | GGTGACGCACCGCGCGGAAACTGATGTGGTGTTAACGTCGGTTCTTCAACGACGCTGGCAACGCCTGCAGCATCGACATCCAGCCGAGGCATCAAATCCTGGGCGATACCGTGACCACGCTGTTGCAAACCAGAGTA > AP009048/2951808‑2951944 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 21 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |