Predicted mutation | ||||||
---|---|---|---|---|---|---|
evidence | seq id | position | mutation | annotation | gene | description |
RA | REL606 | 2,740,354 | G→A | G326G (GGC→GGT) | hycE ← | hydrogenase 3, large subunit |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | REL606 | 2,740,354 | 0 | G | A | 98.2% | 207.6 / ‑4.2 | 57 | G326G (GGC→GGT) | hycE | hydrogenase 3, large subunit |
Reads supporting (aligned to +/- strand): ref base G (0/1); new base A (35/21); total (35/22) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 3.86e-01 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.00e-01 |
GTTGCGCCAGTTGACGGGTCTGGATCATGTCATCTTTCAGCAGATCGCGACGAATCCCGCCGATCAGGTTCAGGCCGTAGGTTTTACGCGCACCGGTAAGGATCTCTGCCATTTTCATGGAGGTTTCACGCACGCGGAAGAACTGC > REL606/2740281‑2740426 | gTTGCGCCAGTTGACGGGTCTGGATCATTTCATCTTTCAGCAGATCGCGACGAATCCCGCCGATCAGGTTCAGAcc < 2:800790/76‑1 (MQ=255) cgcCAGTTGACGGGTCTGGATCATGTCATCTTTCAGCAGATCGCGACGAATCCCGCCGATCAGGTTCAGACCGTAg > 1:1024494/1‑76 (MQ=255) cAGTTGACGGGTCTGGATCATGTCATCTTTCAGCAGATCGCGACGAATCCCGCCGATCAGGTTCAGACCGTAGGtt > 1:1514672/1‑76 (MQ=255) tggcGGGTCTGGATCTTGTCATCTTCCAGCAGACCGCGACGAATCCCGCCGATCAGGTTCAGACCGTAGGTTTTAc < 2:1029079/73‑1 (MQ=37) tGACGGGTCTGGATCATGTCATCTTTCAGCAGATCGCGACGAATCCCGCCGATCAGGTTCAGACCGTAGGTTTTAc > 1:1373067/1‑76 (MQ=255) gACGGGTCTGGATCATGTCATCTTTCAGCAGATCGCGACGAATCCCGCCGATCAGGTTCAGACCGTAGGTTTTAcg > 1:1219284/1‑76 (MQ=255) aCGGGTCTGGATCATGTCATCTTTCAGCAGATCGCGACGAATCCCGCCGATCAGGTTCAGACCGTAGGTTTTAcgc > 1:717448/1‑76 (MQ=255) cGGGTCTGGATCATGTCATCTTTCAGCAGATCGCGACGAATCCCGCCGATCAGGTTCAGACCGTAGGTTTTAcgcg > 1:1180987/1‑76 (MQ=255) gggTCTGGATCATGTCATCTTTCAGCAGATCGCGACGAATCCCGCCGATCAGGTTCAGACCGTAGGTTTTACgcgc > 1:1168078/1‑76 (MQ=255) gTCTGGATCATGTCATCTTTCAGCAGATCGCGACGAATCCCGCCGATCAGGTTCAGACCGTAGGTTTTACGCGCAc > 1:1150484/1‑76 (MQ=255) gTCTGGATCATGTCATCTTTCAGCAGATCGCGACGAATCCCGCCGATCAGGTTCAGACCGTAGGTTTTACGCGCAc > 1:1215315/1‑76 (MQ=255) gTCTGGATCATGTCATCTTTCAGCAGATCGCGACGAATCCCGCCGATCAGGTTCAGACCGTAGGTTTTACGCGCAc > 1:1517414/1‑76 (MQ=255) gTCTGGATCATGTCATCTTTCAGCAGACCGCGACAAACCCCGCCGATCAGGTTCAGACCGTAGGTTTTACGCGCAc < 1:143877/76‑1 (MQ=37) tGGATCATGTCATCTTTCAGCAGATCGCGACGAATCCCGCCGATCAGGTTCAGACCGTAGGTTTTACGCGCACCgg > 1:1852017/1‑76 (MQ=255) ggATCATGTCATCTTTCAGCAGATCGCGACGAATCCCGCCGATCAGGTTCAGACCGTAGGTTTTACGCGCACCGGt > 1:1353845/1‑76 (MQ=255) aTCATTTCATCTTTCAGCAGATCGCGACAAACCCCGCCGATCAGGTTCAGACCGTAGGTTTTACGCGCACCGGTaa < 1:1036373/76‑1 (MQ=25) aTCATGTCATCTTTCAGCAGACCGCGACGAATCCCGCCGATCAGGTTCAGACCGTAGGTTTTACGCGCACCGGTaa < 2:388984/76‑1 (MQ=38) cATGTCATCTTTCAGCAGATCGCGACGAATCCCGCCGATCAGGTTCAGACCGTAGGTTTTACGCGCACCGGTAAgg > 2:1632731/1‑76 (MQ=255) tGTCATCTTTCAGCAGATCGCGACGAATCCCGCCGATCAGGTTCAGACCGTAGGTTTTACGCGCACCGGTAAGGAt < 1:687492/76‑1 (MQ=255) gTCATCTTTCAGCAGATCGCGACGAATCCCGCCGATCAGGTTCAGACCGTAGGTTTTACGCGCACCGGTAAGGAtc > 1:1373715/1‑76 (MQ=255) tCATCTTTCAGCAGATCGCGACGAATCCCGCCGATCAGGTTCAGACCGTAGGTTTTACGCGCACCGGTAAGGAtct > 2:4676/1‑76 (MQ=255) aTCTTTCAGCAGATCGCGACGAATCCCGCCGATCAGGTTCAGACCGTAGGTTTTACGCGCACCGGTAAGGATctct > 1:1084330/1‑76 (MQ=255) aTCTTTCAGCAGATCGCGACGAATCCCGCCGATCAGGTTCAGACCGTAGGTTTTACGCGCACCGGTAAGGATctct > 1:255939/1‑76 (MQ=255) aTCTTTCAGCAGATCGCGACGAATCCCGCCGATCAGGTTCAGACCGTAGGTTTTACGCGCACCGGTAAGGATctct > 1:329539/1‑76 (MQ=255) tCTTTCAGCAGAGCGCGACGAATCCCGCCGATCAGGTTCAGACCGTAGGTTTTACGCGCACCGGTAAGGATCTCTg < 2:78211/76‑1 (MQ=255) cTTTCAGCAGATCGCGACGAATCCCGCCGATCAGGTTCAGACCGTAGGTTTTACGCGCACCGGTAAGGATCTCTGc > 2:530282/1‑76 (MQ=255) tttCAGCAGATCGCGACGAATCCCGCCGATCAGGTTCAGACCGTAGGTTTTACGCGCACCGGTAAGGATCTCTGcc > 2:1479751/1‑76 (MQ=255) tttCAGCAGATCGCGACGAATCCCGCCGATCAGGTTCAGACCGTAGGTTTTACGCGCACCGGTAAGGATCTCTGcc > 1:613074/1‑76 (MQ=255) tccagcagACCGCGACGAACCCCGCCGATCAGGTTCAGACCGTAGGTTTTACGCGCACCGGTAAGGATCTCTGCCa < 1:1245341/74‑1 (MQ=37) agcagATCGCGACGAATCCCGCCGATCAGGTTCAGACCGTAGGTTTTACGCGCACCGGTAAGGATCTCTGCCAttt < 2:1169844/76‑1 (MQ=255) agcagATCGCGACGAATCCCGCCGACCAGGTTCAGACCGTAGGTTTTACGCGCACCGGTAAGGATCTCTGCCAttt > 2:862778/1‑76 (MQ=255) gATCGCGACGAATCCCGCCGATCAGGTTCAGACCGTAGGTTTTACGCGCACCGGTAAGGATCTCTGCCATTTTCAt > 2:1219819/1‑76 (MQ=255) aTCGCGACGAATCCCGCCGATCAGGTTCAGACCGTAGGTTTTACGCGCACCGGTAAGGATCTCTGCCATTTTCATg > 1:1604403/1‑76 (MQ=255) tCGCGACGAATCCCGCCGATCAGGTTCAGACCGTAGGTTTTACGCGCACCGGTAAGGATCTCTGCCATTTTCATgg > 1:1824277/1‑76 (MQ=255) gacgaATCCCGCCGATCAGGTTCAGACCGTAGGTTTTACGCGCACCGGTAAGGATCTCTGCCATTTTCATGGAGGt < 1:296578/76‑1 (MQ=255) ccgaACCCCGCCGATCAGGTTCAGACCGTAGGTTTTACGCGCACCGGTAAGGATCTCTGCCATTTTCATGGAGGtt < 1:1459331/75‑1 (MQ=255) gaATCCCGCCGATCAGGTTCAGACCGTAGGTTTTACGCGCACCGGTAAGGATCTCTGCCATTTTCATGGAGGTTTc > 2:283393/1‑76 (MQ=255) gaATCCCGCCGATCAGGTTCAGACCGTAGGTTTTACGCGCACCGGTAAGGATCTCTGCCATTTTCATGGAGGTTTc > 1:595918/1‑76 (MQ=255) aaTCCCGCCGATCAGGTTCAGACCGTAGGTTTTACGCGCACCGGTAAGGATCTCTGCCATTTTCATGGAGGTTTca < 1:809834/76‑1 (MQ=255) aaTCCCGCCGATCAGGTTCAGACCGTAGGTTTTACGCGCACCGGTAAGGATCTCTGCCATTTTCATGGAGGTTTca > 1:1813570/1‑76 (MQ=255) aaTCCCGCCGATCAGGTTCAGACCGTAGGTTTTACGCGCACCGGTAAGGATCTCTGCCATTTTCATGGAGGTTTca < 2:1367437/76‑1 (MQ=255) aTCCCGCCGATCAGGTTCAGACCGTAGGTTTTACGCGCACCGGTAAGGATCTCTGCCATTTTCATGGAGGTTTcac > 2:1863784/1‑76 (MQ=255) cccGCCGATCAGGTTCAGACCGTAGTTTTTACGCGCACCGGTAAGGATCTCTGCCATTTTCATGGAGGTTTcacgc > 1:1775420/1‑76 (MQ=255) cccGCCGATCAGGTTCAGACCGTAGGTTTTACGCGCACCGGTAAGGATCTCTGCCATTTTCATGGAGGTTTcacgc < 1:284990/76‑1 (MQ=255) cccGCCGATCAGGTTCAGACAGTAGGTTTTACGCGCACCGGTAAGGATCTCTGCCATTTTCATGGAGGTTTcacgc > 2:1523530/1‑76 (MQ=255) cgccgATCAGGTTCAGACCGTAGGTTTTATGCGCACCGGTAAGGATCTCTGCCATTTTCATGGAGGTTTcacgcac < 1:1859415/76‑1 (MQ=255) cgccgATCAGGTTCAGACCGTAGGTTTTACGCGCACCGGTAAGGATCTCTGCCATTTTCATGGAGGTTTcacgcac < 1:1341269/76‑1 (MQ=255) cgATCAGGTTCAGGCCGTAGGTTTTACGCGCACCGGTAAGGATCTCTGCCATTTTCATGGAGGTTTCACGCAcgcg < 1:985968/76‑1 (MQ=255) cgATCAGGTTCAGACCGTAGGTTTTACGCGCACCGGTAAGGATCTCTGCCATTTTCATGGAGGTTTCACGCAcgcg > 1:1325189/1‑76 (MQ=255) aTCAGGTTCAGACCGTAGGTTTTACGCGCACCGGTAAGGATCTCTGCCATTTTCATGGAGGTTTCACGCACGCGga > 2:578347/1‑76 (MQ=255) aTCAGGTTCAGACCGTAGGTTTTACGCGCACCGGTAAGGATCTCTGCCATTTTCATGGAGGTTTCACGCACGCGga < 2:903647/76‑1 (MQ=255) aTCAGGTTCAGACCGTAGGTTTTACGCGCACCGGTAAGGATCTCTGCCATTTTCATGGAGGTTTCACGCACGCGga < 1:1648835/76‑1 (MQ=255) aGGTTCAGACCGTAGGTTTTACGCGCACCGGTAAGGATCTCTGCCATTTTCATGGAGGTTTCACGCACGCGgaaga < 2:262793/76‑1 (MQ=255) aGGTTCAGACCGTAGGTTTTACGCGCACCGGTAAGGATCTCTGCCATTTTCATGGAGGTTTCACGCACGCGgaaga > 2:1815502/1‑76 (MQ=255) tgTTCAGACCGTAGGTTTTACGCGCACCGGTAAGGATCTCTGCCATTTTCATGGAGGTTTCACGCACGCGgaagaa < 2:509803/75‑1 (MQ=255) ggTTCAGCCCGTCGGTTTTACGCGCCCCGGTAAGGATCTCTGCCATTTTCATGGAGGTTTCACGCACGCGgaagaa < 1:944333/76‑1 (MQ=255) ggTTCAGACCGTAGGTTTTGCGCGCACCGGTAAGGATCTCTGCCATTTTCATGGAGGTTTCACGCACGCggaggac > 1:593141/1‑75 (MQ=38) cAGACCGTAGGTTTTACGCGCACCGGTAAGGATCTCTGCCATTTTCATGGAGGTTTCACGCACGCGGAAGAACTGc < 1:505527/76‑1 (MQ=255) | GTTGCGCCAGTTGACGGGTCTGGATCATGTCATCTTTCAGCAGATCGCGACGAATCCCGCCGATCAGGTTCAGGCCGTAGGTTTTACGCGCACCGGTAAGGATCTCTGCCATTTTCATGGAGGTTTCACGCACGCGGAAGAACTGC > REL606/2740281‑2740426 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 23 ≤ ATCG/ATCG < 30 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |