| Predicted mutation | |||||||
|---|---|---|---|---|---|---|---|
| evidence | seq id | position | mutation | freq | annotation | gene | description |
| RA | minE | 2,411,718 | T→C | 10.7% | L330L (TTG→CTG) | hemY → | predicted protoheme IX synthesis protein |
| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | minE | 2,411,718 | 0 | T | C | 10.7% | 156.2 / 9.4 | 56 | L330L (TTG→CTG) | hemY | predicted protoheme IX synthesis protein |
| Reads supporting (aligned to +/- strand): ref base T (34/16); new base C (3/3); total (37/19) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 3.97e-01 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.57e-01 | |||||||||||
CAAACAATCCGGAACAGCTGGAAAAAGTGCTGCGCCAGCAAATCAAAAACGTCGGCGATCGCCCGCTGTTGTGGAGCACACTGGGCCAGTCACTGATGAAGCACGGAGAATGGCAGGAAGCATCGCTCGCCTTCCGCGC > minE/2411650‑2411788 | caaacaaTCCGGAACAGCTGGAAAAAGTGCTGCGCCAGCAAATCAAAAACGTCGGCGATCGCCCGCtgttg < 1:3513659/71‑1 (MQ=255) aacaaTCCGGAACAGCTGGAAAAAGTGCTGCGCCAGCAAATCAAAAACGTCGGCGATCGCCCGCTGTtgtg < 1:581218/71‑1 (MQ=255) acaaTCCGGAACAGCTGGAAAAAGTGCTGCGCCAGCAAATCAAAAACGTCGGCGATCGCCCGCTGTtgtg > 1:3361774/1‑70 (MQ=255) aaTCCGGAACAGCTGGAAAAAGTGCTGCGCCAGCAAATCAAAAACGTCGGCGATCGCCCGCTGTTGTGGAg > 1:3133828/1‑71 (MQ=255) aTCCGGAACAGCTGGAAAAAGTGCTGCGCCAGCAAATCAAAAACGTCGGCGATCGCCCGCTGTTGTGGAg > 1:1880583/1‑70 (MQ=255) ccGGAACAGCTGGAAAAAGTGCTGCGCCAGCAAATCAAAAACGTCGGCGATCGCCCGCtgtt > 1:1867890/1‑62 (MQ=255) aCAGCTGGAAAAAGTGCTGCGCCAGCAAATCAAAAACGTCGGCGATCGCCCGCTGTTGTGGAGCACACTgg < 1:1904804/71‑1 (MQ=255) aCAGCTGGAAAAAGTGCTGCGCCAGCAAATCAAAAACGTCGGCGATCGCCCGCTGTTGTGGAGCACACTgg < 1:3073392/71‑1 (MQ=255) cAGCTGGAAAAAGTGCTGCGCCAGCAAATCAAAAACGTCGGCGATCGCCCGCTGTTGTGGAGCACACTggg > 1:3605242/1‑71 (MQ=255) cTGGAAAAAGTGCTGCGCCAGCAAATCAAAAACGTCGGCGATCGCCCGCTGTTGTGGAGCACACTGGGc > 1:336028/1‑69 (MQ=255) cTGGAAAAAGTGCTGCGCCAGCAAATCAAAAACGTCGGCGATCGCCCGCTGTTGTGGAGCACACTGGGCCa < 1:3518169/71‑1 (MQ=255) tGGAAAAAGTGCTGCGCCAGCAAATCAAAAACGTCGGCGATCGCCCGCtgt > 1:2646774/1‑51 (MQ=255) tGGAAAAAGTGCTGCGCCAGCAAATCAAAAACGTCGGCGATCGCCCGCtgt > 1:1918312/1‑51 (MQ=255) ggAAAAAGTGCTGCGCCAGCAAATCAAAAACGTCGGCGATCGCCCGCTGTtgtg < 1:3671641/54‑1 (MQ=255) ggAAAAAGTGCTGCGCCAGCAAATCAAAAACGTCGGCGATCGCCCGCTGTtgtg < 1:1679766/54‑1 (MQ=255) gAAAAAGTGCTGCGCCAGCAAATCAAAAACGTCGGCGATCGCCCGCTGTTGTGGAGCACACTGGGCCAGTc < 1:1074281/71‑1 (MQ=255) aaaaaGTGCTGCGCCAGCAAATCAAAAACGTCGGCGATCGCCCGCTGTTGTGGAGCACACTGGGCCAGTCa < 1:2686118/71‑1 (MQ=255) aaaGTGCTGCGCCAGCAAATCAAAAACGTCGGCGATCGCCCGCTGTTGTGGAGCACACTGGGCCAGTc > 1:912222/1‑68 (MQ=255) aaaGTGCTGCGCCAGCAAATCAAAAACGTCGGCGATCGCCCGCTGTTGTGGAGCACACTGGGCCAGTCACt > 1:474833/1‑71 (MQ=255) aaGTGCTGCGCCAGCAAATCAAAAACGTCGGCGATCGCCCGCTGTTGTGGAGCACACTGGGCCAGTCACt > 1:973108/1‑70 (MQ=255) aGTGCTGCGCCAGCAAATCAAAAACGTCGGCGATCGCCCGCTGTTGTGGAGCACACTGGGCCAGTCACtga > 1:1401185/1‑71 (MQ=255) gctgcGCCAGCAAATCAAAACCGTCGGCGATCGCCCGCTGTTGTGGAGCACACTGGGCCAGTCACt > 1:1323870/1‑66 (MQ=255) gctgcGCCAGCAAATCAAAAACGTCGGCGATCGCCCGCTGTtgtg < 1:1894841/45‑1 (MQ=255) cgcCAGCAAATCAAAAACGTCGGCGATCGCCCGCTGTTGTGGAGCACACTGGGCCAGTCACTGATGAAg < 1:2142666/69‑1 (MQ=255) cgcCAGCAAATCAAAAACGTCGGCGATCGCCCGCTGTTGTGGAGCACACTGGGCCAGTCACTGATGAAGCa > 1:2050828/1‑71 (MQ=255) gcCAGCAAATCAAAAACGTCGGCGATCGCCCGCtgt > 1:2909586/1‑36 (MQ=255) gcCAGCAAATCAAAAACGTCGGCGATCGCCCGCTGTTGTGGAGCACACTGGGCCAGTCACTGATGAAGCAc > 1:2623543/1‑71 (MQ=255) ccAGCAAATCAAAAACGTCGGCGATCGCCCGCTGTTGTGGAGcaca > 1:1156764/1‑46 (MQ=255) aGCAAATCAAAAACGTCGGCGATCGCCCGCTGTTGTGGAGCACACTGGGCCAGTCACTGATGAAGCAc > 1:2447990/1‑68 (MQ=255) gCAAATCAAAAACGTCGGCGATCGCCCGCTGTTGTGGAGCACACTGGGCCAGTCACTGATGAAg > 1:3117118/1‑64 (MQ=255) gCAAATCAAAAACGTCGGCGATCGCCCGCTGTTGTGGAGCACACTGGGCCAGTCACTGATGAAGCACg > 1:348709/1‑68 (MQ=255) aTCAAAAACGTCGGCGATCGCCCGCTGTTGTGGAGCACACTGGGCCAGTCACtga > 1:714841/1‑55 (MQ=255) aTCAAAAACGTCGGCGATCGCCCGCTGTTGTGGAGCACACTGGGCCAGTCACtga > 1:1230856/1‑55 (MQ=255) aTCAAAAACGTCGGCGATCGCCCGCTGTTGTGGAGCACACTGGGCCAGTCACTGATGAAGCACGGAGAATg > 1:2645248/1‑71 (MQ=255) aTCAAAAACGTCGGCGATCGCCCGCTGTTGTGGAGCACACTGGGCCAGTCACTGATGAAGCACGGAGAATg > 1:3410132/1‑71 (MQ=255) aTCAAAAACGTCGGCGATCGCCCGCTGTTGTGGAGCACACTGGGCCAGTCACTGATGAAGCACGGAGAATg > 1:659351/1‑71 (MQ=255) aTCAAAAACGTCGGCGATCGCCCGCTGTTGTGGAGCACACTGGGCCAGTCACTGATGAAGCACGGAGAATg < 1:1714804/71‑1 (MQ=255) tCAAAAACGTCGGCGATCGCCCGCTGTTGTGGAGCACACTgg < 1:1427471/42‑1 (MQ=255) tCAAAAACGTCGGCGATCGCCCGCTGTTGTGGAGCACACTgg < 1:959998/42‑1 (MQ=255) cAAAAACGTCGGCGATCGCCCGCTGCTGTGGAGCACACTGGGc < 1:3022912/43‑1 (MQ=255) cAAAAACGTCGGCGATCGCCCGCTGCTGTGGAGCACACTGGGc < 1:3489391/43‑1 (MQ=255) cAAAAACGTCGGCGATCGCCCGCTGCTGTGGAGCACACTGGGc < 1:2632948/43‑1 (MQ=255) aaaaaCGTCGGCGATCGCCCGCTGTTGTGGAGCACACTGGGCCAGTCACTGATGAAGCACGGAGAATgg > 1:3172759/1‑69 (MQ=255) aaaaaCGTCGGCGATCGCCCGCTGTTGTGGAGCACACTGGGCCAGTCACTGATGAAGCACGGAGAATGGc < 1:2556074/70‑1 (MQ=255) aaaaCGTCGGCGATCGCCCGCTGTTGTGGAGCACACTGGGCCAGTCACTGATGAACCACGGAGAATGGCAg > 1:2355897/1‑71 (MQ=255) aaaCGTCGGCGATCGCCCGCTGTTGTGGAGCACACTGGGCCAGTCACTGATGAAGCACGgag > 1:625178/1‑62 (MQ=255) aaaCGTCGGCGATCGCCCGCTGTTGTGGAGCACACTGGGCCAGTCACTGATGAAGCACGgag > 1:3487461/1‑62 (MQ=255) aaaCGTCGGCGATCGCCCGCTGTTGTGGAGCACACTGGGCCAGTCACTGATGAAGCACGGAGAATGGCAgg < 1:1081144/71‑1 (MQ=255) aaaCGTCGGCGATCGCCCGCTGTTGTGGAGCACACTGGGCCAGTCACTGATGAAGCACGGAGAATGGCAg > 1:3564143/1‑70 (MQ=255) aaCGTCGGCGATCGCCCGCTGTTGTGGAGCACACTGGGCCAGTCACTGATGAAGCACGGAGAATGGCAgg > 1:573777/1‑70 (MQ=255) aCGTCGGCGATCGCCCGCTGCTGTGGAGCACACTGGGCCAGTCACTGATGAAGCACGGAGAATGGCAgg > 1:1395580/1‑69 (MQ=255) aCGTCGGCGATCGCCCGCTGCTGTGGAGCACACTGGGCCAGTCACTGATGAAGCACGGAGAATGGCAgg > 1:3684795/1‑69 (MQ=255) gCGATCGCCCGCTGTTGTGGAGCACACTGGGCCAGTCACt > 1:439275/1‑40 (MQ=255) cGATCGCCCGCTGTTGTGGAGCACACTGGGCCAGTCACTGATGAAGCACGGAGAATGGCAGGAAGCAtcgc < 1:2266518/71‑1 (MQ=255) tCGCCCGCTGTTGTGGAGCACACTGGGCCAGTCACTGATGAAGCACGGAGAATGGCAGGAAGCa > 1:1707663/1‑64 (MQ=255) cGCCCGCTGTTGTGGAGCACACTGGGCCAGTCACTGATGAAGCACGGAGAATGGCAGGAAGCAtcgctcgc > 1:2727824/1‑71 (MQ=255) cGCCCGCTGTTGTGGAGCACACTGGGCCAGTCACTGATGAAGCACGGAGAATGGCAGGAAGCAtcgctcgc > 1:2070958/1‑71 (MQ=255) cGCCCGCTGCTGTGGAGCACACTGGGCCAGTCACTGATGAAGCACGGAGAATGGCAGGAAGCAtcgctcgc > 1:1464952/1‑71 (MQ=255) cccGCTGTTGTGGAGCACACTGGGCCAGTCACTGATGAAGc > 1:1027414/1‑41 (MQ=255) cccGCTGTTGTGGAGCACACTGGGCCAGTCACTGATGAAGCACGGAGAATGGCAGGAAGCATCGCTCGCCt > 1:874216/1‑71 (MQ=255) cccGCTGTTGTGGAGCACACTGGGCCAGTCACTGATGAAGCACGGAGAATGGCAGGAAGCATCGCTCGCCt > 1:1876233/1‑71 (MQ=255) tgttgtGGAGCACACTGGGCCAGTCACTGATGAAGCACGGAGAATGGCAGGAAGCATCGCTCGCCTTCcgc > 1:1061061/1‑71 (MQ=255) tgttgtGGAGCACACTGGGCCAGTCACTGATGAAGCACGGAGAATGGCAGGAAGCATCGCTCGCCTTCcgc > 1:2945427/1‑71 (MQ=255) ttgtGGAGCACACTGGGCCAGTCACTGATGAAGCACGGAGAATGGCAGGAAGCATCGCTCGcc < 1:1295227/63‑1 (MQ=255) ttgtGGAGCACACTGGGCCAGTCACTGATGAAGCACGGAGAATGGCAGGAAGCATCGCTCGcc < 1:1304644/63‑1 (MQ=255) ttgtGGAGCACACTGGGCCAGTCACTGATGAAGCACGGAGAATGGCAGGAAGCATCGCTCGCCTTCCgcgc < 1:1211084/71‑1 (MQ=255) | CAAACAATCCGGAACAGCTGGAAAAAGTGCTGCGCCAGCAAATCAAAAACGTCGGCGATCGCCCGCTGTTGTGGAGCACACTGGGCCAGTCACTGATGAAGCACGGAGAATGGCAGGAAGCATCGCTCGCCTTCCGCGC > minE/2411650‑2411788 |
| 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: ‑ |