Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | W3110S.gb | 3,348,085 | G→A | 25.6% | E89K (GAA→AAA) | npr → | phosphohistidinoprotein‑hexose phosphotransferase component of N‑regulated PTS system (Npr) |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | W3110S.gb | 3,348,085 | 0 | G | A | 25.6% | 78.0 / 24.5 | 43 | E89K (GAA→AAA) | npr | phosphohistidinoprotein‑hexose phosphotransferase component of N‑regulated PTS system (Npr) |
Reads supporting (aligned to +/- strand): ref base G (19/13); new base A (6/5); total (25/18) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 1.00e+00 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 |
GCGACCGGTCCACAGGAAGAGGAAGCACTGGCCGCCGTTATCGCCCTCTTTAATTCTGGTTTTGATGAAGATTAATCTTCATCACTTTGACATACAACTCCCTTCAAACTCCCCCCGACAATAAGAAAATCA > W3110S.gb/3348019‑3348150 | gCGACCGGTCCACAGGAAGAGGAAGCACTGGCCGCCGTTATCGCCCTCTTTAATTCTGGTTTtgatg > 1:1186854/1‑67 (MQ=255) aCCGGTCCACAGGAAGAGGAAGCACTGGCCGCCGTTATCGCCCTCTTTAATTCTGGTTTTGATGAAGAtt > 1:1798478/1‑70 (MQ=255) aCCGGTCCACAGGAAGAGGAAGCACTGGCCGCCGTTATCGCCCTCTTTAATTCTGGTTTTGATGAAGATTa < 1:1853257/71‑1 (MQ=255) ccGGGCCACAGGAAGAGGAAGCACTGGCCGCCGTTATCGCCCTCTTTAATTCTGGTTTTGATGAAGATTaa < 1:436406/71‑1 (MQ=255) tCCACAGGAAGAGGAAGCACTGGCCGCCGTTATCGCCCTCTTTAATTCTGGTTTTGATGAAGATTAATCtt < 1:2234275/71‑1 (MQ=255) acaGGAAGAGGAAGCACTGGCCGCCGTTATCGCCCTCTTTAATTCTGGTTTTGATGAAGATTAATCTtcat < 1:2761567/71‑1 (MQ=255) agaggaagCACTGGCCGCCGTTATCGCCCTCTTTAATTCTGGTTTTGATGaa < 1:1678953/52‑1 (MQ=255) agCACTGGCCGCCGTTATCGCCCTCTTTAATTCTGGTTTTGATGAAGATTAATCTTCATCACTTTGACATa < 1:1822254/71‑1 (MQ=255) gCACTGGCCGCCGTTATCGCCCTCTTTAATTCTGGTTTTGATGAAGAtt < 1:1783989/49‑1 (MQ=255) tGGCCGCCGTTATCGCCCTCTTTAATTCTGGTTTTGATGAAGATTAATCTTCATCACTTTGACATAc < 1:2175848/67‑1 (MQ=255) tGGCCGCCGTTATCGCCCTCTTTAATTCTGGTTTTGATAAAGATTa < 1:1386574/46‑1 (MQ=255) ggCCGCCGTTATCGCCCTCTTTAATTCTGGTTTTGATGAAGATTAATCTtc > 1:455986/1‑51 (MQ=255) cgccgTTATCGCCCTCTTTAATTCTGGTTTTGATGAAGATTAATCTTCATCACTTTGAc > 1:606357/1‑59 (MQ=255) gTTATCGCCCTCTTTAATTCTGGTTTTGATAAAGATTAATCTTCATCACTTTGACATACAACTCCCTTCa > 1:2279338/1‑70 (MQ=255) ttATCGCCCTCTTTAATTCTGGTTTTGATAAAGATTa < 1:457721/37‑1 (MQ=255) ttATCGCCCTCTTTAATTCTGGTTTTGATAAAGATTAATCTTCATCACTTTGACATAc < 1:66175/58‑1 (MQ=255) ttATCGCCCTCTTTAATTCTGGTTTTGATAAAGATTAATCTTCATCACTTTGACATAc < 1:2237751/58‑1 (MQ=255) ttATCGCACTCTTTAATTCTGGTTTTGATAAAGATTa < 1:1292200/37‑1 (MQ=37) aTCGCCCTCTTTAATTCTGGTTTTGATGAAGATTAATCTTCATCACTTTGACATACAACTCCCTTCAAACt > 1:274040/1‑71 (MQ=255) aTCGCCCTCTTTAATTCTGGTTTTGATAAAGATTAATCTTCATCACTTTGACATACaa > 1:2662207/1‑58 (MQ=255) tctTTAATTCTGGTTTTGATGAAGATTAATCTTCATCACTTTGACATACAACTCCCTTCAAACTc < 1:1622907/65‑1 (MQ=255) ctTTAATTCTGGTTTTGATAAAGATTAATCTTCATCACTTTGACATa > 1:2243322/1‑47 (MQ=255) tttAATTCTGGTTTTGATGAAGATTAATCTTCATCAc > 1:272126/1‑37 (MQ=255) tttAATTCTGGTTTTGATGAAGATTAATCTTCATCAc > 1:2129412/1‑37 (MQ=255) tttAATTCTGGTTTTGATGAAGATTAATCTTCATCACttt < 1:1232403/40‑1 (MQ=255) tttAATTCTGGTTTTGATGAAGATTAATCTTCATCACTTTGACATAc < 1:221756/47‑1 (MQ=255) tttAATTCTGGTTTTGATGAAGATTAATCTTCATCACTTTGACATACAACTCCCTTCAAACTccc > 1:1644324/1‑65 (MQ=255) tttAATTCTGGTTTTGATGAAGATTAATCTTCATCACTTTGACATACAACTCCCTTCAAACTccc > 1:2361700/1‑65 (MQ=255) tttAATTCTGGTTTTGATGAAGATTAATCTTCATCACTTTGACATACAACTCCCTTCAAACTCCCCCCGa > 1:1227194/1‑70 (MQ=255) tttAATTCTGGTTTTGATGAAGATTAATCTTCATCACTTTGACATACAACTCCCTTCAAACTCCCCCCGa > 1:2727912/1‑70 (MQ=255) tAATTCTGGTTTTGATGAAGATTAATCTTCATCACtt > 1:1484642/1‑37 (MQ=255) tAATTCTGGTTTTGATGAAGATTAATCTTCATCACtt > 1:2675082/1‑37 (MQ=255) tAATTCTGGTTTTGATGAAGATTAATCTTCATCACTTTGACATACAACTCCCTTCAAACTCCCCCCGACaa > 1:217902/1‑71 (MQ=255) tAATTCTGGTTTTGATAAAGATTAATCTTCATCACTTTGa > 1:895923/1‑40 (MQ=255) tAATTCTGGTTTTGATAAAGATTAATCTTCATCACTTTGa > 1:1273801/1‑40 (MQ=255) tAATTCTGGTTTTGATAAAGATTAATCTTCATCACTTTGACATACAACTCCCTTCAAACTCCCCCCGACaa > 1:1019832/1‑71 (MQ=255) cTGGTTTTGATGAAGATTAATCTTCATCACTTTGACATACAACTccc > 1:1560461/1‑47 (MQ=255) cTGGTTTTGATGAAGATTAATCTTCATCACTTTGACATACAACTccc > 1:134946/1‑47 (MQ=255) cTGGTTTTGATGAAGATTAATCTTCATCACTTTGACATACAACTccc > 1:1032393/1‑47 (MQ=255) cTGGTTTTGATGAAGATTAATCTTCATCACTTTGACATACAACTCCCTTCAAACTCCCCCCGACAATAAg > 1:1192376/1‑70 (MQ=255) cTGGTTTTGATGAAGATTAATCTTCATCACTTTGACATACAACTCCCTTCAAACTCCCCCCGACAATAAg < 1:113842/70‑1 (MQ=255) tGGTTTTGATGAAGATTAATCTTCATCACTTTGACa > 1:1913415/1‑36 (MQ=255) ttttGATGAAGATTAATCTTCATCACTTTGACATACAACTCCCTTCAAACTCCCCCCGACAATAAGaa < 1:490878/68‑1 (MQ=255) tttGATGAAGATTAATCTTCATCACTTTGACATACAACTCCCTTCAAAc > 1:1488686/1‑49 (MQ=255) tgatgaAGATTAATCTTCATCACTTTGACATACAACTCCCTTCAAACTCCCCCCGACAATAAGAAAATCa < 1:2188290/70‑1 (MQ=255) | GCGACCGGTCCACAGGAAGAGGAAGCACTGGCCGCCGTTATCGCCCTCTTTAATTCTGGTTTTGATGAAGATTAATCTTCATCACTTTGACATACAACTCCCTTCAAACTCCCCCCGACAATAAGAAAATCA > W3110S.gb/3348019‑3348150 |
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: ‑ |