Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | REL606 | 2,515,909 | 0 | A | C | 54.9% | 59.5 / 106.1 | 71 | V260G (GTG→GGG) | ypfI | predicted hydrolase |
Reads supporting (aligned to +/- strand): ref base A (22/10); new base C (1/38); total (23/48) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 9.60e-10 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.52e-10 | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
GTTAACAACGTTCGAGGAAAAC‑GCGATACCAGTTGATGCAACAATGGCGCAGGTATGGCTGCGGCTTCATCGACCACCAGCCAGTC‑GGCTTGCTCA‑T‑CGCTGGCTAACAAGGCAT‑CCGGCGCAATAAAGCGAAACTTCTCGCCCGCAAATT > REL606/2515834‑2515984 | gTTAACAACGTTCGAGGAAAAC‑GCGATACCAGTTGATGCAACAATGGCGCAGGTATGGCTGCGGCTTCATCGacca < 1:2952770/76‑1 (MQ=255) tAACAACGTTCGAGGAAAAC‑GCGATACCAGTTGATGCAACAATGGCGCAGGTATGGCTGCGGCTTCATCGaccacc < 1:1189814/76‑1 (MQ=255) tAACAACGTTCGAGGAAAAC‑GCGATACCAGTTGATGCAACAATGGCGCAGGTATGGCTGCGGCTTCATCGaccacc < 1:517561/76‑1 (MQ=255) aacaacGTTCGAGGAAAAC‑GCGATACCAGTTGATGCAACAATGGCGCAGGTATGGCTGCGGCTTCATCGAccacca < 1:4894199/76‑1 (MQ=255) acaacGTTCGAGGAAAAC‑GCGATACCAGTTGATGCAACAATGGCGCAGGTATGGCTGCGGCTTCATCGACCAccag < 2:3685617/76‑1 (MQ=255) acaacGTTCGAGGAAAAC‑GCGATACCAGTTGATGCAACAATGGCGCAGGTATGGCTGCGGCTTCATCGACCAccag < 2:3915840/76‑1 (MQ=255) acaacGTTCGAGGAAAAC‑GCGATACCAGTTGATGCAACAATGGCGCAGGTATGGCTGCGGCTTCATCGACCAccag < 2:6520452/76‑1 (MQ=255) caacGTTCGAGGAAAAC‑GCGATACCAGTTGATGCAACAATGGCGCAGGTATGGCTGCGGCTTCATCGACCAccagc < 1:2284223/76‑1 (MQ=255) aacGTTCGAGGAAAAC‑GCGATACCAGTTGATGCAACAATGGCGCAGGTATGGCTGCGGCTTCATCGACCAccagcc < 2:4735215/76‑1 (MQ=255) aacGTTCGAGGAAAAC‑GCGATACCAGTTGATGCAACAATGGCGCAGGTATGGCTGCGGCTTCATCGACCAccagcc < 1:573452/76‑1 (MQ=255) acGTTCGAGGAAAAC‑GCGATACCAGTTGATGCAACAATGGCGCAGGTATGGCTGCGGCTTCATCGACCAccagcca < 2:3467830/76‑1 (MQ=255) cGTTCGAGGAAAAC‑GCGATACCAGTTGATGCAACAATGGCGCAGGTATGGCTGCGGCTTCATCGACCAccagccag < 1:2984876/76‑1 (MQ=255) tCGAGGAAAAC‑GCGATACCAGTTGATGCAACAATGGCGCCGGTATGGCTGCGGGTTCATCCAACCCCAGCCAGCC‑g > 1:4730430/1‑76 (MQ=255) ggagaAAAGC‑CCGAAACCGGTAGATCCACCAGTGGCGTTGGTAGGCCTGCGGCTTCATCGCCCCCCAGCCAGTC‑gg < 2:304984/75‑1 (MQ=255) cGAGGAAAAC‑GCGATACCAGTTGATGCAACAATGGCGCAGGTATGGCTGGGGCTTTCTCGGCCCCCCACCAGTC‑gg > 1:3045015/1‑76 (MQ=255) aggggAAGAC‑GCGACACCAGTAGATGCAAAATTGGCGGTGGTAGGGCGGCGTCTTCATCGCCCCCCAGCCAGTC‑gg < 2:6490858/73‑1 (MQ=255) ggagAAAACGGCAATACAAGTGGATGAAACAGGGGGGCAGGTGTGGCTGCGGCTTCATCGCCCCCCAGCCAGTC‑gg < 2:3533895/73‑1 (MQ=255) gagaAAACC‑GCGATACCAGTTGATGCAACAAGGGCGCAGTTAGGGCTGCGGCTTCATCGCCCCCCAGCCAGTC‑GGc < 2:56983/76‑1 (MQ=255) gagaAAAAC‑GCGATCCCTGTTGATGCAACATTGGCGCAGTGAGGGCTGCGGTTTCATCCACCCCCAGCCAGTC‑GGc < 2:4118392/76‑1 (MQ=255) gAGGAAAAC‑GCGATACCAGTTGATGCAACAATGGCGCAGGTATGGCTGGGGCTTTCTCTACCACCCGCCAAGC‑GGc > 2:5101456/1‑76 (MQ=255) ggAAACG‑GCGATACCAGTTGAGGAAAAAAGGGCGCAGGTATGGCGGCGGCTTCATCGACCCCCAGCCAGTC‑GGCtt < 2:650091/76‑1 (MQ=255) aaaaC‑GCGATACCGTTGGATGCAAAATTGCGGAAGGTTTGGTGGCGGCTACATCCCCCCCCAGCCAGTC‑GGCTTGc < 2:6732491/76‑1 (MQ=255) cgATCCGAGTTGATGCAACAATGGCGCAGGTATGCCTGCGGCTTCATCGACCCCCAGCCAGTC‑GGCTTGCTCA‑T‑Cg < 2:2839941/76‑1 (MQ=255) caaTACCGGTTGATGCAACAAGGGCGCAGGTTGCGCTGCGGCTTCATCGCCCCCCAGCCAGTC‑GGCTTGCTCA‑T‑Cg < 2:4957911/74‑1 (MQ=255) tcccAGTGGATGCAACAATGGCGCAGTTAGGTCGGCGGCTTCATCGACCCCCAGCCAGTC‑GGCTTGCTCA‑T‑CGCTg < 2:1629472/74‑1 (MQ=255) aCCAGTTGATGCAACAAGGGCGCAGGTATGGCTGCGTCTTCATCGACCCCCAGCCAGTC‑GGCTTGCTCA‑T‑CGCTgg < 1:4348574/76‑1 (MQ=255) aCCAGTTGAGGCAACAAGGGCGGAGGTAGGGCTGCGGCTTCATCGCCCCCCAGCCAGTC‑GGCTTGCTCA‑T‑CGCTgg < 2:3181454/76‑1 (MQ=255) atttGATGCAACAAGGGCGCAGTTATGGCTGCGGCTTCATCGCCCCCCAGCCAGTC‑GGCTTGCTCA‑T‑CGCTGGCTa < 2:4559594/74‑1 (MQ=255) aGTTGATGCAACAAGGGCGCGGGTATGCCTGCGGCTTCATCGACCCCCAGCCAGTC‑GGCTTGCTCA‑T‑CGCTGGCTa < 2:3947446/76‑1 (MQ=255) gTTGATGCAACAATGGCGCAGGTAGGGCTGCGTTTTCATCGCCCCCCAGCCAGTC‑GGCTTGCTCA‑T‑CGCTGGCTaa < 1:3764371/76‑1 (MQ=255) gTTGAGCCAACAATGGCGCGGGTTGGGCTGCGTCTTCATCGCCCCCCAGCCAGTC‑GGCTTGCTCA‑T‑CGCTGGCTaa < 1:6948083/76‑1 (MQ=255) ttGATGCAACAATGGCGCAGGTATGGCTGCGGCTTTATCGACCACCAGCCAGTC‑GGCCTTCTCC‑C‑AGCTGGCtcaa > 1:5201794/1‑73 (MQ=255) ttGATGCAACAATGGCGCAGGTATGGCTGCGGCTTCATCGACCCCCAGCCAGTC‑GGCTTGCTCA‑T‑CGCTGGCTAAc < 2:6899612/76‑1 (MQ=255) ttGAGGAAAAAAGGGCGGGGGTAGGGCGGCGGCTTCACCGCCCCCCAGCCAGTC‑GGCTTGCTCA‑T‑CGCTGGCTAAc < 1:2223970/76‑1 (MQ=255) gATGCAACAATGGCGCGGGTTTGGCTGCGGCTTCATCGCCCCCCAGCCAGTC‑GGCTTGCTCA‑T‑CGCTGGCTAACaa < 1:1506785/76‑1 (MQ=255) gATGCAACAAGGGCGCGGGTGGGCCTGCGGCTTCATCGCCCCCCAGCCAGTC‑GGCTTGCTCA‑T‑CGCTGGCTAACaa < 1:3862784/76‑1 (MQ=255) ggaaacaaGGGCGCGGGTATGGCTGCGGCTTCACCGCCCCCCAGCCAGTC‑GGCTTGCTCA‑T‑CGCTGGCTAACAAgg < 2:4433833/73‑1 (MQ=255) gCAACAATGGCGCGGGTATGGCGGCGGCTTCATCGCCCCCCAGCCAGTC‑GGCTTGCTCA‑T‑CGCTGGCTAACAAGGc < 2:4252460/76‑1 (MQ=255) gCAACAAGGGCGGGGGTATGGCTGCGCCTTCATCGCCCCCCAGCCAGTC‑GGCTTGCTCA‑T‑CGCTGGCTAACAAGGc < 2:6072447/76‑1 (MQ=255) gCAACAAGGGCGCGGGTTCGCCGGCGGCTTCACCGCCCCCCAGCCAGTC‑GGCTTGCTCA‑T‑CGCTGGCTAACAAGGc < 2:5589307/76‑1 (MQ=255) aaacaaTGGCGCGGTTATGGCTGCGGTTTCATCGCCCCCCAGCCAGTC‑GGCTTGCTCA‑T‑CGCTGGCTAACAAGGCa < 1:432284/75‑1 (MQ=255) aacaaTGGCGCAGGTATGGCTGCGGCTTCATCGACCACCAGCCAGTC‑GGCCTGCTCA‑T‑CGCCGGCTAACAAAGCAt > 1:6828835/1‑76 (MQ=255) aagtGGCGCAGGTGTGCGTGCGTCTTCATCGCCCCCCAGCCAGTC‑GGCTTGCTCA‑T‑CGCTGGCTAACAAGGCAT‑cc < 1:4947464/73‑1 (MQ=255) aagggCGCAGGTGTGGCGGCGTCTTCATCGCCCCCCAGCCAGTC‑GGCTTGCTCA‑T‑CGCTGGCTAACAAGGCAT‑CCg < 1:1051282/73‑1 (MQ=255) aaTGGCGCAGGTATGGCTGCGGCTTCATCGACCACCAGCCAGTC‑GGCTTTCTCC‑TCCCCGGGCTAACAAGGCAA‑cc > 1:1338578/1‑76 (MQ=255) gtGGCGCGGGTAGGGCTGCGGCTTCATCGACCCCCAGCCAGTC‑GGCTTGCTCA‑T‑CGCTGGCTAACAAGGCAT‑CCgg < 2:3447983/75‑1 (MQ=255) aTGGCGCAGTTTTGGCTGCGGCTTCATCGACCCCCAGCCAGTC‑GGCTTGCTCA‑T‑CGCTGGCTAACAAGGCAT‑CCgg < 2:691875/76‑1 (MQ=255) gcgcAGGTATGGCTGGGGCTTCATCCACCCCCCGCCCGTC‑GGCGTTCTCACT‑CCCGGGGTAAAAAGGCAC‑CCgggg > 2:6532344/1‑74 (MQ=255) aGGTATGGCTGCGGCTTTATCCCCCAACCGCCAGCCGGGCTTCCTCA‑C‑CGCCGGTGAACAAGAAAG‑CCGCCGCaga > 2:5724194/1‑74 (MQ=255) gggaTGCCTGCGGTTTCCTCACCCACCAGCCAGTC‑GGCTTGCTCA‑T‑CGCTGGCTAACAAGGCAT‑CCGGCGCAATaa < 2:4921264/73‑1 (MQ=255) ggTATGGCTGCGGCTTTCTCGACCACCCGCCAGGC‑GGCGTTCTCC‑TCCCCGGGCTAAAAAGGACT‑CCGCGGCaaca > 2:1988752/1‑74 (MQ=255) gTATGGCTGGGCCTTCTCCGCCCACCAGCCAGTC‑GGCTTGCTCA‑T‑CGCTGGCTAACAAGGCAT‑CCGGCGCAATaaa < 1:1674248/76‑1 (MQ=255) gTATGGCTGCGGCTTCATCGACCACCAGCCAGTC‑GGCTTGCTCA‑T‑CGCTGGGTAAAAAGGGAT‑CCGGGGCCATaaa > 1:4308160/1‑76 (MQ=255) gTATGGCTGCGGCTTCATCGACCACCAGCCAGTC‑GGCTTGCTCA‑T‑CGCTGGCTAAAAAGGCAT‑CCGGCGCAATaaa > 1:2583743/1‑76 (MQ=255) ttgggCTGCGGCTTCATCCCCCCCCAGCCAGTC‑GGCTTGCTCA‑T‑CGCTGGCTAACAAGGCAT‑CCGGCGCAATAAAg < 1:4023879/73‑1 (MQ=255) gagggcggcggcTTCATCGCCCCCCAGCCAGTC‑GGCTTGCTCA‑T‑CGCTGGCTAACAAGGCAT‑CCGGCGCAATAAAg < 2:1221084/73‑1 (MQ=255) aTGGCTGCGGCTTTATCGACCACCAACCAGTC‑GGCTTGCTCC‑T‑CCCTGGCGAACAAAGGATCCCGGGGCaaaaaaa > 2:6087286/1‑75 (MQ=255) aTGGCTGCGGCTTCATCGACCACCAGCCAGTC‑GGCTTTCTCA‑T‑CGCTGGCTAAAAAAGCAT‑CCGGCGCaaaaaaac > 1:2581980/1‑74 (MQ=255) aTGGCTGCGGCTTCATCGACCACCAGCCAGTC‑GGCTTGCTCC‑T‑CGCTGGCTAACAAAGCCT‑CCCGCGCAATaaaac > 2:3398240/1‑74 (MQ=255) aTGGCTGCGGCTTCATCGACCACCAGCCAGTC‑GGCTTGCTCA‑T‑CGCTGGCTAACAAGGCAT‑CCGGCGCAATaaaac > 1:6550743/1‑74 (MQ=255) aTGGCTGCGGCTTCATCGACCACCAGCCAGTC‑GGCTTGCTCA‑T‑CGCTGGCTAAAAAGGCCT‑CCCGCGCAATAAAGc > 1:5667450/1‑76 (MQ=255) tGGCTGCGGCTTCATCGACCACCAGCCAGTC‑GGCTTGCTCA‑T‑CGCTGGCTAAAAAAGCAC‑CCGGCGCCACAAAACg > 2:3632306/1‑76 (MQ=255) gggcggcGCCTTCACCGCCCCCCAGCCAGTC‑GGCTTGCTCA‑T‑CGCTGGCTAACAAGGCAT‑CCGGCGCAATAAAGCg < 1:3854254/75‑1 (MQ=255) gggCTGCGGCTCCATCGACCACCAGCCAGTC‑GGCTTGCTCA‑T‑CGCTGGCTAACAAGGCAT‑CCGGCGCAATAAAGCg < 1:2138969/75‑1 (MQ=255) ggCTGGGTCTTCATCGCCCACCAGCCAGTC‑GGCTTGCTCA‑T‑CGCTGGCTAACAAGGCAT‑CCGGCGCAATAAAGCGa < 2:854715/76‑1 (MQ=255) gcggcgcCTTCACCGCCCCCCAGCCAGTC‑GGCTTGCTCA‑T‑CGCTGGCTAACAAGGCAT‑CCGGCGCAATAAAGCGaa < 1:4096821/73‑1 (MQ=255) gCTGCGTCTTCACCGCCCCCCAGCCAGTC‑GGCTTGCTCA‑T‑CGCTGGCTAACAAGGCAT‑CCGGCGCAATAAAGCGaa < 2:252596/76‑1 (MQ=255) gCTGCGGCTTCACCGCCCCCCAGCCAGTC‑GGCTTGCTCA‑T‑CGCTGGCTAACAAGGCAT‑CCGGCGCAATAAAGCGaa < 1:5690144/76‑1 (MQ=255) gggCGGCTTCATCCCCCCCCAGCCAGTC‑GGCTTGCTCA‑T‑CGCTGGCTAACAAGGCAT‑CCGGCGCAATAAAGCGaaa < 1:2968717/74‑1 (MQ=255) ccgcgcCTTCATCCACCCCCAGCCAGTC‑GGCTTGCTCA‑T‑CGCTGGCTAACAAGGCAT‑CCGGCGCAATAAAGCGaaa < 1:4136011/74‑1 (MQ=255) cTGCGGCTTCATCGACCCCCAGCCAGTC‑GGCTTGCTCA‑A‑CGCCGGCTAACAAGGCAT‑CCGGCGCAATAAAACGaaa > 2:1409819/1‑76 (MQ=255) cTGCGGCTTCATCGACCACCAGCCAGTC‑GGCTTGCTCA‑T‑CGCTGGCTAACAAGGCAT‑CCGGCGCCATAAAGCGaaa > 1:4119050/1‑76 (MQ=255) tGCGGCTTCATCGACCACCAGCCAGTC‑GGCTTGCTCA‑T‑CGCTGGCTAACAAGGCAT‑CCGGCGCAATAAAGCGAAAc > 1:3434732/1‑76 (MQ=255) gCGGCTTCATCGACCACCAGCCAGTC‑GGCTTGCTCA‑T‑CGCTGGCTAACAAGGCAT‑CCGGCGCAATAAAGCGAAACt > 1:6597948/1‑76 (MQ=255) ggccTTCATCGCCCCCCAGCCAGTC‑GGCTTGCTCA‑T‑CGCTGGCTAACAAGGCAT‑CCGGCGCAATAAAGCGAAACtt < 1:4865551/73‑1 (MQ=255) ggCTTCATCGACCACCAGCCAGTC‑GGCTTGCTCA‑T‑CGCTGGCTAACAAGGCAT‑CCGGCGCAATAAAGCGAAAATtc > 2:2573481/1‑76 (MQ=255) gCTTCCTCGACCCCCAGCCAGTC‑GGCTTGCTCA‑T‑CGCTGGCTAACAAGGCAT‑CCGGCGCAATAAAGCGAAACTtct < 1:2338893/76‑1 (MQ=255) gCTTCATCGACCACCAGCCAGTC‑GGCTTGCTCA‑T‑CGCTGGCTAACAAGGCAT‑CCGGCGCAATAAAGCGAAACTtct > 1:2600815/1‑76 (MQ=255) gCTTCATCGACCACCAGCCAGTC‑GGCTTGCTCA‑T‑CGCTGGCTAACAAGGCAT‑CCGGCGCAATAAAGCGAAACTtcc > 1:2914042/1‑75 (MQ=255) gCTTCATCGACCACCAGCCAGTC‑GGCTTGCTCA‑T‑CGCTGGCTAACAAAGCAT‑CCGGCGCAATAAAGCGAAACtttt > 1:282511/1‑74 (MQ=255) ttCATCGACCACCAGCCAGTC‑GGCTTGCTCA‑T‑CGCTGGCTAACAAGGCAT‑CCGGCGCAATAAAGCGAAACTTCTCg > 1:3542424/1‑76 (MQ=255) ttCATCGACCACCAGCCAGTC‑GGCTTGCTCA‑T‑CGCTGGCTAACAAGGCAT‑CCGGCGCAATAAAGCGAAACTTCTCg > 1:1690569/1‑76 (MQ=255) ttCATCGACCACCAGCCAGTC‑GGCTTGCTCA‑T‑CGCTGGCTAACAAGGCAT‑CCGGCGCAATAAAGCGAAACTTCTCg > 2:562255/1‑76 (MQ=255) ttCATCGACCACCAGCCAGCC‑GGCTTGCTCA‑T‑CGCTGGCTAACAAGGCAT‑CCGGTGCAATAAAGCGAAACTTCTCg > 2:3817125/1‑76 (MQ=255) tCATCGACCACCAGCCAGTC‑GGCTTGCTCA‑T‑CGCTGGCTAACAAGGCAT‑CCGGCGCAATAAAGCGAAACTTCTCGc > 2:1833625/1‑76 (MQ=255) cATCGCCCCCCAGCCAGTC‑GGCTTGCTCA‑T‑CGCTGGCTAACAAGGCAT‑CCGGCGCAATAAAGCGAAACTTCTCGcc < 2:6216216/76‑1 (MQ=255) aTCGACCACCAGCCAGTC‑GGCTTGCTCA‑T‑CGCTGGCTAACAAGGCAT‑CCGGCGCAATAAAGCGAAACTTCTCGccc > 1:5001424/1‑76 (MQ=255) cGACCACCAGCCAGTC‑GGCTTGCTCA‑T‑CGCTGGCTAACAAGGCAT‑CCGGCGCAATAAAGCGAAACTTCTCGCCCGc > 1:5708552/1‑76 (MQ=255) gACCACCAGCCAGTC‑GGCTTGCTCA‑T‑CGCTGGCTAACAAGGCAT‑CCGGCGCAATAAAGCGAAACTTCTCGCCCGCa > 1:5790562/1‑76 (MQ=255) accaccAGCCAGTC‑GGCTTGCTCA‑T‑CGCTGGCTAACAAGGCAT‑CCGGCGCAATAAAGCGAAACTTCTCGCCCGcca > 1:6671342/1‑74 (MQ=255) accaccAGCCAGTC‑GGCTTGCTCA‑T‑CGCTGGCTAACAAGGCAT‑CCGGCGCAATAAAGCGAAACTTCTCGCCCGCaa > 1:5895568/1‑76 (MQ=255) caccaGCCAGTC‑GGCTTGCTCA‑T‑CGCTGGCTAACAAGGCAT‑CCGGCGCAATAAAGCGAAACTTCTCGCCCGCAAAt > 1:3210616/1‑76 (MQ=255) accaGCCAGTC‑GGCTTGCTCA‑T‑CGCTGGCTAACAAGGCAT‑CCGGCGCAATAAAGCGAAACTTCTCGCCCGCCAAtt > 2:1889831/1‑76 (MQ=255) | GTTAACAACGTTCGAGGAAAAC‑GCGATACCAGTTGATGCAACAATGGCGCAGGTATGGCTGCGGCTTCATCGACCACCAGCCAGTC‑GGCTTGCTCA‑T‑CGCTGGCTAACAAGGCAT‑CCGGCGCAATAAAGCGAAACTTCTCGCCCGCAAATT > REL606/2515834‑2515984 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 20 ≤ ATCG/ATCG < 28 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |