Predicted mutation | ||||||
---|---|---|---|---|---|---|
evidence | seq id | position | mutation | annotation | gene | description |
RA | REL606 | 475,640 | T→G | D375A (GAT→GCT) | fsr ← | predicted fosmidomycin efflux system |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | REL606 | 475,640 | 0 | T | G | 91.7% | 274.2 / 7.3 | 72 | D375A (GAT→GCT) | fsr | predicted fosmidomycin efflux system |
Reads supporting (aligned to +/- strand): ref base T (2/4); new base G (27/39); total (29/43) | |||||||||||
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 |
GCAGGAATATGGTCAACATCCCCAATAGTGGCAGGAAAGCACAGATTTTATAGACTAACTCGATGCTGGTGTGATCGGCGATAAGCCCCAGAACTGCCGCTCCCAGACCTCCCATGCCAAAAGCAAAACCGAAAAAGAGTCCAGAAA > REL606/475566‑475712 | gcaggAATATGGTCAACATCCCCAATAGTGGCAGGAAAGCACAGATTTTATAGACTAACTCGATGCTGGTGTGATc < 2:1534961/76‑1 (MQ=255) caggAATATGGTCAACATCCCCAATAGTGGCAGGAAAGCACAGATTTTATAGACTAACTCGATGCTGGTGTGAGCg < 1:731502/76‑1 (MQ=255) caggAATATGGTCAACATCCCCAATAGTGGCACGAAAGCACAGATTTTATAGACTAACTCGATGCTGGTGTGAGCg > 1:235604/1‑76 (MQ=255) aggAATATGGTCAACATCCCCAATAGTGGCAGGAAAGCACAGATTTTATAGACTAACTCGATGCTGGTGTGAGCgg < 2:184206/76‑1 (MQ=255) aggAATATGGTCAACATCCCCAATAGTGGCAGGAAAGCACAGATTTTATAGACTAACTCGATGCTGGTGTGAGCgg < 2:1972699/76‑1 (MQ=255) atGGTCAACATCCCCAATAGTGGCAGGAAAGCACAGATTTTATAGACTAACTCGATGCTGGTGTGAGCGGCGATaa < 1:266815/76‑1 (MQ=255) tGGTCAACATCCCCAATAGTGGCAGGAAAGCACAGATTTTATAGACTAACTCGATGCTGGTGTGAGCGGCGATAAg < 2:251859/76‑1 (MQ=255) tGGTCAACATCCCCAATAGTGGCAGGAAAGCACAGATTTTATAGACTAACTCGATGCTGGTGTGAGCGGCGATAAg < 1:1706017/76‑1 (MQ=255) tGGTCAACATCCCCAATAGTGGCAGGAAAGCACAGATTTTATAGACTAACTCGATGCTGGTGTGAGCGGCGATAAg < 2:366499/76‑1 (MQ=255) ggTCAACATCCCCAATAGTGGCAGGAAAGCACAGATTTTATAGACTAACTCGATGCTGGTGTGAGCGGCGATaaac > 1:493675/1‑74 (MQ=255) ggTCAACATCCCCAATAGTGGCAGGAAAGCACAGATTTTATAGACTAACTCGATGCTGGGGTGGGCGGGGATAAGc > 1:1457082/1‑76 (MQ=255) cAACATCCCCAATAGTGGCAGGAAAGCACAGATTTTATAGACTAACTCGATGCTGGTGTGAGCGGCGATAAGcccc < 2:542165/76‑1 (MQ=255) aaCATCCCCAATAGTGGCAGGAAAGCACAGATTTTATAGACTAACTCGATGCTGGTGTGAGCGGCGATAAGCCCCa < 2:411308/76‑1 (MQ=255) aaCATCCCCAATAGTGGCAGGAAAGCACAGATTTTATAGACTAACTCGATGCTGGTGTGAGCGGCGATAAGCCCCa < 2:933782/76‑1 (MQ=255) aCATCCCCAATAGTGGCAGGAAAGCACAGATTTTATAGACTAACTCGATGCTGGTGTGAGCGGCGATAAGCCCCAg < 1:253482/76‑1 (MQ=255) aCATCCCCAATAGTGGCAGGAAAGCACAGATTTTATAGACTAACTCGATGCTGGTGTGAGCGGCGATAAGCCCCAg < 1:1637120/76‑1 (MQ=255) tCCCCAATAGTGGCAGGAAAGCACAGATTTTATAGACTAACTCGATGCTGGTGTGAGCGGCGATAAGCCCCAGAAc < 2:998962/76‑1 (MQ=255) tCCCCAATAGTGGCAGGAAAGCACAGATTTTATAGACTAACTCGATGCTGGTGTGAGCGGCGATAAGCCCCAGAAc < 2:998952/76‑1 (MQ=255) cccAATAGTGGCAGGAAAGCACAGATTTTATAGACTAACTCGATGCTGGTGTGAGCGGCGATAAGCCCCAGAACTg < 1:1569630/76‑1 (MQ=255) cAATAGTGGCAGGAAAGCACAGATTTTATAGACTAACTCGATGCTGGTGTGAGCGGCGATAAGCCCCAGAACTGcc < 2:79763/76‑1 (MQ=255) cAATAGTGGCAGGAAAGCACAGATTTTATAGACTAACTCGATGCTGGGGTGATCGGCGATAAGCCCCAGAACTGcc > 1:92994/1‑76 (MQ=255) aaTAGTGGCAGGAAAGCACAGATTTTATAGACTAACTCGATGCTGGTGTGAGCGGCGATAAGCCCCAGAACTGCCg < 1:249984/76‑1 (MQ=255) aaTAGTGGCAGGAAAGCACAGATTTTATAGACTAACTCGATGCTGGTGTGAGCGGCGATAAGCCCCAGAACTGCCg < 2:1742829/76‑1 (MQ=255) tGGCAGGAAAGCACAGATTTTATAGACTAACTCGATGCTGGTGTGAGCGGCGATAAGCCCCAGAACTGCCGCTccc < 2:346306/76‑1 (MQ=255) ggCAGGAAAGCACAGATTTTATAGACTAACTCGATGCTGGTGTGAGCGGCGATAAGCCCCAGAACTGCCGCTCCCa < 1:22952/76‑1 (MQ=255) cAGGAAAGCACAGATTTTATAGACTAACTCGATGCTGGTGTGAGCGGCGATAAGCCCCAGAACTGCCGCTCCCAGa < 1:1267769/76‑1 (MQ=255) cAGGAAAGCACAGATTTTATAGACTAACTCGATGCTGGTGTGAGCGGCGATAAGCCCCAGAACTGCCGCTCCCAGa > 2:916707/1‑76 (MQ=255) cAGGAAAGCACAGATTTTATAGACTAACTCGATGCTGGTGTGAGCGGCGATAAGCCCCAGAACTGCCGCTCCCAGa > 1:70833/1‑76 (MQ=255) aGGAAAGCACAGATTTTATAGACTAACTCGATGCTGGTGTGAGCGGCGATAAGCCCCAGAACTGCCGCTCCCAGAc < 2:19838/76‑1 (MQ=255) aGGAAAGCACAGATTTTATAGACTAACTCGATGCTGGTGTGAGCGGCGATAAGCCCCAGAACTGCCGCTCCCAGAc > 2:762119/1‑76 (MQ=255) aGGAAAGCACAGATTTTATAGACTAACTCGATGCTGGTGTGAGCGGCGATAAGCCCCAGAACTGCCGCTCCCAGAc < 2:1202791/76‑1 (MQ=255) aGGAAAGCACAGATTTTATAGACTAACTCGATGCTGGTGTGAGCGGCGATAAGCCCCAGAACTGCCGCTCCCAGAc > 1:86335/1‑76 (MQ=255) gAAAGCACAGATTTTATAGACTAACTCGATGCTGGTGTGAGCGGCGATAAGCCCCAGAACTGCCGCTCCCAGACCt < 2:104136/76‑1 (MQ=255) aaGCACAGATTTTATAGACTAACTCGATGCTGGTGTGAGCGGCGATAAGCCCCAGAACTGCCGCTCCCAGACCTcc < 1:1546412/76‑1 (MQ=255) cacaGATTTTATAGACTAACTCGATGCTGGTGTGAGCGGCGATAAGCCCCAGAACTGCCGCTCCCAGACCTCCCAt < 2:1984064/76‑1 (MQ=255) cacaGATTTTATAGACTAACTCGATGCTGGTGTGAGCGGCGATAAGCCCCAGAACTGCCGCTCCCAGACCTCCCAt < 2:1958453/76‑1 (MQ=255) cacaGATTTTATAGACTAACTCGATGCTGGCGAGAGCGGCGATAAGCCCCAGAACTGCCGCTCCCAGACCTCCCAt > 2:1531285/1‑76 (MQ=255) caGATTTTATAGACTAACTCGATGCTGGTGTGAGCGGCGATAAGCCCCAGAACTGCCGCTCCCAGACCTCCCATGc > 2:1136332/1‑76 (MQ=255) aTTTTATAGACTAACTCGATGCTGGTGTGAGCGGCGATAAGCCCCAGAACTGCCGCTCCCAGACCTCCCATGCCaa > 1:275351/1‑76 (MQ=255) aTTTTATAGACTAACTCGATGCTGGTGTGAGCGGCGATAAGCCCCAGAACTGCCGCTCCCAGACCTCCCATGCCaa < 1:446534/76‑1 (MQ=255) ttttATAGACTAACTCGATGCTGGTGTGAGCGGCGATAAGCCCCAGAACTGCCGCTCCCAGACCTCCCATGCCaaa > 1:1085304/1‑76 (MQ=255) ttttATAGACTAACTCGATGCTGGTGTGAGCGGCGATAAGCCCCAGAACTGCCGCTCCCAGACCTCCCATGCCaaa > 1:1726238/1‑76 (MQ=255) tttATAGACTAACTCGATGCTGGTGTGAGCGGCGATAAGCCCCAGAACTGCCGCTCCCAGACCTCCCATGCCaaaa > 2:492256/1‑76 (MQ=255) ttATAGACTAACTCGATGCTGGTGTGATCGGCGATAAGCCCCAGAACTGCCGCTCCCAGACCTCCCATGCCAAAAg > 2:1481221/1‑76 (MQ=255) ttATAGACTAACTCGATGCTGGTGTGAGCGGCGATAAGCCCCAGAACTGCCGCTCCCAGACCTCCCATGCCAAAAg > 2:1723710/1‑76 (MQ=255) tataGACTAACTCGATGCTGGTGTGATCGGCGATAAGCCCCAGAACTGCCGCTCCCAGACCTCCCATGCCAAAAGc < 1:336616/76‑1 (MQ=255) taGACTAACTCGATGCTGGTGTGAGCGGCGATAAGCCCCAGAACTGCCGCTCCCAGACCTCCCATGCCAAAAGCaa < 1:583016/76‑1 (MQ=255) aGACTAACTCGATGCTGGTGTGATCGGCGATAAGCCCCAGAACTGCCGCTCCCAGACCTCCCATGCCAAAAGCaaa < 1:263892/76‑1 (MQ=255) gACTAACTCGATGCTGGTGTGATCGGCGATAAGCCCCAGAACTGCCGCTCCCAGACCTCCCATGCCAAAAGCaaaa < 2:1175219/76‑1 (MQ=255) gACTAACTCGATGCTGGTGTGAGCGGCGATAAGCCCCAGAACTGCCGCTCCCAGACCTCCCATGCCAAAAGCaaaa < 1:1213692/76‑1 (MQ=255) cTAACTCGATGCTGGTGTGAGCGGCGATAAGCCCCAGAACTGCCGCTCCCAGACCTCCCATGCCAAAAGCAAAAcc > 1:1746649/1‑76 (MQ=255) tAACTCGATGCTGGTGTGAGCGGCGATAAGCCCCAGAACTGCCGCTCCCAGACCTCCCATGCCAAAAGCAAAACCg > 2:908458/1‑76 (MQ=255) aCTCGATGCTGGTGTGAGCGGCGATAAGCCCCAGAACTGCCGCTCCCAGACCTCCCATGCCAAAAGCAAAACCGaa > 2:776703/1‑76 (MQ=255) aCTCGATGCTGGTGTGAGCGGCGATAAGCCCCAGAACTGCCGCTCCCAGACCTCCCATGCCAAAAGCAAAACCGaa < 1:1044678/76‑1 (MQ=255) cTCGATGCTGGTGTGAGCGGCGATAAGCCCCAGAACTGCCGCTCCCAGACCTCCCATGCCAAAATCAAAACCGaaa > 1:1137391/1‑76 (MQ=255) tCGATGCTGGTGTGAGCGGCGATAAGCCCCAGAACTGCCGCTCCCAGACCTCCCATGCCAAAAGCAAAACCGaaaa > 1:356234/1‑76 (MQ=255) cGATGCTGGTGTGAGCGGCGATAAGCCCCAGAACTGCCGCTCCCAGACCTCCCATGCCAAAAGCAAAACCGaaaaa < 1:825377/76‑1 (MQ=255) cGATGCTGGTGTGAGCGGCGATAAGCCCCAGAACTGCCGCTCCCAGACCTCCCATGCCAAAAGCAAAACCGaaaaa > 1:1245658/1‑76 (MQ=255) cGATGCTGGTGTGAGCGGCGATAAGCCCCAGAACTGCCGCTCCCAGACCCCCCATGCCAAAGGCAAAACCGaaaaa > 2:41129/1‑76 (MQ=255) aTGCTGGTGTGAGCGGCGATAAGCCCCAGAACTGCCGCTCCCAGACCTCCCATGCCAAAAGCAAAACCGAAAAaga > 1:1527966/1‑76 (MQ=255) tGCTGGTGTGAGCGGCGATAAGCCCCAGAACTGCCGCTCCCAGACCTCCCATGCCAAAAGCAAAACCGAAAAagag > 1:1512448/1‑76 (MQ=255) gCTGGTGTGAGCGGCGATAAGCCCCAGAACTGCCGCTCCCAGACCTCCCATGCCAAAAGCAAAACCGAAAAAGAGt > 2:1084847/1‑76 (MQ=255) cTGGTGTGAGCGGCGATACGCCCCAGAACTGCCGCTCCCAGACCTCCCATGCCAAAAGCCAAACCGACAAAGAGTc > 2:749035/1‑76 (MQ=255) cTGGTGTGAGCGGCGATAAGCCCCAGAACTGCCGCTCCCAGACCTCCCATGCCAAAAGCAAAACCGAAAAAGAGTc > 2:1691617/1‑76 (MQ=255) tGGTGTGAGCGGCGATAAGCCCCAGAACTGCCGCTCCCAGACCTCCCATGCCAAAAGCAAAACCGAAAAAGAGTcc > 2:1008163/1‑76 (MQ=255) ggTGTGAGCGGCGATAAGCCCCGGAACTGCCGCTCCCAGACCTCCCATGCCAAAAGCAAAACCGAAAAAGAGTCCa < 2:1659607/76‑1 (MQ=255) ggTGTGAGCGGCGATAAGCCCCAGAACTGCCGCTCCCAGACCTCCCATGCCAAAAGCAAAACCGAAAAAGAGTCCa < 2:1812839/76‑1 (MQ=255) ggTGTGAGCGGCGATAAGCCCCAGAACTGCCGCTCCCAGACCTCCCATGCCAAAAGCAAAACCGAAAAAGAGTCCa > 1:1220273/1‑76 (MQ=255) gtgtGAGCGGCGATAAGCCCCAGAACTGCCGCTCCCAGACCTCCCATGCCAAAAGCAAAACCGAAAAAGAGTCCAg < 1:1664854/76‑1 (MQ=255) gtgtGAGCGGCGATAAGCCCCAGAACTGCCGCTCCCAGACCTCCCATGCCAAAAGCAAAACCGAAAAAGAGTCCAg < 2:795813/76‑1 (MQ=255) gtgtGAGCGGCGATAAGCCCCAGAACTGCCGCTCCCAGACCTCCCATGCCAAAAGCAAAACCGAAAAAGAGTCCAg < 2:1121809/76‑1 (MQ=255) gtgAGCGGCGATAAGCCCCAGAACTGCCGCTCCCAGACCTCCCATGCCAAAAGCAAAACCGAAAAAGAGTCCAGaa > 1:816348/1‑76 (MQ=255) gtgAGCGGCGATAAGCCCCAGAACTGCCGCTCCCAGACCTCCCATGCCAAAAGCAAAACCGAAAAAGAGTCCAGaa < 2:131597/76‑1 (MQ=255) tgAGCGGCGATAAGCCCCAGAACTGCCGCTCCCAGACCTCCCATGCCAAAAGCAAAACCGAAAAAGAGTCCAGaaa < 1:113303/76‑1 (MQ=255) | GCAGGAATATGGTCAACATCCCCAATAGTGGCAGGAAAGCACAGATTTTATAGACTAACTCGATGCTGGTGTGATCGGCGATAAGCCCCAGAACTGCCGCTCCCAGACCTCCCATGCCAAAAGCAAAACCGAAAAAGAGTCCAGAAA > REL606/475566‑475712 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 21 ≤ ATCG/ATCG < 29 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |