| Predicted mutation | |||||||
|---|---|---|---|---|---|---|---|
| evidence | seq id | position | mutation | freq | annotation | gene | description |
| RA | minE | 119,298 | A→G | 31.1% | Q57Q (CAA→CAG) | yacL → | conserved hypothetical protein |
| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | minE | 119,298 | 0 | A | G | 31.1% | 51.3 / 28.1 | 61 | Q57Q (CAA→CAG) | yacL | conserved hypothetical protein |
| Reads supporting (aligned to +/- strand): ref base A (24/18); new base G (15/4); total (39/22) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 1.51e-01 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.89e-01 | |||||||||||
TGATGAAGTGGAACAAGCAGCGCACGCACTGAAAGGTAGCGAACGGTCCTGGCAACGGGCAGGGCATGAATACACCCTGTGGATGGACGGTGAAGAGGTGATGGTTCGCGCCAATCA > minE/119244‑119360 | tgatgaAGTGGAACAAGCAGCGCACGCACTGAAAGGTAGCGAACGGTCCTGGCAACGGGCAGGGCATg < 1:2254637/68‑1 (MQ=255)tgatgaAGTGGAACAAGCAGCGCACGCACTGAAAGGTAGCGAACGGTCCTGGCAACGGGCAGGGCATg < 1:2819121/68‑1 (MQ=255) gaAGTGGAAGACGCAGCGCGCACGCTCAAAGGCAGCGAACGATCCTGGCAGCGGGCAGGGCATGAATa < 1:1900362/68‑1 (MQ=255) gaAGTGGAAGACGCAGCGCGCACGCTCAAAGGCAGCGAACGATCCTGGCAGCGGGCAGGGCATGAATa < 1:2544201/68‑1 (MQ=255) gaAGTGGAAGACGCAGCGCGCACGCTCAAAGGCAGCGAACGATCCTGGCAGCGGGCAGGGCATGAATa < 1:1570894/68‑1 (MQ=255) gaAGTGGAAGACGCAGCGCGCACGCTCAAAGGCAGCGAACGATCCTGGCAGCGGGCAGGGCATGAATa < 1:309013/68‑1 (MQ=255) ggAAGACGCAGCGCGCACGCTCAAAGGCAGCGAACGATCCTGGCAGCGGGCAGGGCATGAATacactc > 1:980139/1‑66 (MQ=255) ggAAGACGCAGCGCGCACGCTCAAAGGCAGCGAACGATCCTGGCAGCGGGCAGGGCATGAATACActct > 1:1131717/1‑69 (MQ=255) ggAAGACGCAGCGCGCACGCTCAAAGGCAGCGAACGATCCTGGCAGCGGGCAGGGCATGAATACActct > 1:1913342/1‑69 (MQ=255) ggAAGACGCAGCGCGCACGCTCAAAGGCAGCGAACGATCCTGGCAGCGGGCAGGGCATGAATACActct > 1:3065903/1‑69 (MQ=255) ggAAGACGCAGCGCGCACGCTCAAAGGCAGCGAACGATCCTGGCAGCGGGCAGGGCATGAATACActct > 1:1465438/1‑69 (MQ=255) ggAAGACGCAGCGCGCACGCTCAAAGGCAGCGAACGATCCTGGCAGCGGGCAGGGCATGAATACActct > 1:2976015/1‑69 (MQ=255) ggAAGACGCAGCGCGCACGCTCAAAGGCAGCGAACGATCCTGGCAGCGGGCAGGGCATGAATACActct > 1:1087898/1‑69 (MQ=255) ggAAGACGCAGCGCGCACGCTCAAAGGCAGCGAACGATCCTGGCAGCGGGCAGGGCATGAATACActct > 1:492682/1‑69 (MQ=255) ggAAGACGCAGCGCGCACGCTCAAAGGCAGCGAACGATCCTGGCAGCGGGCAGGGCATGAATACActct > 1:2742559/1‑69 (MQ=255) ggAAGACGCAGCGCGCACGCTCAAAGGCAGCGAACGATCCTGGCAGCGGGCAGGGCATGAATACActct > 1:2696449/1‑69 (MQ=255) ggAAGACGCAGCGCGCACGCTCAAAGGCAGCGAACGATCCTGGCAGCGGGCAGGGCATGAATACActct > 1:1673265/1‑69 (MQ=255) ggAAGACGCAGCGCGCACGCTCAAAGGCAGCGAACGATCCTGGCAGCGGGCAGGGCATGAATACActct > 1:2017031/1‑69 (MQ=255) ggAAGACGCAGCGCGCACGCTCAAAGGCAGCGAACGATCCTGGCAGCGGGCAGGGCATGAATACActct > 1:947099/1‑69 (MQ=255) ggAAGACGCAGCGCGCACGCTCAAAGGCAGCGAACGATCCTGGCAGCGGGCAGGGCATGAATACActct > 1:1013215/1‑69 (MQ=255) ggAAGACGCAGCGCGCACGCTCAAAGGCAGCGAACGATCCTGGCAGCGGGCAGGGCATGAATACActct > 1:2074315/1‑69 (MQ=255) aCAAGCAGCGCACGCACTGAAAGGTAGCGAACGGTCCTGGCAACGGGCAGGGCATGAATACACCCtgtg > 1:2180959/1‑69 (MQ=255) aCAAGCAGCGCACGCACTGAAAGGTAGCGAACGGTCCTGGCAACGGGCAGGGCATGAATACACCCtgtg > 1:2619783/1‑69 (MQ=255) aCAAGCAGCGCACGCACTGAAAGGTAGCGAACGGTCCTGGCAACGGGCAGGGCATGAATACACCCtgtg > 1:363438/1‑69 (MQ=255) aCAAGCAGCGCACGCACTGAAAGGTAGCGAACGGTCCTGGCAACGGGCAGGGCATGAATACACCCtgtg > 1:3155806/1‑69 (MQ=255) aCAAGCAGCGCACGCACTGAAAGGTAGCGAACGGTCCTGGCAACGGGCAGGGCATGAATACACCCtgtg > 1:1891720/1‑69 (MQ=255) aCAAGCAGCGCACGCACTGAAAGGTAGCGAACGGTCCTGGCAACGGGCAGGGCATGAATACACCCtgtg > 1:104995/1‑69 (MQ=255) aCAAGCAGCGCACGCACTGAAAGGTAGCGAACGGTCCTGGCAACGGGCAGGGCATGAATACACCCtgtg > 1:126952/1‑69 (MQ=255) aCAAGCAGCGCACGCACTGAAAGGTAGCGAACGGTCCTGGCAACGGGCAGGGCATGAATACACCCtgtg > 1:1482998/1‑69 (MQ=255) aCAAGCAGCGCACGCACTGAAAGGTAGCGAACGGTCCTGGCAACGGGCAGGGCATGAATACACCCtgtg > 1:1789710/1‑69 (MQ=255) aCAAGCAGCGCACGCACTGAAAGGTAGCGAACGGTCCTGGCAACGGGCAGGGCATGAATACACCCtgtg > 1:1636572/1‑69 (MQ=255) aCAAGCAGCGCACGCACTGAAAGGTAGCGAACGGTCCTGGCAACGGGCAGGGCATGAATACACCCtgtg > 1:1668587/1‑69 (MQ=255) agcagcGCACGCACTGAAAGGTAGCGAACGGTCCTGGCAACGGGCAGGGCATGAATa > 1:2660614/1‑57 (MQ=255) agcagcGCACGCACTGAAAGGTAGCGAACGGTCCTGGCAACGGGCAGGGCATGAATa > 1:1139403/1‑57 (MQ=255) agcagcGCACGCACTGAAAGGTAGCGAACGGTCCTGGCAACGGGCAGGGCATGAATa > 1:3240947/1‑57 (MQ=255) agcagcGCACGCACTGAAAGGTAGCGAACGGTCCTGGCAACGGGCAGGGCATGAATa > 1:463874/1‑57 (MQ=255) agcagcGCACGCACTGAAAGGTAGCGAACGGTCCTGGCAACGGGCAGGGCATGAATa > 1:2813884/1‑57 (MQ=255) agcagcGCACGCACTGAAAGGTAGCGAACGGTCCTGGCAACGGGCAGGGCATGAATa > 1:1060095/1‑57 (MQ=255) agcagcGCACGCACTGAAAGGTAGCGAACGGTCCTGGCAACGGGCAGGGCATGAATa > 1:1759761/1‑57 (MQ=255) agcagcGCACGCACTGAAAGGTAGCGAACGGTCCTGGCAACGGGCAGGGCATGAATa > 1:542744/1‑57 (MQ=255) agcagcGCACGCACTGAAAGGTAGCGAACGGTCCTGGCAACGGGCAGGGCATGAATa > 1:1872224/1‑57 (MQ=255) agcagcGCACGCACTGAAAGGTAGCGAACGGTCCTGGCAACGGGCAGGGCATGAATa > 1:761141/1‑57 (MQ=255) agcagcGCACGCACTGAAAGGTAGCGAACGGTCCTGGCAACGGGCAGGGCATGAATa > 1:1902514/1‑57 (MQ=255) acgcacTGAAAGGTAGCGAACGGTCCTGGCAACGGGCAGGGCATGAATACACCCTGTGGATGGACGGTg < 1:3229804/69‑1 (MQ=255) gcacTGAAAGGTAGCGAACGGTCCTGGCAACGGGCAGGGCATGAATACACCCTGTGGATGGACGGTGaa < 1:1041267/69‑1 (MQ=255) gcacTGAAAGGTAGCGAACGGTCCTGGCAACGGGCAGGGCATGAATACACCCTGTGGATGGACGGTGaa < 1:724980/69‑1 (MQ=255) gcacTGAAAGGTAGCGAACGGTCCTGGCAACGGGCAGGGCATGAATACACCCTGTGGATGGACGGTGaa < 1:3184376/69‑1 (MQ=255) gcacTGAAAGGTAGCGAACGGTCCTGGCAACGGGCAGGGCATGAATACACCCTGTGGATGGACGGTGaa < 1:1134053/69‑1 (MQ=255) gcacTGAAAGGTAGCGAACGGTCCTGGCAACGGGCAGGGCATGAATACACCCTGTGGATGGACGGTGaa < 1:3051024/69‑1 (MQ=255) gcacTGAAAGGTAGCGAACGGTCCTGGCAACGGGCAGGGCATGAATACACCCTGTGGATGGACGGTGaa < 1:1691232/69‑1 (MQ=255) gcacTGAAAGGTAGCGAACGGTCCTGGCAACGGGCAGGGCATGAATACACCCTGTGGATGGACGGTGaa < 1:2299148/69‑1 (MQ=255) gcacTGAAAGGTAGCGAACGGTCCTGGCAACGGGCAGGGCATGAATACACCCTGTGGATGGACGGTGaa < 1:1881490/69‑1 (MQ=255) aGCGAACGGTCCTGGCAACGGGCAGGGCATGAATaca < 1:1099423/37‑1 (MQ=255) aGCGAACGGTCCTGGCAACGGGCAGGGCATGAATaca < 1:649628/37‑1 (MQ=255) aGCGAACGGTCCTGGCAACGGGCAGGGCATGAATaca < 1:703175/37‑1 (MQ=255) aGCGAACGGTCCTGGCAACGGGCAGGGCATGAATaca < 1:1885004/37‑1 (MQ=255) aGCGAACGGTCCTGGCAACGGGCAGGGCATGAATaca < 1:1943392/37‑1 (MQ=255) aGCGAACGGTCCTGGCAACGGGCAGGGCATGAATACACCCTGTGGATGGACGGTGaa < 1:2856109/57‑1 (MQ=255) gCGAACGGTCCTGGCAACGGGCAGGGCATGAATACACCCTGTGGATGGACGGTGAAGAGGTGATGGTTc > 1:2206541/1‑69 (MQ=255) gTCCTGGCAACGGGCAGGGCATGAATACACCCTGTGGATGGACGGTGAAGAGGTGATGGTTCGCGCCa > 1:430825/1‑68 (MQ=255) cTGGCAACGGGCAGGGCATGAATACACCCTGTGGATGGACGGTGAAGAGGTGATGGTTCGCGCCAATCa < 1:1526677/69‑1 (MQ=255) | TGATGAAGTGGAACAAGCAGCGCACGCACTGAAAGGTAGCGAACGGTCCTGGCAACGGGCAGGGCATGAATACACCCTGTGGATGGACGGTGAAGAGGTGATGGTTCGCGCCAATCA > minE/119244‑119360 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |