| Predicted mutation | |||||||
|---|---|---|---|---|---|---|---|
| evidence | seq id | position | mutation | freq | annotation | gene | description |
| RA | minE | 763,302 | G→A | 100% | P133L (CCA→CTA) | potB ← | polyamine transporter subunit |
| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | minE | 763,302 | 0 | G | A | 78.6% | 25.3 / 3.2 | 14 | P133L (CCA→CTA) | potB | polyamine transporter subunit |
| Reads supporting (aligned to +/- strand): ref base G (3/0); new base A (0/11); total (3/11) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 2.75e-03 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 | |||||||||||
| Rejected as polymorphism: Variant not supported by required number of reads on each strand. | |||||||||||
ACCGATAATCACCGCACTGGGCGTGAACATGATACGGATTGGTGTGTCGATAACGCCCAGCCAGAGCAAAAACTCGTTGAGATAGCCTTT > minE/763262‑763351 | aCCGATAATCACCGCACTGGGCGTGAACATGATACGGATTAGTGTGTCGAt < 1:121153/51‑1 (MQ=255)aCCGATAATCACCGCACTGGGCGTGAACATGATACGGATTAGTGTGTCGAt < 1:15160/51‑1 (MQ=255)aCCGATAATCACCGCACTGGGCGTGAACATGATACGGATTAGTGTGTCGAt < 1:2136/51‑1 (MQ=255)aCCGATAATCACCGCACTGGGCGTGAACATGATACGGATTAGTGTGTCGAt < 1:287594/51‑1 (MQ=255)aCCGATAATCACCGCACTGGGCGTGAACATGATACGGATTAGTGTGTCGAt < 1:344146/51‑1 (MQ=255)aCCGATAATCACCGCACTGGGCGTGAACATGATACGGATTAGTGTGTCGAt < 1:397071/51‑1 (MQ=255)aCCGATAATCACCGCACTGGGCGTGAACATGATACGGATTAGTGTGTCGAt < 1:413747/51‑1 (MQ=255)aCCGATAATCACCGCACTGGGCGTGAACATGATACGGATTAGTGTGTCGAt < 1:470970/51‑1 (MQ=255)aCCGATAATCACCGCACTGGGCGTGAACATGATACGGATTAGTGTGTCGAt < 1:48394/51‑1 (MQ=255)aCCGATAATCACCGCACTGGGCGTGAACATGATACGGATTAGTGTGTCGAt < 1:529608/51‑1 (MQ=255)aCCGATAATCACCGCACTGGGCGTGAACATGATACGGATTAGTGTGTCGAt < 1:72231/51‑1 (MQ=255) ggCGTGAACATGATACGGATTGGTGTGTCGATAACGCCCAGCCAGAGCAAAAACTCGTTGAGATAGCCttt > 1:303376/1‑71 (MQ=255) ggCGTGAACATGATACGGATTGGTGTGTCGATAACGCCCAGCCAGAGCAAAAACTCGTTGAGATAGCCttt > 1:345681/1‑71 (MQ=255) ggCGTGAACATGATACGGATTGGTGTGTCGATAACGCCCAGCCAGAGCAAAAACTCGTTGAGATAGCCtt > 1:320536/1‑70 (MQ=255) | ACCGATAATCACCGCACTGGGCGTGAACATGATACGGATTGGTGTGTCGATAACGCCCAGCCAGAGCAAAAACTCGTTGAGATAGCCTTT > minE/763262‑763351 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |