Predicted mutation
evidence seq id position mutation annotation gene description
MC JC REL606 4,202,621 Δ579 bp [metH] [metH]

Missing coverage evidence...
   seq id start end size ←reads reads→ gene description
* * ÷ REL606 4202621 4203199 579 44 [1] [1] 43 [metH] [metH]

New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? REL606 = 42026201 (0.020)42 (0.790) 38/150 0.4 97.7% intergenic (‑61/‑139) iclR/metH DNA‑binding transcriptional repressor/B12‑dependent methionine synthase
?REL606 4203200 = 1 (0.020)coding (442/3684 nt) metH B12‑dependent methionine synthase

AATGGGTGCGACCATGACAGTCTCCTTTTTCTGTATCGTGGAAATCATTTTCATTTTTATTGTTAGCTAATGCAATAGTTGCTGAACTGATCCGATGAGTTAATGTTGAACAAATCTCATGTTGCGTGGTGGTCGCTTTTACCACAGAT  >  REL606/4202545‑4202693
                                                                           |                                                                         
aaTGGGTGCGACCATGACAGTCTCCTTTTTCTGTATCGTGGAAATCATTTTCATTTTTATTGTTAGCTAATGCAAt                                                                           <  1:505553/76‑1 (MQ=255)
       gCGACCATGACAGTCTCCTTTTTCTGTATCGTGGAAATCATTTTCATTTTTATTGTTAGCTAATGCAATgatccgg                                                                    <  2:290977‑M1/76‑8 (MQ=255)
           ccATGACAGTCTCCTTTTTCTGTATCGTGGAAATCATTTTCATTTTTATTGTTAGCTAATGCCATgatccggcatt                                                                >  2:420518‑M1/1‑65 (MQ=255)
             aTGACAGTCTCCTTTTTCTGTATCGTGGAAATCATTTTCATTTTTATTGTTAGCTAATGCAATgatccggcatttc                                                              >  2:1386204‑M1/1‑63 (MQ=255)
             aTGACAGTCTCCTTTTTCTGTATCGTGGAAATCATTTTCATTTTTATTGTTAGCTAATGCAATgatccggcatttc                                                              <  1:721025‑M1/76‑14 (MQ=255)
               gACAGTCTCCTTTTTCTGTATCGTGGAAATCATTTTCATTTTTATTGTTAGCTAATGCAATgatccggcatttcgt                                                            <  1:1343100‑M1/76‑16 (MQ=255)
                 cAGTCTCCTTTTTCTGTATCGTGGAAATCATTTTCATTTTTATTGTTAGCTAATGCAATgatccggcatttcgtaa                                                          <  2:354326‑M1/76‑18 (MQ=255)
                 cAGTCTCCTTTTTCTGTATCGTGGAAATCATTTTCATTTTTATTGTTAGCTAATGCAATgatccggcatttcgtaa                                                          <  2:1592240‑M1/76‑18 (MQ=255)
                   gTCTCCTTTTTCTGTATCGTGGAAATCATTTTCATTTTTATTGTTAGCTAATGCAATgatccggcatttcgtaata                                                        <  1:140290‑M1/76‑20 (MQ=255)
                    tctcCTTTTTCTGTATCGTGGAAATCATTTTCATTTTTATTGTTAGCTAATGCAATgatccggcatttcgtaatat                                                       <  1:1314004‑M1/76‑21 (MQ=255)
                    tctcCTTTTTCTGTATCGTGGAAATCATTTTCATTTTTATTGTTAGCTAATGCAATgatccggcatttcgtaatat                                                       >  1:786987‑M1/1‑56 (MQ=255)
                      tcCTTTTTCTGTATCGTGGAAATCATTTTCATTTTTATTGTTAGCTAATGCAATgatccggcatttcgtaatatca                                                     <  1:494169‑M1/76‑23 (MQ=255)
                       ccTTTTTCTGTATCGTGGAAATCATTTTCATTTTTATTGTTAGCTAATGCAATgatccggcatttcgtaatatcac                                                    >  2:905108‑M1/1‑53 (MQ=255)
                          ttttCTGTATCGTGGAAATCATTTTCATTTTTATTGTTAGCTAATGCAATgatccggcatttcgtaatatcacttt                                                 >  1:784762‑M1/1‑50 (MQ=255)
                          ttttCTGTATCGTGGAAATCATTTTCATTTTTATTGTTAGCTAATGCAATgatccggcatttcgtaatatcacttt                                                 >  1:1045847‑M1/1‑50 (MQ=255)
                          ttttCTGTATCGTGGAAATCATTTTCATTTTTATTGTTAGCTAATGCAATgatccggcatttcgtaatatcacttt                                                 <  1:1570418‑M1/76‑27 (MQ=255)
                           tttCTGTATCGTGGAAATCATTTTCATTTTTATTGTTAGCTAATGCAATgatccggcatttcgtaatatcactttt                                                <  2:652614‑M1/76‑28 (MQ=255)
                           tttCTGTATCGTGGAAATCATTTTCATTTTTATTGTTAGCTAATGCAATgatccggcatttcgtaatatcactttt                                                <  1:674689‑M1/76‑28 (MQ=255)
                            ttCTGTATCGTGGAAATCATTTTCATTTTTATTGTTAGCTAATGCAATgatccggcatttcgtaatatcacttttg                                               >  1:914213‑M1/1‑48 (MQ=255)
                                gTATCGTGGAAATCATTTTCATTTTTATTGTTAGCTAATGCAATgatccggcatttcgtaatatcacttttgacgg                                           >  2:416724‑M1/1‑44 (MQ=255)
                                 tATCGTGGAAATCATTTTCATTTTTATTGTTAGCTAATGCAATgatccggcatttcgtaatatcacttttgacggg                                          >  1:1475696‑M1/1‑43 (MQ=255)
                                    cGTGGAAATCATTTTCATTTTTATTGTTAGCTAATGCAATgatccggcatttcgtaatatcacttttgacgggctg                                       <  2:112883‑M1/76‑37 (MQ=255)
                                       ggAAATCATTTTCATTTTTATTGTTAGCTAATGCAATgatccggcatttcgtaatatcacttttgacgggctggtg                                    <  1:267029‑M1/76‑40 (MQ=255)
                                        gAAATCATTTTCATTTTTATTGTTAGCTAATGCAATgatccggcatttcgtaatatcacttttgacgggctggtgg                                   >  1:217377‑M1/1‑36 (MQ=255)
                                             cattttcattttTATTGTTAGCTAATGCAATgatccggcatttcgtaatatcacttttgacgggctggtggcggcc                              >  1:519275‑M1/1‑31 (MQ=255)
                                              attttcattttTATTGTTAGCTGATGCGATggtccggcgtttcgtgatatcacttttgacgggctggtggcggcct                             <  2:814186‑M1/76‑47 (MQ=255)
                                                tttcattttTATTGTTAGCTAATGCAATgatccggcatttcgtaatatcacttttgacgggctggtggcggcctat                           >  1:862252‑M1/1‑28 (MQ=255)
                                                 ttcattttTATTGTTAGCTAATGCAATgatccggcatttcgtaatatcacttttgacgggctggtggcggcctatc                          >  1:1614565‑M1/1‑27 (MQ=255)
                                                 ttcattttTATTGTTAGCTAATGCAATgatccggcatttcgtaatatcacttttgacgggctggtggcggcctatc                          <  2:590279‑M1/76‑50 (MQ=255)
                                                    attttTATTGTTAGCTAATGCAATgatccggcatttcgtaatatcacttttgacgggctggtggcggcctatcgag                       <  2:1624268‑M1/76‑53 (MQ=255)
                                                       tttATTGTTAGCTAATGCAATgatccggcatttcgtaatatcacttttgacgggctggtggcggcctatcgagagt                    >  1:169853‑M1/1‑21 (MQ=255)
                                                         tATTGTTAGCTAATGCAATgatccggcatttcgtaatatcacttttgacgggctggtggcggcctatcgagagtcc                  >  2:711168‑M1/1‑19 (MQ=255)
                                                         tATTGTTAGCTAATGCAATgatccggcatttcgtaatatcacttttgacgggctggtggcggcctatcgagagtcc                  >  1:1386536‑M1/1‑19 (MQ=255)
                                                           ttGTTAGCTAATGCAATgatccggcatttcgtaatatcacttttgacgggctgggggcggcctatcgagagtccac                >  1:1368087‑M1/1‑17 (MQ=255)
                                                            tGTTAGCTAATGCAATgatccggcatttcgtaatatcacttttgacgggctggtggcggcctatcgagagtccacc               >  1:94069‑M1/1‑16 (MQ=255)
                                                             gTTAGCTAATGCAATAGTTGCTGAACTGATCCGATGAGTTAATGTTGAACAAATCTCATGTTGCGTGGTGGTCGCt              >  1:1540725/1‑76 (MQ=255)
                                                              ttAGCTAATGCAATgatccggcatttcgtaatatcacttttgacgggctggtggcggcctatcgagagtccaccaa             <  1:192075‑M1/76‑63 (MQ=255)
                                                              ttAGCTAATGCAATgatccggcatttcgtaatatcacttttgacgggctggtggcggcctatcgagagtccaccaa             >  1:907268‑M1/1‑14 (MQ=255)
                                                                aGCTAATGCAATgatccggcatttcgtaatatcacttttgacgggctggtggcggcctatcgagagtccaccaaag           >  1:434962‑M1/1‑12 (MQ=255)
                                                                  cTAAGGCAATgatccggcattgcgtaatatcacttttgacgggctggtggcggcctatcgagagtccaccaaagcg         <  2:621819‑M1/76‑67 (MQ=255)
                                                                   tAATGCAATgatccggcatttcgtaatatcacttttgacgggctggtggcggcctatcgagagtccaccaaagcgc        >  2:1524889‑M1/1‑9 (MQ=255)
                                                                    aaTGCAATgatccggcatttcgtaatatcacttttgacgggctggtggcggcctatcgagagtccaccaaagcgct       >  1:1155686‑M1/1‑8 (MQ=255)
                                                                        cAATgatccggcatttcgtaatatcacttttgacgggctggtggcggcctatcgagagtccaccaaagcgctggtg   >  1:1143448‑M1/1‑4 (MQ=255)
                                                                         aaTgatccggcatttcgtaatatcacttttgacgggctggtggcggcctatcgagagtccaccaaagcgctggtgg  <  2:1157285‑M1/76‑74 (MQ=255)
                                                                           |                                                                         
AATGGGTGCGACCATGACAGTCTCCTTTTTCTGTATCGTGGAAATCATTTTCATTTTTATTGTTAGCTAATGCAATAGTTGCTGAACTGATCCGATGAGTTAATGTTGAACAAATCTCATGTTGCGTGGTGGTCGCTTTTACCACAGAT  >  REL606/4202545‑4202693

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 25 ≤ ATCG/ATCG < 31 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: 
Reads not counted as support for junction
read_name Not counted due to insufficient overlap past the breakpoint.
read_name Not counted due to not crossing MOB target site duplication.