Predicted mutation | ||||||
---|---|---|---|---|---|---|
evidence | seq id | position | mutation | annotation | gene | description |
RA | REL606 | 4,273,299 | A→T | intergenic (+227/‑115) | nrfG → / → gltP | heme lyase (NrfEFG) for insertion of heme into c552, subunit NrfG/glutamate/aspartate:proton symporter |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | REL606 | 4,273,299 | 0 | A | T | 96.3% | 213.8 / ‑1.1 | 54 | intergenic (+227/‑115) | nrfG/gltP | heme lyase (NrfEFG) for insertion of heme into c552, subunit NrfG/glutamate/aspartate:proton symporter |
Reads supporting (aligned to +/- strand): ref base A (0/2); new base T (18/34); total (18/36) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 5.47e-01 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-01 |
GTTATTATGCATTAATGCAGCGAAAAGCTCTGTTGTTAAAGGGTTGCGCAACATACCGCGCAAATGATACTGATCATAAGCGTTAAAAAAATCTACAAACCAACGCAACACAATTCATGCCCTGGCAGTATGTCACGTTCTCGC > REL606/4273227‑4273370 | gTTATTATGCATTAATGCAGCGAAAAGCTCTGTTGTTAAAGGGTTGCGCAACATACCGCGCAAATGATACTGTTCa < 1:1436231/76‑1 (MQ=255) gTTATTATGCATTAATGCAGCGAAAAGCTCTGTTGTTAAAGGGTTGCGCAACATACCGCGCAAATGATACTGTTCa < 1:1385366/76‑1 (MQ=255) tattatGCATTAATGCAGCGAAAAGCTCTGTTGTTAAAGGGTTGCGCAACATACCGCGCAAATGATACTGTTCATa < 2:676898/76‑1 (MQ=255) attatGCATTAATGCAGCGAAAAGCTCTGTTGTTAAAGGGTTGCGCAACATACCGCGCAAATGATACTGTTCATaa < 1:1260007/76‑1 (MQ=255) tatGCATTAATGCAGCGAAAAGCTCTGTTGTTAAAGGGTTGCGCAACATACCGCGCAAATGATACTGTTCATAAGc > 1:228051/1‑76 (MQ=255) tatGCATTAATGCAGCGAAAAGCTCTGTTGTTAAAGGGTTGCGCAACATACCGCGCAAATGATACTGTTCATAAGc > 1:494692/1‑76 (MQ=255) tGCATTAATGCAGCGAAAAGCTCTGTTGTTAAAGGGTTGCGCAACATACCGCGCAAATGATACTGTTCATAAGCGt < 1:341313/76‑1 (MQ=255) tGCATTAATGCAGCGAAAAGCTCTGTTGTTAAAGGGTTGCGCAACATACCGCGCAAATGATACTGTTCATAAGCGt < 1:760388/76‑1 (MQ=255) tGCATTAATGCAGCGAAAAGCTCTGTTGTTAAAGGGTTGCGCAACATACCGCGCAAATGATACTGTTCATAAGCGt < 1:1301484/76‑1 (MQ=255) gCATTAATGCAGCGAAAAGCTCTGTTGTTAAAGGGTTGCGCAACATACCGCGCAAATGATACTGTTCATAAGCGtt < 2:1069699/76‑1 (MQ=255) gCATTAATGCAGCGAAAAGCTCTGTTGTTAAAGGGTTGCGCAACATACCGCGCAAATGATACTGTTCATAAGCGtt < 2:619962/76‑1 (MQ=255) aTTAATGCAGCGAAAAGCTCTGTTGTTAAAGGGTTGCGCAACATACCGCGCAAATGATACTGTTCATAAGCGTTaa > 1:206338/1‑76 (MQ=255) aTTAATGCAGCGAAAAGCTCTGTTGTTAAAGGGTTGCGCAACATACCGCGCAAATGATACTGTTCATAAGCGTTaa < 1:948761/76‑1 (MQ=255) tAATGCAGCGAAAAGCTCTGTTGTTAAAGGGTTGCGCAACATACCGCGCAAATGATACTGTTCATAAGCGTTaaaa < 2:154609/76‑1 (MQ=255) tAATGCAGCGAAAAGCTCTGTTGTTAAAGGGTTGCGCAACATACCGCGCAAATGATACTGTTCATAAGCGTTaaaa < 2:39103/76‑1 (MQ=255) aTGCAGCGAAAAGCTCTGTTGTTAAAGGGTTGCGCAACATACCGCGCAAATGATACTGTTCATAAGCGTTaaaaan < 1:235457/76‑2 (MQ=255) tGCAGCGAAAAGCTCTGTTGTTAAAGGGTTGCGCAACATACCGCGCAAATGATACTGTTCATAAGCGTTaaaaaaa > 1:1288162/1‑76 (MQ=255) gCAGCGAAAAGCTCTGTTGTTAAAGGGTTGCGCAACATACCGCGCAAATGATACTGTTCATAAGCGTTAAAAAAAt < 2:823731/76‑1 (MQ=255) cGAAAAGCTCTGTTGTTAAAGGGTTGCGCAACATACCGCGCAAATGATACTGTTCATAAGCGTTAAAAAAATCTAc < 2:1294737/76‑1 (MQ=255) cGAAAAGCTCTGTTGTTAAAGGGTTGCGCAACATACCGCGCAAATGATACTGTTCATAAGCGTTAAAAAAATCTAc < 1:254643/76‑1 (MQ=255) cGAAAAGCTCTGTTGTTAAAGGGTTGCGCAACATACCGCGCAAATGATACTGTTCATAAGCGTTAAAAAAATCTAc < 1:39070/76‑1 (MQ=255) cGAAAAGCTCTGTTGTTAAAGGGTTGCGCAACATACCGCGCAAATGATACTGTTCATAAGCGTTAAAAAAATCTAc < 1:1124149/76‑1 (MQ=255) cGAAAAGCTCTGTTGTTAAAGGGTTGCGCAACATACCGCGCAAATGATACTGTTCATAAGCGTTAAAAAAATCTAc < 1:1124132/76‑1 (MQ=255) gAAAAGCTCTGTTGTTAAAGGGTTGCGCAACATACCGCGCAAATGATACTGTTCATAAGCGTTAAAAAAATCTACa < 1:1034404/76‑1 (MQ=255) aaaaGCTCTGTTGTTAAAGGGTTGCGCAACATACCGCGCAAATGATACTGTTCATAAGCGTTAAAAAAATCTACaa < 1:65974/76‑1 (MQ=255) aaaGCTCTGTTGTTAAAGGGTTGCGCAACATACCGCGCAAATGATACTGTTCATAAGCGTTAAAAAAATCTACaaa < 2:1021286/76‑1 (MQ=255) aGCTCTGTTGTTAAAGGGTTGCGCAACATACCGCGCAAATGATACTGTTCATAAGCGTTAAAAAAATCTACAAAcc < 1:1102898/76‑1 (MQ=255) aGCTCTGTTGTTAAAGGGTTGCGCAACATACCGCGCAAATGATACTGTTCATAAGCGTTAAAAAAATCTACAAAcc > 2:1069844/1‑76 (MQ=255) ctGTTGTTAAAGGGTTGCGTAACATACCGCGCAAATGATACTGTTCATAAGCGTTAAAAAAATCTACAAACCAACg > 2:620734/1‑76 (MQ=255) ctGTTGTTAAAGGGTTGCGCAACATACCGCGCAAATGATACTGTTCATAAGCGTTAAAAAAATCTACAAACCAACg > 2:274723/1‑76 (MQ=255) tgttgtTAAAGGGTTGCGCAACATACCGCGCAAATGATACTGTTCATAAGCGTTAAAAAAATCTACAAACCAACGc > 1:89654/1‑76 (MQ=255) tgttgtTAAAGGGTTGCGCAACATACCGCGCAAATGATACTGTTCATAAGCGTTAAAAAAATCTACAAACCAACGc < 1:1042770/76‑1 (MQ=255) gttgttAAAGGGTTGCGCAACATACCGCGCAAATGATACTGTTCATAAGCGTTAAAAAAATCTACAAACCAACGCa < 1:841691/76‑1 (MQ=255) gttAAAGGGTTGCGCAACATACCGCGCAAATGATACTGTTCATAAGCGTTAAAAAAATCTACAAACCAACGCAaca < 1:295795/76‑1 (MQ=255) ttAAAGGGTTGCGCAACATACCGCGCAAATGATACTGTTCATAAGCGTTAAAAAAATCTACAAACCAACGCAacac > 1:1471951/1‑76 (MQ=255) ttAAAGGGTTGCGCAACATACCGCGCAAATGATACTGTTCATAAGCGTTAAAAAAATCTACAAACCAACGCAacac > 2:1074124/1‑76 (MQ=255) ttAAAGGGTTGCGCAACATACCGCGCAAATGATACTGTTCATAAGCGTTAAAAAAATCTACAAACCAACGCAacac > 2:1146727/1‑76 (MQ=255) tAAAGGGTTGCGCAACATACCGCGCAAATGATACTGTTCATAAGCGTTAAAAAAATCTACAAACCAACGCAAcaca < 2:1236839/76‑1 (MQ=255) tAAAGGGTTGCGCAACATACCGCGCAAATGATACTGATCATAAGCGTTAAAAAAATCTACAAACCAACGCAAcaca < 1:499345/76‑1 (MQ=255) aaaGGGTTGCGCAACATACCGCGCAAATGATACTGTTCATAAGCGTTAAAAAAATCTACAAACCAACGCAACACaa > 2:1313041/1‑76 (MQ=255) aaGGGTTGCGCAACATACCGCGCAAATGATACTGTTCATAAGCGTTAAAAAAATCTACAAACCAACGCAACACAAt > 2:217767/1‑76 (MQ=255) aGGGTTGCGCAACATACCGCGCAAATGATACTGTTCATAAGCGTTAAAAAAATCTACAAACCAACGCAACACAAtt > 2:1243467/1‑76 (MQ=255) aGGGTTGCGCAACATACCGCGCAAATGATACTGTTCATAAGCGTTAAAAAAATCTACAAACCAACGCAACACAAtt < 2:266014/76‑1 (MQ=255) ggTTGCGCAACATACCGCGCAAATGATACTGTTCATAAGCGTTAAAAAAATCTACAAACCAACGCAACACAATTCa > 2:922825/1‑76 (MQ=255) gTTGCGCAACATACCGCGCAAATGATACTGATCATAAGCGTTAAAAAAATCTACAAACCAACGCAACACAATTCAt < 1:409702/76‑1 (MQ=255) aCATACCGCGCAAATGATACTGTTCATAAGCGTTAAAAAAATCTACAAACCAACGCAACACAATTCATGCCCTGGc < 2:1446422/76‑1 (MQ=255) aTACCGCGCAAATGATACTGTTCATAAGCGTTAAAAAAATCTACAAACCAACGCAACACAATTCATGCCCTGGCAg < 2:903668/76‑1 (MQ=255) tACCGCGCAAATGATACTGTTCATAAGCGTTAAAAAAATCTACAAACCAACGCAACACAATTCATGCCCTGGCAGt < 2:222785/76‑1 (MQ=255) aCCGCGCAAATGATACTGTTCATAAGCGTTAAAAAAATCTACAAACCAACGCAACACAATTCATGCCCTGGCAGTa > 1:362006/1‑76 (MQ=255) aCCGCGCAAATGATACTGTTCATAAGCGTTAAAAAAATCTACAAACCAACGCAACACAATTCATGCCCTGGCAGTa < 1:1157853/76‑1 (MQ=255) cAAATGATACTGTTCATAAGCGTTAAAAAAATCTACAAACCAACGCAACACAATTCATGCCCTGGCAGTATGTCAc > 1:470714/1‑76 (MQ=255) aaaTGATACTGTTCATAAGCGTTAAAAAAATCTACAAACCAACGCAACACAATTCATGCCCTGGCAGTATGTCACg < 1:1285968/76‑1 (MQ=255) aTACTGTTCATAAGCGTTAAAAAAATCTACAAACCAACGCAACACAATTCATGCCCTGGCAGTATGTCACGTtctc < 2:1066368/76‑1 (MQ=255) aCTGTTCATAAGCGTTAAAAAAATCTACAAACCAACGCAACACAATTCATGCCCTGGCAGTATGTCACGTTCTcgc > 1:964224/1‑76 (MQ=255) | GTTATTATGCATTAATGCAGCGAAAAGCTCTGTTGTTAAAGGGTTGCGCAACATACCGCGCAAATGATACTGATCATAAGCGTTAAAAAAATCTACAAACCAACGCAACACAATTCATGCCCTGGCAGTATGTCACGTTCTCGC > REL606/4273227‑4273370 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 26 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |