| Predicted mutation | ||||||
|---|---|---|---|---|---|---|
| evidence | seq id | position | mutation | annotation | gene | description |
| RA | REL606 | 3,928,260 | T→A | I15N (ATC→AAC) | rho → | transcription termination factor Rho |
| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | REL606 | 3,928,260 | 0 | T | A | 93.9% | 182.1 / 0.6 | 49 | I15N (ATC→AAC) | rho | transcription termination factor Rho |
| Reads supporting (aligned to +/- strand): ref base T (0/3); new base A (24/22); total (24/25) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 2.35e-01 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 | |||||||||||
TTTACCCCAAGTTTAAGAACTCACACCACTATGAATCTTACCGAATTAAAGAATACGCCGGTTTCTGAGCTGATCACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTATGCGTAAGCAGGACATTATTTTTGCCATCCT > REL606/3928187‑3928332 | tttACCCCAATTGTAAGAACGCGCACCGCTATGAAACTTACCGAAATAAATAAAACGCCGGATTCTGAGCTGATCa < 2:54710/76‑1 (MQ=255) ccccAAGTTTAAGAACTCACACCACTATGAATCTTACCGAATTAAAGAATACGCCGGTTTCTGAGCTGAACActct < 1:747937/76‑1 (MQ=255) cccAAGTTTAAGAACTCACACCACTATGAATCTTACCGAATTAAAGAATACGCCGGTTTCTGAGCTGAACACtctc > 2:861078/1‑76 (MQ=255) cccAAGTTTAAGAACTCACACCACTATGAATCTTACCGAATTAAAGAATACGCCGGTTTCTGAGCTGAACACtctc > 2:535937/1‑76 (MQ=255) ccAAGTTTAAGAACTCACACCACTATGAATCTTACCGAATTAAAGAATACGCCGGTTTCTGAGCTGAACACTCTCg < 1:239561/76‑1 (MQ=255) aGTTTAAGAACTCACACCACTATGAATCTTTCCGAATTAAAGAATACGCCGGTTTCTGAGCTGATCACTCTCGGCg < 2:605604/76‑1 (MQ=255) gTTTAAGAACTCACACCACTATGAATCTTACCGAATTAAAGAATACGCCGGTTTCTGAGCTGAACACTCTCGGCGa < 2:1443807/76‑1 (MQ=255) ttAAGAACTCACACCACTATGAATCTTACCGAATTAAAGAATACGCCGGTTTCTGAGCTGAACACTCTCGGCGaaa > 1:851119/1‑76 (MQ=255) aGAACTCACACCACTATGAATCTTACCGAATTAAAGAATACGCCGGTTTCTGAGCTGAACACTCTCGGCGAAAata < 1:367219/76‑1 (MQ=255) cTCACACCACTATGAATCTTACCGAATTAAAGAATACGCCGGTTTCTGAGCTGAACACTCTCGGCGAAAATATggg > 1:332634/1‑76 (MQ=255) tCACACCACTATGAATCTTACCGAATTAAAGAATACGCCGGTTTCTGAGCTGGACACTCTCGGCGAAAATATgggg > 2:744904/1‑76 (MQ=255) tCACACCACTATGAATCTTACCGAATTAAAGAATACGCCGGTTTCTGAGCTGAACACTCTCGGCGCAAATATgggg > 2:744899/1‑76 (MQ=255) acacCACTATGAATCTTACCGAATTAAAGAATACGCCGGTTTCTGAGCTGAACACTCTCGGCGAAAATATGGGGCt > 1:507998/1‑76 (MQ=255) caccacTATGAATCTTACCGAATTAAAGAATACGCCGGTTTCTGAGCTGAACACTCTCGGCGAAAATATGGGGCTg < 1:1302102/76‑1 (MQ=255) ccacTATGAATCTTACCGAATTAAAGAATACGCCGGTTTCTGAGCTGAACACTCTCGGCGAAAATATGGGGCTGGa < 2:566244/76‑1 (MQ=255) ccacTATGAATCTTACCGAATTAAAGAATACGCCGGTTTCTGAGCTGAACACTCTCGGCGAAAATATGGGGCTGGa > 1:592710/1‑76 (MQ=255) cacTATGAATCTTACCGAATTAAAGAATACGCCGGTTTCTGAGCTGAACACTCTCGGCGAAAATATGGGGCTGGaa > 2:438489/1‑76 (MQ=255) cacTATGAATCTTACCGAATTAAAGAATACGCCGGTTTCTGAGCTGAACACTCTCGGCGAAAATATGGGGCTGGaa < 2:923094/76‑1 (MQ=255) acTATGAATCTTACCGAATTAAAGAATACGCCGGTTTCTGAGCTGAACACTCTCGGCGAAACTATGGGGCTGGaaa > 2:1412753/1‑76 (MQ=255) cTATGAATCTTACCGAATTAAAGAATACGCCGGTTTCTGAGCTGAACACTCTCGGCGAAAATATGGGGCTGGaaaa < 2:1452231/76‑1 (MQ=255) tATGAATCTTACCGAATTAAAGAATACGCCGGTTTCTGAGCTGAACACTCTCGGCGAAAATATGGGGCTGGAAAAc < 1:983311/76‑1 (MQ=255) ttACCGAATTAAAGAATACGCCGGTTTCTGAGCTGAACACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCg > 1:245688/1‑76 (MQ=255) tACCGAATTAAAGAATACGCCGGTTTCTGAGCTGAACACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGt > 1:1078081/1‑76 (MQ=255) tACCGAATTAAAGAATACGCCGGTTTCTGAGCTGAACACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGt > 2:427384/1‑76 (MQ=255) ccGAATTAAAGAATACGCCGGTTTCTGAGCTGAACACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTAt < 1:869365/76‑1 (MQ=255) gAATTAAAGAATACGCCGGTTTCTGAGCTGAACACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTATGc > 2:75719/1‑76 (MQ=255) aaTTAAAGAATACGCCGGTTTCTGAGCTGAACACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTATGCg > 1:76794/1‑76 (MQ=255) tAAAGAATACGCCGGTTTCTGAGCTGAACACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTATGCGTaa > 1:447634/1‑76 (MQ=255) tAAAGAATACGCCGGTTTCTGAGCTGAACACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTATGCGTaa > 1:1038594/1‑76 (MQ=255) tAAAGAATACGCCGGTTTCTGAGCTGAACACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTATGCGTaa < 2:1164053/76‑1 (MQ=255) tAAAGAATACGCCGGTTTCTGAGCTGAACACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTATGCGTaa < 2:849259/76‑1 (MQ=255) aaaGGATACGCCGGTTTCTGAGCTGAACACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTATGCGTAAg < 1:226262/76‑1 (MQ=255) aaaGAATACGCCGGTTTCTGAGCTGAACACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTATGCGTAAg < 1:787113/76‑1 (MQ=255) aaaGAATACGCCGGTTTCTGAGCTGAACACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTATGCGTAAg > 1:1292173/1‑76 (MQ=255) aaaGAATACGCCGGTTTCTGAGCTGAACACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTATGCGTAAg < 2:1274840/76‑1 (MQ=255) aGAATACGCCGGTTTCTGAGCTGATCACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTATGCGTAAGCa < 1:786082/76‑1 (MQ=255) aGAATACGCCGGTTTCTGAGCTGAACACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTATGCGGAAGCa > 1:595704/1‑76 (MQ=255) aTACGCCGGTTTCTGAGCTGAACACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTATGCGTAAGCAGGa < 2:83324/76‑1 (MQ=255) gCCGGTTTCTGAGCTGAACACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTATGCGTAAGCAGGACatt < 2:1170976/76‑1 (MQ=255) ccGGTTTCTGAGCTGAACACTCTCGGCGAAAATATAGGGCTGGAAAACCTGGCTCGTATGCGTAAGCAGGACatta < 1:907188/76‑1 (MQ=255) cGGTTTCTGAGCTGAACACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTATGCGGAAGCAGGACattat < 1:801682/76‑1 (MQ=255) gTTTCTGAGCTGAACACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTATGCGTAAGCAGGACATTAttt > 1:1272718/1‑76 (MQ=255) tttCTGAGCTGAACACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTATGCGTAAGCAGGACATTAtttt > 1:1200102/1‑76 (MQ=255) tCTGAGCTGAACACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTATGCGTAAGCAGGACATTATTTTTg > 1:685889/1‑76 (MQ=255) tGAGCTGAACACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTATGCGTAAGCAGGACATTATTTTTGcc < 1:1094739/76‑1 (MQ=255) gAGCTGAACACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTATGCGTAAGCAGGACATTATTTTTGCCa < 1:1461515/76‑1 (MQ=255) aGCTGAACACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTATGCGTAAGCAGGACATTATTTTTGCCAt < 2:1382960/76‑1 (MQ=255) gCTGAACACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTATGCGTAAGCAGGACATTATTTTTGCCATc < 1:632266/76‑1 (MQ=255) tGATCACTCTCGGCGAAAATATGGGGCTGGAAAGCCTGGCTCGTATGCGTAAGCAGGACATTATTTTTGCCATCCt < 1:879015/76‑1 (MQ=255) tGAACACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTATGCGTAAGCAGGACATTATTTTTGCCATCCt > 1:1350786/1‑76 (MQ=255) tGAACACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTATGCGTAAGCAGGACATTATTTTTGCCATCCt > 1:631259/1‑76 (MQ=255) | TTTACCCCAAGTTTAAGAACTCACACCACTATGAATCTTACCGAATTAAAGAATACGCCGGTTTCTGAGCTGATCACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTATGCGTAAGCAGGACATTATTTTTGCCATCCT > REL606/3928187‑3928332 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 25 ≤ ATCG/ATCG < 31 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |