| Predicted mutation | |||||||
|---|---|---|---|---|---|---|---|
| evidence | seq id | position | mutation | freq | annotation | gene | description |
| RA | W3110S.gb | 4,295,216 | G→A | 54.4% | G38D (GGC→GAC) | nrfE → | heme lyase (NrfEFG) for insertion of heme into c552, subunit NrfE |
| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | W3110S.gb | 4,295,216 | 0 | G | A | 54.4% | 9.0 / 78.6 | 57 | G38D (GGC→GAC) | nrfE | heme lyase (NrfEFG) for insertion of heme into c552, subunit NrfE |
| Reads supporting (aligned to +/- strand): ref base G (5/21); new base A (10/21); total (15/42) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 3.68e-01 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 | |||||||||||
CCATGATGCGACTCACTTGCATCGGCATTCTGGCGCAGTTCGCGCTCCTGCTGCTCGCCTTTGGCGTACTGACGTATTGTTTTCTCATCAGCGATTTCTCGGTCATTTATGTCGCGCAACATAGCTATAGCCTG > W3110S.gb/4295153‑4295286 | ccATGATGCGACTCACTTGCATCGGCATTCTGGCGCAGTTCGCGCTCCTGCTGCTCGCCTTTGACGTACt < 1:44988/70‑1 (MQ=255) atgCGACTCACTTGCATCGGCATTCTGGCGCAGTTCGCGCTCCTGCTGCTCGCCTTTGACGTACTGACGTa > 1:1037973/1‑71 (MQ=255) atgCGACTCACTTGCATCGGCATTCTGGCGCAGTTCGCGCTCCTGCTGCTCGCCTTTGACGTACTGACGTa > 1:1922880/1‑71 (MQ=255) atgCGACTCACTTGCATCGGCATTCTGGCGCAGTTCGCGCTCCTGCTGCTCGCCTTTGACGTACTGACGTa > 1:1934449/1‑71 (MQ=255) atgCGACTCACTTGCATCGGCATTCTGGCGCAGTTCGCGCTCCTGCTGCTCGCCTTTGACGTACTGACGTa > 1:2261188/1‑71 (MQ=255) atgCGACTCACTTGCATCGGCATTCTGGCGCAGTTCGCGCTCCTGCTGCTCGCCTTTGACGTACTGACGTa > 1:2644795/1‑71 (MQ=255) atgCGACTCACTTGCATCGGCATTCTGGCGCAGTTCGCGCTCCTGCTGCTCGCCTTTGACGTACTGACGTa > 1:1281299/1‑71 (MQ=255) atgCGACTCACTTGCATCGGCATTCTGGCGCAGTTCGCGCTCCTGCTGCTCGCCTTTGACGTACTGACGTa > 1:497384/1‑71 (MQ=255) tgCGACTCACTTGCATCGGCATTCTGGCGCAGTTCGCGCTCCTGCTGCTCGCCTTTGGCGTACTGACGTAt < 1:108895/71‑1 (MQ=255) gCGACTCACTTGCATCGGCATTCTGGCGCAGTTCGCGCTGCTGCTGCTCGCCTTTGACGTACTGACGTAtt > 1:3250222/1‑71 (MQ=255) aCTCACTTGCATCGGCATTCTGGCGCAGTTCGCGCTCCTGCTGCTCGCCTTTGACGTACt < 1:332202/60‑1 (MQ=255) aCTCACTTGCATCGGCATTCTGGCGCAGTTCGCGCTCCTGCTGCTCGCCTTTGACGTACt < 1:3000004/60‑1 (MQ=255) aCTCACTTGCATCGGCATTCTGGCGCAGTTCGCGCTCCTGCTGCTCGCCTTTGACGTACt < 1:875671/60‑1 (MQ=255) tCACTTGCATCGGCATTCTGGCGCAGTTCGCGCTCCTGCTGCTCGCCTTTGGCGTACTGACGTATTGtttt > 1:3615544/1‑71 (MQ=255) tCACTTGCATCGGCATTCTGGCGCAGTTCGCGCTCCTGCTGCTCGCCTTTGACGTACTGACGTATTGttt > 1:3110012/1‑70 (MQ=255) cACTTGCATCGGCATTCTGGCGCAGTTCGCGCTCCTGCTGCTCGCCTTTGACGTACt < 1:471858/57‑1 (MQ=255) gCATCGGCATTCTGGCGCAGTTCGCGCTCCTGCTGCTCGCCTTTGGCGTACTGACGTATTGtttt < 1:4084394/65‑1 (MQ=255) tCGGCATTCTGGCGCAGTTCGCGCTCCTGCTGCTCGCCTTTGACGTACTGACGTATTGTTTTCTCATCAGc < 1:3153604/71‑1 (MQ=255) tCGGCATTCTGGCGCAGTTCGCGCTCCTGCTGCTCGCCTTTGACGTACTGACGTATTGTTTTCTCATCAGc < 1:104318/71‑1 (MQ=255) ttCTGGCGCAGTTCGCGCTCCTGCTGCTCGCCTTTGGCGTACTGACGTATTGtttt < 1:850687/56‑1 (MQ=255) ttCTGGCGCAGTTCGCGCTCCTGCTGCTCGCCTTTGGCGTACTGACGTATTGttt < 1:3878407/55‑1 (MQ=255) ttCTGGCGCAGTTCGCGCTCCTGCTGCTCGCCTTTGGCGTACTGACGTATTGttt < 1:4035562/55‑1 (MQ=255) cTGGCGCAGTTCGCGCTCCTGCTGCTCGCCTTTGGCGTACTGACGTATTGTTTTCTCATCAGCGATTtctc > 1:2688593/1‑71 (MQ=255) gcgcAGTTCGCGCTCCTGCTGCTCGCCTTTGGCGTACTGACGTATTGTTTTCTCATCAGCGATTTCTCGGt < 1:3549277/71‑1 (MQ=255) gcgcAGTTCGCGCTCCTGCTGCTCGCCTTTGACGTACTGACGTATTGTTTTCTCATCAGCGATTTCTCGGt < 1:2704358/71‑1 (MQ=255) gcgcAGTTCGCGCTCCTGCTGCTCGCCTTTGACGTACTGACGTATTGTTTTCTCATCAGCGATTTCTCGGt < 1:1372649/71‑1 (MQ=255) gcgcAGTTCGCGCTCCTGCTGCTCGCCTTTGACGTACTGACGTATTGTTTTCTCATCAGCGATTTCTCGGt < 1:1893159/71‑1 (MQ=255) aGTTCGCGCTCCTGCTGCTCGCCTTTGGCGTACTGACGTATTGTTTTCTCATCAGCGATTTCTCGGTCAtt < 1:1941919/71‑1 (MQ=255) aGTTCGCGCTCCTGCTGCTCGCCTTTGGCGTACTGACGTATTGTTTTCTCATCAGCGATTTCTCGGTCAtt < 1:1777496/71‑1 (MQ=255) aGTTCGCGCTCCTGCTGCTCGCCTTTGGCGTACTGACGTATTGTTTTCTCATCAGCGATTTCTCGGTCAtt < 1:3803554/71‑1 (MQ=255) aGTTCGCGCTCCTGCTGCTCGCCTTTGACGTAATGACGTATTGTTTTCTCATCAGCGATTTCTCGGTCAtt < 1:2944302/71‑1 (MQ=255) aGTGCGCGCTCCTGCTGCTCGCCTTTGGCGTACTGACGTAGTGTTTTCTCATCAGCGATTTCTCGGTCAtt < 1:1133719/71‑1 (MQ=255) ttCGCGCTCCTGCTGCTCGCCTTTGGCGTACTGACGTATTGTTTTCTCATCAGCGATTTCTCGGTCATTTa < 1:460759/71‑1 (MQ=255) ttCGCGCTCCTGCTGCTCGCCTTTGGCGTACTGACGTATTGTTTTCTCATCAGCGATTTCTCGGTCATTTa < 1:1903169/71‑1 (MQ=255) tCGCGCTCCTGCTGCTCGCCTTTGACGTACTGACGTATTGTTTTCTCATCAGCGATTTCTCGGTCATTTAt > 1:1346518/1‑71 (MQ=255) cgcgCTCCTGCTGCTCGCCTTTGGCGTACTGACGTAtt > 1:2993414/1‑38 (MQ=255) cgcgCTCCTGCTGCTCGCCTTTGGCGTACTGACGTAtt > 1:1395250/1‑38 (MQ=255) cgcgCTCCTGCTGCTCGCCTTTGGCGTACTGACGTATTGTTTTCtca > 1:1912495/1‑47 (MQ=255) cgcgCTCCTGCTGCTCGCCTTTGACGTACTGACGTATTGTTTTCTCATCAGCGATTTCTCGGTCATTTAt < 1:1642573/70‑1 (MQ=255) cgcgCTCCTGCTGCTCGCCTTTGACGTACTGACCTATTGTTTTCTCATCAGCGATTTCTCGGTCATTTAt < 1:1519642/70‑1 (MQ=255) gcgcTCCTGCTGCTCGCCTTTGGCGTACTGACGTATTGttt < 1:3581551/41‑1 (MQ=255) cgcTCCTGCTGCTCGCCTTTGGCGTACTGACGTATTGTTTTCTCATCAGCGATTTCTCGGTCATTTAt < 1:3145802/68‑1 (MQ=255) cgcTCCTGCTGCTCGCCTTTGGCGTACTGACGTATTGTTTTCTCATCAGCGATTTCTCGGTCATTTAt < 1:690066/68‑1 (MQ=255) tgctgcTCGCCTTTGGCGTACTGACGTATTGTTTTCTCATCAGCGATTTCTCGGTCAtt < 1:340541/59‑1 (MQ=255) tgctgcTCGCCTTTGGCGTACTGACGTATTGTTTTCTCATCAGCGATTTCTCGGTCAtt < 1:178661/59‑1 (MQ=255) gCCTTTGGCGTACTGACGTATTGTTTTCTCATCAGCGATTTCTCGGTCATTTATGTCGCGCAACATAGc < 1:3352006/69‑1 (MQ=255) gCCTTTGGCGTACTGACGTATTGTTTTCTCATCAGCGATTTCTCGGTCATTTATGTCGCGCAACATAGc < 1:1260740/69‑1 (MQ=255) gCCTTTGACGTACTGACGTATTGTTTTCTCATCAGCGATTTCTCGGTCATTTATGTCGCGCAACATAGc < 1:3995694/69‑1 (MQ=255) gCCTTTGACGTACTGACGTATTGTTTTCTCATCAGCGATTTCTCGGTCATTTATGTCGCGCAACATAGCt < 1:3373352/70‑1 (MQ=255) gCCTTTGACGTACTGACGTATTGTTTTCTCATCAGCGATTTCTCGGTCATTTATGTCGCGCAACATAGCt < 1:152306/70‑1 (MQ=255) gCCTTTGACGTACTGACGTATTGTTTTCTCATCAGCGATTTCTCGGTCATTTATGTCGCGCAACATAGCt < 1:1553868/70‑1 (MQ=255) ccTTTGGCGTACTGACGTATTGTTTTCTCATCAGCGATTTCTCGGTCAtt < 1:634896/50‑1 (MQ=255) cTTTGGCGTACTGACGTATTGTTTTCTCATCAGCGATTTCTCGGTCATTTATGTCGCGCAACATAGCtata < 1:1868457/71‑1 (MQ=255) cTTTGACGTACTGACGTATTGTTTTCTCATCAGCGATTTCTCGGTCATTTATGTCGCGCAACATAGCtata < 1:471228/71‑1 (MQ=255) cTTTGACGTACTGACGTATTGTTTTCTCATCAGCGATTTCTCGGTCATTTATGTCGCGCAACATAGCtata < 1:697107/71‑1 (MQ=255) cTTTGACGTACTGACGTATTGTTTTCTCATCAGCGATTTCTCGGTCATTTATGTCGCGCAACATAGCtata < 1:704314/71‑1 (MQ=255) cTTTGACGTACTGACGTATTGTTTTCTCATCAGCGATTTCTCGGTCATTTATGTCGCGCAACATAGCtata < 1:2196775/71‑1 (MQ=255) gCGTACTGACGTATTGTTTTCTCATCAGCGATTTCTCGGTCATTTATGTCGCGCAACATAGCTATAGCctg < 1:231821/71‑1 (MQ=255) gCGTACTGACGTATTGTTTTCTCATCAGCGATTTCTCGGTCATTTATGTCGCGCAACATAGCTATAGCctg < 1:1653875/71‑1 (MQ=255) gCGTACTGACGTATTGTTTTCTCATCAGCGATTTCTCGGTCATTTATGTCGCGCAACATAGCTATAGCctg < 1:1445012/71‑1 (MQ=255) | CCATGATGCGACTCACTTGCATCGGCATTCTGGCGCAGTTCGCGCTCCTGCTGCTCGCCTTTGGCGTACTGACGTATTGTTTTCTCATCAGCGATTTCTCGGTCATTTATGTCGCGCAACATAGCTATAGCCTG > W3110S.gb/4295153‑4295286 |
| 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: ‑ |