| Predicted mutation | |||||||
|---|---|---|---|---|---|---|---|
| evidence | seq id | position | mutation | freq | annotation | gene | description |
| RA | AP009048 | 3,639,730 | T→C | 6.7% | K222R (AAA→AGA) | xerC ← | site‑specific tyrosine recombinase |
| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | AP009048 | 3,639,730 | 0 | T | C | 6.7% | 203.8 / 6.2 | 74 | K222R (AAA→AGA) | xerC | site‑specific tyrosine recombinase |
| Reads supporting (aligned to +/- strand): ref base T (34/35); new base C (2/3); total (36/38) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 1.00e+00 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.32e-01 | |||||||||||
GAGTGACGTAATTTATGCGGATGAACGTGATTATTCAGCCCTTGTTTTATGCCCCATTCGGCAAAGCGTTTCTGCACATTACGCGCGGAGATACGCTTGCCCAGTTTCGACAGAAAAAGCGCGTCGTCTTCGCTA > AP009048/3639661‑3639795 | gAGTGACGTAATTTATGCGGATGAACGTGATTATTCAGCCCTTGTTTTATGCCCCATTCGGCAAAGCGttt < 1:981759/71‑1 (MQ=255)gAGTGACGTAATTTATGCGGATGAACGTGATTATTCAGCCCTTGTTTTATGCCCCATTCGGCAAAGCGttt > 1:1509841/1‑71 (MQ=255) tGACGTAATTTATGCGGATGAACGTGATTATTCAGCCCTTGTTTTATGCCCCATTCGGCAAAGCGTTTCTg > 1:23952/1‑71 (MQ=255) gACGTAATTTATGCGGATGAACGTGATTATTCAGCCCTTGTTTTATGCCCCATTCGGCAAAGCGttt < 1:3419889/67‑1 (MQ=255) gACGTAATTTATGCGGATGAACGTGATTATTCAGCCCTTGTTTTATGCCCCATTCGGCAAAGCGTTTCTg > 1:3591317/1‑70 (MQ=255) aCGTAATTTATGCGGATGAACGTGATTATTCAGCCCTTGTTTTATGCCCCATTCGGCAAAGCGTTTCTGca > 1:3165646/1‑71 (MQ=255) aCGTAATTTATGCGGATGAACGTGATTATTCAGCCCTTGTTTTATGCCCCATTCGGCAAAGCGTTTCTGc < 1:1108589/70‑1 (MQ=255) cGTAATTTATGCGGATGAACGTGATTATTCAGCCCTTGTTTTATGCCCCATTCGGCAAAGCGTTTCTGcac < 1:1791419/71‑1 (MQ=255) gTAATTTATGCGGATGAACGTGATTATTCAGCCCTTGTTTTATGCCCCATTCGGCAAAGCGTTTCt > 1:2196831/1‑66 (MQ=255) gTAATTTATGCGGATGAACGTGATTATTCAGCCCTTGTTTTATGCCCCATTCGGCAAAGCGTTTCTGcaca > 1:1936386/1‑71 (MQ=255) gTAATTTATGCGGATGAACGTGATTATTCAGCACTTGTTTTATGCCCCATTCGGCAAAGCGTTTCt > 1:1024315/1‑66 (MQ=255) tAATTTATGCGGATGAACGTGATTATTCAGCCCTTGTTTTATGCCCCATTCGGCAAAGCGTTTCTGCACAt < 1:1864357/71‑1 (MQ=255) aaTTTATGCGGATGAACGTGATTATTCAGCCCTTGTTTTATGCCCCATTCGGCAAAGCGTTTCTGc > 1:2526644/1‑66 (MQ=255) aTTTATGCGGATGAACGTGATTATTCAGCCCTTGTTTTATGCCCCATTCGGCAAAGCGTTTCTGCACAtt < 1:2362615/70‑1 (MQ=255) aTTTATGCGGATGAACGTGATTATTCAGCCCTTGTTTTATGCCCCATTCGGCAAAGCGTTTCTGCACAtt < 1:3976190/70‑1 (MQ=255) aTGCGGATGAACGTGATTATTCAGCCCTTGTTTTATGCCCCATTCGGCAAAGCGTTTCTGCACATTAcgc < 1:2915125/70‑1 (MQ=255) tGCGGATGAACGTGATTATTCAGCCCTTGTTTTATGCCCCATTCGGCAAAGCGTTTCTGCACATTAcgcgc < 1:2091549/71‑1 (MQ=255) tGCGGATGAACGTGATTATTCAGCCCTTGTTTTATGCCCCATTCGGCAAAGCGTTTCTGCACATTAcgcgc < 1:1669825/71‑1 (MQ=255) ggATGAACGTGATTATTCAGCCCTTGTTTTATGCCCCATTCGGCAAAGCGTTTCTGCACATTAcgcgcg < 1:1307645/69‑1 (MQ=255) gATGAACGTGATTATTCAGCCCTTGTTTTATGCCCCATTCGGCAAAGCGTTTCTGCACATTAcgcgcg > 1:3824723/1‑68 (MQ=255) gATGAACGTGATTATTCAGCCCTTGTTTTATGCCCCATTCGGCAAAGCGTTTCTGCACATTACGCGCGgag < 1:4326147/71‑1 (MQ=255) aTGAACGTGATTATTCAGCCCTTGTTTTATGCCCCATTCGGCAAAGCGTTTCTGcaca < 1:3911018/58‑1 (MQ=255) aTGAACGTGATTATTCAGCCCTTGTTTTATGCCCCATTCGGCAAAGCGTTTCTGCACATTACGCGCGgaga < 1:755931/71‑1 (MQ=255) aTGAACGTGATTATTCAGCCCTTGTTTTATGCCCCATTCGGCAAAGCGTCTCTGCACATTAcgcgcg > 1:2204287/1‑67 (MQ=255) gAACGTGATTATTCAGCCCTTGTTTTATGCCCCATTCGGCAAAGCGTTTCTGcac > 1:3944171/1‑55 (MQ=255) gAACGTGATTATTCAGCCCTTGTTTTATGCCCCATTCGGCAAAGCGTTTCTGCACAtt > 1:600488/1‑58 (MQ=255) cGTGATTATTCAGCCCTTGTTTTATGCCCCATTCGGCAAAGCGTTTCTGCACATTACGCGCGGAGATACGc < 1:2435203/71‑1 (MQ=255) cGTGATTATTCAGCCCTTGTTTTATGCCCCATTCGGCAAAGCGTTTCTGCACATTACGCGCGGAGATACGc > 1:3866998/1‑71 (MQ=255) tGATTATTCAGCCCTTGTTTTATGCCCCATTCGGCAAAGCGTTTCTGCACATTACGCGCGGAGATACGCtt < 1:429678/71‑1 (MQ=255) gATTATTCAGCCCTTGTTTTATGCCCCATTCGGCAAAGCGTCTCTGCACATTACGCGCGGAGATACGCTTg > 1:932595/1‑71 (MQ=255) ttattCAGCCCTTGTTTTATGCCCCATTCGGCAAAGCGTTTCTGCACATTACGCGCGGAGATACGCTTGcc > 1:3903523/1‑71 (MQ=255) attCAGCCCTTGTTTTATGCCCCATTCGGCAAAGCGTTTCTGCACATTACGCGCGGAGATACGCTTGcc > 1:2999960/1‑69 (MQ=255) cAGCCCTTGTTTTATGCCCCATTCGGCAAAGCGTCTCTGCACATTACGCGCGGAGATACGCTTGccc < 1:1138333/67‑1 (MQ=255) aGCCCTTGTTTTATGCCCCATTCGGCAAAGCGTTTCTGCACATTa > 1:3092371/1‑45 (MQ=255) aGCCCTTGTTTTATGCCCCATTCGGCAAAGCGTTTCTGCACATTACGCGCGGAGATACGCTTGCCCAGtt < 1:4606657/70‑1 (MQ=255) aGCCCTTGTTTTATGCCCCATTCGGCAAAGCGTTTCTGCACATTACGCGCGGAGATACGCTTGCCCAGtt < 1:4018130/70‑1 (MQ=255) gCCCTTGTTTTATGCCCCATTCGGCAAAGCGTTTCTg < 1:2924519/37‑1 (MQ=255) gCCCTTGTTTTATGCCCCATTCGGCAAAGCGTTTCTGCACATTACGCGCGGAGATACGCTTGCCCAGTTTc < 1:1608436/71‑1 (MQ=255) cccTTGTTTTATGCCCCATTCGGCAAAGCGTTTCTTCACATTACGCGCGGAGATACGCTTGccc < 1:242187/64‑1 (MQ=255) cccTTGTTTTATGCCCCATTCGGCAAAGCGTCTCTGCACATTACGCGCGGAGATACGCTTGCCCAGTTTCg < 1:4569797/71‑1 (MQ=255) tGTTTTATGCCCCATTCGGCAAAGCGTTTCTGCACATTACGCGCGGAGATACGCTTGCCCAGTTTCGa > 1:4265120/1‑68 (MQ=255) tGTTTTATGCCCCATTCGGCAAAGCGTTTCTGCACATTACGCGCGGAGATACGCTTGCCCAGTTTCGa > 1:3710152/1‑68 (MQ=255) tGTTTTATGCCCCATTCGGCAAAGCGTCTCTGCACATTACGCGCGga < 1:430919/47‑1 (MQ=255) gTTTTATGCCCCATTCGGCAAAGCGTTTCTGcaca > 1:1885166/1‑35 (MQ=255) gTTTTATGCCCCATTCGGCAAAGCGTTTCTGCACATTACGCGCTGAGATACGCTTGCCCAGTTTCGACAg < 1:3125695/70‑1 (MQ=255) gTTTTATGCCCCATTCGGCAAAGCGTTTCTGCACATTACGCGCGGAGATACGCTTGCCCAGTTTCGACAg < 1:1145920/70‑1 (MQ=255) ttttATGCCCCATTCGGCAAAGCGTTTCTGCACAtt < 1:4579432/36‑1 (MQ=255) ttttATGCCCCATTCGGCAAAGCGTTTCTGCACATTACGCGCGGAGATACGCTTg < 1:3403042/55‑1 (MQ=255) tttATGCCCCATTCGGCAAAGCGTTTCTGCACATTACGCGCGGAGATACGCTTGCCCAGTTTCGACAGaa < 1:491162/70‑1 (MQ=255) ttATGCCCCATTCGGCAAAGCGTTTCTGCACATTACGCGCGGAGATACGc > 1:2180563/1‑50 (MQ=255) tATGCCCCATTCGGCAAAGCGTTTCTGCACATTACGCGCGGAGATACGCTTGCCCAGTTTCGACAGaaaa < 1:513036/70‑1 (MQ=255) aTGCCCCATTCGGCAAAGCGTTTCTGCACATTACGCGCGga > 1:2578392/1‑41 (MQ=255) aTGCCCCATTCGGCAAAGCGTTTCTGCACATTACGCGCGGAGATACGCTTGCCCAGTTTCGACAGaaaaa < 1:2578621/70‑1 (MQ=255) tGCCCCATTCGGCAAAGCGTTTCTGCACATTACGCGCGGAGATACGCTTGCCCAGTTTCGACAGaaaaa < 1:3494814/69‑1 (MQ=255) tGCCCCATTCGGCAAAGCGTTTCTGCACATTACGCGCGGAGATACGCTTGCCCAGTTTCGACAGaaaa > 1:696440/1‑68 (MQ=255) gCCCCATTCGGCAAAGCGTTTCTGCACATTACGCGCGGAGATACGCTTGCCCAGTTTCGACAGAAAAAgcg > 1:3844155/1‑71 (MQ=255) ccccATTCGGCAAAGCGTTTCTGCACATTACGCGCGGAGATACGCTTGCCCAGTTTCGACAGAAAAAgcgc > 1:3581873/1‑71 (MQ=255) cccATTCGGCAAAGCGTTTCTGCACATTACGCGCGGAGATACGCTTGCCCAGTTTCGACAGAAAAAGcgcg < 1:2113388/71‑1 (MQ=255) cccATTCGGCAAAGCGTTTCTGCACATTACGCGCGGAGATACGCTTGCCCAGTTTCGACAGAAAAAGcgcg < 1:3408014/71‑1 (MQ=255) ccATTCGGCAAAGCGTTTCTGCACATTACGCGCGGAGATACGCTTGCCCAGTTTCGACAGAAAAAGCGcgt > 1:370509/1‑71 (MQ=255) ccATTCGGCAAAGCGTTTCTGCACATTACGCGCGGAGATACGCTTGCCCAGTTTCGACAGAAAAAGCGcgt > 1:2690617/1‑71 (MQ=255) cATTCGGCAAAGCGTTTCTGCACATTACGCGCGGAGATACGCTTGCCCAGTTTCGa > 1:3268347/1‑56 (MQ=255) cATTCGGCAAAGCGTTTCTGCACATTACGCGCGGAGATACGCTTGCCCAGTTTCGACAGAAAAAGCGcgtc > 1:2107393/1‑71 (MQ=255) cATTCGGCAAAGCGTTTCTGCACATTACGCGCGGAGATACGCTTGCCCAGTTTCGACAGAAAAAGCGcgtc > 1:290306/1‑71 (MQ=255) aTTCGGCAAAGCGTTTCTGCACATTACGCGCGGAGATACGCTTGCCCAGTTTCGACAGAAAAAGCGcgtcg < 1:2546202/71‑1 (MQ=255) ttCGGCAAAGCGTTTCTGCACATTACGCGCGGAGATACGCTTGCCCAGTTTCGACAGAAAAAGCGcgtcgt > 1:410311/1‑71 (MQ=255) tCGGCAAAGCGTTTCTGCACATTACGCGCGGAGATACGCTTGCCCAGTTTCGACAGAAAAAGCGcgtcgt > 1:1231468/1‑70 (MQ=255) tCGGCAAAGCGTTTCTGCACATTACGCGCGGAGATACGCTTGCCCAGTTTCGACAGAAAAAGCGCgtcgtc > 1:3320927/1‑71 (MQ=255) ggCAAAGCGTTTCTGCACATTACGCGCGGAGATACGCTTGCCCAGTTTCGACAGAAAAAGCGCGTCGTCt < 1:2302631/70‑1 (MQ=255) ggCAAAGCGTTTCTGCACATTACGCGCGGAGATACGCTTGCCCAGTTTCGACAGAAAAAGCGCGTCGTCt < 1:2302568/70‑1 (MQ=255) ggCAAAGCGTTTCTGCACATTACGCGCGGAGATACGCTTGCCCAGTTTCGACAGAAAAAGCGCGTCGTCt < 1:2102521/70‑1 (MQ=255) ggCAAAGCGTTTCTGCACATTACGCGCGGAGATACGCTTGCACAGTTTCGACAGAAAAAGCGCGTCGTCtt > 1:1613422/1‑71 (MQ=255) aaaGCGTTTCTGCACATTACGCGCGGAGATACGCTTGCCCAGTTTCGACAGAAAAAGCGCGTCGTCTTCg > 1:3011758/1‑70 (MQ=255) aaaGCGTTTCTGCACATTACGCGCGGAGATACGCTTGCCCAGTTTCGACAGAAAAAGCGCGTCGTCTTCGc < 1:2056293/71‑1 (MQ=255) aaGCGTTTCTGCACATTACGCGCGGAGATACGCTTGCCCAGTTTCGACAGAAAAAGCGCGTCGTCTTCGCt < 1:4329877/71‑1 (MQ=255) aGCGTTTCTGCACATTACGCGCGGAGATACGCTTGCCCAGTTTCGACAGAAAAAGCGCGTCGTCTTCGCTa > 1:902410/1‑71 (MQ=255) gCGTTTCTGCACATTACGCGCGGAGATACGCTTGCCCAGTTTCGACAGAAAAAGCGcgtcgt > 1:823274/1‑62 (MQ=255) | GAGTGACGTAATTTATGCGGATGAACGTGATTATTCAGCCCTTGTTTTATGCCCCATTCGGCAAAGCGTTTCTGCACATTACGCGCGGAGATACGCTTGCCCAGTTTCGACAGAAAAAGCGCGTCGTCTTCGCTA > AP009048/3639661‑3639795 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |