Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | minE | 585,755 | 0 | G | T | 36.0% | 78.0 / 44.1 | 50 | A47E (GCG→GAG) | cydC | fused cysteine transporter subunits and membrane component and ATP‑binding component of ABC superfamily |
Reads supporting (aligned to +/- strand): ref base G (0/32); new base T (18/0); total (18/32) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 5.54e-14 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.98e-01 | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
AGTCCGGCAACCCCCGCAACCGCTGAGGCCGAGAGGAACCAGCCGGAAAGTGTCAACAGACCGATACTGGCGAGCAGCGTCA > minE/585734‑585815 | aGTCCGGCAACCCCCGCAACCTCTGAGGCCGAGAGGCACCAGCCGGAAAGTGTCAACAGAcc > 1:207829/1‑62 (MQ=255) aGTCCGGCAACCCCCGCAACCTCTGAGGCCGAGAGGAACCAGCCGGAAAGTGTCAACAGa > 1:108139/1‑60 (MQ=255) aGTCCGGCAACCCCCGCAACCTCTGAGGCCGAGAGGAACCAGCCGGAAAGTGTCAACAGAcc > 1:795036/1‑62 (MQ=255) aGTCCGGCAACCCCCGCAACCTCTGAGGCCGAGAGGAACCAGCCGGAAAGTGTCAACAGAcc > 1:553211/1‑62 (MQ=255) aGTCCGGCAACCCCCGCAACCTCTGAGGCCGAGAGGAACCAGCCGGAAAGTGTCAACAGAcc > 1:41995/1‑62 (MQ=255) aGTCCGGCAACCCCCGCAACCTCTGAGGCCGAGAGGAACCAGCCGGAAAGTGTCAACAGAcc > 1:660088/1‑62 (MQ=255) aGTCCGGCAACCCCCGCAACCTCTGAGGCCGAGAGGAACCAGCCGGAAAGTGTCAACAGAcc > 1:352129/1‑62 (MQ=255) aGTCCGGCAACCCCCGCAACCTCTGAGGCCGAGAGGAACCAGCCGGAAAGTGTCAACAGAcc > 1:350673/1‑62 (MQ=255) aGTCCGGCAACCCCCGCAACCTCTGAGGCCGAGAGGAACCAGCCGGAAAGTGTCAACAGAcc > 1:257895/1‑62 (MQ=255) aGTCCGGCAACCCCCGCAACCTCTGAGGCCGAGAGGAACCAGCCGGAAAGTGTCAACAGAcc > 1:217210/1‑62 (MQ=255) aGTCCGGCAACCCCCGCAACCTCTGAGGCCGAGAGGAACCAGCCGGAAAGTGTCAACAGAcc > 1:174722/1‑62 (MQ=255) aGTCCGGCAACCCCCGCAACCTCTGAGGCCGAGAGGAACCAGCCGGAAAGTGTCAACAGAcc > 1:793408/1‑62 (MQ=255) aGTCCGGCAACCCCCGCAACCTCTGAGGCCGAGAGGAACCAGCCGGAAAGTGTCAACAGAcc > 1:847616/1‑62 (MQ=255) aGTCCGGCAACCCCCGCAACCTCTGAGGCCGAGAGGAACCAGCCGGAAAGTGTCAACAGAcc > 1:101520/1‑62 (MQ=255) aGTCCGGCAACCCCCGCAACCTCTGAGGCCGAGAGGAACCAGCCGGAAAGTGTCAACAGAcc > 1:1025716/1‑62 (MQ=255) aGTCCGGCAACCCCCGCAACCTCTGAGGCCGAGAGGAACCAGCCGGAAAGTGTCAACAGAcc > 1:1048941/1‑62 (MQ=255) aGTCCGGCAACCCCCGCAACCTCTGAGGCCGAGAGGAACCAGCCGGAAAGTGTCAACAGAc > 1:857509/1‑61 (MQ=255) aGTCCGGCAACCCCCGCAACCTCTGAGGCCGAGAGGAACCAGCCGGAAAGTGTCAACAGAc > 1:1045184/1‑61 (MQ=255) gCAACCGCTGAGGCCGAGAGGAACCAGCCGGAAAGTGTCAACAGACCGATACTGGCGagcag < 1:846875/62‑1 (MQ=255) cGCTGAGGCCGAGAGGAACCAGCCGGAAAGTGTCAACAGACCGATACTGGCGAGCAGCGTca < 1:1000482/62‑1 (MQ=255) cGCTGAGGCCGAGAGGAACCAGCCGGAAAGTGTCAACAGACCGATACTGGCGAGCAGCGTca < 1:635455/62‑1 (MQ=255) cGCTGAGGCCGAGAGGAACCAGCCGGAAAGTGTCAACAGACCGATACTGGCGAGCAGCGTca < 1:637368/62‑1 (MQ=255) cGCTGAGGCCGAGAGGAACCAGCCGGAAAGTGTCAACAGACCGATACTGGCGAGCAGCGTca < 1:683503/62‑1 (MQ=255) cGCTGAGGCCGAGAGGAACCAGCCGGAAAGTGTCAACAGACCGATACTGGCGAGCAGCGTca < 1:955450/62‑1 (MQ=255) cGCTGAGGCCGAGAGGAACCAGCCGGAAAGTGTCAACAGACCGATACTGGCGAGCAGCGTca < 1:724962/62‑1 (MQ=255) cGCTGAGGCCGAGAGGAACCAGCCGGAAAGTGTCAACAGACCGATACTGGCGAGCAGCGTca < 1:776321/62‑1 (MQ=255) cGCTGAGGCCGAGAGGAACCAGCCGGAAAGTGTCAACAGACCGATACTGGCGAGCAGCGTca < 1:948663/62‑1 (MQ=255) cGCTGAGGCCGAGAGGAACCAGCCGGAAAGTGTCAACAGACCGATACTGGCGAGCAGCGTca < 1:933612/62‑1 (MQ=255) cGCTGAGGCCGAGAGGAACCAGCCGGAAAGTGTCAACAGACCGATACTGGCGAGCAGCGTca < 1:893960/62‑1 (MQ=255) cGCTGAGGCCGAGAGGAACCAGCCGGAAAGTGTCAACAGACCGATACTGGCGAGCAGCGTca < 1:963429/62‑1 (MQ=255) cGCTGAGGCCGAGAGGAACCAGCCGGAAAGTGTCAACAGACCGATACTGGCGAGCAGCGTca < 1:795075/62‑1 (MQ=255) cGCTGAGGCCGAGAGGAACCAGCCGGAAAGTGTCAACAGACCGATACTGGCGAGCAGCGTca < 1:81447/62‑1 (MQ=255) cGCTGAGGCCGAGAGGAACCAGCCGGAAAGTGTCAACAGACCGATACTGGCGAGCAGCGTca < 1:158263/62‑1 (MQ=255) cGCTGAGGCCGAGAGGAACCAGCCGGAAAGTGTCAACAGACCGATACTGGCGAGCAGCGTca < 1:1013805/62‑1 (MQ=255) cGCTGAGGCCGAGAGGAACCAGCCGGAAAGTGTCAACAGACCGATACTGGCGAGCAGCGTca < 1:1034000/62‑1 (MQ=255) cGCTGAGGCCGAGAGGAACCAGCCGGAAAGTGTCAACAGACCGATACTGGCGAGCAGCGTca < 1:1089719/62‑1 (MQ=255) cGCTGAGGCCGAGAGGAACCAGCCGGAAAGTGTCAACAGACCGATACTGGCGAGCAGCGTca < 1:1126619/62‑1 (MQ=255) cGCTGAGGCCGAGAGGAACCAGCCGGAAAGTGTCAACAGACCGATACTGGCGAGCAGCGTca < 1:1140101/62‑1 (MQ=255) cGCTGAGGCCGAGAGGAACCAGCCGGAAAGTGTCAACAGACCGATACTGGCGAGCAGCGTca < 1:1158165/62‑1 (MQ=255) cGCTGAGGCCGAGAGGAACCAGCCGGAAAGTGTCAACAGACCGATACTGGCGAGCAGCGTca < 1:1178452/62‑1 (MQ=255) cGCTGAGGCCGAGAGGAACCAGCCGGAAAGTGTCAACAGACCGATACTGGCGAGCAGCGTca < 1:129512/62‑1 (MQ=255) cGCTGAGGCCGAGAGGAACCAGCCGGAAAGTGTCAACAGACCGATACTGGCGAGCAGCGTca < 1:546705/62‑1 (MQ=255) cGCTGAGGCCGAGAGGAACCAGCCGGAAAGTGTCAACAGACCGATACTGGCGAGCAGCGTca < 1:159740/62‑1 (MQ=255) cGCTGAGGCCGAGAGGAACCAGCCGGAAAGTGTCAACAGACCGATACTGGCGAGCAGCGTca < 1:173974/62‑1 (MQ=255) cGCTGAGGCCGAGAGGAACCAGCCGGAAAGTGTCAACAGACCGATACTGGCGAGCAGCGTca < 1:263945/62‑1 (MQ=255) cGCTGAGGCCGAGAGGAACCAGCCGGAAAGTGTCAACAGACCGATACTGGCGAGCAGCGTca < 1:348550/62‑1 (MQ=255) cGCTGAGGCCGAGAGGAACCAGCCGGAAAGTGTCAACAGACCGATACTGGCGAGCAGCGTca < 1:394690/62‑1 (MQ=255) cGCTGAGGCCGAGAGGAACCAGCCGGAAAGTGTCAACAGACCGATACTGGCGAGCAGCGTca < 1:465601/62‑1 (MQ=255) cGCTGAGGCCGAGAGGAACCAGCCGGAAAGTGTCAACAGACAGATACTGGCGAGCAGCGTca < 1:666590/62‑1 (MQ=255) cGCTGAGGCCGAAAGGAACCAGCCGGAAAGTGTCAACAGACCGATACTGGCGAGCAGCGTca < 1:3680/62‑1 (MQ=255) gCTGAGGCCGAGAGGAACCAGCCGGAAAGTGTCAACAGACCGATACTGGCGAGCAGCGTca < 1:13636/61‑1 (MQ=255) | AGTCCGGCAACCCCCGCAACCGCTGAGGCCGAGAGGAACCAGCCGGAAAGTGTCAACAGACCGATACTGGCGAGCAGCGTCA > minE/585734‑585815 |
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: ‑ |