Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	REL606	4,203,168	+A	coding (410/3684 nt)	metH →	B12‑dependent methionine synthase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	REL606	4,203,166	1	.	A	100.0%	147.9 / NA	39	T136T (ACC→ACA)	metH	B12‑dependent methionine synthase
Reads supporting (aligned to +/- strand):  ref base . (0/0);  new base A (21/18);  total (21/18)

GCGCCTGTGCTGACGAGTGGACCGCGCGCACGCCGGAGAAACCGCGCTACGTTGCCGGTGTTCTCGGCCCAACC‑AACCGCACGGCGTCTATTTCTCCGGACGTCAACGATCCGGCATTTCGTAATATCACTTTTGACGGGCTGG  >  REL606/4203093‑4203236                                                                           |                                                                       gcgcCTGTGCTGACGAGTGGACCGCGCGCACGCCGGAGAAACCGCGCTACGTTGCCGGTGTTCTCGGCccaacc‑aa                                                                      >  1:1585987/1‑76 (MQ=255)   gcCTGTGCTGACGAGTGGACCGCGCGCACGCCGGAGAAACCGCGCTACGTTGCCGGTGTTCTCGGCccaacc‑aaac                                                                    <  2:1686735/76‑3 (MQ=255)    ccTGTGCTGACGAGTGGACCGCGCGCACGCCGGAGAAACCGCGCTACGTTGCCGGTGTTCTCGGCccaacc‑aaacc                                                                   >  2:1425134/1‑73 (MQ=255)    ccTGTGCTGACGAGTGGACCGCGCGCACGCCGGAGAAACCGCGCTACGTTGCCGGTGTTCTCGGCccaacc‑aaacc                                                                   <  1:348160/76‑4 (MQ=255)        ggCTGACGAGTGGACCGCGCGCACGCCGGAGAAACCGCGCTACGTTGCCGGTGTTCTCGGCCCAACCAAACCGCAc                                                                <  1:697522/75‑1 (MQ=255)          cTGACGAGTGGACCGCGCGCACGCCGGAGAAACCGCGCTACGTTGCCGGTGTTCTCGGCCCAACCAAACCGCACgg                                                              >  1:320909/1‑76 (MQ=255)          cTGACGAGTGGACCGCGCGCACGCCGGAGAAACCGCGCTACGTTGCCGGTGTTCTCGGCCCAACCAAACCGCACgg                                                              >  1:833890/1‑76 (MQ=255)            gACGAGTGGACCGCGCGCACGCCGGAGAAACCGCGCTACGTTGCCGGTGTTCTCGGCCCAACCAAACCGCACGGCg                                                            <  1:696068/76‑1 (MQ=255)               gAGTGGACCGCGCGCACGCCGGAGAAACCGCGCTACGTTGCCGGTGTTCTCGGCCCAACCAAACCGCACGGCGTCt                                                         <  1:1496504/76‑1 (MQ=255)                  tGGACCGCGCGCACGCCGGAGAAACCGCGCTACGTTGCCGGTGTTCTCGGCCCAACCAAACCGCACGGCGTCTAtt                                                      <  2:854498/76‑1 (MQ=255)                      ccGCGCGCACGGCGGAGATACCGCGCTACGTTGCCGGTGTTCTCGGCCCAACCAAACCGCACGGCGTCTATTtctc                                                  <  1:1010872/76‑1 (MQ=255)                      ccGCGCGCACGCCGGAGAAACCGCGCTACGTTGCCGGTGTTCTCGGCCCAACCAAACCGCACGGCGTCTATTtctc                                                  >  2:352209/1‑76 (MQ=255)                       cgcgcgCACGCCGGAGAAACCGCGCTACGTTGCCGGTGTTCTCGGCCCAACCAAACCGCACGGCGTCTATTTCTcc                                                 <  1:136743/76‑1 (MQ=255)                         cgcgcACGCCGGAGAAACCGCGCTACGTTGCCGGTGTTCTCGGCCCAACCAAACCCCACGGCGTCTATTTCTCCgg                                               >  2:1442813/1‑76 (MQ=255)                             cACGCCGGAGAAACCGCGCTACGTTGCCGGTGTTCTCGGCCCAACCAAACCGCACGGCGTCTATTTCTCCGGACGt                                           >  1:200921/1‑76 (MQ=255)                             cACGCCGGAGAAACCGCGCTACGTTGCCGGTGTTCTCGGCCCAACCAAACCGCACGGCGTCTATTTCTCCGGACGt                                           >  2:1688161/1‑76 (MQ=255)                             cACGCCGGAGAAACCGCGCTACGTTGCCGGTGTTCTCGGCCCAACCAAACCGCACGGCGTCTATTTCTCCGGACGt                                           >  2:1660386/1‑76 (MQ=255)                             cACGCCGGAGAAACCGCGCTACGTTGCCGGTGTTCTCGGCCCAACCAAACCGCACGGCGTCTATTTCTCCGGACGt                                           >  2:1402492/1‑76 (MQ=255)                              aCGCCGGAGAAACCGCGCTACGTTGCCGGTGTTCTCGGCGCAACCAAACCGCACGGCGTCTATTTCTCCGGACGTc                                          <  2:978908/76‑1 (MQ=255)                               cGCCGGAGACACCGCGCTACGTTGCCGGTGTTCTCGGCCCAACCAAACCGCACGGCGTCTATTTCTCCGGACGTCa                                         <  2:144091/76‑1 (MQ=255)                                gCCGGAGAAACCGCGCTACGTTGCCGGTGTTCTCGGCCCAACCAAACCGCACGGCGTCTATTTCTCCGGACGTCaa                                        <  1:1037306/76‑1 (MQ=255)                                  cGGAGAAACCGCGCTACGTTGCCGGTGTTCTCGGCCCAACCAAACCGCACGGCGTCTATTTCTCCGGACGTCAACg                                      >  1:464483/1‑76 (MQ=255)                                     agaAACCGCGCTACGTTGCCGGTGTTCTCGGCCCAACCAAACCGCACGGCGTCTATTTCTCCGGACGTCAACGATc                                   <  2:1092914/76‑1 (MQ=255)                                      gaAACCGCGCTACGTTGCCGGTGTTCTCGGCCCAACCAAACCGCACGGCGTCTATTTCTCCGGACGTCAACGATcc                                  <  1:1679078/76‑1 (MQ=255)                                        aaCCGCGCTACGTTGCCGGTGTTCTCGGCCCAACCAAACCGCACGGCGTCTATTTCTCCGGACGTCAACGATCCgg                                <  2:808608/76‑1 (MQ=255)                                         aCCGCGCTACGTTGCCGGTGTTCTCGGCCCAACCAAACCGCACGGCGTCTATTTCTCCGGACGTCAACGATCCGGc                               <  2:1570600/76‑1 (MQ=255)                                               cTACGTTGCCGGTGTTCTCGGCCCAACCAAACCGCACGGCGTCTATTTCTCCGGACGTCAACGATCCGGCATTTCg                         >  2:725584/1‑76 (MQ=255)                                               cTACGTTGCCGGTGTTCTCGGCCCAACCAAACCGCACGGCGTCTATTTCTCCGGACGTCAACGATCCGGCATTTCg                         >  2:1806325/1‑76 (MQ=255)                                                tACGTTGCCGGTGTTCTCGGCCCAACCAAACCGCACGGCGTCTATTTCTCCGGACGTCAACGATCCGGCATTTCGt                        <  2:626134/76‑1 (MQ=255)                                                 aCGTTGCCGGTGTTCTCGGCCCAACCAAACCGCACGGCGTCTATTTCTCCGGACGTCAACGATCCGGCATTTCGTa                       >  2:5678/1‑76 (MQ=255)                                                       ccGGTGTTCTCGGCCCAACCAAACCGCACGGCGTCTATTTCTCCGGACGTCAACGATCCGGCATTTCGTAATATCa                 <  2:1710542/76‑1 (MQ=255)                                                         ggTGTTCTCGGCCCAACCAAACCGCACGGCGTCTATTTCTCCGGACGTCAACGATCCGGCATTTCGTAATATCACt               <  1:1790972/76‑1 (MQ=255)                                                         ggTGTTCGCGGCCCAACCAAACCGCACGGCGTCTATTTCTCCGGACGTGAACGATCCGGCATTTCGTAATATCACt               <  1:1508659/76‑1 (MQ=255)                                                            gtTCTCGGCCCAACCAAACCGCACGGCGTCTATTTCTCCGGACGTCAACGATCCGGCATTTCGTAATATCACtttt            >  1:885334/1‑76 (MQ=255)                                                              tctcGGCCCAACCAAACCGCACGGCGTCTATTTCTCCGGACGTCAACGATCCGGCATTTCGTAATATCACTTTTGa          >  2:639642/1‑76 (MQ=255)                                                              tctcGGCCCAACCAAACCGCACGGCGTCTATTTCTCCGGACGTCAACGATCCGGCATTTCGTAATATCACTTTTGa          >  1:473150/1‑76 (MQ=255)                                                               ctcGGCCCAACCAAACCGCACGGCGTCTATTTCTCCGGACGTCAACGATCCGGCATTTCGTAATATCACTTTTGAc         >  1:717555/1‑76 (MQ=255)                                                                tcGGCCCAACCAAACCGCACGGCGTCTATTTCTCCGGACGTCAACGATCCGGCATTTCGTAATATCACTTTTGACg        >  2:1327990/1‑76 (MQ=255)                                                                tcGGCCCAACCAAACCGCACGGCGTCTATTTCTCCGGACGTCAACGATCCGGCATTTCGTAATATCACTTTTGACg        <  2:1301271/76‑1 (MQ=255)                                                                 cGGCCCAACCAAACCGCACGGCGTCTATTTCTCCGGACGTCAACGATCCGGCATTTCGTAATATCACTTTTGACgg       >  1:624650/1‑76 (MQ=255)                                                                  ggCCCAACCAAACCGCACGGCGTCTATTTCTCCGGACGTCAACGATCCGGCATTTCGTAATATCACTTTTGACggg      >  1:681933/1‑76 (MQ=255)                                                                    cccAACCAAACCGCACGGCGTCTATTTCTCCGGACGTCAACGATCCGGCATTTCGTAATATCACTTTTGACGGGCt    >  1:233117/1‑76 (MQ=255)                                                                    cccAACCAAACCGCACGGCGTCTATTTCTCCGGACGTCAACGATCCGGCATTTCGTAATATCACTTTTGACGGGCt    >  1:761630/1‑76 (MQ=255)                                                                      caaccaaaCCGCACGGCGTCTATTTCTCCGGACGTCAACGATCCGGCATTTCGTAATATCACTTTTGACGGGCtgg  >  2:1553691/1‑76 (MQ=255)                                                                           |                                                                       GCGCCTGTGCTGACGAGTGGACCGCGCGCACGCCGGAGAAACCGCGCTACGTTGCCGGTGTTCTCGGCCCAACC‑AACCGCACGGCGTCTATTTCTCCGGACGTCAACGATCCGGCATTTCGTAATATCACTTTTGACGGGCTGG  >  REL606/4203093‑4203236

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 20 ≤ ATCG/ATCG < 28 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑