Predicted mutation | ||||||
---|---|---|---|---|---|---|
evidence | seq id | position | mutation | annotation | gene | description |
RA | REL606 | 3,264,050 | (A)7→6 | intergenic (‑101/‑147) | greA ← / → dacB | transcription elongation factor GreA/D‑alanyl‑D‑alanine carboxypeptidase |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | REL606 | 3,264,044 | 0 | A | . | 88.0% | 86.0 / 7.5 | 25 | intergenic (‑95/‑153) | greA/dacB | transcription elongation factor GreA/D‑alanyl‑D‑alanine carboxypeptidase |
Reads supporting (aligned to +/- strand): ref base A (0/3); new base . (10/12); total (10/15) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 2.50e-01 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.25e-01 |
TGATGGGGCAAGTCTTCACCCACCCCGGCTGTTCGGCCTGCCCGGCATAACTCCTCACCGGGACCAGAAGCAAAAAAATACCGACCCGGGTACAAGTCCCAGGTCAGCTACAATTCACATTTTGATAGTCATTTTA > REL606/3263973‑3264108 | tGATGGGGCAAGTCTTAACCCACCCCGGCTGTTCGGCCTGCCCGGCATAACTCCTAACCGGGACCAGAAGCaaaaa > 1:1759608/1‑76 (MQ=255) gATGGGGCAAGTCTTCACCCACCCCGGCTGTTCGGCCTGCCCGGCATAACTCCTCACCGGGACCAGAAGCaaaaaa > 1:782447/1‑76 (MQ=255) gATGGGGCAAGTCTTCACCCACCCCGGCTGTTCGGCCTGCCCGGCATAACTCCTCACCGGGACCAGAAGCaaaaaa < 1:975699/76‑1 (MQ=255) tGGGGCAAGTCTTCACCCACCCCGGCTGTTCGGCCTGCCCGGCATAACTCCTCACCGGGACCAGAAGCAAAAAAAt < 2:829607/76‑1 (MQ=255) ggggCAAGTCTTCACCCACCCCGGCTGTTCGGCCTGCCCGGCATAACTCCTCACCGGGCCCAGAAGCaaaaaatac > 1:1397532/1‑73 (MQ=255) ggCAAGTCTTCACCCACCCCGGCTGTTCGGCCTGCCCGGCATAACTCCTCACCGGGACCAGAAGC‑AAAAAATACCg > 1:316404/1‑76 (MQ=255) aaGTCTTCACCCACCCCGGCTGTTCGGCCTGCCCGGCATAACTCCTCACCGGGACCAGAAGC‑AAAAAATACCGAcc > 2:447779/1‑76 (MQ=255) tCTTCACCCCCCCCGGCTGTTCGGCCTGCCCGGCATAACTCCTCACCGGGACCAGAAGC‑AAAAAATACCGACCCgg < 2:955496/76‑1 (MQ=255) tcccccacccCGGCTGTTCGGCCTGCCCGGCATAACTCCTCACCGGGACCAGAAGC‑AAAAAATACCGACCCGGGTa < 2:1092194/73‑1 (MQ=255) cacccaccCCGGCTGTTCGGCCTGCCCGGCATACCTCCTCACCGGGACCAGAAGC‑AAAAAATACCGACCCGGGTAc < 1:1030436/76‑1 (MQ=255) ccacccCGGCTGTTCGGCCTGCCCGGCATAACTCCTCACCGGGACCAGACGC‑AAAAAATACCGACCCGGGTACAAg > 2:1141442/1‑76 (MQ=255) acccCGGCTGTTCGGCCCGCCCGGCATAACTCCTCACCGNGACCATACGC‑AAAAANTACCGATCCGGGTACAAGTc > 2:1747524/1‑76 (MQ=255) cccGGCTGTTCGGCCTGCCCGGCATAACTCCTCACCGGGACCAGAAGC‑AAAAAATACCGACCCGGGTACAAGTccc > 1:1574343/1‑76 (MQ=255) cccGGCTGTTCGGCCTGCCCGGCATAACTCCTCACCGGGACCAGAAGC‑AAAAAATACCGACCCGGGTACAAGTccc < 2:1530855/76‑1 (MQ=255) ccGCCTGTTCGGCCTGCCCGGCATACCCCCTCACCGGGACCAGAAGC‑AAAAAATACCGACCCGGGTACAAGTCCCa < 2:640522/76‑1 (MQ=255) cTGTTCGGCCTGCCCGGCATAACTCCTCACCGGGACCAGAAGC‑AAAAAATACCGACCCGGGTACAAGTCCCAGGTc > 1:787839/1‑76 (MQ=255) cTGTTCGGCCTGCCCGGCATAACTCCTCACCGGGACCAGAAGC‑AAAAAATACCGACCCGGGTACAAGTCCCAGGTc < 1:854996/76‑1 (MQ=255) gTTCGGCCTGCCCGGCATAACTCCTCACCGGGACCAGAAGC‑AAAAAATACCGACCCGGGTACAAGTCCCAGGTCAg > 2:1207016/1‑76 (MQ=255) gTTCGGCCTGCCCGGCATAACTCCTCACCGGGACCAGAAGC‑AAAAAATACCGACCCGGGTACAAGTCCCAGGTCAg > 2:1352749/1‑76 (MQ=255) tCGGCCTGCCCGGCATACCTCCTCACCGGGACCAGAAGC‑AAAAAATACCGACCCGGGTACAAGTCCCAGGTCAGCt < 1:1293912/76‑1 (MQ=255) cggccggcCCGGCATACCTCCTCCCCGGGACCAGAAGC‑AAGAAATACCGACCCGGGTACAAGTCCCAGGTCAGCTa < 2:1097860/76‑1 (MQ=255) cGGCCTGCCCGGCATAACTCCTCACCGGGACCAGAAGCAAAAAAATACCGACCCGGGTACAAGTCCCAGGTCAGCt < 2:1462279/76‑1 (MQ=255) gCCTGCCCNGCAGAACTCCTCACCGGGCCCAGATGC‑AAAAAATACCGACCCGGGTACAAGTCCCAGGTCAGCTACa < 2:292460/76‑1 (MQ=255) cTGCCCGGCATAACTCCTCACCGGGACCAGAAG‑AAAAAAATACCGACCCGGGTACAAGTCCCAGGTCAGCTACAAt < 1:1564157/76‑1 (MQ=255) cGGCATAACTCCTCACCGGGACCAGAAGC‑AAAAAATACCGACCCGGGTACAAGTCCCAGGTCAGCTACAATTcaca > 2:1221452/1‑76 (MQ=255) gCATAACTCCTCACCGGGACCAGAAGC‑AAAAAATACCGACCCGGGTACAAGTCCCAGGTCAGCTACAATTCACAtt < 2:441561/76‑1 (MQ=255) aTAACTCCTCACCGGGACCAGAAGC‑AAAAAATACCGACCCGGGTACAAGTCCCAGGTCAGCTACAATTCACAtttt > 1:499289/1‑76 (MQ=255) cctcACCGGGACCAGAAGC‑AAAAGATACCAACCCGGTAACACGTCCCCGGGCAGTTACATCTCACCTTTTGATAGt > 2:12425/1‑76 (MQ=255) cccccGGGACCAGAAGC‑AAAAAATACCGACCCGGGTACAAGTCCCAGGTCAGCTACAATTCACATTTTGATAGTCa < 1:1275197/73‑1 (MQ=255) ccGGGCCCAGAAGCAAAAANTNCTCGACCNGTGTTCAAGTCCCGGGTCAGCNACNAGTCACATTTTGATAGTCAtt < 2:262010/76‑1 (MQ=255) gggACCAGAAGC‑AAAAAATACCGACCCGGGTACAAGTCCCAGGTCAGCTACAATTCACATTTTGATAGTCATTTTa < 2:487563/76‑1 (MQ=255) gggACCAGAAGC‑AAAAAATACCGACCCGGGTACAAGTCCCAGGTCAGCTACAATTCACATTTTGATAGTCATTTTa < 2:1033290/76‑1 (MQ=255) | TGATGGGGCAAGTCTTCACCCACCCCGGCTGTTCGGCCTGCCCGGCATAACTCCTCACCGGGACCAGAAGCAAAAAAATACCGACCCGGGTACAAGTCCCAGGTCAGCTACAATTCACATTTTGATAGTCATTTTA > REL606/3263973‑3264108 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 18 ≤ ATCG/ATCG < 26 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |