| Predicted mutation | |||||||
|---|---|---|---|---|---|---|---|
| evidence | seq id | position | mutation | freq | annotation | gene | description |
| RA | AP009048 | 4,083,310 | G→A | 12.3% | G325E (GGG→GAG) | rtcB → | conserved hypothetical protein |
| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | AP009048 | 4,083,310 | 0 | G | A | 12.3% | 163.3 / 11.8 | 57 | G325E (GGG→GAG) | rtcB | conserved hypothetical protein |
| Reads supporting (aligned to +/- strand): ref base G (17/33); new base A (5/2); total (22/35) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 9.52e-02 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.24e-01 | |||||||||||
TCTGCGCGTGCTGGTCAATATGGAATTATTCCCGGTTCGATGGGAGCAAAAAGCTTTATCGTCCGTGGGCTGGGAAATGAAGAGTCGTTCTGTTCGTGCAGCCACGGTGCCGGGCGGGTAATGAGCCGAACTAAAGCG > AP009048/4083243‑4083380 | tCTGCGCGTGCTGGTCAATATGGAATTATTCCCGGTTCGATGGGAGCAAAAAGCTTTATCGTCCGTGGGc < 1:3358569/70‑1 (MQ=255) cgcgTGCTGGTCAATATGGAATTATTCCCGGTTCGATGGGAGCAAAAAGCTTTATCGTCCGTGGGCTg < 1:2156037/68‑1 (MQ=255) gTGCTGGTCAATATGGAATTATTCCCGGTTCGATGGGAGCAAAAAGCTTTATCGTCCGTGGGCTGGGAAAt < 1:754418/71‑1 (MQ=255) tGCTGGTCAATATGGAATTATTCCCGGTTCGATGGGAGCAAAAAGCTTTATCGTCCGTGGGCTGGGAAATg < 1:3833147/71‑1 (MQ=255) tGCTGGTCAATATGGAATTATTCCCGGTTCGATGGGAGCAAAAAGCTTTATCGTCCGTGGGCTGGGAAATg < 1:218885/71‑1 (MQ=255) cTGGTCAATATGGAATTATTCCCGGTTCGATGGGAGCAAAAAGCTTTATCGTCCGTGGGCTGGGAAATGaa < 1:4057105/71‑1 (MQ=255) tCAATATGGAATTATTCCCGGTTCGATGGGAGCAAAAAGCTTTATCGTCCGTGGGCTGGGAAATGAAGAGt < 1:3608614/71‑1 (MQ=255) aaTATGGAATTATTCCCGGTTCGATGGGAGCAAAAAGCTTTATCGTCCGTGGGCTGGGAAATGAAGAGTc < 1:4218013/70‑1 (MQ=255) aaTATGGAATTATTCCCGGTTCGATGGGAGCAAAAAGCTTTATCGTCCGTGGGCTGGGAAATGAAGAGTCg > 1:1331771/1‑71 (MQ=255) tatGGAATTATTCCCGGTTCGATGGGAGCAAAAAGCTTTATCGTCCGTGGGCTGGGAAATGAAGAGTCGtt < 1:4492400/71‑1 (MQ=255) tatGGAATTATTCCCGGTTCGATGGGAGCAAAAAGCTTTATCGTCCGTGAGCTGGGAAATGAAGAGTCGt > 1:248188/1‑70 (MQ=255) atGGAATTATTCCCGGTTCGATGGGAGCAAAAAGCTTTATCGTCCGTGGGCTGGGAAATGAAGAGTCGtt < 1:4497670/70‑1 (MQ=255) tGGAATTATTCCCGGTTCGATGGGAGCAAAAAGCTTTATCGTCCGTGGGCTggg > 1:3978433/1‑54 (MQ=255) gAATTATTCCCGGTTCGATGGGAGCAAAAAGCTTTATCGTCCGTGGGCTGGGAAATGAAGAGTCGTTCTg > 1:4176012/1‑70 (MQ=255) gAATTATTCCCGGTTCGATGGGAGCAAAAAGCTTTATCGTCCGTGAGCTGGGAAATGAAGAGTCGTTCTGt > 1:417638/1‑71 (MQ=255) aaTTATTCCCGGTTCGATGGGAGCAAAAAGCTTTATCGTCCGTGGGCTGGGAAATGAAGAGTCGTTCTGtt < 1:1846679/71‑1 (MQ=255) aaTTATTCCCGGTTCGATGCGAGCAAAAAGCTTTATCGTCCGTGGGCTggg > 1:2194682/1‑51 (MQ=255) ttattCCCGGTTCGATGGGAGCAAAAAGCTTTATCGTCCGTGGGCTGGGAAATGAAGAGTCg < 1:2478668/62‑1 (MQ=255) attCCCGGTTCGATGGGAGCAAAAAGCTTTATCGTCCGTgg > 1:1727973/1‑41 (MQ=255) attCCCGGTTCGATGGGAGCAAAAAGCTTTATCGTCCGTgg > 1:699777/1‑41 (MQ=255) attCCCGGTTCGATGGGAGCAAAAAGCTTTATCGTCCGTGGTCTGGGAAATGAAGAGTCGTTCTGTTCGt < 1:652732/70‑1 (MQ=255) attCCCGGTTCGATGGGAGCAAAAAGCTTTATCGTCCGTGGGCTGGGAAATGAAGAGTCGTTCTGTTCGt < 1:3384823/70‑1 (MQ=255) attCCCGGTTCGATGGGAGCAAAAAGCTTTATCGTCCGTGGGCTGGGAAATGAAGAGTCGTTCTGTTCGTg > 1:1394229/1‑71 (MQ=255) ttCCCGGTTCGATGGGAGCAAAAAGCTTTATCGTCCGTGGGCTGGGAAATGaa < 1:3496611/53‑1 (MQ=255) ttCCCGGTTCGATGGGAGCAAAAAGCTTTATCGTCCGTGGGCTGGGAAATGAag > 1:1529088/1‑54 (MQ=255) attccGGTTCGATGGGAGCAAAAAGCTTTATCGTCCGTGGGCTGGGAAATGAAGAGTCGTTCTGTTCGTg > 1:2332864/4‑70 (MQ=255) cccGGTTCGATGGGAGCAAAAAGCTTTATCGTCCGTGGGCTGGGAAATGAAGAGTCGTTCTGTTCGTGCAg > 1:1937086/1‑71 (MQ=255) ccGGTTCGATGGGAGCAAAAAGCTTTATCGTCCGTGGGCTGGGAAATGAAGAGTCGTTCTGTTCGTGc < 1:3363190/68‑1 (MQ=255) ggTTCGATGGGAGCAAAAAGCTTTATCGTCCGTGGGCTGGGAAATGAAGAGTCGTTCTGTTCGTGCAGCCa < 1:326788/71‑1 (MQ=255) ggTTCGATGGGAGCAAAAAGCTTTATCGTCCGTGGGCTGGGAAATGAAGAGTCGTTCTGTTCGTGCAGCCa < 1:2138170/71‑1 (MQ=255) ggTTCGATGGGAGCAAAAAGCTTTATCGTCCGTGAGCTGGGAAATGAAGAGTCGTTCTGTTCGTGCAGCCa > 1:3855970/1‑71 (MQ=255) gTTCGATGGGAGCAAAAAGCTTTATCGTCCGTGGGCTGGGAAATGAAGAGTCGTTCTGTTCGTGCAGCCAc < 1:1106931/71‑1 (MQ=255) gTTCGATGGGAGCAAAAAGCTTTATCGTCCGTGGGCTGGGAAATGAAGAGTCGTTCTGTTCGTGCAGCCAc < 1:2020348/71‑1 (MQ=255) tCGATGGGAGCAAAAAGCTTTATCGTCCGTGAGCTTGGAAATGa > 1:2366826/1‑44 (MQ=38) cGATGGGAGCAAAAAGCTTTATCGTCCGTGGGCTGGGAAATGAAGAGTCGTTCTGTTCGt > 1:2610400/1‑60 (MQ=255) gATGGGAGCAAAAAGCTTTATCGTCCGTGGGCTGGGAAATGAAGAGTCGTTCTGTTc < 1:1627950/57‑1 (MQ=255) gATGGGAGCAAAAAGCTTTATCGTCCGTGGGCTGGGAAATGAAGAGTCGTTCTGTTCGTGCAGc < 1:3030926/64‑1 (MQ=255) aTGGGAGCAAAAAGCTTTATCGTCCGTGAGCTGGGAAATGAAGAGTCGTTCTGTTCGTGCAGCCACGGTGc < 1:3746883/71‑1 (MQ=255) tGGGAGCAAAAAGATTTATCGTCCGTGAGCTGGGAAATGAAGAGTCGTTCTGTTCGTGCAGCCACGGTGcc > 1:1731608/1‑71 (MQ=255) ggAGCAAAAAGCTTTATCGTCCGTGGGCTGGGAAATGAAGAGTCGTTCTGTTCGTGCAGCCACGGTGCcgg < 1:4311830/71‑1 (MQ=255) ggAGCAAAAAGCTTTATCGTCCGTGGGCTGGGAAATGAAGAGTCGTTCTGTTCGTGCAGCCACGGTGCcgg > 1:173958/1‑71 (MQ=255) ggAGCAAAAAGCTTTATCGTCCGTGGGCTGGGAAATGAAGAGTCGTTCTGTTCGTGCAGCCACGGTGCcgg > 1:637433/1‑71 (MQ=255) aGCAAAAAGCTTTATCGTCCGTGGGCTGGGAAATGAAGAGTCGTTCTGTTCGTGCAGCCACGGTGCcgggc > 1:1038872/1‑71 (MQ=255) gCAAAAAGCTTTATCGTCCGTGGGCTGGGAAATGAAGAGTCGTTCTGTTCGTGCAGCCACGGTGCcgggcg > 1:279972/1‑71 (MQ=255) aaaaaGCTTTATCGTCCGTGGGCTGGGAAATGAAGAGTCGTTCTGTTCGTGCAGCCACGGTGCcgggcgg < 1:3438424/70‑1 (MQ=255) aaaaGCTTTATCGTCCGTGGGCTGGGAAATGAAGAGTCGTTCTGTTCGTGCAGCCACGGTGCcgggcggg > 1:3798578/1‑70 (MQ=255) aaaaGCTTTATCGTCCGTGGGCTGGGAAATGAAGAGTCGTTCTGTTCGTGCAGCCACGGTGCCGGGCGGGt < 1:2683885/71‑1 (MQ=255) aaaaGCTTTATCGTCCGTGGGCTGGGAAATGAAGAGTCGTTCTGTTCGTGCAGCCACGGTGCCGGGCGGGt < 1:1752232/71‑1 (MQ=255) aaaaGCTTTATCGTCCGTGGGCTGGGAAATGAAGAGTCGTGCTGTTCGTGCAGCCACGGTGCCGGGCGGGt < 1:2689315/71‑1 (MQ=255) aaaaGCTTTATCGTCCGTGAGCTGGGAAATGAAGAGTCGTTCTGTTCGTGCAGCCACGGTGCcgg < 1:3063136/65‑1 (MQ=255) aaaaCTTTATCGTCCGTGGGCTGGGAAATGAAGAGTCGTTCTGTTCGTGCAGCCACGGTGCCGGGCGGGt > 1:4185992/1‑70 (MQ=255) aaaGCTTTATCGTCCGTGGGCTGGGAAATGAAGAGTCGTTCTGTTCGTg < 1:4548841/49‑1 (MQ=255) aaaGCTTTATCGTCCGTGGGCTGGGAAATGAAGAGTCGTTCTGTTCGTGCAGCCACGGTGCCGGGCGGGTa < 1:1334683/71‑1 (MQ=255) aaaGCTTTATCGTCCGTGGGCTGGGAAATGAAGAGTCGTTCTGTTCGTGCAGCCACGGTGCCAGGCGGGTa < 1:3806061/71‑1 (MQ=255) aaGCTTTATCGTCCGTGGGCTGGGAAATGAAGAGTCGTTCTGTTCGTGCAGCCACGGTGCcggg > 1:2185653/1‑64 (MQ=255) aGCTTTATCGTCCGTGGGCTGGGAAATGAAGAGTCGTTCTGTTCGTGCAGCCACGGTGCCGGGCGGGTa > 1:203339/1‑69 (MQ=255) aGCTTTATCGTCCGTGGGCTGGGAAATGAAGAGTCGTTCTGTTCGTGCAGCCACGGTGCCGGGCGGGTAAt < 1:3091469/71‑1 (MQ=255) tATCGTCCGTGGGCTGGGAAATGAAGAGTCGTTCTGt < 1:3519576/37‑1 (MQ=255) tCGTCCGTGGGCTGGGAAATGAAGAGTCGTTCTGTTCGTGCAGCCACGGTGCcgggcggg < 1:1770192/60‑1 (MQ=255) ggCTGGGAAATGAAGAGTCGTTCTGTTCGTGCAGCCACGGTGCCGGGCGGGTAATGAGCCGAACTAAAGCg < 1:4320116/71‑1 (MQ=255) | TCTGCGCGTGCTGGTCAATATGGAATTATTCCCGGTTCGATGGGAGCAAAAAGCTTTATCGTCCGTGGGCTGGGAAATGAAGAGTCGTTCTGTTCGTGCAGCCACGGTGCCGGGCGGGTAATGAGCCGAACTAAAGCG > AP009048/4083243‑4083380 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |