Predicted mutation | ||||||
---|---|---|---|---|---|---|
evidence | seq id | position | mutation | annotation | gene | description |
RA | NC_000913 | 4,268,487 | (T)7→8 | intergenic (+180/+322) | tyrB → / ← yjbS | tyrosine aminotransferase, tyrosine‑repressible, PLP‑dependent/uncharacterized protein |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_000913 | 4,268,480 | 1 | . | T | 97.2% | 112.0 / ‑3.4 | 36 | intergenic (+173/+329) | tyrB/yjbS | tyrosine aminotransferase, tyrosine‑repressible, PLP‑dependent/uncharacterized protein |
Reads supporting (aligned to +/- strand): ref base . (1/0); new base T (18/17); total (19/17) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 1.00e+00 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.86e-01 |
TGAAATTAAACTGTCGTCGCTTTCACTCTTTCTTTATAGATGATTTTTTTGATGCCATCGTTCTACGTGAGAGATAATAAACGTTGTTAGTTCTTTTATTGTTAAGCTTATCCCAATTATCTGGAATTCCTTATCCTG‑TTTTTTTGGGTGGAGTAATAATTTCTGCACATCTAACCATTATTAAATTTAATAATAAGTATCGGGTTTGTTACGATAGAATATTTCTATATCCTCAAGAATATGCTTTGCCTTATTGCCTGTGACTTTCTGA > NC_000913/4268343‑4268613 | tGAAATTAAACTGTCGTCGCTTTCACTCTTTCTTTATAGATGATTTTTTTGATGCCATCGTTCTACGTGAGAGATAATAAACGTTGTTAGTTCTTTTATTGTTAAGCTTATCCCAATTATCTGGAATTCCTTATCCTG‑t < 2:137080/139‑1 (MQ=255) cTGTCGTCGCTTTCACTCTTTCTTTATAGATGATTTTTTTGATGCCATCGTTCTACGTGAGAGATAATAAACGTTGTTAGTTCTTTTATTGTTAAGCTTATCCCAATTATCTGGAATTCCTTATCCTG‑ttttttttggg < 2:598022/139‑5 (MQ=255) gtcgCTTTCACTCTTTCTTTATAGATGATTTTTTTGATGCCATCGTTCTACGTGAGAGATAATAAACGTTGTTAGTTCTTTTATTGTTAAGCTTATCCCAATTATCTGGAATTCCTTATCCTGTTTTTTTTGGGTGGAg > 2:150221/1‑139 (MQ=255) cACTCTTTCTTTATAGATGATTTTTTTGATGCCATCGTTCTACGTGAGAGATAATAAACGTTGTTAGTTCTTTTATTGTTAAGCTTATCCCAATTATCTGGAATTCCTTATCCTGTTTTTTTTGGGTGGAGtaat > 1:941709/1‑135 (MQ=255) cACTCTTTCTTTATAGATGATTTTTTTGATGCCATCGTTCTACGTGAGAGATAATAAACGTTGTTAGTTCTTTTATTGTTAAGCTTATCCCAATTATCTGGAATTCCTTATCCTGTTTTTTTTGGGTGGAGtaat < 2:941709/135‑1 (MQ=255) ttctttATAGATGATTTTTTTGATGCCATCGTTCTACGTGAGAGATAATAAACGTTGTTAGTTCTTTTATTGTTAAGCTTATCCCAATTATCTGGAATTCCTTATCCTGTTTTTTTTGGGTGGAGTAATAATTTCTGca < 2:398629/139‑1 (MQ=255) tataGATGATTTTTTTGATGCCATCGTTCTACGTGAGAGATAATAAACGTTGTTAGTTCTTTTATTGTTAAGCTTATCCCAATTATCTGGAATTCCTTATCCTGTTTTTTTTGGGTGGAGTAATAATTTCTGca > 1:1205662/1‑134 (MQ=255) tataGATGATTTTTTTGATGCCATCGTTCTACGTGAGAGATAATAAACGTTGTTAGTTCTTTTATTGTTAAGCTTATCCCAATTATCTGGAATTCCTTATCCTGTTTTTTTTGGGTGGAGTAATAATTTCTGca < 2:1205662/134‑1 (MQ=255) taGATGATTTTTTTGATGCCATCGTTCTACGTGAGAGATAATAAACGTTGTTAGTTCTTTTATTGTTAAGCTTATCCCAATTATCTGGAATTCCTTATCCTG‑ttttttttggg < 1:1192697/113‑5 (MQ=255) taGATGATTTTTTTGATGCCATCGTTCTACGTGAGAGATAATAAACGTTGTTAGTTCTTTTATTGTTAAGCTTATCCCAATTATCTGGAATTCCTTATCCTG‑ttttttttggg > 2:1192697/1‑109 (MQ=255) gatTTTTTTGATGCCATCGTTCTACGTGAGAGATAATAAACGTTGTTAGTTCTTTTATTGTTAAGCTTATCCCAATTATCTGGAATTCCTTATCCTGTTTTTTTTGGGTGGAGTAATAATTTCTg < 2:376284/125‑1 (MQ=255) gatTTTTTTGATGCCATCGTTCTACGTGAGAGATAATAAACGTTGTTAGTTCTTTTATTGTTAAGCTTATCCCAATTATCTGGAATTCCTTATCCTGTTTTTTTTGGGTGGAGTAATAATTTCTg > 1:376284/1‑125 (MQ=255) gatTTTTTTGATGCCATCGTTCTACGTGAGAGATAATAAACGTTGTTAGTTCTTTTATTGTTAAGCTTATCCCAATTATCTGGAATTCCTTATCCTGTTTTTTTTGGGTGGAGTAATAATTTCTGCACATCTAACCatt > 1:377272/1‑139 (MQ=255) cATCGTTCTACGTGAGAGATAATAAACGTTGTTAGTTCTTTTATTGTTAAGCTTATCCCAATTATCTGGAATTCCTTATCCTGTTTTTTTTGGGTGGAGTAATAATTTCTGCACATCTAACCATTATTAAATTtaataa < 1:152514/139‑1 (MQ=255) agaTAATAAACGTTGTTAGTTCTTTTATTGTTATGCTTAGCCCAATTATCTGGAATTCCTTATCCTG‑ttttttttgg < 2:914606/77‑4 (MQ=255) agaTAATAAACGTTGTTAGTTCTTTTATTGTTAAGCTTATCCCAATTATCTGGAATTCCTTATCCTGTTTTTTTTGGGTGGAGTAATAATTTCTGc < 1:382195/96‑1 (MQ=255) agaTAATAAACGTTGTTAGTTCTTTTATTGTTAAGCTTATCCCAATTATCTGGAATTCCTTATCCTGTTTTTTTTGGGTGGAGTAATAATTTCTGc > 2:382195/1‑96 (MQ=255) agaTAATAAACGTTGCTAGTTCTTTTATTGTTAAGCTTATCCCAATTATCTGGAATTCCTTATCCTG‑TTTTTTTtgg > 1:914606/1‑77 (MQ=255) aataaACGTTGTTAGTTCTTTTATTGTTAAGCTTATCCCAATTTTCTGGAATTCCTTATCCTGTTTTTTTTGGGTGGAGTAATAATTTCTGCACATCTAACCATTATTAAATTTAATAATAAGTATCGGGTTTGt < 2:459389/135‑1 (MQ=255) aataaACGTTGTTAGTTCTTTTATTGTTAAGCTTATCCCAATTTTCTGGAATTCCTTATCCTGTTTTTTTTGGGTGGAGTAATAATTTCTGCACATCTAACCATTATTAAATTTAATAATAAGTATCGGGTTTGt > 1:459389/1‑135 (MQ=255) aataaACGTTGTTAGTTCTTTTATTGTTAAGCTTATCCCAATTATCTGGAATTCCTTATCCTGTTTTTTTTGGGTGGAGTAATAATTTCTGCACATCTAACCATTATTAAATTtaat < 1:921895/117‑1 (MQ=255) aataaACGTTGTTAGTTCTTTTATTGTTAAGCTTATCCCAATTATCTGGAATTCCTTATCCTGTTTTTTTTGGGTGGAGTAATAATTTCTGCACATCTAACCATTATTAAATTtaat > 2:921895/1‑117 (MQ=255) cGTTGTTAGTTCTTTTATTGTTAAGCTTATCCCAATTATCTGGAATTCCTTATCCTGTTTTTTTTGGGTGGAGtaat > 1:827049/1‑77 (MQ=255) cGTTGTTAGTTCTTTTATTGTTAAGCTTATCCCAATTATCTGGAATTCCTTATCCTGTTTTTTTTGGGTGGAGtaat < 2:827049/77‑1 (MQ=255) cTTTTATTGTTAAGCTTATCCCAATTATCTGGAATTCCTTATCCTGTTTTTTTTGGGTGGAGTAATAATTTCTGCACATCTAACCa > 1:827778/1‑86 (MQ=255) cTTTTATTGTTAAGCTTATCCCAATTATCTGGAATTCCTTATCCTGTTTTTTTTGGGTGGAGTAATAATTTCTGCACATCTAACCa < 2:827778/86‑1 (MQ=255) tttATTGTTAAGCTTATCCCAATTATCTGGAATTCCTTATCCTGTTTTTTTTGGGTGGAGTAATAATTTCTGCACATCTAACCATTATTAAATTTAATAATAAGTATCGGGTTTGTTACGATAGAATATTTCTATATcc > 1:576255/1‑139 (MQ=255) ttATTGTTCAGTTTATCCCAATTATCTGGAATTCCTTATCCTGTTTTTTTTGGGTGGAGTAATAATTTCTGCACATCTAACCATTATTAAATTTAATAATAAGTATCGGGTTTGTTACGATAGAATATTTCTATATCCt < 1:831385/139‑1 (MQ=255) ttATTGTTAAGCTTATCCCAATTATCTGGAATTCCTTATCCTGTTTTTTTTGGGTGGAGTAATAATTTCTGCACATCTAACCATTATTAAATTTAATAATAAGTATCGGGTTTGTTACGATAGAATATTTCTATATCCt < 2:576255/139‑1 (MQ=255) aaGCTTATCCCAATTATCTGGAATTCCTTATCCTGTTTTTTTTGGGTGGAGTAATAATTTCTGCACATCTAACCATTATTAAATTTAATAATAAGTATCGGGTTTGTTACGATAGAATATTTCTATATCCTCAAGAata > 1:503997/1‑139 (MQ=255) aTCCCAATTATCTGGAATTCCTTATCCTGTTTTTTTTGGGTGGAGTAATAATTTCTGCACATCTAACCATTATTAAATTtaataa > 2:418446/1‑85 (MQ=255) aTCCCAATTATCTGGAATTCCTTATCCTGTTTTTTTTGGGTGGAGTAATAATTTCTGCACATCTAACCATTATTAAATTtaataa < 1:418446/85‑1 (MQ=255) aTTATCTGGAATTCCTTATCCTGTTTTTTTTGGGTGGAGTAATAATTTCTg < 2:313942/51‑1 (MQ=255) aTTATCTGGAATTCCTTATCCTGTTTTTTTTGGGTGGAGTAATAATTTCTg > 1:313942/1‑51 (MQ=255) ttATCTGGAATTCCTTATCCTGTTTTTTTTGGGTGGAGTaataat > 2:1329087/1‑45 (MQ=255) ttATCTGGAATTCCTTATCCTGTTTTTTTTGGGTGGAGTaataat < 1:1329087/45‑1 (MQ=255) tGGAATTCCTTATCCTGTTTTTTTTGGGTGGAGTAATAATTTCTGCACATCt > 2:301496/1‑52 (MQ=39) tGGAATTCCTTATCCTGTTTTTTTTGGGTGGAGTAATAATTTCTGCACATCt < 1:301496/52‑1 (MQ=39) cTTATCCTGTTTTTTTTGGGTGGAGTAATAATTTCTg < 2:177035/37‑1 (MQ=255) cTTATCCTGTTTTTTTTGGGTGGAGTAATAATTTCTg > 1:177035/1‑37 (MQ=255) tCCTGTTTTTTTTGGGTGGAGTAATAATTTCTGCACATCTAACCATTATTAAATTTAATAATAAGTATCGGGTTTGTTACGATAGAATATTTCTATATCCTCAAGAATATGCTTTGCCTTATTGCCTGTGACTTTCTga > 1:197515/1‑139 (MQ=255) | TGAAATTAAACTGTCGTCGCTTTCACTCTTTCTTTATAGATGATTTTTTTGATGCCATCGTTCTACGTGAGAGATAATAAACGTTGTTAGTTCTTTTATTGTTAAGCTTATCCCAATTATCTGGAATTCCTTATCCTG‑TTTTTTTGGGTGGAGTAATAATTTCTGCACATCTAACCATTATTAAATTTAATAATAAGTATCGGGTTTGTTACGATAGAATATTTCTATATCCTCAAGAATATGCTTTGCCTTATTGCCTGTGACTTTCTGA > NC_000913/4268343‑4268613 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |