| Predicted mutation | |||||||
|---|---|---|---|---|---|---|---|
| evidence | seq id | position | mutation | freq | annotation | gene | description |
| RA | Exported | 2,823,871 | A→G | 9.2% | intergenic (+616/+342) | nrfG → / ← pheU | heme lyase (NrfEFG) for insertion of heme into c552, subunit NrfG/tRNA‑Phe |
| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | Exported | 2,823,871 | 0 | A | G | 9.2% | 152.4 / 6.7 | 54 | intergenic (+616/+342) | nrfG/pheU | heme lyase (NrfEFG) for insertion of heme into c552, subunit NrfG/tRNA‑Phe |
| Reads supporting (aligned to +/- strand): ref base A (24/25); new base G (2/3); total (26/28) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 1.00e+00 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.85e-01 | |||||||||||
TGCCAAAGAAAAATGTTCAGAAAATAATCCATGGAAAATTGTCCGGAGCACTTACTATTTTAATGGATTATTAGTCTTTGCATGAGCAAGCGAACTGATACATTTCTCTTTGTTCTCATTCAGAAAATCTCATCAGT > Exported/2823802‑2823938 | tGCCAAAGAAAAATGTTCAGAAAATAATCCATGGAAAATTGTCCGGAGCACTTACTATTTTAATGGatta < 1:1616382/70‑1 (MQ=255) gAAAAATGTTCAGAAAATAATCCATGGAAAATTGTCCGGAGCACTTACTATTTTAATGGATTATTAGTCtt > 1:603626/1‑71 (MQ=255) gTTCAGAAAATAATCCATGGAAAATTGTCCGGAGCACTTACTATTTTAATGGATTATTAGTCTTTGCAt > 1:715542/1‑69 (MQ=255) ttCAGAAAATAATCCATGGAAAATTGTCCGGAGCACTTACTATTTTAATGGATTATTAGTCTTTGCATGAg < 1:1092205/71‑1 (MQ=255) ttCAGAAAATAATCCATGGAAAATTGTCCGGAGCACTTACTATTTTAATGGATTATTAGTCTTTGCATGAg < 1:2000238/71‑1 (MQ=255) aGAAAATAATCCATGGAAAATTGTCCGGAGCACTTACTATTTTAATGGATTATTAGTCTTTGCATGAGCa < 1:39662/70‑1 (MQ=255) gAAAATAATCCATGGAAAATTGTCCGGAGCACTTACTATTTTAATGGATTATTAGTCTTTGCATGAGCa < 1:3501584/69‑1 (MQ=255) aaaTAATCCATGGAAAATTGTCCGGAGCACTTACTATTTTAATGGATTATTAGTCTTTGCATGAGCAAGCg < 1:549668/71‑1 (MQ=255) aataatCCATGGAAAATTGTCCGGAGCACTTACTATTTTAATGGATTATTAGTCTTTGCATGAGCAAGCGa < 1:1023178/71‑1 (MQ=255) ataatCCATGGAAAATTGTCCGGAGCACTTACTATTTTAATGGATTATTAGTCTTTGCATGAGCAAGCGaa < 1:2426717/71‑1 (MQ=255) ataatCCATGGAAAATTGTCCGGAGCACTTACTATTTTAATGGATTATTAGTCTTTGCATGAGCAAGCGaa < 1:648093/71‑1 (MQ=255) atCCATGGAAAATTGTCCGGAGCACTTACTATTTTAATGGATTATTAGTCTTTGCATGAGCAAGCGAACTg < 1:2364189/71‑1 (MQ=255) ccATGGAAAATTGTCCGGAGCACTTACTATTTTAATGGATTATTAGTCTTTGCATGAGCAAGCGAACTGa > 1:1837396/1‑70 (MQ=255) gAAAATTGTCCGGAGCACTTACTATTTTAATGGATTGTTAGTc < 1:3477032/43‑1 (MQ=255) gAAAATTGTCCGGAGCACTTACTATTTTAATGGATTATTAGTc > 1:3042167/1‑43 (MQ=255) gAAAATTGTCCGGAGCACTTACTATTTTAATGGATTATTAGTCTTTGCATGAGCAAGCGa < 1:272795/60‑1 (MQ=255) aaaTTGTCCGGAGCACTTACTATTTTAATGGATTATTAGTCTTTGCATGAGCAAGCGAAc > 1:917854/1‑60 (MQ=255) aaTTGTCCGGAGCACTTACTATTTTAATGGATTATTAGTCTTTGCAt > 1:47448/1‑47 (MQ=255) aaTTGTCCGGAGCACTTACTATTTTAATGGATTATTAGTCTTTGCATGAGCAAGCGAACTg < 1:4813/61‑1 (MQ=255) aaTTGTCCGGAGCACTTACTATTTTAATGGATTATTAGTCTTTGCATGAGCAAGCGAACTGATACATTtct < 1:2197948/71‑1 (MQ=255) aaTTGTCCGGAGCACTTACTATTTTAATGGATTATTAGTCTTTGCATGAGCAAGCGAACTGATACATTtc > 1:2144430/1‑70 (MQ=255) ttGTCCGGAGCACTTACTATTTTAATGGATTATTAGTCTTTGCATGAGCAAGCg < 1:2720594/54‑1 (MQ=255) ttGTCCGGAGCACTTACTATTTTAATGGATTATTAGTCTTTGCATGAGCAAGCGAACTGATACATTTctct < 1:469026/71‑1 (MQ=255) gTCCGGAGCACTTACTATTTTAATGGATTATTAGTCTTTGCATGAGCAAGCGaa > 1:3443946/1‑54 (MQ=255) cGGAGCACTTACTATTTTAATGGATTGTTAGTCTTTGCATGAGCAAGCGAACTGATACATTTCTCTTTGt > 1:1221707/1‑70 (MQ=255) ggAGCACTTACTATTTTAATGGATTATTAGTCTTTGCATGAGCAAGCGAACTGATACATTtctc > 1:168609/1‑64 (MQ=255) gCACTTACTATTTTAATGGATTATTAGTCTTTGCATGAGCAAGCGAACTGATACATTTCTCttt < 1:751095/64‑1 (MQ=255) cTTACTATTTTAATGGATTATTAGTCTTTGCATGAGCAAGCg < 1:2947295/42‑1 (MQ=255) cTTACTATTTTAATGGATTATTAGTCTTTGCATGAGCAAGCGAACTg < 1:640988/47‑1 (MQ=255) cTTACTATTTTAATGGATTATTAGTCTTTGCATGAGCAAGCGAACTg < 1:1635917/47‑1 (MQ=255) ttACTATTTTAATGGATTGTTAGTCTTTGCATGAGCAAGCGAACTTAt < 1:994949/48‑1 (MQ=255) ttACTATTTTAATGGATTATTAGTCTTTGCATGAGCAAGCg < 1:528674/41‑1 (MQ=255) ttACTATTTTAATGCATTGTTAGTCTTTGCATGAGCAAGCGAACTGAt < 1:517769/48‑1 (MQ=255) tACTATTTTAATGGATTATTAGTCTTTGCATGAGCAAGCGaa > 1:2264549/1‑42 (MQ=255) tACTATTTTAATGGATTATTAGTCTTTGCATGAGCAAGCGaa > 1:2507491/1‑42 (MQ=255) tACTATTTTAATGGATTATTAGTCTTTGCATGAGCAAGCGAAc > 1:2024131/1‑43 (MQ=255) tACTATTTTAATGGATTATTAGTCTTTGCATGAGCAAGCGAAc > 1:433843/1‑43 (MQ=255) tACTATTTTAATGGATTATTAGTCTTTGCATGAGCAAGCGAACTGAt > 1:1151220/1‑47 (MQ=255) tACTATTTTAATGGATTATTAGTCTTTGCATGAGCAAGCGAACTGAt > 1:3494341/1‑47 (MQ=255) tACTATTTTAATGGATTATTAGTCTTTGCATGAGCAAGCGAACTGAt > 1:2577231/1‑47 (MQ=255) tACTATTTTAATGGAATATTAGTCTTTGCATGAGCAAGCGAACTGATACATTTCTCTTTGTTCTCAtt > 1:283601/1‑68 (MQ=255) aCTATTTTAATGGATTATTAGTCTTTGCATGAGCAAGCGAACTGATACATTTCTCTTTGTtc < 1:2212647/62‑1 (MQ=255) cTATTTTAATGGATTGTTAGTCTTTGCATGAGCAAGCGAACTGATACATTTCTCTTTGTTCTCATTCAGa > 1:1646200/1‑70 (MQ=255) cTATTTTAATGGATTATTAGTCTTTGCATGAGCAAGCGAACTGATACATTTCTCTTTGTTCTCATTCAGa > 1:3130664/1‑70 (MQ=255) tATTTTAATGGATTATTAGTCTTTGCATGAGCAAGCGAACTGATACATTTCTCTTTGTtct < 1:3295082/61‑1 (MQ=255) tATTTTAATGGATTATTAGTCTTTGCATGAGCAAGCGAACTGATACATTTCTCTTTGTTCTCa > 1:3605607/1‑63 (MQ=255) tATTTTAATGGATTATTAGTCTTTGCATGAGCAAGCGAACTGATACATTTCTCTTTGTTCTCa > 1:3435820/1‑63 (MQ=255) aTTTTAATGGATTATTAGTCTTTGCATGAGCAAGCGaa > 1:1098113/1‑38 (MQ=255) ttttAATGGATTATTAGTCTTTGCATGAGCAAGCGAACTGATACATTTCTCTTTGTTCTCa > 1:2972981/1‑61 (MQ=255) ttttAATGGATTATTAGTCTTTGCATGAGCAAGCGAACTGATACATTTCTCTTTGTTCTCa > 1:3093234/1‑61 (MQ=255) tttAATGGATTATTAGTCTTTGCATGAGCAAGCGAACTGATACATTTCTCTTTGTTCTCATTc < 1:3373701/63‑1 (MQ=255) aTGGATTATTAGTCTTTGCATGAGCAAGCGAACTGATACATTTCTCTTTGTTCTCATTCAGAAAAtctc < 1:3703379/69‑1 (MQ=255) tGGATTATTAGTCTTTGCATGAGCAAGCGAACTGATACATTTCTCTTTGTTCTCATTCAGAAAATCtcatc < 1:2692518/71‑1 (MQ=255) tGGATTATTAGTCTTTGCATGAGCAAGCGAACTGATACATTTCTCTTTGTTCTCATTCAGAAAATCtcatc < 1:1754082/71‑1 (MQ=255) gATTATTAGTCTTTGCATGAGCAGGCGAACTGATACATTTCTCTTTGTTCTCATTCAGAAAATCTCATCAg > 1:1792504/1‑71 (MQ=255) attattAGTCTTTGCATGAGCAAGCGAACTGATACATTTCTCTTTGTTCTCATTCAGAAAATCtcatca < 1:1915435/69‑1 (MQ=255) attattAGTCTTTGCATGAGCAAGCGAACTGATACATTTCTCTTTGTTCTCATTCAGAAAATCTCATCAg > 1:875258/1‑70 (MQ=255) attattAGTCTTTGCATGAGCAAGCGAACTGATACATTTCTCTTTGTTCTCATTCAGAAAATCTCATCAGt > 1:2202377/1‑71 (MQ=255) | TGCCAAAGAAAAATGTTCAGAAAATAATCCATGGAAAATTGTCCGGAGCACTTACTATTTTAATGGATTATTAGTCTTTGCATGAGCAAGCGAACTGATACATTTCTCTTTGTTCTCATTCAGAAAATCTCATCAGT > Exported/2823802‑2823938 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |