Predicted mutation | ||||||
---|---|---|---|---|---|---|
evidence | seq id | position | mutation | annotation | gene | description |
RA | REL606 | 3,271,926 | T→C | D38G (GAC→GGC) | yrbA ← | predicted DNA‑binding transcriptional regulator |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | REL606 | 3,271,926 | 0 | T | C | 98.2% | 212.8 / ‑3.9 | 58 | D38G (GAC→GGC) | yrbA | predicted DNA‑binding transcriptional regulator |
Reads supporting (aligned to +/- strand): ref base T (1/0); new base C (36/20); total (38/20) | |||||||||||
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.44e-01 |
TATACGCTTTGATCGACACAGCATGAATGCGGTTATCCGCAATATATTCCATCAGCGGACCATAGACCGTCTGCTGTTTTTTAACCCGACTCATGCCGTCAAACAACTCACCCACGGCAATAACCTGAAAGTGGCTGCCATCGCCGGAAACGTGGACTTCCTGGAGGGAGAGAGCGTTCATCAACACGCTCTGA > REL606/3271828‑3272021 | tataCGCTTTGATCGACACAGCATGAATGCGGTTATCCGCAATATATTCCATCAGCGGACCATAGACCGTCTGCTGTTTTTTAACCCGACTCATGCCGCCa > 1:213429/1‑101 (MQ=255) ataCGCTTTGATCGACACAGCATGAATGCGGTTATCCGCAATATATTCCATCAGCGGACCATAGACCGTCTGCTGTTTTTTAACCCGACTCATGCCGCCaa > 2:752146/1‑101 (MQ=255) aCGCTTTGATCGACACAGCATGAATGCGGTTATCCGCAATATATTCCATCAGCGGACCATAGACCGTCTGCTGTTTTTTAACCCGACTCATGCCGCCAaac < 1:1671397/101‑1 (MQ=255) cTTTGATCGACACAGCATGAATGCGGTTATCCGCAATATATTCCATCAGCGGACCATAGACCGTCTGCTGTTTTTTAACCCGACTCATGCCGCCAaacaac > 1:765891/1‑101 (MQ=255) cTTTGATCGACACAGCATGAATGCGGTTATCCGCAATATATTCCATCAGCGGACCATAGACCGTCTGCTGTTTTTTAACCCGACTCATGCCGCCAaacaac < 1:1315015/101‑1 (MQ=255) ttGATCGACACAGCATGAATGCGGTTATCCGCAATATATTCCATCAGCGGACCATAGACCGTCTGCTGTTTTTTAACCCGACTCATGCCGCCAAACAACTc > 2:783539/1‑101 (MQ=255) ttGATCGACACAGCATGAATGCGGTTATCCGCAATATATTCCATCAGCGGACCATAGACCGTCTGCTGTTTTTTAACCCGACTCATGCCGCCAAACAACTc > 1:1410459/1‑101 (MQ=255) gATCGACACAGCATGAATGCGGTTATCCGCAATATATTCCATCAGCGGACCATAGACCGTCTGCTGTTTTTTAACCCGACTCATGCCGCCAAACAACTCAc < 2:1193095/101‑1 (MQ=255) tCGACACAGCATGAATGCGGTTATCCGCAATATATTCCATCAGCGGACCATAGACCGTCTGCTGTTTTTTAACCCGACTCATGCCGCCAAACAACTCAccc < 2:1160668/101‑1 (MQ=255) tCGACACAGCATGAATGCGGTTATCCGCAATATATTCCATCAGCGGACCATAGACCGTCTGCTGTTTTTTAACCCGACTCATGCCGCCAAACAACTCAccc < 2:721891/101‑1 (MQ=255) cacaGCATGAATGCGGTTATCCGCAATATATTCCATCAGCGGACCATAGACCGTCTGCTGTTTTTTAACCCGACTCATGCCGCCAAACAACTCACCCACgg > 2:113211/1‑101 (MQ=255) gCATGAATGCGGTTATCCGCAATATATTCCATCAGCGGACCATAGACCGTCTGCTGTTTTTTAACCCGACTCATGCCGCCAAACAACTCACCCACGGCAAt < 2:993406/101‑1 (MQ=255) cATGAATGCGGTTATCCGCAATATATTCCATCAGCGGACCATAGACCGTCTGCTGTTTTTTAACCCGACTCATGCCGCCAAACAACTCACCCACGGCAATa < 1:1026985/101‑1 (MQ=255) cATGAATGCGGTTATCCGCAATATATTCCATCAGCGGACCATAGACCGTCTGCTGTTTTTTAACCCGACTCATGCCGCCAAACAACTCACCCACGGCAATa > 2:941254/1‑101 (MQ=255) aTGAATGCGGTTATCCGCAATATATTCCATCAGCGGACCATAGACCGTCTGCTGTTTTTTAACCCGACTCATGCCGCCAAACAACTCACCCACGGCAATaa > 2:1634393/1‑101 (MQ=255) aTGAATGCGGTTATCCGCAATATATTCCATCAGCGGACCATAGACCGTCTGCTGTTTTTTAACCCGACTCATGCCGACAAACAACTCACCCACGGCAATaa > 2:1545367/1‑101 (MQ=255) gAATGCGGTTATCCGCAATATATTCCATCAGCGGACCATAGACCGTCTGCTGTTTTTTAACCCGACTCATGCCGCCAAACAACTCACCCACGGCAATAAcc < 1:942770/101‑1 (MQ=255) aaTGCGGTTATCCGCAATATATTCCATCAGCGGACCATAGACCGTCTGCTGTTTTTTAACCCGACTCATGCCGCCAAACAACTCACCCACGGCAATAACCt > 2:910851/1‑101 (MQ=255) aTGCGGTTATCCGCAATATATTCCATCAGCGGACCATAGACCGTCTGCTGTTTTTTAACCCGACTCATGCCGCCAAACAACTCACCCACGGCAATAACCTg < 2:528698/101‑1 (MQ=255) ggTTATCCGCAATATATTCCATCAGCGGACCATAGACCGTCTGCTGTTTTTTAACCCGACTCATGCCGCCAAACAACTCACCCACGGCAATAACCTGAAAg > 1:285534/1‑101 (MQ=255) gTTATCCGCAATATATTCCATCAGCGGACCATAGACCGTCTGCTGTTTTTTAACCCGACTCATGCCGCCAAACAACTCACCCACGGCAATAACCTGAAAGt > 1:467185/1‑101 (MQ=255) ccGCAATATATTCCATCAGCGGACCATAGACCGTCTGCTGTTTTTTAACCCGACTCATGCCGCCAAACAACTCACCCACGGCAATAACCTGAAAGTGGCTg > 2:233539/1‑101 (MQ=255) cGCAATATATTCCATCAGCGGACCATAGACCGTCTGCTGTTTTTTAACCCGACTCATGCCGCCAAACAGCTCGCCCACGGCAATAACCTGAAAGTGGCTGc < 1:752151/101‑1 (MQ=255) atatatTCCATCAGCGGACCATAGACCGTCTGCTGTTTTTTAACCCGACTCATGCCGCCAAACAACTCACCCACGGCAATAACCTGAAAGTGGCTGCCATc > 2:316814/1‑101 (MQ=255) atatatTCCATCAGCGGACCATAGACCGTCTGCTGTTTTTTAACCCGACTCATGCCGCCAAACAACTCACCCACGGCAATAACCTGAAAGTGGCTGCCATc < 1:674201/101‑1 (MQ=255) tatTCCATCAGCGGACCATAGACCGTCTGCTGTTTTTTAACCCGACTCATGCCGCCAAACAACTCACCCACGGCAATAACCTGAAAGTGGCTGCCATCGcc > 1:454832/1‑101 (MQ=255) tatTCCATCAGCGGACCATAGACCGTCTGCTGTTTTTTAACCCGACTCATGCCGCCAAACAACTCACCCACGGCAATAACCTGAAAGTGGCTGCCATCGcc < 1:225958/101‑1 (MQ=255) ttCCATCAGCGGACCATAGACCGTCTGCTGTTTTTTAACCCGACTCATGCCGCCAAACAACTCACCCACGGCAATAACCTGAAAGTGGCTGCCATCGCCgg > 1:752470/1‑101 (MQ=255) ccATCAGCGGACCATAGACCGTCTGCTGTTTTTTAACCCGACTCATGCCGCCAAACAACTCACCCACGGCAATAACCTGAAAGTGGCTGGCATCGCCGGaa > 1:1641122/1‑101 (MQ=255) ccATCAGCGGACCATAGACCGTCTGCTGTTTTTTAACCCGACTCATGCCGCCAAACAACTCACCCACGGCAATAACCTGAAAGTGGCTGCCATCGCCGGaa > 2:1137613/1‑101 (MQ=255) cATCAGCGGACCATAGACCGTCTGCTGTTTTTTAACCCGACTCATGCCGCCAAACAACTCACCCACGGCAATAACCTGAAAGTGTCTGCCATCGCCGGaaa > 2:1192941/1‑101 (MQ=255) cATCAGCGGACCATAGACCGTCTGCTGTTTTTTAACCCGACTCATGCCGCCAAACAACTCA‑CCACGGCAATAACCTGAAAGTGGCTGCCATCGCCGGAAAc < 1:596050/101‑1 (MQ=255) cATCAGCGAACCATAGACCGTCTGCTGTTTTTTAACCCGACTCATGCCGCCAAACAACTCACCCACGGCAATAACCTGAAAGTGGCTGCCATCGCCGGaaa > 2:1535345/1‑101 (MQ=255) aTCAGCGGACCATAGACCGTCTGCTGTTTTTTAACCCGACTCATGCCGCCAAACAACTCACCCACGGCAATAACCTGAAAGTGGCTGCCATCGCCGGAAAc > 1:397581/1‑101 (MQ=255) aGCGGACCATAGACCGTCTGCTGTTTTTTAACCCGACTCATGCCGCCAAACAACTCACCCACGGCAATAACCTGAAAGTGGCTGCCATCGCCGGAAACGTg < 1:484037/101‑1 (MQ=255) aGCGGACCATAGACCGTCTGCTGTTTTTTAACCCGACTCATGCCGCCAAACAACTCACCCACGGCAATAACCTGAAAGTGGCTGCCATCGCCGGAAACGTg > 1:193507/1‑101 (MQ=255) cGGACCATAGACCGTCTGCTGTTTTTTAACCCGACTCATGCCGCCAAACAACTCACCCACGGCAATAACCTGAAAGTGGCTGCCATCGCCGGAAACGTGGa > 1:700383/1‑101 (MQ=255) ggACCATAGACCGTCTGCTGTTTTTTAACCCGACTCATGCCGCCAAACAACTCACCCACGGCAATAACCTGAAAGTGGCTGCCATCGCCGGAAACGTGGAc < 1:1634399/101‑1 (MQ=255) gACCATAGACCGTCTGCTGTTTTTTAACCCGACTCATGCCGCCAAACAACTCACCCACGGCAATAACCTGAAAGTGGCTGCCATCGCCGGAAACGTGGACt < 2:1531194/101‑1 (MQ=255) gACCATAGACCGTCTGCTGTTTTTTAACCCGACTCATGCCGCCAAACAACTCACCCACGGCAATAACCTGAAAGTGGCTGCCATCGCCGGAAACGTGGACt > 2:1492648/1‑101 (MQ=255) tAGACCGTCTGCTGTTTTTTAACCCGACTCATGCCGCCAAACAACTCACCCACGGCAATAACCTGAAAGTGGCTGCCATCGCCGGAAACGTGGACTTCCTg < 2:1410453/101‑1 (MQ=255) aCCGTCTGCTGTTTTTTAACCCGACTCATGCCGCCAAACAACTCACCCACGGCAATAACCTGAAAGTGGCTGCCATCGCCGGAAACGTGGACTTCCTggag < 2:710731/101‑1 (MQ=255) cGTCTGCTGTTTTTTAACCCGACTCATGCCGCCAAACAACTCACCCACGGCAATAACCTGAAAGTGGCTGCCATCGCCGGAAACGTGGACTTCCTggaggg > 2:26465/1‑101 (MQ=255) tCTGCTGTTTGTTAACCCGACTCATGCCGCCAAACAACTCACCCACGGCAATAACCTGAAAGTGGCTGCCATCGCCGGAAACGTGGACTTCCTggagggag > 2:1135775/1‑101 (MQ=255) ctgctgTTTTTTAACCCGACTCATGCCGCCAAACAACTCACCCACGGCAATAACCTGAAAGTGGCTGCCATCGCCGGAAACGTGGACTTCCTGGAGGgaga > 2:1628132/1‑101 (MQ=255) tgctgTTTTTTAACCCGACTCATGCCGCCAAACAACTCACCCACGGCAATAACCTGAAAGTGGCTGCCATCGCCGGAAACGTGGACTTCCTGGAGGgagag > 2:1096519/1‑101 (MQ=255) gctgTTTTTTAACCCGACTCATGCCGCCAAACAACTCACCCACGGCAATAACCTGAAAGTGGCTGCCATCGCCGGAAACGTGGACTTCCTGGAGGgagaga > 1:920079/1‑101 (MQ=255) gTTTTTTAACCCGACTCATGCCGCCAAACAACTCACCCACGGCAATAACCTGAAAGTGGCTGCCATCGCCGGAAACGTGGACTTCCTGGAGGGAGAGAGCg > 1:257318/1‑101 (MQ=255) ttttAACCCGACTCATGCCGCCAAACAACTCACCCACGGCAATAACCTGAAAGTGGCTGCCATCGCCGGAAACGTGGACTTCCTGGAGGGAGAGAGCGTtc < 1:1635849/101‑1 (MQ=255) tAACCCGACTCATGCCGCCAAACAACTCACCCACGGCAATAACCTGAAAGTGGCTGCCATCGCCGGAAACGTGGACTTCCTGGAGGGAGAGAGCGTtcatc > 1:1468263/1‑101 (MQ=255) tAACCCGACTCATGCCGCCAAACAACTCACCCACGGCAATAACCTGAAAGTGGCTGCCATCGCCGGAAACGTGGACTTCCTGGAGGGAGAGAGCGTtaatc > 1:1667872/1‑101 (MQ=255) cGACTCATGCCGCCAAACAACTCACCCACGGCAATAACCTGAAAGTGGCTGCCATCGCCGGAAACGTGGACTTCCTGGAGGGAGAGAGCGTTCATCAacac > 1:1421197/1‑101 (MQ=255) cGACGCATGCCGCCAAACAACTCACCCACGGCAATAACCTGAAAGTGGCTGCCATCGCCGGAAACGTGGACTTCCTGGAGGGAGAGAGCGTTCATCAacac > 2:1718365/1‑101 (MQ=255) cTCATGCCGCCAAACAACTCACCCACGGCAATAACCTGAAAGTGGCTGCCATCGCCGGAAACGTGGACTTCCTGGAGGGAGAGAGCGTTCATCAACACGct < 1:650253/101‑1 (MQ=255) tCATGCCGCCAAACAACTCACCCACGGCAATAACCTGAAAGTGGCTGCCATCGCCGGAAACGTGGACTTCCTGGAGGGAGAGAGCGTTCATCAACACGctc > 1:1151934/1‑101 (MQ=255) aTGCCGCCAAACAACTCACCCACGGCAATAACCTGAAAGTGGCTGCCATCGCCGGAAACGTGGACTTCCTGGAGGGAGAGAGCGTTCATCAACACGCTCTg > 2:455911/1‑101 (MQ=255) aTGCCGCCAAACAACTCACCCACGGCAATAACCTGAAAGTGGCTGCCATCGCCGGAAACGTGGACTTCCTGGAGGGAGAGAGCGTTCATCAACACGCTCTg < 1:316814/101‑1 (MQ=255) tGCCGTCAAACAACTCACCCACGGCAATAACCTGAAAGTGGCTGCCATCGCCGGAAACGTGGACTTCCTGGAGGGAGAGAGCGTTCATCAACACGCTCTGa > 1:47578/1‑101 (MQ=255) tGCCGCCAAACAACTCACCCACGGCAATAACCTGAAAGTGGCTGCCATCGCCGGAAACGTGGACTTCCTGGAGGGAGAGAGCGTTCATCAACACGCTCTGa > 2:517604/1‑101 (MQ=255) | TATACGCTTTGATCGACACAGCATGAATGCGGTTATCCGCAATATATTCCATCAGCGGACCATAGACCGTCTGCTGTTTTTTAACCCGACTCATGCCGTCAAACAACTCACCCACGGCAATAACCTGAAAGTGGCTGCCATCGCCGGAAACGTGGACTTCCTGGAGGGAGAGAGCGTTCATCAACACGCTCTGA > REL606/3271828‑3272021 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 22 ≤ ATCG/ATCG < 28 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |