| Predicted mutation | |||||||
|---|---|---|---|---|---|---|---|
| evidence | seq id | position | mutation | freq | annotation | gene | description |
| RA | minE | 2,224,699 | T→A | 100% | P116P (CCT→CCA) | prmA → | methylase for 50S ribosomal subunit protein L11 |
| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | minE | 2,224,699 | 0 | T | A | 88.2% | 124.5 / 8.4 | 59 | P116P (CCT→CCA) | prmA | methylase for 50S ribosomal subunit protein L11 |
| Reads supporting (aligned to +/- strand): ref base T (2/5); new base A (0/52); total (2/57) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 1.23e-02 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.95e-01 | |||||||||||
| Rejected as polymorphism: Variant not supported by required number of reads on each strand. | |||||||||||
AATTTCCACCCGATGCGCTTTGGTGAACGACTGTGGATCTGCCCTAGCTGGCGTGATGTGCCGGACGAAAACGCCGTCAACGTGATGTTAGATCC > minE/2224655‑2224749 | aaTTTCCACCCGATGCGCTTTGGTGAACGACTGTGGATCTGCCCTAGCTGGCGTGATGTGCCGGACGaa < 1:1857715/69‑1 (MQ=255) gATGCGCTTTGGTGAACGACTGTGGATCTGCCCAAGCTGGCGTGATGTGCCGGACGAAAACGCCGTCaa < 1:1158837/69‑1 (MQ=255) gATGCGCTTTGGTGAACGACTGTGGATCTGCCCAAGCTGGCGTGATGTGCCGGACGAAAACGCCGTCaa < 1:3334165/69‑1 (MQ=255) gATGCGCTTTGGTGAACGACTGTGGATCTGCCCAAGCTGGCGTGATGTGCCGGACGAAAACGCCGTCaa < 1:3090125/69‑1 (MQ=255) gATGCGCTTTGGTGAACGACTGTGGATCTGCCCAAGCTGGCGTGATGTGCCGGACGAAAACGCCGTCaa < 1:28959/69‑1 (MQ=255) gATGCGCTTTGGTGAACGACTGTGGATCTGCCCAAGCTGGCGTGATGTGCCGGACGAAAACGCCGTCaa < 1:282655/69‑1 (MQ=255) gATGCGCTTTGGTGAACGACTGTGGATCTGCCCAAGCTGGCGTGATGTGCCGGACGAAAACGCCGTCaa < 1:2776444/69‑1 (MQ=255) gATGCGCTTTGGTGAACGACTGTGGATCTGCCCAAGCTGGCGTGATGTGCCGGACGAAAACGCCGTCaa < 1:2636053/69‑1 (MQ=255) gATGCGCTTTGGTGAACGACTGTGGATCTGCCCAAGCTGGCGTGATGTGCCGGACGAAAACGCCGTCaa < 1:2445282/69‑1 (MQ=255) gATGCGCTTTGGTGAACGACTGTGGATCTGCCCAAGCTGGCGTGATGTGCCGGACGAAAACGCCGTCaa < 1:2219195/69‑1 (MQ=255) gATGCGCTTTGGTGAACGACTGTGGATCTGCCCAAGCTGGCGTGATGTGCCGGACGAAAACGCCGTCaa < 1:2205261/69‑1 (MQ=255) gATGCGCTTTGGTGAACGACTGTGGATCTGCCCAAGCTGGCGTGATGTGCCGGACGAAAACGCCGTCaa < 1:206187/69‑1 (MQ=255) gATGCGCTTTGGTGAACGACTGTGGATCTGCCCAAGCTGGCGTGATGTGCCGGACGAAAACGCCGTCaa < 1:1020217/69‑1 (MQ=255) gATGCGCTTTGGTGAACGACTGTGGATCTGCCCAAGCTGGCGTGATGTGCCGGACGAAAACGCCGTCaa < 1:125209/69‑1 (MQ=255) gATGCGCTTTGGTGAACGACTGTGGATCTGCCCAAGCTGGCGTGATGTGCCGGACGAAAACGCCGTCaa < 1:1491385/69‑1 (MQ=255) gATGCGCTTTGGTGAACGACTGTGGATCTGCCCAAGCTGGCGTGATGTGCCGGACGAAAACGCCGTCaa < 1:124426/69‑1 (MQ=255) gATGCGCTTTGGTGAACGACTGTGGATCTGCACAAGCTGGCGTGATGTGCCGGACGAAAACGCCGTCaa < 1:39177/69‑1 (MQ=255) cgcTTTGGTGAACGACTGTGGATCTGCCCAAGCTGGCGTGATGTGCCGGACGAAAACGCCGTCAACGTg < 1:2695140/69‑1 (MQ=255) cgcTTTGGTGAACGACTGTGGATCTGCCCAAGCTGGCGTGATGTGCCGGACGAAAACGCCGTCAACGTg < 1:157974/69‑1 (MQ=255) cgcTTTGGTGAACGACTGTGGATCTGCCCAAGCTGGCGTGATGTGCCGGACGAAAACGCCGTCAACGTg < 1:2488789/69‑1 (MQ=255) cgcTTTGGTGAACGACTGTGGATCTGCCCAAGCTGGCGTGATGTGCCGGACGAAAACGCCGTCAACGTg < 1:410378/69‑1 (MQ=255) cgcTTTGGTGAACGACTGTGGATCTGCCCAAGCTGGCGTGATGTGCCGGACGAAAACGCCGTCAACGTg < 1:1863547/69‑1 (MQ=255) cgcTTTGGTGAACGACTGTGGATCTGCCCAAGCTGGCGTGATGTGCCGGACGAAAACGCCGTCAACGTg < 1:2265176/69‑1 (MQ=255) tttGGTGAACGACTGTGGATCTGCCCAAGCTGGCGTGATGTGCCGGACGAAAACGCCGTCAACGTGATg < 1:1331709/69‑1 (MQ=255) tttGGTGAACGACTGTGGATCTGCCCAAGCTGGCGTGATGTGCCGGACGAAAACGCCGTCAACGTGATg < 1:1513520/69‑1 (MQ=255) tttGGTGAACGACTGTGGATCTGCCCAAGCTGGCGTGATGTGCCGGACGAAAACGCCGTCAACGTGATg < 1:1274490/69‑1 (MQ=255) tttGGTGAACGACTGTGGATCTGCCCAAGCTGGCGTGATGTGCCGGACGAAAACGCCGTCAACGTGATg < 1:2784622/69‑1 (MQ=255) tttGGTGAACGACTGTGGATCTGCCCAAGCTGGCGTGATGTGCCGGACGAAAACGCCGTCAACGTGATg < 1:3240701/69‑1 (MQ=255) tttGGTGAACGACTGTGGATCTGCCCAAGCTGGCGTGATGTGCCGGACGAAAACGCCGTCAACGTGATg < 1:366385/69‑1 (MQ=255) tttGGTGAACGACTGTGGATCTGCCCAAGCTGGCGTGATGTGCCGGACGAAAACGCCGTCAACGTGATg < 1:396905/69‑1 (MQ=255) tttGGTGAACGACTGTGGATCTGCCCAAGCTGGCGTGATGTGCCGGACGAAAACGCCGTCAACGTGATg < 1:486531/69‑1 (MQ=255) tttGGTGAACGACTGTGGATCTGCCCAAGCTGGCGTGATGTGCCGGACGAAAACGCCGTCAACGTGATg < 1:594065/69‑1 (MQ=255) tttGGTGAACGACTGTGGATCTGCCCAAGCTGGCGTGATGTGCCGGACGAAAACGCCGTCAACGTGATg < 1:690611/69‑1 (MQ=255) tttGGTGAACGACTGTGGATCTGCCCAAGCTGGCGTGATGTGCCGGACGAAAACGCCGTCAACGTGATg < 1:2522201/69‑1 (MQ=255) tttGGTGAACGACTGTGGATCTGCCCAAGCTGGCGTGATGTGCCGGACGAAAACGCCGTCAACGTGATg < 1:1796994/69‑1 (MQ=255) tttGGTGAACGACTGTGGATCTGCCCAAGCTGGCGTGATGTGCCGGACGAAAACGCCGTCAACGTGATg < 1:1980674/69‑1 (MQ=255) tttGGTGAACGACTGTGGATCTGCCCAAGCTGGCGTGATGTGCCGGACGAAAACGCCGTCAACGTGATg < 1:2556843/69‑1 (MQ=255) tttGGTGAACGACTGTGGATCTCCCCAAGCTGGCGTGATGTGCCGGACGAAAACGCCGTCAACGTGATg < 1:1506724/69‑1 (MQ=255) ttGGTGAACGACTGTGGATCTGCCCAAGCTGGCGTGATGTGCCGGACGAAAACGCCGTCAACGTGATg < 1:1442473/68‑1 (MQ=255) ttGGTGAACGACTGTGGATCTGCCCAAGCTGGCGTGATGTGCCGGACGAAAACGCCGTCAAAGTGATg < 1:442094/68‑1 (MQ=255) gTGAACGACTGTGGATCTGCCCTAGCTGGCGTGATGTTCCGGACGAAAACGCCGTCAACGTGAt > 1:1497035/1‑64 (MQ=255) gTGAACGACTGTGGATCTGCCCTAGCTGGCGTGATGTGCCGGACGAAAACGCCGTCAACGTGAt > 1:403692/1‑64 (MQ=255) gAACGACTGTGGATCTGCCCTAGCTGGCGTGATGTGCCGGACGAAAACGCCGTCAACGTGATGTTAGa < 1:2284293/68‑1 (MQ=255) gAACGACTGTGGATCTGCCCTAGCTGGCGTGATGTGCCGGACGAAAACGCCGTCAACGTGATGTTAGa < 1:5858/68‑1 (MQ=255) gAACGACTGTGGATCTGCCCTAGCTGGCGTGATGTGCCGGACGAAAACGCCGTCAACGTGATGTTAGa < 1:2355334/68‑1 (MQ=255) gAACGACTGTGGATCTGCCCTAGCTGGCGTGATGTGCCGGACGAAAACGCCGTCAACGTGATGTTAGa < 1:484033/68‑1 (MQ=255) aaCGACTGTGGATCTGCCCAAGCTGGCGTGATGTGCCGGACGAAAACGCCGTCAACGTGATg < 1:1487230/62‑1 (MQ=255) cGACTGTGGATCTGCCCAAGCTGGCGTGATGTGCCGGACGAAAACGCCGTCAACGTGATGTTAGAccc < 1:1781854/68‑1 (MQ=255) cGACTGTGGATCTGCCCAAGCTGGCGTGATGTGCCGGACGAAAACGCCGTCAACGTGATGTTAGAccc < 1:702906/68‑1 (MQ=255) cGACTGTGGATCTGCCCAAGCTGGCGTGATGTGCCGGACGAAAACGCCGTCAACGTGATGTTAGAccc < 1:1551476/68‑1 (MQ=255) cGACTGTGGATCTGCCCAAGCTGGCGTGATGTGCCGGACGAAAACGCCGTCAACGTGATGTTAGAccc < 1:495677/68‑1 (MQ=255) cGACTGTGGATCTGCCCAAGCTGGCGTGATGTGCCGGACGAAAACGCCGTCAACGTGATGTTAGAccc < 1:1696206/68‑1 (MQ=255) cGACTGTGGATCTGCCCAAGCTGGCGTGATGTGCCGGACGAAAACGCCGTCAACGTGATGTTAGAccc < 1:175003/68‑1 (MQ=255) cGACTGTGGATCTGCCCAAGCTGGCGTGATGTGCCGGACGAAAACGCCGTCAACGTGATGTTAGAccc < 1:3111451/68‑1 (MQ=255) cGACTGTGGATCTGCCCAAGCTGGCGTGATGTGCCGGACGAAAACGCCGTCAACGTGATGTTAGAccc < 1:1233768/68‑1 (MQ=255) cGACTGTGGATCTGCCCAAGCTGGCGTGATGTGCCGGACGAAAACGCCGTCAACGTGATGTTAGAccc < 1:3057200/68‑1 (MQ=255) cGACTGTGGATCTGCCCAAGCTGGCGTGATGTGCCGGACGAAAACGCCGTCAACGTGATGTTAGAccc < 1:2952278/68‑1 (MQ=255) cGACTGTGGAGCTGCCCAAGCTGGCGTGATGTGCCGGACGAAAACGCCGTCAACGTGATGTTAGAccc < 1:1491862/68‑1 (MQ=255) aTCTGCCCAAGCTGGCGTGATGTGCCGGACGAAAACGCCGTCAACGTg < 1:968515/48‑1 (MQ=255) | AATTTCCACCCGATGCGCTTTGGTGAACGACTGTGGATCTGCCCTAGCTGGCGTGATGTGCCGGACGAAAACGCCGTCAACGTGATGTTAGATCC > minE/2224655‑2224749 |
| 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: ‑ |