| Predicted mutation | |||||||
|---|---|---|---|---|---|---|---|
| evidence | seq id | position | mutation | freq | annotation | gene | description |
| RA | minE | 798,949 | C→T | 24.4% | M214I (ATG→ATA) | treA ← | periplasmic trehalase |
| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | minE | 798,949 | 0 | C | T | 24.4% | 67.2 / 19.2 | 37 | M214I (ATG→ATA) | treA | periplasmic trehalase |
| Reads supporting (aligned to +/- strand): ref base C (11/16); new base T (4/5); total (16/21) | |||||||||||
| 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 | |||||||||||
TTCTTTTTGCATTTGCGGCAGGTATTGCTTCAACGCGGCATCGCCTTCATGCTGCGCCAGTAACTCTACCATCAGGGCAAAGAAGGGCGGTTGCGAGCGGCTTAAATAGTAACTGCGGTTGCCGTTGGGAATATGA > minE/798880‑799015 | ttCTTTTTGCATTTGCGGCAGGTATTGCTTCAACGCGGCATCGCCTTCATGCTGCGCCAGTAACTCTACCa > 1:1563644/1‑71 (MQ=255)ttCTTTTTGCATTTGCGGCAGGTATTGCTTCAACGCGGCATCGCCTTCATGCTGCGCCAGTAACTCTACCa > 1:1995240/1‑71 (MQ=255) tCTTTTTGCATTTGCGGGAGGTATTGCTTCAACGCGGCATCGCCTTCATGCTGCGCCAGTAACTCTACAAt > 1:1939147/1‑71 (MQ=255) tCTTTTTGCATTTGCGGCAGGTATTGCTTCAACGCGGCATCGCCTTCATGCTGCGCCAGTAACTCTACTAt > 1:297066/1‑71 (MQ=255) tttGCATTTGCGGCAGGTATTGCTTCAACGCGGCATCGCCTTCATGCTGCGCCAGTAACTCTACTATCAgg > 1:504275/1‑71 (MQ=255) cATTTGCGGCAGGTATTGCTTCAACGCGGCATCGCCTTCATGCTGCGCCAGTAACTCTACTATCAGGGCaa < 1:492201/71‑1 (MQ=255) tttGCGGCAGGTATTGCTTCAACGCGGCATCGCCTTCATGCTGCGCCAGTAACTCTACCATCAGGGCAaag < 1:2188831/71‑1 (MQ=255) tGCGGCAGGTATTGCTTCAACGCGGCATCGCCTTCATGCTGCGCCAGTAACTCTACTATCAGGGCAaag > 1:1054637/1‑69 (MQ=255) tGCGGCAGGTATTGCTTCAACGCGGCATCGCCTTCATGCTGCGCCAGTAACTCTACCATCAGGGCAaag > 1:1008200/1‑69 (MQ=255) gCGGCAGGTATTGCTTCAACGCGGCATCGCCTTCATGCTGCGCCAGTAACTCTACCATCAGGGCAaagaa < 1:959434/70‑1 (MQ=255) gCGGCAGGTATTGCTTCAACGCGGCATCGCCTTCATGCTGCGCCAGTAACTCTACCATCAGGGCAaagaa < 1:468418/70‑1 (MQ=255) cAGGTATTGCTTCAACGCGGCATCGCCTTCATGCTGCGCCAGTAACTCTACCATCAGGGCAAAGAAGGGCg < 1:1732606/71‑1 (MQ=255) aTTGCTTCAACGCGGCATCGCCTTCATGCTGCGCCAGTAACTCTACCATCAGGGCAaag < 1:757275/59‑1 (MQ=255) ttGCTTCAACGCGGCATCGCCTTCATGCTGCTCCAGTAACTCTACCATCAGGGCAaag < 1:681816/58‑1 (MQ=255) ttGCTTCAACGCGGCATCGCCTTCATGCTGCGCCAGTAACTCTACCATCAGGGCAaaga < 1:683113/59‑1 (MQ=255) ttGCTTCAACGCGGCATCGCCTTCATGCTGCGCCAGTAACTCTACCATCAGGGCAaaga < 1:1707719/59‑1 (MQ=255) ttGCTTCAACGCGGCATCGCCTTCATGCTGCGCCAGTAACTCTACCATCAGGGCAaag < 1:352131/58‑1 (MQ=255) ttGCTTCAACGCGGCATCGCCTTCATGCTGCGCCAGTAACTCTACCATCAGGGCAAAGAAGGGCGGTTGCg < 1:156055/71‑1 (MQ=255) cTTCAACGCGGCATCGCCTTCATGCTGCGCCAGTAACTCTACCATCAGGGCAAAGAAGGGCGGTTGCGAg < 1:2146679/70‑1 (MQ=255) aaCGCGGCATCGCCTTCATGCTGCGCCAGTAACTCTACCATCAGGGCAAAGAAGGGCGGt > 1:766112/1‑60 (MQ=255) cgcgGCATCGCCTTCATGCTGCGCCAGTAACTCTACCATCAGGGCAAAGAAGGGCGGTTGCGAGCGGCTTa < 1:42490/71‑1 (MQ=255) gCATCGCCTTCATGCTGCGCCAGTAACTCTACCATCAGGGCAAAGAAGGGCGGTTGCGAGCGGCTTAAATa < 1:1354229/71‑1 (MQ=255) aTCGCCTTCATGCTGCGCCAGTAACTCTACCATCAGGGCAAAGAAGGGCGGTTGCGAGCGGCTTAAATAGt < 1:775418/71‑1 (MQ=255) cTTCATGCTGCGCCAGTAACTCTACCATCAGGGCAAAGAAGGGCGGTTGCGAGCGGCTTAAATAGt > 1:1453068/1‑66 (MQ=255) cATGCTTCTCCAGTAACTCTACCATCAGGGCAAAGAAgg > 1:1145470/1‑39 (MQ=38) cATGCTGCGCCAGTAACTCTACCATCAGGGCAAAGAAGGGCGGTTGCGAGCGGCTTAAATAg > 1:27109/1‑62 (MQ=255) aGTAACTCTACCATCAGGGCAAAGAAGGGCGGTTGCGAGCGGCTTAAATAGTAACTGCGGTTGc > 1:895564/1‑64 (MQ=255) aGTAACTCTACCATCAGGGCAAAGAAGGGCGGTTGCGAGCGGCTTAAATAGTAACTGCGGTTGCCGTTggg < 1:114943/71‑1 (MQ=255) gTAACTCTACCATCAGGGCAAAGAAGGGCGGTTGCGAGCGGCTTAAATAGTAACTGCGGTTGCCGTTggg > 1:786989/1‑70 (MQ=255) tAACTCTACTATCAGGGCAAAGAAGGGCGGTTGCGAGCGGCTTaaa < 1:1340324/46‑1 (MQ=255) tAACTCTACTATCAGGGCAAAGAAGGGCGGTTGCGAGCGGCTTaaa < 1:224670/46‑1 (MQ=255) tAACTCTACTATCAGGGCAAAGAAGGGCGGTTGCGAGCGGCTTAAATAGTAACTGCGGTTGCCGTTGGGaa > 1:1875956/1‑71 (MQ=255) tAACTCTACCATCAGGGCAAAGAAGGGCGGTTGCGAGCGGCTTaaa < 1:395751/46‑1 (MQ=255) tAACTCTACCATCAGGGCAAAGAAGGGCGGTTGCGAGCGGCTTAAATAGTAACTGCGGTTGCCGTTggg > 1:1872024/1‑69 (MQ=255) aaCTCTACCATCAGGGCAAAGAAGGGCGGTTGCGAGCGGCTTAAATAGTAACTGCGGTTGCCGTTGGGaa > 1:367698/1‑70 (MQ=255) aCTCTACTATCAGGGCAAAGAAGGGCGGTTGCGAGCGGCTTAAATAGTAACTGCGGTTGCCGTTGGGAat < 1:296054/70‑1 (MQ=255) tctACTATCAGGGCAAAGAAGGGCGGTTGCGAGCGGCTTAAATAGTAACTGCGGTTGCCGTTGGGaa < 1:1752988/67‑1 (MQ=255) ctactaTCAGGGCAAAGAAGGGCGGTTGCGAGCGGCTTAAATAGTAACTGCGGTTGCCGTTGGGAATATGa > 1:855438/1‑71 (MQ=255) | TTCTTTTTGCATTTGCGGCAGGTATTGCTTCAACGCGGCATCGCCTTCATGCTGCGCCAGTAACTCTACCATCAGGGCAAAGAAGGGCGGTTGCGAGCGGCTTAAATAGTAACTGCGGTTGCCGTTGGGAATATGA > minE/798880‑799015 |
| 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: ‑ |