Predicted mutation
evidence seq id position mutation annotation gene description
MC JC REL606 4,205,069 Δ5 bp coding (2311‑2315/3684 nt) metH → B12‑dependent methionine synthase

Missing coverage evidence...
   seq id start end size ←reads reads→ gene description
* * ÷ REL606 4205069 4205073 5 87 [0] [0] 88 metH B12‑dependent methionine synthase

New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? REL606 = 42050680 (0.000)79 (1.030) 55/96 0.1 100% coding (2310/3684 nt) metH B12‑dependent methionine synthase
?REL606 4205074 = 0 (0.000)coding (2316/3684 nt) metH B12‑dependent methionine synthase

GTGAAGGGCGACGTCCACGACATCGGTAAAAATATCGTTGGTGTGGTGCTGCAATGTAACAACTACGAAATTGTCGATCTCGGCGTTATGGTGCCTGCGGA  >  REL606/4205018‑4205118
                                                  |                                                  
gTGAAGGGCGACGTCCACGACATCGGTAAAAATATCGTTGGTGTGGTGCTG                                                    >  1:604‑M1/1‑51 (MQ=255)
 tGAAGGGCGACGTCCACGACATCGGTAAAAATATCGTTGGTGTGGTGCTGt                                                   >  1:1561482‑M1/1‑50 (MQ=255)
 tGAAGGGCGACGTCCACGACATCGGTAAAAATATCGTTGGTGTGGTGCTGt                                                   >  1:4636019‑M1/1‑50 (MQ=255)
 tGAAGGGCGACGTCCACGACATCGGTAAAAATATCGTTGGTGTGGTGCTGt                                                   <  1:6022572‑M1/51‑2 (MQ=255)
  gAAGGGCGACGTCCACGACATCGGTAAAAATATCGTTGGTGTGGTGCTGta                                                  <  1:762320‑M1/51‑3 (MQ=255)
  gAAGGGCGACGTCCACGACATCGGTAAAAATATCGTTGGTGTGGTGCTGta                                                  <  1:1962508‑M1/51‑3 (MQ=255)
  gAAGGGCGACGTCCACGACATCGGTAAAAATATCGTTGGTGTGGTGCTGta                                                  >  1:5939149‑M1/1‑49 (MQ=255)
  gAAGGGCGACGTCCACGACATCGGTAAAAATATCGTTGGTGTGGTGCTGta                                                  >  1:5246396‑M1/1‑49 (MQ=255)
    aGGGCGACGTCCACGACATCGGTAAAAATATCGTTGGTGTGGTGCTGtaac                                                <  1:2588665‑M1/51‑5 (MQ=255)
    aGGGCGACGTCCACGACATCGGTAAAAATATCGTTGGTGTGGTGCTGtaac                                                >  1:956415‑M1/1‑47 (MQ=255)
     gggCGACGTCCACGACATCGGTAAAAATATCGTTGGTGTGGTGCTGtaaca                                               >  1:3612646‑M1/1‑46 (MQ=255)
      ggCGACGTCCACGACATCGGTAAAAATATCGTTGGTGTGGTGCTGtaacaa                                              <  1:5903713‑M1/51‑7 (MQ=255)
       gCGACGTCCACGACATCGGTAAAAATATCGTTGGTGTGGTGCTGtaacaac                                             <  1:6716749‑M1/51‑8 (MQ=255)
       gCGACGTCCACGACATCGGTAAAAATATCGTTGGTGTGGTGCTGtaacaac                                             <  1:2374485‑M1/51‑8 (MQ=255)
        cGACGTCCACGACATCGGTAAAAATATCGTTGGTGTGGTGCTGtaacaact                                            <  1:6160930‑M1/51‑9 (MQ=255)
        cGACGTCCACGACATCGGTAAAAATATCGTTGGTGTGGTGCTGtaacaact                                            <  1:3429069‑M1/51‑9 (MQ=255)
          aCGTCCACGACATGGGTAAAAATATCGTTGGTGTGGTGCTGt                                                   <  1:6374668‑M1/42‑2 (MQ=255)
          aCGTCCACGACATCGGTAAAAATATCGTTGGTGTGGTGCTGtaacaactac                                          >  1:1373335‑M1/1‑41 (MQ=255)
          aCGTCCACGACATCGGTAAAAATATCGTTGGTGTGGTGCTGtaacaactac                                          <  1:1501316‑M1/51‑11 (MQ=255)
           cGTCCACGACATCGGTAAAAATATCGTTGGTGTGGTGCTGtaacaactacg                                         >  1:486075‑M1/1‑40 (MQ=255)
           cGTCCACGACATCGGTAAAAATATCGTTGGTGTGGTGCTGtaacaactacg                                         >  1:186396‑M1/1‑40 (MQ=255)
           cGTCCACGACATCGGTAAAAATATCGTTGGTGTGGTGCTGtaacaactacg                                         <  1:3560341‑M1/51‑12 (MQ=255)
           cGTCCACGACATCGGTAAAAATATCGTTGGTGTGGTGCTGtaa                                                 <  1:4236305‑M1/43‑4 (MQ=255)
            gTCCACGACATCGGTAAAAATATCGTTGGTGTGGTGCTGttacaactacga                                        >  1:6593047‑M1/1‑39 (MQ=255)
            gTCCACGACATCGGTAAAAATATCGTTGGTGTGGTGCTGtaacaactacga                                        <  1:5333591‑M1/51‑13 (MQ=255)
            gTCCACGACATCGGTAAAAATATCGTTGGTGTGGTGCTGtaacaactacga                                        <  1:605549‑M1/51‑13 (MQ=255)
            gTCCACGACATCGGTAAAAATATCGTTGGTGTGGTGCTGtaacaactacga                                        >  1:4553038‑M1/1‑39 (MQ=255)
            gTCCACGACATCGGTAAAAATATCGTTGGTGTGGTGCTGtaacaactacga                                        <  1:6046052‑M1/51‑13 (MQ=255)
             tCCACGACATCGGTAAAAATATCGTTGGTGTGGTGCTGtaacaactacgaa                                       <  1:1335677‑M1/51‑14 (MQ=255)
              ccACGACATCGGTAAAAATATCGTTGGTGTGGTGCTGtaacaactacgaaa                                      >  1:5215757‑M1/1‑37 (MQ=255)
               cACGACATCGGTAAAAATATCGTTGGTGTGGTGCTGtaacaactacgaaat                                     <  1:5810641‑M1/51‑16 (MQ=255)
                aCGACATCGGTAAAAATATCGTTGGTGTGGTGCTGtaacaactacgaaatt                                    <  1:4739221‑M1/51‑17 (MQ=255)
                 cGACATCGGTAAAAATATCGTTGGTGTGGTGCTGtaacaactacgaaattg                                   <  1:2649487‑M1/51‑18 (MQ=255)
                 cGACATCGGTAAAAATATCGTTGGTGTGGTGCTGtaacaactacgaaattg                                   >  1:4322337‑M1/1‑34 (MQ=255)
                 cGACATCGGTAAAAATATCGTTGGTGTGGTGCTGtaacaactacgaa                                       <  1:2566473‑M1/47‑14 (MQ=255)
                   aCATCGGTAAAAATATCGTTGGTGTGGTGCTGtaacaactacgaaattgtc                                 >  1:1990309‑M1/1‑32 (MQ=255)
                   aCATCGGTAAAAATATCGTTGGTGTGGTGCTGtaacaactacgaaattgtc                                 <  1:3093725‑M1/51‑20 (MQ=255)
                    cATCGGTAAAAATATCGTTGGTGTGGTGCTGtaacaactacgaaattgtcg                                <  1:435042‑M1/51‑21 (MQ=255)
                     aTCGGTAAAAATATCGTTGGTGTGGTGCTGtaacaactacgaaattgtcga                               >  1:2698980‑M1/1‑30 (MQ=255)
                     aTCGGTAAAAATATCGTTGGTGTGGTGCTGtaacaactacgaaattgtcga                               >  1:2584785‑M1/1‑30 (MQ=255)
                     aTCGGTAAAAATATCGTTGGTGTGGTGCTGtaacaactacgaaattgtcga                               >  1:2518651‑M1/1‑30 (MQ=255)
                      tCGGTAAAAATATCGTTGGTGTGGTGCTGtaacaactacgaaattgtcgat                              <  1:323279‑M1/51‑23 (MQ=255)
                      tCGGTAAAAATATCGTTGGTGTGGTGCTGtaacaactacgaaattgtcgat                              >  1:198381‑M1/1‑29 (MQ=255)
                        ggTAAAAATATCGTTGGTGTGGTGCTGtaacaactacgaaattgtcgatct                            >  1:1896616‑M1/1‑27 (MQ=255)
                        ggTAAAAATATCGTTGGTGTGGTGCTGtaacaactacgaaattgtcgatct                            >  1:4687148‑M1/1‑27 (MQ=255)
                         gTAAAAATATCGTTGGTGTGGTGCTGtaacaactacgaaattgtcgatctc                           >  1:6418414‑M1/1‑26 (MQ=255)
                         gTAAAAATATCGTTGGTGTGGTGCTGtaacaactacgaaattgtcgatctc                           >  1:2458588‑M1/1‑26 (MQ=255)
                          tAAAAATATCGTTGGTGTGGTGCTGtaacaactacgaaattgtcgatctcg                          <  1:7012600‑M1/51‑27 (MQ=255)
                          tAAAAATATCGTTGGTGTGGTGCTGtaacaactacgaaattgtcgatctcg                          <  1:6979705‑M1/51‑27 (MQ=255)
                          tAAAAATATCGTTGGTGTGGTGCTGtaacaactacgaaattgtcg                                >  1:753479‑M1/1‑25 (MQ=255)
                           aaaaaTATCGTTGGTGTGGTGCTGtaacaactacgaaattgtcgatctcgg                         >  1:3676363‑M1/1‑24 (MQ=255)
                           aaaaaTATCGTTGGTGTGGTGCTGtaacaactacgaaattgtcgatctcgg                         <  1:561602‑M1/51‑28 (MQ=255)
                           aaaaaTATCGTTGGTGTGGTGCTGtaacaactacgaaattgtcgatctcgg                         <  1:4927100‑M1/51‑28 (MQ=255)
                            aaaaTATCGTTGGTGTGGTGCTGtaacaactacgaaattgtcgatctcggc                        <  1:4985300‑M1/51‑29 (MQ=255)
                             aaaTATCGTTGGTGTGGTGCTGtaacaactacgaaattgtcgatctcggcg                       >  1:3311502‑M1/1‑22 (MQ=255)
                              aaTATCGTTGGTGTGGTGCTGtaacaactacgaaattgtcgatctcggcgt                      >  1:6700715‑M1/1‑21 (MQ=255)
                               atatCGTTGGTGTGGTGCTGtaacaactacgaaattgtcgatctcggcgtt                     <  1:4444899‑M1/51‑32 (MQ=255)
                               atatCGTTGGTGTGGTGCTGtaacaactacgaaattgtcgatctcggcgtt                     <  1:6339234‑M1/51‑32 (MQ=255)
                                  tCGTTGGTGTGGTGCTGtaacaactacgaaattgtcgatctcggcgttatg                  <  1:337530‑M1/51‑35 (MQ=255)
                                   cGTTGGTGTGGTGCTGtaacaactacgaaattgtcgatctcggcgttatgg                 >  1:6443899‑M1/1‑16 (MQ=255)
                                   cGTTGGTGTGGTGCTGtaacaactacgaaattgtcgatctcggcgttatgg                 >  1:4683197‑M1/1‑16 (MQ=255)
                                    gTTGGTGTGGTGCTGtaacaactacgaaattgtcgatctcggcgttatggt                <  1:940135‑M1/51‑37 (MQ=255)
                                     ttGGTGTGGTGCTGtaacaactacgaaattgtcgatctcggcgttatggtg               <  1:1028747‑M1/51‑38 (MQ=255)
                                     ttGGTGTGGTGCTGtaacaactacgaaattgtcgatctcggcgtta                    >  1:333278‑M1/1‑14 (MQ=255)
                                      tggtgtggtgCTGtaacaactacgaaattgtcgatctcggcgttatggtgc              >  1:7053918‑M1/1‑13 (MQ=255)
                                       ggtgtggtgCTGtaacaactacgaaattgtcgatctcggcgttatggtgcc             >  1:1141004‑M1/1‑12 (MQ=255)
                                        gtgtggtgCTGtaacaactacgaaattgtcgatctcggcgttatggtgcct            <  1:4441460‑M1/51‑41 (MQ=255)
                                        gtgtggtgCTGtaacaactacgaaattgtcgatctcggcgttatggtgcct            <  1:5306763‑M1/51‑41 (MQ=255)
                                        gtgtggtgCTGtaacaactacgaaattgtcgatctcggcgttatggtgcct            >  1:6263258‑M1/1‑11 (MQ=255)
                                         tgtggtgCTGtaacaactacgaaattgtcgatctcggcgttatggtgcctg           >  1:1877426‑M1/1‑10 (MQ=255)
                                         tgtggtgCTGtaacaactacgaaattgtcgatctcggcgttatggtgcctg           >  1:6232795‑M1/1‑10 (MQ=255)
                                         tgtggtgCTGtaacaactacgaaattgtcgatctcggcgttatggtgcctg           <  1:7296591‑M1/51‑42 (MQ=255)
                                         tgtggtgCTGtaacaactacgaaattgtcgatctcggcgttatggtgcctg           <  1:4117505‑M1/51‑42 (MQ=255)
                                          gtggtgCTGtaacaactacgaaattgtcgatctcggcgttatggtgcctgc          <  1:1179346‑M1/51‑43 (MQ=255)
                                          gtggtgCTGtaacaactacgaaattgtcgatctcggcgttatggtgcctgc          >  1:2763003‑M1/1‑9 (MQ=255)
                                          gtggtgCTGtaacaactacgaaattgtcgatctcggcgttatggtgcctgc          >  1:1473446‑M1/1‑9 (MQ=255)
                                          gtggtgCTGtaacaactacgaaattgtcgatctcggcgttatggtgcctgc          >  1:1696213‑M1/1‑9 (MQ=255)
                                            ggtgCTGtaacaactacgaaattgtcgatctcggcgttatggtgcctgcgg        >  1:5533124‑M1/1‑7 (MQ=255)
                                             gtgCTGtaacaactacgaaattgtcgatctcggcgttatggtgcctgcgga       <  1:4553195‑M1/51‑46 (MQ=255)
                                               gCTGtaacaactacgaaattgtcgatctcggcgttatggtgcctgcggaaa     >  1:5991646‑M1/1‑4 (MQ=255)
                                                 tGtaacaactacgaaattgtcgatctcggcgttatggtgcctgcggaaaaa   >  1:3278916‑M1/1‑2 (MQ=255)
                                                 tGtaacaactacgaaattgtcgatctcggcgttatggtgcctgcggaaaaa   <  1:1037955‑M1/51‑50 (MQ=255)
                                                 tGtaacaactacgaaattgtcgatctcggcgttatggtgcctgcggaaaaa   <  1:4514334‑M1/51‑50 (MQ=255)
                                                  gtaacaactacgaaattgtcgatctcggcgttatggtgcctgcggaaaaaa  <  1:312336‑M1/51‑51 (MQ=255)
                                                  gtaacaactacgaaattgtcgatctcggcgttatggtgcctgcggaaaaaa  <  1:5385309‑M1/51‑51 (MQ=255)
                                                  gtaacaactacgaaattgtcgatctcggcgttatggtgcctgcggaaaaaa  >  1:7291005‑M1/1‑1 (MQ=255)
                                                  gtaacaactacgaaattgtcgatctcggcgttatggtgcctgcggaaaaaa  >  1:2513170‑M1/1‑1 (MQ=255)
                                                  |                                                  
GTGAAGGGCGACGTCCACGACATCGGTAAAAATATCGTTGGTGTGGTGCTGCAATGTAACAACTACGAAATTGTCGATCTCGGCGTTATGGTGCCTGCGGA  >  REL606/4205018‑4205118

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 30 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: