| Predicted mutation | ||||||
|---|---|---|---|---|---|---|
| evidence | seq id | position | mutation | annotation | gene | description |
| RA | REL606 | 1,886,177 | A→C | intergenic (‑7/+662) | yobF ← / ← yebO | hypothetical protein/hypothetical protein |
| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | REL606 | 1,886,177 | 0 | A | C | 97.4% | 151.8 / ‑5.4 | 39 | intergenic (‑7/+662) | yobF/yebO | hypothetical protein/hypothetical protein |
| Reads supporting (aligned to +/- strand): ref base A (1/0); new base C (22/16); total (23/16) | |||||||||||
| 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.88e-01 | |||||||||||
CGGCAGAGAAGCGGTATTCAACGTCAACGTGTTTACTCAGGACTTCTTTACTGAAAATGCCACACATAAACAGAACTGTACCTCGTTTAACCCGAAATCTGTTATCACATACAACGTTAATTATGGCAAGCCATTTTTAAACA > REL606/1886104‑1886246 | cGGCAGAGAAGCGGTATTCAACGTCAACGTGTTTACTCAGGACTTCTTTACTGAAAATGCCACACATAAACAGCAc < 2:843187/76‑1 (MQ=255) agagAAGCGGTATTCAACGTCAACGTGTTTACTCAGGACTTCTTTACTGAAAATGCCACACATAAACAGCACTGTa > 2:984351/1‑76 (MQ=255) agagAAGCGGTATTCAACGTCAACGTGTTTACTCAGGACTTCTTTACTGAAAATGCCACACATAAACAGCACTGTa > 1:559480/1‑76 (MQ=255) aaGCGGTATTCAACGTCAACGTGTTTACTCAGGACTTCTTTACTGAAAATGCCACACATAAACAGCACTGTACCTc > 2:681955/1‑76 (MQ=255) aaGCGGTATTCAACGTCAACGTGTTTACTCAGGACTTCTTTACTGAAAATGCCACACATAAACAGCACTGTACCTc > 2:681931/1‑76 (MQ=255) aGCGGTATTCAACGTCAACGTGGTTACTCAGGACTTCTTTACTGAAAATGCCACACATAAACAGCACTGTACCTCg < 2:40217/76‑1 (MQ=255) tATTCAACGTCAACGTGTTTACTCAGGACTTCTTTACTGAAAATGCCACACATAAACAGCACTGTACCTCGTTTaa > 1:1174348/1‑76 (MQ=255) aacgtcaacgtGTTTACTCAGGACTTCTTTACTGAAAATGCCACACATAAACAGCACTGTACCTCGTTTAACCCGa > 2:349469/1‑76 (MQ=255) tcaacgtGTTTACTCAGGACTTCTTTACTGAAAATGCCACACATAAACAGCACTGTACCTCGTTTAGCACGCAATc > 2:1307675/1‑76 (MQ=255) tcaacgtGTTTACTCAGGACTTCTTTACTGAAAATGCCACACATAAACAGCACTGTACCTCGTTTAACCCGAAATc < 1:51207/76‑1 (MQ=255) caacgtGTTTACTCAGGACTTCTTTACTGAAAATGCCACACATAAACAGCACTGTACCTCGTTTAACCCGAAATCt < 2:1161853/76‑1 (MQ=255) caacgtGTTTACTCAGGACTTCTTTACTGAAAATGCCACACATAAACAGCACTGTACCTCGTTTAACCCGAAATCt < 2:41252/76‑1 (MQ=255) caacgtGTTTACTCAGGACTTCTTTACTGAAAATGCCACACATAAACAGCACTGTACCTCGTTTAACCCGAAATCt > 2:327174/1‑76 (MQ=255) caacgtGTTTACTCAGGACTTCTTTACTGAAAATGCCACACATAAACAGCACTGTACCTCGTTTAACCCGAAATCt < 1:231129/76‑1 (MQ=255) gtgtTTACTCAGGACTTCTTTACTGAAGATGCCACACATAAACAGCACTGTACCTCGTTTAACCCGAAAGCTGTTa > 1:1297722/1‑76 (MQ=255) gtgtTTACTCAGGACTTCTTTACTGAAAATGCCACACATAAACAGCACTGTACCTCGTTTAACCCGAAATCTGTTa > 2:8570/1‑76 (MQ=255) tgtTTACTCAGGACTTCTTTACTGAAAATGCCACACATAAACAGCACTGTACCTCGTTTAACCCGAAATCTGTTAt < 1:960123/76‑1 (MQ=255) tgtTTACTCAGGACNTCTTTACTGAAAATGCCACACATAAACAGCACTGTACCTCGTTTAACCCGAAATCTGTTan < 2:178738/76‑2 (MQ=255) gtTTACTCAGGACTTCTTTACTGAAAATGCCACACATAAACAGAACTGTACCTCGTTTAACCCGAAATCTGTTATc > 2:480855/1‑76 (MQ=255) ttACTCAGGACTTCTTTACTGTAACTGCCACACATAAACAGCACTGTACCTCGTTTAACCCGAAATCTGTTATcac < 1:1241180/76‑1 (MQ=255) aCTCAGGACTTCTTTACTGAAAATGCCACACATAAACAGCACTGTACCTCGTTTAACCCGAAATCTGTTATCACAt > 2:330456/1‑76 (MQ=255) cAGGACTTCTTTACTGAAAATGCCACACATAAACAGCACTGTACCTCGTTTAACCCGAAATCTGTTATCACATACa > 2:830845/1‑76 (MQ=255) aGGACTTCTTTACTGAAAATGCCACACATAAACAGCACTGTACCTCGTTTAACCCGAAATCTGTTATCACATACaa > 2:806110/1‑76 (MQ=255) aGGACTTCTTTACTGAAAATGCCACACATAAACAGCACTGTACCTCGTTTAACCCGAAATCTGTTATCACATACaa < 1:1281528/76‑1 (MQ=255) tttACTGAAAATGCCACACATAAACAGCACTGTACCTCGTTTAACCCGAAATCTGTTATCACATACAACGTTAAtt > 2:741776/1‑76 (MQ=255) tACTGAAAATGCCACACATAAACAGCACTGTACCTCGTTTAACCCGAAATCTGTTATCACATACAACGTTAATTAt < 1:457338/76‑1 (MQ=255) gAAAATGCCACACATAAACAGCACTGTACCTCGTTTAACCCGACNTCNGTTATCACATACAANGTTAATTATGGCa > 2:1190258/1‑76 (MQ=255) gAAAATGCCACACATAAACAGCACTGTACCTCGTTTAACCCGAAATCTGTTATCACATACAACGTTAATTATGGCa < 1:1045861/76‑1 (MQ=255) gAAAATGCCACACATAAACAGCACTGTACCTCGTTTAACCCGAAATCTGTTATCACATACAACGTTAATTATGGCa < 2:468388/76‑1 (MQ=255) aaaaTGCCACACATAAACAGCACTGTACCTCGTTTAACCCGAAATCTGTTATCACATACAACGTTAATTATGGCaa > 1:551471/1‑76 (MQ=255) aaaTGCCACACATAAACAGCACTGTACCTCGTTTAACCCGAAATCTGTTATCACATACAACGTTAATTATGGCAAg > 1:685863/1‑76 (MQ=255) aaTGCCACACATAAACAGCACTGTACCTCGTTTAACCCGAAATCTGTTATCACATACAACGTTAATTATGGCAAGc > 2:367696/1‑76 (MQ=255) aTGCCACACATAAACAGCACTGTACCTCGTTTAACCCGAAATCTGTTATCACATACAACGTTAATTATGGCAAGcc < 1:327174/76‑1 (MQ=255) ccACACATAAACAGCACTGTACCTCGTTTAACCCGAAATCTGTTATCACATACAACGTTAATTATGGCAAGCCAtt > 2:943032/1‑76 (MQ=255) acaTAAACAGCACTGTACCTCGTTTAACCCGAAATCTGTTATCACATACAACGTTAATTATGGCAAGCCATTTTTa < 1:147543/76‑1 (MQ=255) acaTAAACAGCACTGTACCTCGTTTAACCCGAAATCTGTTATCACATACAACGTTAATTATGCCAAGCCATTTTTa < 1:84441/76‑1 (MQ=255) tAAACAGCACTGTACCTAGTTTAACCCGAAATCTGTTATCACATACAACGTTAATTATGGCAAGCCATTTTTAAAc < 1:147811/76‑1 (MQ=255) aaaCAGCACTGTACCTCGTTTAACCCGCAATCTGTTATCACATCNAANGTTAATTATGGCAANCCATTTTTAAACa > 2:1184737/1‑76 (MQ=255) aaaCAGCACTGTACCTCGTTTAACCCGAAATCTGTTATCACATACAACGTTAATTATGGCAAGCCATTTTTAAACa > 2:606590/1‑76 (MQ=255) aaaCAGCACTGTACCTCGTTTAACCCGAAATCTGTTATCACATACAACGTTAATTATGGCAAGCCATTTTTAAACa > 2:1031209/1‑76 (MQ=255) | CGGCAGAGAAGCGGTATTCAACGTCAACGTGTTTACTCAGGACTTCTTTACTGAAAATGCCACACATAAACAGAACTGTACCTCGTTTAACCCGAAATCTGTTATCACATACAACGTTAATTATGGCAAGCCATTTTTAAACA > REL606/1886104‑1886246 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 24 ≤ ATCG/ATCG < 30 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |