| Predicted mutation | |||||||
|---|---|---|---|---|---|---|---|
| evidence | seq id | position | mutation | freq | annotation | gene | description |
| RA | minE | 584,816 | G→T | 60.4% | S360Y (TCT→TAT) | cydC ← | fused cysteine transporter subunits and membrane component and ATP‑binding component of ABC superfamily |
| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | minE | 584,816 | 0 | G | T | 60.4% | 16.2 / 43.3 | 48 | S360Y (TCT→TAT) | cydC | fused cysteine transporter subunits and membrane component and ATP‑binding component of ABC superfamily |
| Reads supporting (aligned to +/- strand): ref base G (11/8); new base T (13/16); total (24/24) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 5.56e-01 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.94e-01 | |||||||||||
AGTGTTGATTTGCCGCATCCGGTTCGCCCGAGAATCGCTATATGTTCCCCGGCGTTTACCTGAAGAGAAATCCCTTTAAGTGCCTGTTGAGATTGCTCCGGATAAGTGAACTGAACATCCCGTAACGTCAGCGA > minE/584750‑584883 | aGTGTTGATTTGCCGCATCCGGTTCGCCCGAGAATCGCTATATGTTCCCCGGCGTTTACCTGAAGAGaa < 1:1100124/69‑1 (MQ=255)aGTGTTGATTTGCCGCATCCGGTTCGCCCGAGAATCGCTATATGTTCCCCGGCGTTTACCTGAAGAGaa < 1:2023650/69‑1 (MQ=255)aGTGTTGATTTGCCGCATCCGGTTCGCCCGAGAATCGCTATATGTTCCCCGGCGTTTACCTGAAGAGaa < 1:1286643/69‑1 (MQ=255) gATTTGCCGCATCCGGTTCGCCCGAGAATCGCTATATGTTCCCCGGCGTTTACCTGAAGAGAAATccc > 1:949302/1‑68 (MQ=255) tGCCGCATCCGGTTCGCCCGAGAATCGCTATATGTTCCCCGGCGTTTACCTGAagag > 1:966876/1‑57 (MQ=255) tGCCGCATCCGGTTCGCCCGAGAATCGCTATATGTTCCCCGGCGTTTACCTGAAgaga > 1:1235738/1‑58 (MQ=255) tGCCGCATCCGGTTCGCCCGAGAATCGCTATATGTTCCCCGGCGTTTACCTGAAgaga > 1:2865872/1‑58 (MQ=255) tGCCGCATCCGGTTCGCCCGAGAATCGCTATATGTTCCCCGGCGTTTACCTGAAgaga > 1:2131580/1‑58 (MQ=255) tGCCGCATCCGGTTCGCCCGAGAATCGCTATATGTTCCCCGGCGTTTACCTGAAgaga > 1:3161298/1‑58 (MQ=255) tGCCGCATCCGGTTCGCCCGAGAATCGCTATATGTTCCCCGGCGTTTACCTGAAgaga > 1:1164450/1‑58 (MQ=255) tGCCGCATCCGGTTCGCCCGAGAATCGCTATATGTTCCCCGGCGTTTACCTGAAgaga > 1:1941431/1‑58 (MQ=255) tGCCGCATCCGGTTCGCCCGAGAATCGCTATATGTTCCCCGGCGTTTACCTGAAgaga > 1:1489610/1‑58 (MQ=255) tGCCGCATCCGGTTCGCCCGAGAATCGCTATATGTTCCCCGGCGTTTACCTGAAgaga > 1:962986/1‑58 (MQ=255) cGGTTCGCCCGAGAATCGCTATATGTTCCCCGGCGTTTACCTGAAGAGAAATCCCTTTAAGTGCCTg > 1:1682363/1‑67 (MQ=255) cGGTTCGCCCGAGAATCGCTATATGTTCCCCGGCGTTTACCTGAAGAGAAATCCCTTTAAGTGCCTg > 1:866304/1‑67 (MQ=255) cGGTTCGCCCGAGAATCGCTATATGTTCCCCGGCGTTTACCTGAAGAGAAATCCCTTTAAGTGCCTg > 1:2075250/1‑67 (MQ=255) cGGTTCGCCCGAGAATCGCTATATGTTCCCCGGCGTTTACCTGAAGAGAAATCCCTTTAAGTGCCTg > 1:2914413/1‑67 (MQ=255) ggTTCGCCCGAGAATCGCAATATGTTCCCCGGCGTTTACCTGAAGATAAATCCCttt > 1:688586/1‑57 (MQ=255) ggTTCGCCCGAGAATCGCAATATGTTCCCCGGCGTTTACCTGAAGATAAATCCCttt > 1:113975/1‑57 (MQ=255) ggTTCGCCCGAGAATCGCAATATGTTCCCCGGCGTTTACCTGAAGATAAATCCCttt > 1:61813/1‑57 (MQ=255) ggTTCGCCCGAGAATCGCAATATGTTCCCCGGCGTTTACCTGAAGATAAATCCCttt > 1:2932022/1‑57 (MQ=255) ggTTCGCCCGAGAATCGCAATATGTTCCCCGGCGTTTACCTGAAGATAAATCCCttt > 1:1226913/1‑57 (MQ=255) ggTTCGCCCGAGAATCGCAATATGTTCCCCGGCGTTTACCTGAAGATAAATCCCttt > 1:2899566/1‑57 (MQ=255) ggTTCGCCCGAGAATCGCAATATGTTCCCCGGCGTTTACCTGAAGATAAATCCCttt > 1:1046937/1‑57 (MQ=255) ggTTCGCCCGAGAATCGCAATATGTTCCCCGGCGTTTACCTGAAGATAAATCCCttt > 1:2743267/1‑57 (MQ=255) ggTTCGCCCGAGAATCGCAATATGTTCCCCGGCGTTTACCTGAAGATAAATCCCttt > 1:2588093/1‑57 (MQ=255) ggTTCGCCCGAGAATCGCAATATGTTCCCCGGCGTTTACCTGAAGATAAATCCCttt > 1:2515739/1‑57 (MQ=255) ggTTCGCCCGAGAATCGCAATATGTTCCCCGGCGTTTACCTGAAGATAAATCCCttt > 1:1845131/1‑57 (MQ=255) ggTTCGCCCGAGAATCGCAATATGTTCCCCGGCGTTTACCTGAAGATAAATCCCttt > 1:2151806/1‑57 (MQ=255) ggTTCGCCCGAGAATCGCAATATGTTCCCCGGCGTTTACCTGAAGATAAATCCCttt > 1:2093628/1‑57 (MQ=255) gagaATCGCAATATGTTCCCCGGCGTTTACCTGAAGATAAATCCCTTTAAGTGCCTGTTGGGATTGCTc < 1:1299601/69‑1 (MQ=255) gagaATCGCAATATGTTCCCCGGCGTTTACCTGAAGATAAATCCCTTTAAGTGCCTGTTGGGATTGCTc < 1:2464840/69‑1 (MQ=255) gagaATCGCAATATGTTCCCCGGCGTTTACCTGAAGATAAATCCCTTTAAGTGCCTGTTGGGATTGCTc < 1:1539882/69‑1 (MQ=255) gagaATCGCAATATGTTCCCCGGCGTTTACCTGAAGATAAATCCCTTTAAGTGCCTGTTGGGATTGCTc < 1:2716144/69‑1 (MQ=255) gagaATCGCAATATGTTCCCCGGCGTTTACCTGAAGATAAATCCCTTTAAGTGCCTGTTGGGATTGCTc < 1:220931/69‑1 (MQ=255) gagaATCGCAATATGTTCCCCGGCGTTTACCTGAAGATAAATCCCTTTAAGTGCCTGTTGGGATTGCTc < 1:2036439/69‑1 (MQ=255) gagaATCGCAATATGTTCCCCGGCGTTTACCTGAAGATAAATCCCTTTAAGTGCCTGTTGGGATTGCTc < 1:482215/69‑1 (MQ=255) gagaATCGCAATATGTTCCCCGGCGTTTACCTGAAGATAAATCCCTTTAAGTGCCTGTTGGGATTGCTc < 1:324056/69‑1 (MQ=255) gcaatatTTTCCCCGGCGTTTACCTGAAGATAAATCCCTTTAAGTGCCTGTTGGGATTGCTc < 1:3115626/59‑1 (MQ=255) tGTTCCCCGGCGTTTACCTGAAGAGAAATCCCTTTAAGTGCCTGTTGAGATTGCTCCGGATAAGtgaa < 1:1208402/68‑1 (MQ=255) tGTTCCCCGGCGTTTACCTGAAGAGAAATCCCTTTAAGTGCCTGTTGAGATTGCTCCGGATAAGtgaa < 1:263834/68‑1 (MQ=255) ggCGTTTACCTGAAGATAAATCCCTTTAAGTGCCTGTTGGGATTGCTCCGGATAAGTGAACTGAACATc < 1:3163524/69‑1 (MQ=255) ggCGTTTACCTGAAGATAAATCCCTTTAAGTGCCTGTTGGGATTGCTCCGGATAAGTGAACTGAACATc < 1:2039717/69‑1 (MQ=255) ggCGTTTACCTGAAGATAAATCCCTTTAAGTGCCTGTTGGGATTGCTCCGGATAAGTGAACTGAACATc < 1:404235/69‑1 (MQ=255) ggCGTTTACCTGAAGATAAATCCCTTTAAGTGCCTGTTGGGATTGCTCCGGATAAGTGAACTGAACATc < 1:518216/69‑1 (MQ=255) ggCGTTTACCTGAAGATAAATCCCTTTAAGTGCCTGTTGGGATTGCTCCGGATAAGTGAACTGAACATc < 1:2097561/69‑1 (MQ=255) ggCGTTTACCTGAAGATAAATCCCTTTAAGTGCCTGTTGGGATTGCTCCGGATAAGTGAACTGAACATc < 1:227732/69‑1 (MQ=255) ggCGTTTACCTGAAGATAAATCCCTTTAAGTGCCTGTTGGGATTGCTCCGGATAAGTGAACTGAACATc < 1:724596/69‑1 (MQ=255) tACCTGAAGAGAAATCCCTTTAAGTGCCTGTTGAGATTGCTCCGGATAAGTGAACTGAACATCCCGTaa > 1:2550794/1‑69 (MQ=255) tACCTGAAGAGAAATCCCTTTAAGTGCCTGTTGAGATTGCTCCGGATAAGTGAACTGAACATCCCGTaa > 1:2550959/1‑69 (MQ=255) tACCTGAAGAGAAATCCCTTTAAGTGCCTGTTGAGATTGCTCCGGATAAGTGAACTGAACATCCCGTaa > 1:941322/1‑69 (MQ=255) tACCTGAAGAGAAATCCCTTTAAGTGCCTGTTGAGATTGCTCCGGATAAGTGAACTGAACATCCCGTaa > 1:895444/1‑69 (MQ=255) tACCTGAAGAGAAATCCCTTTAAGTGCCTGTTGAGATTGCTCCGGATAAGTGAACTGAACATCCCGTaa < 1:1737159/69‑1 (MQ=255) tACCTGAAGAGAAATCCCTTTAAGTGCCTGTTGAGATTGCTCCGGATAAGTGAACTGAACATCCCGTaa < 1:61750/69‑1 (MQ=255) tACCTGAAGAGAAATCCCTTTAAGTGCCTGTTGAGATTGCTCCGGATAAGTGAACTGAACATCCCGTaa > 1:3124747/1‑69 (MQ=255) tACCTGAAGAGAAATCCCTTTAAGTGCCTGTTGAGATTGCTCCGGATAAGTGAACTGAACATCCCGTaa > 1:3309389/1‑69 (MQ=255) tACCTGAAGAGAAATCCCTTTAAGTGCCTGTTGAGATTGCTCCGGATAAGTGAACTGAACATCCCGTaa < 1:3265829/69‑1 (MQ=255) agaAATCCCTTTAAGTGCCTGTTGAGATTGCTCCGGATAAGTGAACTGAACATCCCGTAACGTCAGCGa < 1:3238394/69‑1 (MQ=255) | AGTGTTGATTTGCCGCATCCGGTTCGCCCGAGAATCGCTATATGTTCCCCGGCGTTTACCTGAAGAGAAATCCCTTTAAGTGCCTGTTGAGATTGCTCCGGATAAGTGAACTGAACATCCCGTAACGTCAGCGA > minE/584750‑584883 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 21 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |