Predicted mutation | ||||||
---|---|---|---|---|---|---|
evidence | seq id | position | mutation | annotation | gene | description |
RA | REL606 | 2,749,029 | G→A | V292M (GTG→ATG) | hypE → | carbamoyl phosphate phosphatase, hydrogenase 3 maturation protein |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | REL606 | 2,749,029 | 0 | G | A | 97.9% | 173.9 / ‑3.1 | 48 | V292M (GTG→ATG) | hypE | carbamoyl phosphate phosphatase, hydrogenase 3 maturation protein |
Reads supporting (aligned to +/- strand): ref base G (1/0); new base A (26/21); total (27/21) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 1.00e+00 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.25e-01 |
CAACGCGGCAGAGCAAGTGCTGGCAGCGTTACATTCCCATCCACTGGGGAAAGACGCGGCGCTGATTGGTGAAGTGGTGGAACGTAAAGGTGTTCGTCTTGCCGGTCTGTATGGCGTGAAACGAACCCTCGATTTACCACACGCC > REL606/2748956‑2749100 | caacgcGGCAGAGCAAGTGCTGGCAGCGTTACATTCCCATCCACTGGGGAAAGACGCGGCGCTGATTGGTGAAatg > 1:1491915/1‑76 (MQ=255) aacgcGGCAGAGCAAGTGCTGGCAGCGTTACATTCCCATCCACTGGGGAAAGACGCGGCGCTGATTGGTGAAatgg > 1:1278274/1‑76 (MQ=255) aacgcGGCAGAGCAAGTGCTGGCAGCGTTACATTCCCATCCACTGGGGAAAGACGCGGCGCTGATTGGTGAAatgg > 2:224956/1‑76 (MQ=255) aacgcGGCAGAGCAAGTGCTGGCAGCGTTACATTCCCATCCACTGGGGAAAGACGCGGCGCTGATTGGTGAAatgg > 2:1424833/1‑76 (MQ=255) gcgGCAGAGCAAGTGCTGGCAGCGTTACATTCCCATCCACTGGGGAAAGACGCGGCGCTGATTGGTGAAAtggtgg > 2:232556/1‑76 (MQ=255) cgGCAGAGCAAGTGCTGGCAGCGTTACATTCCCATCCACTGGGGAAAGACGCGGCGCTGATTGGTGAAATGGTGGa > 1:1679015/1‑76 (MQ=255) ggCAGAGCAAGTGCTGGCAGCGTTACATTCCCATCCACTGGGGAAAGACGCGGCGCTGATTGGTGAAATGGTGGaa < 2:988730/76‑1 (MQ=255) ggCAGAGCAAGTGCTGGCAGCGTTACATTCCCATCCACTGGGGAAAGACGCGGCGCTGATTGGTGAAATGGTGGaa > 2:1713352/1‑76 (MQ=255) agagCAAGTGCTGGCAGCGTTACATTCCCATCCACTGGGGAAAGACGCGGCGCTGATTGGTGAAATGGTGGAACGt > 1:932958/1‑76 (MQ=255) gagCAAGTGCTGGCAGCGTTACATTCCCATCCACTGGGGAAAGACGCGGCGCTGATTGGTGAAATGGTGGAACGTa > 1:116170/1‑76 (MQ=255) gagCAAGTGCTGGCAGCGTTACATTCCCATCCACTGGGGAAAGACGCGGCGCTGATTGGTGAAATGGTGGAACGTa > 2:1136209/1‑76 (MQ=255) caGTGCTGGCAGCGTTACATTCCCATCCACTGGGGAAAGACGCGGCGCTGATTGGTGAAATGGTGGAACGTAAAgg < 1:210592/75‑1 (MQ=255) aaGTGCTGGCAGCGTTACATTCCCATCCACTGGGGAAAGACGCGGCGCTGATTGGTGAAATGGTGGAACGTAAAgg > 1:1701017/1‑76 (MQ=255) aaGTGCTGGCAGCGTTACATTCCCATCCACTGGGGAAAGACGCGGCGCTGATTGGTGAAATGGTGGAACGTAAAgg > 2:1345377/1‑76 (MQ=255) aGTGCTGGCAGCGTTACATTCCCATCCACTGGGGAAAGACGCGGCGCTGATTGGTGAAATGGTGGAACGTAAAGgt < 2:176261/76‑1 (MQ=255) tGCTGGCAGCGTTACATTCCCATCCACTGGGGAAAGACGCGGCGCTGATTGGTGAAATGGTGGAACGTAAAGgtgt < 1:1401654/76‑1 (MQ=255) ggCAGCGTTACATTCCCATCCACTGGGGAAAGACGCGGCGCTGATTGGTGAAATGGTGGAACGTAAAGGTGTTcac < 2:1396502/76‑3 (MQ=255) cAGCGTTACATTCCCATCCACTGGGGAAAGACGCGGCGCTGATTGGTGAAATGGTGGAACGTAAAGGTGTTCGTCt < 2:470005/76‑1 (MQ=255) cAGCGTTACATTCCCATCCACTGGGGAAAGACGCGGCGCTGATTGGTGAAATGGTGGAACGTAAAGGTGTTCGTCt > 1:1367315/1‑76 (MQ=255) ggCGTTACATTCCCATCCACTGGGGAAAGACGCGGCGCTGATTGGTGAAATGGTGGAACGTAAAGGTGTTCGTCtt < 2:856945/75‑1 (MQ=255) cGTTACATTCCCATCCACTGGGGAAAGACGCGGCGCTGATTGGTGAAATGGTGGAACGTAAAGGTGTTCGTCTTGc < 2:1545913/76‑1 (MQ=255) cGTTACATTCCCATCCACTGGGGAAAGACGCGGCGCTGATTGGTGAAATGGTGGAACGTAAAGGTGTTCGTCTTGc > 1:545120/1‑76 (MQ=255) cGTTACATTCCCATCCACTGGGGAAAGACGCGGCGCTGATTGGTGAAATGGTGGAACGTAAAGGTGTTCGTCTTGc < 2:1179158/76‑1 (MQ=255) gTTACATTCCCATCCACTGGGGAAAGACGCGGCGCTGATTGGTGAAATGGTGGAACGTAAAGGTGTTCGTCTTGcc > 1:1346567/1‑76 (MQ=255) gTTACATTCCCATCCACTGGGGAAAGACGCGGCGCTGATTGGTGAAATGGTGGAACGTAAAGGTGTTCGTCTTGcc > 1:1451582/1‑76 (MQ=255) ttACATTCCCATCCACTGGGGAAAGACGCGGCGCTGATTGGTGAAATGGTGGAACGTAAAGGTGTTCGTCTTGCCg < 1:1713352/76‑1 (MQ=255) ttACATTCCCATCCACTGGGGAAAGACGCGGCGCTGATTGGTGAAATGGTGGAACGTAAAGGTGTTCGTCTTGCCg > 2:136215/1‑76 (MQ=255) aCATTCCCATCCACTGGGGAAAGACGCGGCGCTGATTGGTGAAATGGTGGAACGTAAAGGTGTTCGTCTTGCCGGt > 2:1349832/1‑76 (MQ=255) aTTCCCATCCACTGGGGAAAGACGCGGCGCTGATTGGTGAAATGGTGGAACGTAAAGGTGTTCGTCTTGCCGGTCt > 1:1826425/1‑76 (MQ=255) aTTCCCATCCACTGGGGAAAGACGCGGCGCTGATTGGTGAAATGGTGGAACGTAAAGGTGTTCGTCTTGCCGGTCt > 1:312612/1‑76 (MQ=255) tCCCATCCACTGGGGAAAGACGCGGCGCTGATTGGTGAAATGGTGGAACGTAAAGGTGTTCGTCTTGCCGGTCTGt > 1:1021840/1‑76 (MQ=255) cATCCACTGGGGAAAGACGCGGCGCTGATTGGTGAAATGGTGGAACGTAAAGGTGTTCGTCTTGCCGGTCTGTATg < 1:1371756/76‑1 (MQ=255) cATCCACTGGGGAAAGACGCGGCGCTGATTGGTGAAATGGTGGAACGTAAAGGTGTTCGTCTTGCCGGTCTGTATg < 2:1056926/76‑1 (MQ=255) gtCCACTGGGGAAAGACGCGGCGCTGATTGGTGAAATGGTGGAACGTAAAGGTGTTAGTCTTGCCGGCCGGTATgg > 1:50311/2‑76 (MQ=255) aTCCACTGGGGAAAGACGCGGCGCTGATTGGTGAAATGGTGGAACGTAAAGGTGTTCGTCTTGCCGGTCTGTATgg > 1:764827/1‑76 (MQ=255) aTCCACTGGGGAAAGACGCGGCGCTGATTGGTGAAATGGTGGAACGTAAAGGTGTTCGTCTTGCCGGTCTGTATgg < 2:1712366/76‑1 (MQ=255) ccACTGGGGAAAGACGCGGCGCTGATTGGTGAAATGGTGGAACGTAAAGGTGTTCGTCTTGCCGGTCTGTATGGCg < 2:951470/76‑1 (MQ=255) cTGGGGAAAGACGCGGCGCTGATTGGTGAAATGGTGGAACGTAAAGGTGTTCGTCTTGCCGGTCTGTATGGCGTgg > 1:973321/1‑75 (MQ=255) cTGGGGAAAGACGCGGCGCTGATTGGTGAAATGGTGGAACGTAAAGGTGTTCGTCTTGCCGGTCTGTATGGCGTGa > 2:1452155/1‑76 (MQ=255) tGGGGAAAGACGCGGCGCTGATTGGTGAAATGGTGGAACGTAAAGGTGTTCGTCTTGCCGGTCTGTATGGCGTGaa > 1:1732215/1‑76 (MQ=255) aaaGACGCGGCGCTGATTGGTGAAATGGTGGAACGTAAAGGTGTTCGTCTTGCCGGTCTGTATGGCGTGAAACGaa > 2:257765/1‑76 (MQ=255) aaGACGCGGCGCTGATTGGGGAAAGGGTGGAACGTAAAGGTGTTCGTCTTGCTGGTCTGTATGGCGTGAAACGAAc > 2:918551/1‑76 (MQ=255) cgcgGCGCTGATTGGTGAAATGGTGGAACGTAAAGGTGTTCGTCTTGCCGGTCTGTATGGCGTGAAACGAACCCTc > 2:603302/1‑76 (MQ=255) ggcgCTGATTGGTGAAATGGTGGAACGTAAAGGTGTTCGTCTTGCCGGTCTGTATGGCGTGAAACGAACCCTCGAt < 1:379524/76‑1 (MQ=255) gcTGATTGGTGAAATGGTGGAACGTAAAGGTGTTCGTCTTGCCGGTCTGTATGGCGTGAAACGAACCCTCGATTTa < 2:87239/76‑1 (MQ=255) gcTGATTGGTGAAATGGTGGAACGTAAAGGTGTTCGTCTTGCCGGTCTGTATGGCGTGAAACGAACCCTCGATTTa < 2:1293057/76‑1 (MQ=255) tGATTGGTGAAATGGTGGAACGTAAAGGTGTTCGTCTTGCCGGTCTGTATGGCGTGAAACGAACCCTCGATTTAcc < 1:317984/76‑1 (MQ=255) aTTGGTGAAGTGGTGGAACGTAAAGGTGTTCGTCTTGCCGGTCTGTATGGCGTGAAACGAACCCTCGATTTACcac > 1:196307/1‑76 (MQ=255) aTTGGTGAAATGGTGGAACGTAAAGGTGTTCGTCTTGCCGGTCTGTATGGCGTGAAACGAACCCTCGATTTACcac < 1:1476606/76‑1 (MQ=255) aTTGGTGAAATGGTGGAACGTAAAGGTGTTCGTCTTGCCGGTCTGTATGGCGTGAAACGAACCCTCGATTTACcac < 1:1575413/76‑1 (MQ=255) ggTGAAATGGTGGAACGTAAAGGTGTTCGTCTTGCCGGTCTGTATGGCGTGAAACGAACCCTCGATTTACCACACg < 2:495705/76‑1 (MQ=255) tGAAATGGTGGAACGTAAAGGTGTTCGTCTTGCCGGTCTGTATGGCGTGAAACGAACCCTCGATTTACCACACGcc > 1:55441/1‑76 (MQ=255) | CAACGCGGCAGAGCAAGTGCTGGCAGCGTTACATTCCCATCCACTGGGGAAAGACGCGGCGCTGATTGGTGAAGTGGTGGAACGTAAAGGTGTTCGTCTTGCCGGTCTGTATGGCGTGAAACGAACCCTCGATTTACCACACGCC > REL606/2748956‑2749100 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 23 ≤ ATCG/ATCG < 28 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |