Predicted mutation
evidence seq id position mutation annotation gene description
MC JC REL606 4,097,564 Δ4 bp coding (6‑9/846 nt) glpF ← glycerol facilitator

Missing coverage evidence...
   seq id start end size ←reads reads→ gene description
* * ÷ REL606 4097564 4097567 4 50 [0] [0] 51 glpF glycerol facilitator

New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? REL606 = 40975630 (0.000)49 (0.780) 34/150 1.0 100% coding (10/846 nt) glpF glycerol facilitator
?REL606 4097568 = 0 (0.000)coding (5/846 nt) glpF glycerol facilitator

CGCAACCCACACCGAAGAAAATCAACAACCCGGTACCGAGGAATTCAGCAATGCACTGGCCTTTCAAGGTTGATGTTTGACTCATAATCGGATCCTGAAGAGTTAATGTTTGTTGTATGCGTGAAAGTCACGGACCTCCACGATGCT  >  REL606/4097488‑4097634
                                                                           |                                                                       
cgcaacCCACACCGAAGAAAATCAACAACCCGGTACCGAGGAATTCAGCAATGCACTGGCCTTTCAAGGTTGATGt                                                                         >  1:1259777/1‑76 (MQ=255)
 gCAACCCACACCGAAGAAAATCAACAACCCGGTACCGAGGAATTCAGCAATGCACTGGCCTTTCAAGGTTGATGTc                                                                        >  2:435715‑M1/1‑75 (MQ=255)
  cAACCCACACCGAAGAAAATCAACAACCCGGTACCGAGGAATTCAGCAATGCACTGGCCTTTCAAGGTTGATGTct                                                                       >  2:886940‑M1/1‑74 (MQ=255)
  cAACCCACACCGAAGAAAATCAACAACCCGGTACCGAGGAATTCAGCAATGCACTGGCCTTTCAAGGTTGATGTct                                                                       <  1:1340384‑M1/76‑3 (MQ=255)
  cAACCCACACCGAAGAAAATCAACAACCCGGTACCGAGGAATTCAGCAATGCACTGGCCTTTCAAGGTTGATGTct                                                                       >  2:1788100‑M1/1‑74 (MQ=255)
  cAACCCACACCGAAGAAAATCAACAACCCGGAACCGAGGAATTCAGCAATGCACTGGCCTTTCAAGGTTGATTTct                                                                       >  2:886953‑M1/1‑74 (MQ=255)
   aaCCCACACCGAAGAAAATCAACAACCCGGTACCGAGGAATTCAGCAATGCACTGGCCTTTCAAGGTTGATGTctc                                                                      >  2:1033795‑M1/1‑73 (MQ=255)
    aCCCACACCGAAGAAAATCAACAACCCGGTACCGAGGAATTCAGCAATGCACTGGCCTTTCAAGGTTGATGTctca                                                                     <  1:601043‑M1/76‑5 (MQ=255)
    aCCCACACCGAAGAAAATCAACAACCCGGTACCGAGGAATTCAGCAATGCACTGGCCTTTCAAGGTTGATGTctca                                                                     >  1:284951‑M1/1‑72 (MQ=255)
       cacaCCGAAGAAATTCAACAACCCGGTACCGAGGAATTCAGCAATGCACTGGCCTTTCAAGGTTGATGTctcataa                                                                  <  2:939459‑M1/76‑8 (MQ=255)
       cacaCCGAAGAAAATCAACAACCCGGTACCGAGGAATTCAGCAATGCACTGGCCTTTCAAGGTTGATGTctcataa                                                                  <  1:1666369‑M1/76‑8 (MQ=255)
        acacCGAAGAAAATCAACAACCCGGTACCGAGGAATTCAGCAATGCACTGGCCTTTCAAGGTTGATGTctcataat                                                                 >  2:540476‑M1/1‑68 (MQ=255)
         cacCGAAGAAAATCAACAACCCGGTACCGAGGAATTCAGCAATGCACTGGCCTTTCAAGGTTGATGTctcataatc                                                                <  1:1870099‑M1/76‑10 (MQ=255)
         cacCGAAGAAAATCAACAACCCGGTACCGAGGAATTCAGCAATGCACTGGCCTTTCAAGGTTGATGTctcataatc                                                                <  1:1183520‑M1/76‑10 (MQ=255)
           ccGAAGAAAATCAACAACCCGGTACCGAGGAATTCAGCAATGCACTGGCCTTTCAAGGTTGATGTctcataatcgg                                                              <  1:1531313‑M1/76‑12 (MQ=255)
               agaaAATCAACAACCCGGTACCGAGGAATTCAGCAATGCACTGGCCTTTCAAGGTTGATGTctcataatcggatcc                                                          <  1:1846278‑M1/76‑16 (MQ=255)
                 aaaaTCAACAACCCGGTACCGAGGAATTCAGCAATGCACTGGCCTTTCAGGGTTGATGTctcataatcggatcctg                                                        >  1:543197‑M1/1‑59 (MQ=255)
                 aaaaTCAACAACCCGGTACCGAGGAATTCAGCAATGCACTGGCCTTTCAAGGTTGATGTctcataatcggatcctg                                                        <  1:1998001‑M1/76‑18 (MQ=255)
                 aaaaTCAACAACCCGGTACCGAGGAATTCAGCAATGCACTGGCCTTTCAAGGTTGATGTctcataatcggatcctg                                                        <  2:168977‑M1/76‑18 (MQ=255)
                 aaaaTCAACAACCCGGTACCGAGGAATTCAGCAATGCACTGGCCTTTCAAGGTTGATGTctcataatcggatcctg                                                        <  2:1388638‑M1/76‑18 (MQ=255)
                         caacCCGGTACCGAGGAATTCAGCAATGCACTGGCCTTTCAAGGTTGATGTctcataatcggatcctgaagagtta                                                >  2:1438334‑M1/1‑51 (MQ=255)
                         caacCCGGTACCGAGGAATTCAGCAATGCACTGGCCTTTCAAGGTTGATGTctcataatcggatcctgaagagtta                                                >  1:1406809‑M1/1‑51 (MQ=255)
                           acCCGGTACCGAGGAATTCAGCAATGCACTGGCCTTTCAAGGTTGATGTctcataatcggatcctgaagagttaat                                              <  2:70381‑M1/76‑28 (MQ=255)
                               ggTACCGAGGAATTCAGCAATGCACTGGCCTTTCAAGGTTGATGTctcataatcggatcctgaagagttaatgttt                                          <  2:1084770‑M1/76‑32 (MQ=255)
                                   ccGAGGAATTCAGCAATGCACTGGCCTTTCAAGGTTGATGTctcataatcggatcctgaagagttaatgtttgttg                                      >  1:1299249‑M1/1‑41 (MQ=255)
                                    cGAGGAATTCAGCAATGCACTGGCCTTTCAAGGTTGATGTctcataatcggatcctgaagagttaatgtttgttgt                                     >  1:1821815‑M1/1‑40 (MQ=255)
                                     gAGGAATTCAGCAATGCACTGGCCTTTCAAGGTTGATGTctcataatcggatcctgaagagttaatgtttgttgta                                    >  2:1375455‑M1/1‑39 (MQ=255)
                                        gAATTCAGCAATGCACTGGCCTTTCAAGGTTGATGTctcataatcggatcctgaagagttaatgtttgttgtatgc                                 <  2:1116403‑M1/76‑41 (MQ=255)
                                           ttCAGCAATGCACTGGCCTTTCAAGGTTGATGTctcataatcggatcctgaagagttaatgtttgttgtatgcgtg                              <  1:167333‑M1/76‑44 (MQ=255)
                                           ttCAGCAATGCACTGGCCTTTCAAGGTTGATGTctcataatcggatcctgaagagttaatgtttgttgtatgcgtg                              <  1:1170557‑M1/76‑44 (MQ=255)
                                              aGCAATGCACTGGCCTTTCAAGGTTGATGTctcataatcggatcctgaagagttaatgtttgttgtatgcgtgaaa                           >  2:1758607‑M1/1‑30 (MQ=255)
                                                cAATGCACTGGCCTTTCAAGGTTGATGTctcataatcggatcctgaagagttaatgtttgttgtatgcgtgaaagt                         <  2:227478‑M1/76‑49 (MQ=255)
                                                cAATGCACTGGCCTTTCAAGGTTGATGTctcataatcggatcctgaagagttaatgtttgttgtatgcgtgaaagt                         <  1:1033796‑M1/76‑49 (MQ=255)
                                                    gCACTGGCCTTTCAAGGTTGATGTctcataatcggatcctgaagagttaatgtttgttgtatgcgtgaaagtcacg                     <  1:599290‑M1/76‑53 (MQ=255)
                                                    gCACTGGCCTTTCAAGGTTGATGTctcataatcggatcctgaagagttaatgtttgttgtatgcgtgaaagtcacg                     >  2:252781‑M1/1‑24 (MQ=255)
                                                     cACTGGCCTTTCAAGGTTGATGTctcataatcggatcctgaagagttaatgtttgttgtatgcgtgaaagtcacgg                    <  1:933257‑M1/76‑54 (MQ=255)
                                                     cACTGGCCTTTCAAGGTTGATGTctcataatcggatcctgaagagttaatgtttgttgtatgcgtgaaagtcacgg                    <  1:875150‑M1/76‑54 (MQ=255)
                                                     cACTGGCCTTTCAAGGTTGATGTctcataatcggatcctgaagagttaatgtttgttgtatgcgtgaaagtcacgg                    <  1:292497‑M1/76‑54 (MQ=255)
                                                        tGGCCTTTCAAGGTTGATGTctcataatcggatcctgaagagttaatgtttgttgtatgcgtgaaagtcacggacc                 >  2:1303774‑M1/1‑20 (MQ=255)
                                                        tGGCCTTTCAAGGTTGATGTctcataatcggatcctgaagagttaatgtttgttgtatgcgtgaaagtcacggacc                 >  1:889704‑M1/1‑20 (MQ=255)
                                                         ggCCTTTCAAGGTTGATGTctcataatcggatcctgaagagttaatgtttgttgtatgcgtaaaagtcacggacct                <  1:1421497‑M1/76‑58 (MQ=255)
                                                            cTTTCAAGGTTGATGTctcataatcggatcctgaagagttaatgtttgttgtatgcgtgaaagtcacggacctcca             <  2:1484309‑M1/76‑61 (MQ=255)
                                                               tCAAGGTTGATGTctcataatcggatcctgaagagttaatgtttgttgtatgcgtgaaagtcacggacctccacga          >  1:234478‑M1/1‑13 (MQ=255)
                                                                cAAGGTTGATGTctcataatcggatcctgaagagttaatgtttgttgtatgcgtgaaagtcacggacctccacgat         >  1:1497614‑M1/1‑12 (MQ=255)
                                                                cAAGGTTGATGTctcataatcggatcctgaagagttaatgtttgttgtatgcgtgaaagtcacggacctccacgat         >  1:1438668‑M1/1‑12 (MQ=255)
                                                                cAAGGTTGATGTctcataatcggatcctgaagagttaatgtttgttgtatgcgtgaaagtcacggacctccacgat         >  1:332766‑M1/1‑12 (MQ=255)
                                                                cAAGGTTGATGTctcataatcggatcctgaagagttaatgtttgttgtatgcgtgaaagtcacggacctccacgat         >  1:1432641‑M1/1‑12 (MQ=255)
                                                                 aaGGTTGATGTctcataatcggatcctgaagagttaatgtttgttgtatgcgtgaaagtcacggacctccacgatg        >  2:122865‑M1/1‑11 (MQ=255)
                                                                    gTTGATGTctcataatcggatcctgaagagttaatgtttgttgtatgcgtgaaagtcacggacctccacgatgctt     <  2:451508‑M1/76‑69 (MQ=255)
                                                                       gATGTctcataatcggatcctgaagagttaatgtttgttgtatgcgtgaaagtcacggacctccacgatgcttgta  <  1:1502715‑M1/76‑72 (MQ=255)
                                                                           |                                                                       
CGCAACCCACACCGAAGAAAATCAACAACCCGGTACCGAGGAATTCAGCAATGCACTGGCCTTTCAAGGTTGATGTTTGACTCATAATCGGATCCTGAAGAGTTAATGTTTGTTGTATGCGTGAAAGTCACGGACCTCCACGATGCT  >  REL606/4097488‑4097634

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.