Predicted mutation | ||||||
---|---|---|---|---|---|---|
evidence | seq id | position | mutation | annotation | gene | description |
RA | REL606 | 2,775,879 | Δ1 bp | intergenic (‑49/+200) | insJ‑3 ← / ← cysH | IS150 hypothetical protein/phosphoadenosine phosphosulfate reductase |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | REL606 | 2,775,878 | 0 | G | . | 98.3% | 277.2 / ‑1.5 | 59 | intergenic (‑48/+201) | insJ‑3/cysH | IS150 hypothetical protein/phosphoadenosine phosphosulfate reductase |
Reads supporting (aligned to +/- strand): ref base G (0/0); major base . (30/28); minor base A (0/1); total (30/29) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 4.92e-01 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.50e-01 |
GGGCTATTCCATTTCATCGTCCAACAAAATGGGTGCAGTACAGGTTGATAGAAATATCGCCTGGGGCTTTTGTCCATCTGGAACCTCGC > REL606/2775836‑2775924 | gggCTATTCCATTTCATCGTCCAACAAAATGGGTGCAGTACA‑GTTGATAGa < 1:1760044/51‑1 (MQ=255) ggCTATTCCATTTCATCGTCCAACAAAATGGGTGCAGTACA‑GTTGATAGaa < 1:6226600/51‑1 (MQ=255) ggCTATTCCATTTCATCGTCCAACAAAATGGGTGCAGTACA‑GTTGATAGaa < 1:2934350/51‑1 (MQ=255) gCTATTCCATTTCATCGTCCAACAAAATGGGTGCAGTACA‑GTTGATAGaaa > 1:574784/1‑51 (MQ=255) gCTATTCCATTTCATCGTCCAACAAAATGGGTGCAGTACA‑GTTGATAGaaa > 1:2693330/1‑51 (MQ=255) gCTATTCCATTTCATCGTCCAACAAAATGGGTGCAGTACA‑GTTGATAGaaa < 1:250300/51‑1 (MQ=255) cTATTCCATTTCATCGTCCAACAAAATGGGTGCAGTACA‑GTTGATAGAAat < 1:2342345/51‑1 (MQ=255) cTATTCCATTTCATCGTCCAACAAAATGGGTGCAGTACA‑GTTGATAGAAat < 1:3826564/51‑1 (MQ=255) cTATTCCATTTCATCGTCCAACAAAATGGGTGCAGTACA‑GTTGATAGAAat < 1:838421/51‑1 (MQ=255) cTATTCCATTTCATCGTCCAACAAAATGGGTGCAGTACA‑GTTGATAGAAat > 1:2393717/1‑51 (MQ=255) cTATTCCATTTCATCGTCCAACAAAATGGGTGCAGTACA‑GTTGATAGAAat < 1:1232261/51‑1 (MQ=255) tATTCCATTTCATCGTCCAACAAAATGGGTGCAGTACA‑GTTGATAGAAata > 1:2700231/1‑51 (MQ=255) aTTCCATTTCATCGTCCAACAAAATGGGTGCAGTACA‑GTTGATAGAAatat < 1:5198260/51‑1 (MQ=255) ttCCATTTCATCGTCCAACAAAATGGGTGCAGTACA‑GTTGATAGAAATATc > 1:3428358/1‑51 (MQ=255) tCCATTTCATCGTCCAACAAAATGGGTGCAGTACA‑GTTGATAGAAATATCg > 1:5669934/1‑51 (MQ=255) tCCATTTCATCGTCCAACAAAATGGGTGCAGTACA‑GTTGATAGAAATATCg < 1:4638243/51‑1 (MQ=255) ccATTTCATCGTCCAACAAAATGGGTGCAGTACA‑GTTGATAGAAATATCGc > 1:6568564/1‑51 (MQ=255) aTTTCATCGTCCAACAAAATGGGTGCAGTACA‑GTTGATAGAAATATCGCCt > 1:586746/1‑51 (MQ=255) aTTTCATCGTCCAACAAAATGGGTGCAGTACA‑GTTGATAGAAATATCGCCt < 1:3294411/51‑1 (MQ=255) tttCATCGTCCAACAAAATGGGTGCAGTACA‑GTTGATAGAAATATCGCCTg < 1:1358552/51‑1 (MQ=255) ttCATCGTCCAACAAAATGGGTGCAGTACA‑GTTGATAGAAATATCGCCTgg > 1:2790543/1‑51 (MQ=255) ttCATCGTCCAACAAAATGGGTGCAGTACA‑GTTGATAGAAATATCGCCTgg > 1:1171428/1‑51 (MQ=255) ttCATCGTCCAACAAAATGGGTGCAGTACA‑GTTGATAGAAATATCGCCTgg < 1:2481093/51‑1 (MQ=255) tCATCGTCCAACAAAATGGGTGCAGTACA‑GTTGATAGAAATATCGCCTggg > 1:183625/1‑51 (MQ=255) cATCGTCCAACAAAATGGGTGCAGTACA‑GTTGATAGAAATATCGCCTgggg > 1:3762330/1‑51 (MQ=255) cGTCCAACAAAATGGGTGCAGTACA‑GTTGATAGAAATATCGCCTGGGGCtt > 1:1920997/1‑51 (MQ=255) cGTCCAAC‑AAATGGGTGCAGTACA‑GTTGATAGAAATATCGCCTGGGGCtt < 1:3124501/50‑1 (MQ=25) gTCCAACAAAATGGGTGCAGTACA‑GTTGATAGAAATATCGCCTGGGGCttt < 1:1974002/51‑1 (MQ=255) gTCCAACAAAATGGGTGCAGTACA‑GTTGATAGAAATATCGCCTGGGGCttt > 1:1586727/1‑51 (MQ=255) gTCCAACAAAATGGGTGCAGTACA‑GTTGATAGAAATATCGCCTGGGGCttt > 1:6204760/1‑51 (MQ=255) gTCCAACAAAATGGGTGCAGTACA‑GTTGATAGAAATATCGCCTGGGGCttt < 1:6017716/51‑1 (MQ=255) tCCAACAAAATGGGTGCAGTACA‑GTTGATAGAAATATCGCCTGGGGCtttt > 1:6553812/1‑51 (MQ=255) tCCAACAAAATGGGTGCAGTACA‑GTTGATAGAAATATCGCCTGGGGCtttt < 1:7280371/51‑1 (MQ=255) tCCAACAAAATGGGTGCAGTACA‑GTTGATAGAAATATCGCCTGGGGCtttt < 1:2852418/51‑1 (MQ=255) tCCAACAAAATGGGTGCAGTACA‑GTTGATAGAAATATCGCCTGGGGCtttt > 1:6002299/1‑51 (MQ=255) tCCAACAAAATGGGTGCAGTACA‑GTTGATAGAAATATCGCCTGGGGCtttt > 1:6384997/1‑51 (MQ=255) caaAATGGGTGCAGTACA‑GTTGATAGAAATATCGCCTGGGGCTTTTGTCCa < 1:5259812/51‑1 (MQ=255) caaAATGGGTGCAGTACA‑GTTGATAGAAATATCGCCTGGGGCTTTTGTCCa < 1:4624461/51‑1 (MQ=255) aaaaTGGGTGCAGTACA‑GTTGATAGAAATATCGCCTGGGGCTTTTGTCCAt > 1:5655888/1‑51 (MQ=255) aaaTGGGTGCAGTACA‑GTTGATAGAAATATCGCCTGGGGCTTTTGTCCATc < 1:2934459/51‑1 (MQ=255) aaTGGGTGCAGTACA‑GTTGATAGAAATATCGCCTGGGGCTTTTGt > 1:3755827/1‑45 (MQ=38) aTGGGTGCAGTACA‑GTTGATAGAAATATCGCCTGGGGCTTTTGTCCATCTg > 1:1008149/1‑51 (MQ=255) aTGGGTGCAGTACA‑GTTGATAGAAATATCGCCTGGGGCTTTTGTCCATCTg < 1:2382496/51‑1 (MQ=255) aTGGGTGCAGTACA‑GTTGATAGAAATATCGCCTGGGGCTTTTGTCCATCTg > 1:2196692/1‑51 (MQ=255) aTGGGTGCAGTACA‑GTTGATAGAAATATCGCCTGGGGCTTTTGTCCATCTg < 1:813699/51‑1 (MQ=255) ggTGCAGTA‑AAGTTGATAGAAATATCGCCTGGGGCTTTTGTCCATCTGGaa < 1:1396976/51‑1 (MQ=38) gTGCAGTACA‑GTTGATAGAAATATCGCCTGGGGCTTTTGTCCATCTGGAAc > 1:890665/1‑51 (MQ=255) gTGCAGTACA‑GTTGATAGAAATATCGCCTGGGGCTTTTGTCCATCTGGAAc > 1:5432522/1‑51 (MQ=255) tGCAGTACA‑GTTGATAGAAATATCGCCTGGGGCTTTTGTCCATCTGGAAcc > 1:4539840/1‑51 (MQ=255) tGCAGTACA‑GTTGATAGAAATATCGCCTGGGGCTTTTGTCCATCTGGAAcc < 1:5757739/51‑1 (MQ=255) tGCAGTACA‑GTTGATAGAAATATCGCCTGGGGCTTTTGTCCATCTGGAAcc > 1:1171085/1‑51 (MQ=255) tGCAGTACA‑GTTGATAGAAATATCGCCTGGGGCTTTTGTCCATCTGGAAcc < 1:6535078/51‑1 (MQ=255) gCAGTACA‑GTTGATAGAAATATCGCCTGGGGCTTTTGTCCATCTGGAACCt > 1:2878867/1‑51 (MQ=255) cAGTACA‑GTTGATAGAAATATCGCCTGGGGCTTTTGTCCATCTGGAACCTc > 1:324949/1‑51 (MQ=255) cAGTACA‑GTTGATAGAAATATCGCCTGGGGCTTTTGTCCATCTGGAACCTc < 1:806368/51‑1 (MQ=255) aGTACA‑GTTGATAGAAATATCGCCTGGGGCTTTTGTCCATCTGGAACCTcg > 1:4874273/1‑51 (MQ=255) aGTACA‑GTTGATAGAAATATCGCCTGGGGCTTTTGTCCATCTGGAACCTcg < 1:5839725/51‑1 (MQ=255) aGTACA‑GTTGATAGAAATATCGCCTGGGGCTTTTGTCCATCTGGAACCTcg < 1:5787157/51‑1 (MQ=255) gTACA‑GTTGATAGAAATATCGCCTGGGGCTTTTGTCCATCTGGAACCTcgc > 1:6077982/1‑51 (MQ=255) | GGGCTATTCCATTTCATCGTCCAACAAAATGGGTGCAGTACAGGTTGATAGAAATATCGCCTGGGGCTTTTGTCCATCTGGAACCTCGC > REL606/2775836‑2775924 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 31 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |