Predicted mutation | ||||||
---|---|---|---|---|---|---|
evidence | seq id | position | mutation | annotation | gene | description |
RA | REL606 | 3,911,997 | (A)6→7 | intergenic (‑198/‑125) | yifB ← / → ilvL | predicted bifunctional enzyme and transcriptional regulator/ilvG operon leader peptide |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | REL606 | 3,911,991 | 1 | . | A | 80.0% | 92.0 / 17.5 | 30 | intergenic (‑192/‑131) | yifB/ilvL | predicted bifunctional enzyme and transcriptional regulator/ilvG operon leader peptide |
Reads supporting (aligned to +/- strand): ref base . (2/4); new base A (10/14); total (12/18) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 1.00e+00 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 |
AAAGGATGAACTACGAGGAAGGGAACAACATTCATACTGAAATTGAATTTTTTTCACTCACTATTTTAT‑T‑T‑TT‑AAAA‑A‑ACAACAATTTATATTGAAATTATTAAACGCATCATAAAAATCGGCCAAAAAATATCTTGTACTATT > REL606/3911919‑3912062 | aaaGGATGAACTACGAGGAAGGGAACAACATTCATACTGAAATTGAATTTTTTTCACTCACTATTTTAT‑T‑T‑TT‑aaa > 2:270083/1‑76 (MQ=255) aaGGATGAACTACGAGGAAGGGAACAACATTCATACTGAAATTGAATTTTTTTCACTCACTATTTTAT‑T‑T‑TT‑aaaa > 2:2645/1‑76 (MQ=255) ngATGAACTACGAGGAAGGGAACAACATTCATACTGAAATTGAATTTTTTTCACTCACTATTTTAT‑T‑T‑TT‑aaaa‑a‑a < 1:1004409/75‑1 (MQ=255) ggATGAACTACGAGGAAGGGAACAACATTCATACTGAAATTGAATTTTTTTCACTCACTATTTTAT‑T‑T‑TT‑aaaa‑a‑a < 1:635835/76‑1 (MQ=255) gATGAACTACGAGGTAGGGAACAACATTCATACTGAAATTGAATTTTTTTCACTCACTATTTTAT‑T‑T‑TT‑aaaa‑a‑aa < 2:747464/76‑2 (MQ=255) tACGAGGGAGGGAGCAACATTCATACTGATATTGCATTTTTTTCACTCACTATTTTAT‑T‑T‑TT‑AAAA‑ATGCAACAAt < 2:1032078/76‑1 (MQ=255) tACGAGGAAGGGAACAACATTCATACTGAAATTGAATTATTTTCACTCACTATTTTAT‑T‑T‑TTAAAAA‑A‑ACAACAAt > 1:840060/1‑76 (MQ=255) cGAGGAAGGGAACAACATTCATACTGAAATTGAATTTTTTTCACTCACTATTTTAT‑TCT‑TA‑AAAA‑A‑acaaca < 1:31175/72‑1 (MQ=255) cGAGGAAGGGAACAACATTCATACTGAAATTGAATTTTTTTCACTCACTATTTTAT‑T‑T‑TTAAAAA‑A‑acaaca > 2:31175/1‑72 (MQ=255) cGAGGAAGGGAACAACATTCATACTGAAATTGAATTTTTTTCACTCACTATTTTAT‑T‑T‑TTAAAAA‑A‑ACAACAAttt < 1:1131028/76‑1 (MQ=255) aaGGGAACAACATTCATACTGAAATTGAATTTTTTTCACTCACTATTTTAT‑T‑T‑TTAAAAA‑A‑AAAAAAATTTTTAtt > 2:714258/1‑76 (MQ=255) ggAACAACATTCATACTGAAATTGAATTTTTTTCACTCACTATTTTAT‑T‑T‑TTAAAAA‑A‑ACAACAATTTATATTGaa < 1:200711/76‑1 (MQ=255) gAACAACACTCATACTGAAATTGAATTTTTTTCATTCACTATTTTAT‑TCT‑TA‑AAAA‑A‑ACAACGATTTATATTGaaa > 2:1030757/1‑76 (MQ=255) aacaTTCATACTGAAATTGAATTTTTTTCACTCACTATTTTAT‑T‑T‑TTAAAAA‑A‑ACAACAATTTATATTGAAattat < 2:507663/76‑1 (MQ=255) caTTCATACTGAAATTGAATTTTTTTCACTCACTATTTTAT‑T‑T‑TTAAAAA‑A‑ACAACAATTTATATTGAAAttatta > 1:1128498/1‑76 (MQ=255) aTTCATACTGAAATTGAATTTTTTTCACTCACTATTTTAT‑T‑T‑TTAAAAA‑A‑ACAACAATTTATATTGAAATTATTaa < 2:1148780/76‑1 (MQ=255) aTTCATACTGAAATTGAATTTTTTTCACTCACTATTTTAT‑T‑T‑TTAAAAA‑A‑ACAACAATTTATATTGAAATTATTaa < 1:437857/76‑1 (MQ=255) ttCATACTGAAATTGAATTTTTTTCACTCACTATTTTAT‑T‑T‑TTAAAAA‑A‑ACAACAATTTATATTGAAATTATTaaa > 1:617932/1‑76 (MQ=255) tACTGAAATTGAATTTTTTTCACTCACTATTTTAT‑T‑T‑TTAAAAA‑A‑ACAACAATTTATATTGAAATTATTAAACGca > 2:844429/1‑76 (MQ=255) cTGAAATTGAATTTTTTTCACTCACTATTTTAT‑T‑T‑TTAAAAA‑A‑ACAACAATTTATATTGAAATTATTAAACGcatc < 2:518153/76‑1 (MQ=255) tGAAATTGAATTTTTTTCACTCACTATTTTAT‑T‑T‑TTAAAAA‑A‑ACAACAATTTATATTGAAATTATTAAACGcatca < 1:710758/76‑1 (MQ=255) tGAAATTGAATTTTTTTCACCGACTATTTTAT‑T‑T‑TTAAAAA‑A‑ACAACAATTTATATTGAAATTATTAAACGcatca < 2:973930/76‑1 (MQ=255) aaaTTGAATTTTTTTCACTCACTATTTTAT‑T‑T‑TTAAAAA‑A‑ACAACAATTTATATTGAAATTATTAAACGCATCATa < 2:332400/76‑1 (MQ=255) aTTGAATTTTTTTCACTCACTATTTTAT‑T‑T‑TTAAAAA‑A‑ACAACAATTTATATTGAAATTATTAAACGCATCATaaa < 2:668277/76‑1 (MQ=255) gAATTTTTTTCACTCACTATTTTAT‑T‑T‑TTAAAAA‑A‑ACAACAATTTATATTGAAATTATTAAACGCATCATAAAAAt > 1:684227/1‑76 (MQ=255) ttttCACTCACTATTTTAT‑T‑T‑TTAAAAA‑A‑ACAACAATTTATATTGAAATTATTAAACGCATCATAAAAATCGGCCa < 1:218631/76‑1 (MQ=255) tattCACTCACTATTTTAT‑T‑TATA‑AAAA‑A‑ACATCAATTTATATTGAAATTATTAATCGCATCATAAAAATCGGCCa < 1:140512/74‑1 (MQ=255) tcactcacTATTTTATCT‑T‑TA‑AAAA‑A‑ACAACAATTTATATTGAAATTATTAAACGCATCATAAAAATCGGCCaaaa < 2:486802/76‑1 (MQ=255) ccctcactATTTTAT‑T‑T‑TTAAAAA‑A‑ACAACAATTTATATTGCAATTATTAAAAGCATCATAAAAATCGGCCaaaaa > 2:403353/3‑76 (MQ=255) cactcactATTTTAT‑T‑T‑TTAAAAA‑A‑ACCACAATTTATATTGAAATTATTAAACGCATCATAAAAATCCTCCaaaaa < 2:1036171/76‑1 (MQ=255) cactcactATTTTAT‑T‑T‑TTAAAAA‑A‑ACAACAATTTATATTGAAATTATTAAACGCATCATAAAAATCGGCCaaaaa > 2:66131/1‑76 (MQ=255) actATTTTAT‑T‑T‑TT‑AAAACA‑ACAACAATTTATATTGAAATTATTAAACGCATCATAAAAATCGGCCAAAAAATATc > 2:1048875/1‑76 (MQ=255) tattttat‑t‑t‑TTAAAAA‑A‑ACAACAATTTATATTGAAATTATTAAACGCATCATAAAAATCGGCCAAAAAATATCtt > 1:634018/1‑76 (MQ=255) attttat‑t‑t‑tTAAAAA‑A‑ACAACAATTTATATTGAAATTATTAAACGCATCATAAAAATCGGCCAAAAAATATCTTg > 2:360167/1‑76 (MQ=255) ttttat‑t‑t‑tTAAAAA‑A‑ACAACAATTTATATTGAAATTATTAAACGCATCATAAAAATCGGCCAAAAAATATCTTGt < 2:1147380/76‑1 (MQ=255) ttttat‑t‑t‑tTAAAAA‑A‑ACAACAATTTATATTGAAATTATTAAACGCATCATAAAAATCGGCCAAAAAATATCTTGt < 2:730109/76‑1 (MQ=255) att‑t‑t‑ta‑AAAA‑A‑ACAACAATTTATATTGAAATTATTAAACGCATCATAAAAATCGGCCAAAAAATATCTTGTACTa > 1:435277/3‑76 (MQ=255) att‑t‑t‑ta‑AAAA‑A‑ACAACAATTTATATTGAAATTATTAAACGCATCATAAAAATCGGCCAAAAAATATCTTGTACTa > 1:824685/3‑76 (MQ=255) nt‑t‑t‑ta‑AAAA‑A‑ACAACAATTTATATTGAAATTATTAAACGCATCATAAAAATCGGCCAAAAAATATCTTGTACTAt < 1:1002455/75‑1 (MQ=255) nt‑t‑t‑ta‑AAAA‑A‑ACAACAATTTATATTGAAATTATTAAACGCATCATAAAAATCGGCCAAAAAATATCTTGTACTAt < 1:1002450/75‑1 (MQ=255) t‑t‑t‑ta‑AAAA‑A‑ACAACAATTTATATTGAAATTATTAAACGCATCATAAAAATCGGCCAAAAAATATCTTGTACTAtt > 1:537519/1‑76 (MQ=255) | AAAGGATGAACTACGAGGAAGGGAACAACATTCATACTGAAATTGAATTTTTTTCACTCACTATTTTAT‑T‑T‑TT‑AAAA‑A‑ACAACAATTTATATTGAAATTATTAAACGCATCATAAAAATCGGCCAAAAAATATCTTGTACTATT > REL606/3911919‑3912062 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 26 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |