| Predicted mutation | ||||||
|---|---|---|---|---|---|---|
| evidence | seq id | position | mutation | annotation | gene | description |
| RA | NC_000913 | 1,799,179 | Δ1 bp | intergenic (‑237/+47) | pheS ← / ← pheM | phenylalanine‑‑tRNA ligase subunit alpha/pheST‑ihfA operon leader peptide |
| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | NC_000913 | 1,799,177 | 0 | T | . | 93.3% | 54.4 / ‑1.0 | 15 | intergenic (‑235/+49) | pheS/pheM | phenylalanine‑‑tRNA ligase subunit alpha/pheST‑ihfA operon leader peptide |
| Reads supporting (aligned to +/- strand): ref base T (1/0); new base . (8/6); total (9/6) | |||||||||||
| 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 | |||||||||||
GCTGCGACAATTACACCTGATGAGACAGGCTTTTTATTTTTCAAAACGCGCATACAAAAAAAGCCTCCACTGGGAGGCTTTCAGGCGCTGTTTTCCGTTTCTCTTCTCACGCGCTAGCCTCCTGGATTCAGGTGCTAAAGTAAAAAAAGAAGC > NC_000913/1799099‑1799251 | gCTGCGACAATTACACCTGATGAGACAGGCTTTTTATTTTTCAAAACGCGCATACAAAAAAAGCCTCCACTGGGAGGCtt < 1:323467/80‑1 (MQ=255) cTGCGACAATTACACCTGATGAGACAGGCTTTTTATTTTTCAAAACGCGCATACAAAAAAAGCCTCCACTGGGAGGCttca < 2:323468/81‑3 (MQ=255) tGCGACAATTACACCTGATGAGACAGGCTTTTTATTTTTCAAAACGCGCATACAAAAAAAGCCTCCACTGGGAGGCttcag > 1:323498/1‑78 (MQ=255) atgaGACAGGCTTTTTATTTTTCAAAACGCGCATACAAAAAAAGCCTCCACTGGGAGGC‑TTCAGGCGCTGTTTTCCGttt < 1:323469/80‑1 (MQ=255) atgaGACAGGCTTTTTATTTTTCAAAACGCGCATACAAAAAAAGCCTCCACTGGGAGGC‑TTCAGGCGCTGTTTTCCGTttc > 1:323504/1‑81 (MQ=255) ttttATTTTTCAAAACGCGCATAC‑AAAAAAGCCTCCACTGGGAGGCTTT‑AGGCGCCGGTTTCCGTTTCCCTTCTCAcgcgg > 2:323493/1‑80 (MQ=255) tttATTTTTCAAAACGCGCATACAAAAAAAGCCTCCACTGGGAGGC‑TTCAGGCGCTGTTTTCCGTTTCTCTTCTCAcgc < 1:323470/79‑1 (MQ=255) aTTTTTCAAAACGCGCATACAAAAAAAGCCTCCACTGGGAGGC‑TTCAGGCGCTGTTTTCCGTTTCTCTTCTCACGCGCTAg > 1:323502/1‑81 (MQ=255) tCAAAACGCGCATACAAAAAAAGCCTCCACTGGGAGGC‑TTCAGGCGCTGTTTTCCGTTTCTCTTCTCACGCGCTAGcctc < 2:323471/80‑1 (MQ=255) aaaaCGCGCATACAAAAAAAGCCTCCACTGGGAGGC‑TTCAGGCGCTGTTTTCCGTTTCTCTTCTCACGCGCTAGcctcct < 2:323472/80‑1 (MQ=255) aaaaCGCGCATACAAAAAAAGCCTCCACTGGGAGGC‑TTCAGGCGCTGTTTTCCGTTTCTCTTCTCACGCGCTAACCTCCTg > 2:323507/1‑81 (MQ=255) cgcATACAAAAAAAGCCTCCACTGGGAGGC‑TTCAGGCGCTGTTTTCCGTTTCTCTTCTCACGCGCTAGCCTCCTGGATTCa < 1:323473/81‑1 (MQ=255) cAAAAAAAGCCTCCACTGGGAGGC‑TTCAGGCGCTGTTTTCCGTTTCTCTTCTCACGCGCTAGCCTCCTGGATTCAGGTGc > 2:323503/1‑80 (MQ=255) aGCCTCCACTGGGAGGC‑TTCAGGCGCTGTTTTCCGTTTCTCTTCTCACGCGCTAGCCTCCTGGATTCAGGTGCTAAAGTaa < 2:323474/81‑1 (MQ=255) tCCACTGGGAGGC‑TTCAGGCGCTGTTTTCCGTTTCTCTTCTCACGCGCTAGCCTCCTGGATTCAGGTGCTAAAGTaaaaa > 2:323501/1‑80 (MQ=255) ccACTGGGAGGC‑TTCAGGCGCTGTTTTCCGTTTCTCTTCTCACGCGCTAGCCTCCTGGATTCAGGTGCTAAAGTaaaaa > 2:323490/1‑79 (MQ=255) aCTGGGAGGC‑TTCAGGCGCTGTTTTCCGTTTCTCTTCTCACGCGCTAGCCTCCTGGATTCAGGTGCTAAAGTAAAAAaag > 1:323499/1‑80 (MQ=255) gggAGGC‑TTCAGGCGCTGTTTTCCGTTTCTCTTCTCACGCGCTAGCCTCCTGGATTCAGGTGCTAAAGTAAAAAAAGAAGc > 2:323518/1‑81 (MQ=255) | GCTGCGACAATTACACCTGATGAGACAGGCTTTTTATTTTTCAAAACGCGCATACAAAAAAAGCCTCCACTGGGAGGCTTTCAGGCGCTGTTTTCCGTTTCTCTTCTCACGCGCTAGCCTCCTGGATTCAGGTGCTAAAGTAAAAAAAGAAGC > NC_000913/1799099‑1799251 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 31 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |