Predicted mutation | ||||||
---|---|---|---|---|---|---|
evidence | seq id | position | mutation | annotation | gene | description |
RA | REL606 | 4,179,757 | T→G | intergenic (+51/‑136) | yjaG → / → hupA | hypothetical protein/HU, DNA‑binding transcriptional regulator, alpha subunit |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | REL606 | 4,179,757 | 0 | T | G | 95.8% | 185.3 / ‑3.9 | 48 | intergenic (+51/‑136) | yjaG/hupA | hypothetical protein/HU, DNA‑binding transcriptional regulator, alpha subunit |
Reads supporting (aligned to +/- strand): ref base T (0/2); new base G (26/20); total (26/22) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 2.05e-01 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.91e-01 |
TGGTATAATTTTTCAGCAATAAGACCAGAAAACGTGATTTAACGCCTGATTTGTCGTACCTGGAGTCTTCCCTTTCGCCCCCCGTCTGGTCTACATTTGGGGGGCGAAAAAAAGTGGCTATCGGTGCGTGTATGCAGGAGAGTG > REL606/4179685‑4179828 | tGGTATAATTTTTCAGCAATAAGACCAGAAAACGTGATTTAACTCCTGATTTGTCGTACCTGGAGTCTTCCCGTTc < 2:672598/76‑1 (MQ=255) tGGTATAATTTTTCAGCAATAAGACCAGAAAACGTGATTTAACGCCTGATTTGTCGTACCTGGAGTCTTCCCTTTc < 2:572493/76‑1 (MQ=255) tgTATAATTTTTCAGCAATAAGACCAGAAAACGTGATTTAACGCCTGATTTGTCGTACCTGGAGTCTTCCCGTTCg > 2:709150/2‑76 (MQ=255) tgTATAATTTTTCAGCAATAAGACCAGAAAACGTGATTTAACGCCTGATTTGTCGTACCTAGAGTCTTCCCGTTCg > 2:841717/2‑76 (MQ=255) gTATAATTTTTCAGCAATAAGACCAGAAAACGTGATTTAACGCCTGATTTGTCGTACCTGGAGTCTTCCCGTTCGc > 2:1564301/1‑76 (MQ=255) ataATTTTTCAGCAATAAGACCAGAAAACGTGATTTAACGCCTGATTTGTCGTACCTGGAGTCTTCCCGTTCGccc > 1:1037448/1‑76 (MQ=255) aagttttCAGCAATAAGACCAGAAAACGTGATGTAACGCCTGATTTGTCGTACCTGGAGTCTTCCCGTTCGccccc < 2:1260642/73‑1 (MQ=255) aTTTTTCAGCAATAAGACCAGAAAACGTGATTTAACGCCTGATTTGTCGTACCTGGAGTCTTCCCGTTCGcccccc > 1:1071106/1‑76 (MQ=255) tttttCAGCAATAANACCAGAAAAAGTGATTTAACGCCTGATTTGTCGTACCTGGAGTCTTCCCGTTCGCCCCCCg > 1:238/1‑76 (MQ=255) tttttCAGCAATAAGACCAGAAAACGTGATTTAACGCCTGATTTGTCGTACCTGGAGTCTTCCCGTTCGCCCCCCg < 2:892442/76‑1 (MQ=255) cAGCAATAAGACCAGAAAACGTGATTTAACGCCTGATTTGTCGTACCTGGAGTCTTCCCGTTCGCCCCCCGCCTgg > 1:843703/1‑76 (MQ=255) cAATAAGACCAGAAAACGTGATTTAACGCCTGATTTGTCGTACCTGGAGTCTTCCCGTTCGCCCCCCGTCTGGTCt > 2:803379/1‑76 (MQ=255) cAATAAGACCAGAAAACGTGATTTAACGCCTGATTTGTCGTACCTGGAGTCTTCCCGTTCGCCCCCCGTCTGGTCt > 2:444524/1‑76 (MQ=255) cAATAAGACCAGAAAACGTGATTTAACGCCTGATTTGTCGTACCTGGAGTCTTCCCGTTCGCCCCCCGCCTGGTCt > 2:1179685/1‑76 (MQ=255) aaTAAGACCAGAAAACGTGATTTAACGCCTGATTTGTCGTACCTGGAGTCTTCCCGTTCGCCCCCCGTCTGGTCTa < 1:1267844/76‑1 (MQ=255) aaTAAGACCAGAAAACGTGATTTAACGCCTGATTTGTCGTACCTGGAGTCTTCCCGTTCGCCCCCCGTCTGGTCTa > 2:1604013/1‑76 (MQ=255) aaTAAGACCAGAAAACGTGATTTAACGCCTGATTTGTCGTACCTGGAGTCTTCCCGTTCGCCCCCCGTCTGGTCTa > 2:1542727/1‑76 (MQ=255) gACCAGAAAACGTGATTTAACGCCTGATTTGTCGTACCTGGAGTCTTCCCGTTCGCCCCCCGTCTGGTCTACAttt < 1:928971/76‑1 (MQ=255) gACCAGAAAACGTGATTTAACGCCTGATTTGTCGTACCTGGAGTCTTCCCGTTCGCCCCCCGTCTGGTCTACAttt < 2:641566/76‑1 (MQ=255) ccAGAAAACGTGATTTAACGCCTGATTTGTCGTACCTGGAGTCTTCCCGTTCGCCCCCCGTCTGGTCTACATTTgg > 2:578278/1‑76 (MQ=255) aaaaCGTGATTTAACGCCTGATTTGTCGTACCTGGAGTCTTCCCGTTCGCCCCCCGTCTGGTCTACATTTgggggg < 2:1544118/76‑1 (MQ=255) aaCGTGATTTAACGCCTGATTTGTCGTACCTGGAGTCTTCCCGTTCGCCCCCCGTCTGGTCTACATTTGGGGGGCg > 2:329953/1‑76 (MQ=255) ttAACGCCTGATTTGTCGTACCTGGAGTCTTCCCGTTCGCCCCCCGTCTGGTCTACATTTGGGGGGCGAAAAAAAg < 2:1015829/76‑1 (MQ=255) gCCTGATTTGTCGTACCTGGAGTCTTCCCGTTCGCCCCCCGTCTGGTCTACATTTGGGGGGCGAAAAAAAGTGGCt < 2:585638/76‑1 (MQ=255) cttaTTTGTCGTCCCTGGACTCTCCCCTTCCCCCCCCCGTCTGGTCTACATTTGGGGGGCGAAAAAAAGTGGCTAt < 1:775343/73‑1 (MQ=255) cTGATTTGTCGTACCTGGAGTCTTCCCGTTCGCCCCCCGTCTGGTCTACATTTGGGGGGCGAAAAAAAGTGGCTAt > 2:166050/1‑76 (MQ=255) cTGATTTGTCGTACCTGGAGTCTTCCCGTTCGCCCCCCGTCTGGTCTACATTTGGGGGGCGAAAAAAAGTGGCTAt < 1:1632251/76‑1 (MQ=255) cTGATTTGTCGTACATGCACTCATCCCGTTCGCCCCCCGTCTGGTCTACATTTGGTGGGCGAAAAAAAGTGGCTAt < 1:86267/76‑1 (MQ=255) gATTTGTCGTACCTGGAGTCTTCCCGTTCGCCCCCCGTCTGGTCTACATTTGGGGGGCGAAAAAAAGTGGCTATCg > 2:1151806/1‑76 (MQ=255) gATTTGTCGTACCTGGAGTCTTCCCGTTCGCCCCCCGTCTGGTCTACATTTGGGGGGCGAAAAAAAGTGGCTATCg > 2:571435/1‑76 (MQ=255) aTTTGTCGTACCTGGAGTCTTCCCGTTCGCCCCCCGTCTGGTCTACATTTGGGGGGCGAAAAAAAGTGGCTATCgg < 1:1634684/76‑1 (MQ=255) ttGTCGTACCTGGAGTCTTCCCGTTCGCCCCCCGTCTGGTCTACATTTGGGGGGCGAAAAAAAGTGGCTATCGGTg > 1:1201814/1‑76 (MQ=255) gTACCTGGAGTCTTCCCGTTCGCCCCCCGTCTGGTCTACATTTGGGGGGCGAAAAAAAGTGGCTATCGGTGCgtgt > 1:1027441/1‑76 (MQ=255) tACCTGGAGTCTTCCCGTTCGCCCCCCGTCTGGTCTACATTTGGGGGGCGAAAAAAAGTGGCTATCGGTGCGTGTa < 2:1326326/76‑1 (MQ=255) tACCTGGAGTCTTCCCGTTCGCCCCCCGTCTGGTCTACATTTGGGGGGCGAAAAAAAGTGGCTATCGGTGCGTGTa < 2:1071391/76‑1 (MQ=255) aCCTGGAGTCTTCCCGTTCGCCCCCCGTCTGGTCTACATTTGGGGGGCGAAAAAAAGTGGCTATCGGTGCGTGTAt < 2:298297/76‑1 (MQ=255) aCCTGGAGTCTTCCCGTTCGCCCCCCGTCTGGTCTACATTTGGGGGGCGAAAAAAAGTGGCTATCGGTGCGTGTAt < 2:1085644/76‑1 (MQ=255) cTGGAGTCTTCCCGTTCGCCCCCCGTCTGGTCTACATTTGGGGGGCGAAAAAAAGTGGCTATCGGTGCGTGTATGc < 2:434466/76‑1 (MQ=255) cTGGAGTCTTCCCGTTCGCCCCCCGACTGATCTACTTTTGCCCGGCGATAAGAAGTCACTTTTCTTGNGTGTATGc > 2:279669/1‑76 (MQ=255) tGGNGTCTTCCCGTTCGCCCCCCGTCTGGTCTTCATTTGGGGGGCGAAAAAAAGTGGCTATCGGGGCGTGTATGCa > 1:1532631/1‑76 (MQ=255) tGGAGTCTTCCCGTTCGCCCCCCGTCTGGTCTACATTTGGGGGGCGAAAAAAAGTGGCTATCGGTGCGTGTATGCa < 1:1204729/76‑1 (MQ=255) tGGAGTCTTCCCGTTCGCCCCCCGTCTGGTCTACATTTGGGGGGCGAAAAAAAGTGGCTATCGGTGCGTGTATGCa > 1:749321/1‑76 (MQ=255) gAGTCTTCCCTTTCGCCCCCCGTCTGGTCTACATTTGGGGGGCGAAAAAAAGTGGCTATCGGTGCGTGTATGCAgg < 2:1555443/76‑1 (MQ=255) gAGTCTCCCCGTTCGCCCCCCGTCTGGTCTACATTTGGGGGGCGAAAAAAAGTGGCTATCGGTGCGTGTATGCAgg < 2:832913/76‑1 (MQ=255) gTCTTCCCGTTCGCCCCCCGTCTGGTCTACATTTGGGGGGCGAAAAAAAGTGGCTATCGGTGCGTGTATGCAGgag < 2:201271/76‑1 (MQ=255) gTCTTCCCGTTCGCCCCCCGTCTGGTCTACATTTGGGGGGCGAAAAAAAGTGGCTATCGGTGCGTGTATGCAGgag < 1:983880/76‑1 (MQ=255) cTTCCCGTTCGCCCCCCGTCTGGTCTACATTTGGGGGGCGAAAAAAAGTGGCTATCGGTGCGTGTATGCAGGagag > 2:121139/1‑76 (MQ=255) cTTCCCGTTCGCCCCCCGTCTGGTCTACATTTGGGGGGCGAAAAAAAGTGGCTATCGGTGCGTGTATGCAGGagag > 2:582459/1‑76 (MQ=255) ttCCCGTTCGCCCCCCGTCTGGTCTACATTTGGGGGGCGAAAAAAAGTGGCTATCGGTGCGTGTATGCAGGAGAGt < 2:849638/76‑1 (MQ=255) tCCCGTTCGCCCCCCGTCTGGTCTACATTTGGGGGGCGAAAAAAAGTGGCTATCGGTGCGTGTATGCAGGagaggg > 2:324834/1‑74 (MQ=255) tCCCGTTCGCCCCCCGTCTGGTCTACATTTGGGGGGCGAAAAAAAGCGGCTATCGGTGCGTGTATGCAGGAGGGTg > 2:1621575/1‑76 (MQ=255) | TGGTATAATTTTTCAGCAATAAGACCAGAAAACGTGATTTAACGCCTGATTTGTCGTACCTGGAGTCTTCCCTTTCGCCCCCCGTCTGGTCTACATTTGGGGGGCGAAAAAAAGTGGCTATCGGTGCGTGTATGCAGGAGAGTG > REL606/4179685‑4179828 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 20 ≤ ATCG/ATCG < 29 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |