Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | minE | 1,467,905:2 | +C | 29.1% | intergenic (+64/+64) | ubiG → / ← yfaL | bifunctional 3‑demethylubiquinone‑9 3‑methyltransferase and 2‑octaprenyl‑6‑hydroxy phenol methylase/adhesin |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | minE | 1,467,905 | 2 | . | C | 29.1% | 162.6 / 81.0 | 55 | intergenic (+64/+64) | ubiG/yfaL | bifunctional 3‑demethylubiquinone‑9 3‑methyltransferase and 2‑octaprenyl‑6‑hydroxy phenol methylase/adhesin |
Reads supporting (aligned to +/- strand): ref base . (26/13); new base C (8/8); total (34/21) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 3.60e-01 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.64e-01 |
GTGAGGTTGATGTTTGGCCGCGCCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACAAAT‑‑‑‑GCTCC‑‑CCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAACCGTGCCGACTAGACAGTGATTACCAT > minE/1467838‑1467975 | gTGAGGTTGATGTTTGGCCGCGCCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACAAAT‑‑‑‑GCt < 1:1940781/71‑1 (MQ=255) gTTGATGTTTGGCCGCGCCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACAAAT‑‑‑‑GCTcc‑‑cc < 1:843942/70‑1 (MQ=255) tGATGTTTGGCCGCGCCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACAAAT‑‑‑‑GCTCC‑‑CCGt > 1:2447871/1‑70 (MQ=255) tGTTTGGCCGCGCCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACAAAT‑‑‑‑GCTCC‑‑CCGTAGGc > 1:3144410/1‑71 (MQ=255) tGTTTGGCCGCGCCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACAAAT‑‑‑‑GCTCC‑‑CCGTAGGc < 1:2023096/71‑1 (MQ=255) tGGCCGCGCCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACAAAT‑‑‑‑GCTCC‑‑CCGTAGGCCGGa > 1:3473915/1‑71 (MQ=255) cgcgCCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACAAAT‑‑‑‑GCTCC‑‑CCGTAGGCCGGATAAg < 1:3733454/71‑1 (MQ=255) cgcgCCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACAAAT‑‑‑‑GCTCC‑‑CCGTAGGCCGGATAAg < 1:468940/71‑1 (MQ=255) cgcCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACAAAT‑‑‑‑GCTCC‑‑CCGTAGGCCGGATAAGGc < 1:150269/71‑1 (MQ=255) ccAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACAAAT‑‑‑‑GCTCC‑‑CCGTAGGCCGGATAAGGCGt > 1:1534616/1‑71 (MQ=255) ccAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACAAAT‑‑‑‑GCTCC‑‑CCGTAGGCCGGATAAGGCGt > 1:5322616/1‑71 (MQ=255) cAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACAAAT‑‑‑‑GCTCC‑‑CCGTAGGCCGGATAAGGCGtt < 1:1838129/71‑1 (MQ=255) cAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACAAAT‑‑‑‑GCTCC‑‑CCGTAGGCCGGATAAGGCGtt < 1:647468/71‑1 (MQ=255) aaTGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACAAAT‑‑‑‑GCTCC‑‑CCGTAGGCCGGATAAGGCGtt > 1:125423/1‑70 (MQ=255) ccTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACAAAT‑‑‑‑GCTCC‑‑CCGt > 1:5569488/1‑49 (MQ=255) tGATGCGACGCTTGCCGCGTCTTATCAGGCCTACAAAT‑‑‑‑GCTCC‑‑CCg < 1:2097987/46‑1 (MQ=255) aTGCGACGCTTGCCGCGTCTTATCAGGCCTACAAAT‑‑‑‑GCTCC‑‑CCGTAGGCCGGATAAGGCGTTTACGCCGCa > 1:2453231/1‑71 (MQ=255) aTGCGACGCTTG‑CGCGTCTTATCAGGCCTACAAATTCCCGCACCCTCCGTAGGCCGGATAAGGCGTTTAcg > 1:4890722/1‑71 (MQ=0) aTGCGACGCTTG‑CGCGTCTTATCAGGCCTACAAATTCCCGCACCCTCCGTAGGCCGGATAAGGCGTTTAcg > 1:4963371/1‑71 (MQ=0) aTGCGACGCTTG‑CGCGTCTTATCAGGCCTACAAATTCCCGCACCCTCCGTAGGCCGGATAAGGCGTTTAcg < 1:5374231/71‑1 (MQ=0) aTGCGACGCTTG‑CGCGTCTTATCAGGCCTACAAATTCCCGCACCCTCCGTAGGCCGGATAAGGCGTTTAc < 1:4719260/70‑1 (MQ=0) tGCGACGCTTGCCGCGTCTTATCAGGCCTACAAAT‑‑‑‑GCTCC‑‑CCg > 1:1700515/1‑43 (MQ=255) tGCGACGCTTGCCGCGTCTTATCAGGCCTACAAAT‑‑‑‑GCTCC‑‑CCg > 1:37603/1‑43 (MQ=255) tGCGACGCTTGCCGCGTCTTATCAGGCCTACAAAT‑‑‑‑GCTCC‑‑CCGTAGGCCGGATAAGGCGTTTACGCCGCa < 1:4427823/70‑1 (MQ=255) tGCGACGCTTG‑CGCGTCTTATCAGGCCTACAAATTCCCGCACCCTCCGTAGGCCGGATAAGGCGTTTAcgc > 1:1345902/1‑71 (MQ=0) gCGACGCTTG‑CGCGTCTTATCAGGCCTACAAATTCCCGCACCATCCGTAGGCCGGATAAGGCGTTTAcgcc > 1:29656/1‑71 (MQ=0) gCGACGCTTG‑CGCGTCTTATCAGGCCTACAAATTCCCGCACCATCCGTAGGCCGGATAAGGCGTTTAcgcc > 1:4485858/1‑71 (MQ=0) cGACGCTTG‑CGCGTCTTATCAGGCCTACAAATTCCCGCACCATCCGTAGGCCGGATAAGGCGTTTAcgcc < 1:4567717/70‑1 (MQ=0) aCGCTTGCCGCGTCTTATCAGGCCTACAAAT‑‑‑‑GCTCC‑‑CCGTAGGCCGGATAAGGCGTTTACGCCGCATCCgg > 1:1858644/1‑71 (MQ=255) aCGCTTG‑CGCGTCTTATCAGGCCTACAAATTCCCGCACCCTCCGTAGGCCGGATAAGGCGTTTACGCCGCa > 1:1204526/1‑71 (MQ=0) aCGCTTG‑CGCGTCTTATCAGGCCTACAAATTCCCGCACCCTCCGTAGGCCGGATAAGGCGTTTACGCCGCa > 1:3608400/1‑71 (MQ=0) aCGCTTG‑CGCGTCTTATCAGGCCTACAAATTCCCGCACCCTCCGTAGGCCGGATAAGGCGTTTACGCCGCa < 1:3533806/71‑1 (MQ=0) cGCTTG‑CGCGTCTTATCAGGCCTACAAATTCCCGCACCCTCCGTAGGCCGGATAAGGCGTTTACGCCGCAt < 1:5289091/71‑1 (MQ=0) cGCTTG‑CGCGTCTTATCAGGCCTACAAATTCCCGCACCCTCCGTAGGCCGGATAAGGCGTTTACGCCGCAt < 1:4720923/71‑1 (MQ=0) gCTTGCCGCGTCTTATCAGGCCTACAAAT‑‑‑‑GCTCC‑‑CCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGc > 1:3375370/1‑70 (MQ=255) gCTTG‑CGCGTCTTATCAGGCCTACAAATTCCCGCACCATCCGTAGGCCGGATAAGGCGTTTACGCCGCATc < 1:445947/71‑1 (MQ=0) tgcgcgTCTTATCAGGCCTACAAATTCCCGCACCCTCCGTAGGCCGGATAAGGCGTTTACGCCGCATCCgg > 1:3564206/3‑71 (MQ=2) acgcgTCTTATCAGGCCTACAAATTCCCGCGCCATCCGTAGGCCGGATAAGGCGTTTACGCCGCATCCg < 1:2023343/68‑1 (MQ=0) gcgTCTTATCAGGCCTACAAAT‑‑‑‑GCTCC‑‑CCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAACCGTGc > 1:4659247/1‑71 (MQ=255) gcgTCTTATCAGGCCTACAAAT‑‑‑‑GCTCC‑‑CCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAACCGTGc > 1:338519/1‑71 (MQ=255) gcgTCTTATCAGGCCTACAAAT‑‑‑‑GCTCC‑‑CCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAACCGTGc > 1:5225554/1‑71 (MQ=255) gcgTCTTATCAGGCCTACAAAT‑‑‑‑GCTCC‑‑CCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAACCGTGc > 1:4990484/1‑71 (MQ=255) cgTCTTATCAGGCCTACAAAT‑‑‑‑GCTCC‑‑CCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAACCGTGcc < 1:3636583/71‑1 (MQ=255) ttATCAGGCCTACAAAT‑‑‑‑GCTCC‑‑CCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAACCGTGCCGa > 1:4731841/1‑69 (MQ=255) ttATCAGGCCTACAAAT‑‑‑‑GCTCC‑‑CCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAACCGTGCCGAc > 1:5660861/1‑70 (MQ=255) ttATCAGGCCTACAAAT‑‑‑‑GCTCC‑‑CCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAACCGTGCCGAc > 1:2112006/1‑70 (MQ=255) tACAAAT‑‑‑‑GCTCC‑‑CCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAACCGTGCCGACTAGACAGTGAt > 1:545890/1‑71 (MQ=255) cAAAT‑‑‑‑GCTCC‑‑CCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAAc > 1:2670849/1‑49 (MQ=255) cAAAT‑‑‑‑GCTCC‑‑CCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAACCGTGCCGACTAGACAGTGATTa > 1:1377418/1‑71 (MQ=255) cAAAT‑‑‑‑GCTCC‑‑CCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAACCGTGCCGACTAGACAGTGATTa > 1:1430356/1‑71 (MQ=255) cAAAT‑‑‑‑GCTCC‑‑CCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAACCGTGCCGACTAGACAGTGATTa > 1:5600342/1‑71 (MQ=255) aaaT‑‑‑‑GCTCC‑‑CCGTAGGCCGGATAAGGCGTTTACGCCGCa < 1:108964/39‑1 (MQ=255) aaaT‑‑‑‑GCTCC‑‑CCGTAGGCCGGATAAGGCGTTTACGCCGCa < 1:3969461/39‑1 (MQ=255) aaaT‑‑‑‑GCTCC‑‑CCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAACCGTGCCGACTAGACAGTGATTAc > 1:695619/1‑71 (MQ=255) aT‑‑‑‑GCTCC‑‑CCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAAACGTGCCGACTAGACAGTGATTACCa > 1:5178548/1‑71 (MQ=255) t‑‑‑‑GCTCC‑‑CCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAACCGTGCCGACTAGACAGTGATTACCAt < 1:2224276/71‑1 (MQ=255) | GTGAGGTTGATGTTTGGCCGCGCCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACAAAT‑‑‑‑GCTCC‑‑CCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAACCGTGCCGACTAGACAGTGATTACCAT > minE/1467838‑1467975 |
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: ‑ |