Predicted mutation
evidence seq id position mutation annotation gene description
JC REL606 2,362,391 +GGA coding (378/543 nt) yfcD ← predicted NUDIX hydrolase

New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? REL606 = 23623881 (0.020)40 (0.850)
+GGA
30/138 0.4 97.6% coding (381/543 nt) yfcD predicted NUDIX hydrolase
?REL606 2362389 = 1 (0.020)coding (380/543 nt) yfcD predicted NUDIX hydrolase

CATCGCAGCGTGCGGTGATTTCTTCCGGCGTCAGCCAGCAAACTTCACTGACTTCATCTTCCTGTAGGGCGAAGGGACCGTGAGAGACGCAGCTGAACAACGCGCCCCAGACACGGCAATTTTTATCTTCGAAATAGAACTGCCCGT  >  REL606/2362315‑2362461
                                                                         |                                                                         
cATCGCCGCGCTCCGAGATTTCTTCCGGCGTCAGCCAGCAACCTTCACTGACTTCATCTTCCTTTAGTGCGAAggg                                                                         >  2:332432/1‑76 (MQ=255)
cATCGCAGCGTGCGGTGATTTCTTCCGGCGTCAGCCAGCAAACTTCACTGACTTCATCTTCCTGTAGGGCGAAGgg                                                                         >  2:332420‑M1/1‑74 (MQ=255)
cATCGCAGCGTGCGGTGATTTCTTCCGGCGTCAGCCAGCAAACTTCACTGACTTCATCTTCCTGTAGGGCGAAGgg                                                                         >  2:1208931‑M1/1‑74 (MQ=255)
cATCGCAGCGTGCGGTGATTTCTTCCGGCGTCAGCCAGCAAACTTCACTGACTTCATCTTCCTGTAGGGCGAAGgg                                                                         <  1:1545629‑M1/76‑3 (MQ=255)
  tCGCAGCGTGCGGTGATTTCTTCCGGCGTCAGCCAGCAAACTTCACTGACTTCATCTTCCTGTAGGGCGAAGggag                                                                       <  1:84903‑M1/76‑5 (MQ=255)
   cGCAGCGTGCGGTGATTTCTTCCGGCGTCAGCCAGCAAACTTCACTGACTTCATCTTCCTGTAGGGCGAAGggagg                                                                      >  1:1163826‑M1/1‑71 (MQ=255)
   cGCAGCGTGCGGTGATTTCTTCCGGCGTCAGCCAGCAAACTTCACTGACTTCATCTTCCTGTAGGGCGAAGggagg                                                                      >  2:465017‑M1/1‑71 (MQ=255)
   cGCAGCGTGCGGTGATTTCTTCCGGCGTCAGCCAGCAAACTTCACTGACTTCATCTTCCTGTAGGGCGAAGggagg                                                                      >  1:1141995‑M1/1‑71 (MQ=255)
   cGCAGCGTGCGGTGATTTCTTCCGGCGTCAGCCAGCAAACTTCACTGACTTCATCTTCCTGTAGGGCGAAGggagg                                                                      >  1:1322399‑M1/1‑71 (MQ=255)
     cAGCGTGCGGTGATTTCTTCCGGCGTCAGCCAGCAAACTTCACTGACTTCATCTTCCTGTAGGGCGAAGggaggac                                                                    >  1:509521‑M1/1‑69 (MQ=255)
      aGCGTGCGGTGATTTCTTCCGGCGTCAGCCAGCAAACTTCACTGACTTCATCTTCCTGTAGGGCGAAGggaggacc                                                                   >  2:1518745‑M1/1‑68 (MQ=255)
          tGCGGTGATTTCTTCCGGCGTCAGCCAGCAAACTTCACTGACTTCATCTTCCTGTAGGGCGAAGggaggaccgtga                                                               >  2:1229245‑M1/1‑64 (MQ=255)
            cGGTGATTTCTTCCGGCGTCAGCCAGCAAACTTCACTGACTTCATCTTCCTGTAGGGCGAAGggaggaccgtgaga                                                             >  1:1271404‑M1/1‑62 (MQ=255)
                 aTTTCTTCCGGCGTCAGCCAGCAAACTTCACTGACTTCATCTTCCTGTAGGGCGAAGggaggaccgtgagagacgc                                                        >  1:1225589‑M1/1‑57 (MQ=255)
                 aTTTCTTCCGGCGTCAGCCAGCAAACTTCACTGACTTCATCTTCCTGTAGGGCGAAGggaggaccgtgagagacgc                                                        <  2:1130604‑M1/76‑20 (MQ=255)
                        ccGGCGTCAGCCAGCAAACTTCACTGACTTCATCTTCCTGTAGGGCGAAGggaggaccgtgagagacgcagctgaa                                                 <  1:471325‑M1/76‑27 (MQ=255)
                          ggCGTCAGCCAGCAAACTTCACTGACTTCATCTTCCTGTAGGGCGAAGggaggaccgtgagagacgcagctgaaca                                               >  2:437046‑M1/1‑48 (MQ=255)
                           gCGTCAGCCAGCAAACTTCACTGACTTCATCTTCCTGTAGGGCGAAGggaggaccgtgagagacgcagctgaacaa                                              <  2:720864‑M1/76‑30 (MQ=255)
                            cGTCAGCCAGCAAACTTCACTGACTTCATCTTCCTGTAGGGCGAAGggaggaccgtgagagacgcagctgaacaac                                             >  1:1165940‑M1/1‑46 (MQ=255)
                            cGTCAGCCAGCAAACTTCACTGACTTCATCTTCCTGTAGGGCGAAGggaggaccgtgagagacgcagctgaacaac                                             <  2:1301963‑M1/76‑31 (MQ=255)
                            cGTCAGCCAGCAAACTTCACTGACTTCATCTTCCTGTAGGGCGAAGggaggaccgtgagagacgcagctgaacaaa                                             >  2:313613‑M1/1‑46 (MQ=255)
                              tCAGCCAGCAAACTTCACTGACTTCATCTTCCTGTAGGGCGAAGggaggaccgtgagagacgcagctgaacaacgc                                           <  1:1532003‑M1/76‑33 (MQ=255)
                                    agcAAACTTCACTGACTTCATCTTCCTGTAGGGCGAAGggaggaccgtgagagacgcagctgaacaacgcgcccca                                     <  2:1521917‑M1/76‑39 (MQ=255)
                                    agcAAACTTCACTGACTTCATCTTCCTGTAGGGCGAAGggaggaccgtgagagacgcagctgaacaacgcgcccca                                     <  1:572520‑M1/76‑39 (MQ=255)
                                       aaaCTTCACTGACTTCATCTTCCTGTAGGGCGAAGggaggaccgtgagagacgcagctgaacaacgcgccccagac                                  <  1:1313542‑M1/76‑42 (MQ=255)
                                       aaaCTTCACTGACTTCATCTTCCTGTAGGGCGAAGggaggaccgtgagagacgcagctgaacaacgcgccccagac                                  <  1:1357137‑M1/76‑42 (MQ=255)
                                       aaaCTTCACTGACTTCATCTTCCTGTAGGGCGAAGggaggaccgtgagagacgcagctgaacaacgcgccccagac                                  <  1:144829‑M1/76‑42 (MQ=255)
                                           ttCACTGACTTCATCTTCCTGTAGGGCGAAGggaggaccgtgagagatgcagctgaacaacgcgccccagacacgg                              >  2:421866‑M1/1‑31 (MQ=255)
                                               cTGACTTCATCTTCCTGTAGGGCGAAGggaggaccgtgagagacgcagctgaacaacgcgccccagacacggcaat                          >  2:202725‑M1/1‑27 (MQ=255)
                                                tGACTTCATCTTCCTGTAGGGCGAAGggaggaccgtgagagacgcagctgaacaacgcgccccagacacggcaatt                         <  1:907451‑M1/76‑51 (MQ=255)
                                                tGACTTCATCTTCCTGTAGGGCGAAGggaggaccgtgagagacgcagctgaacaacgcgccccagacacggcaatt                         <  1:1471860‑M1/76‑51 (MQ=255)
                                                 gACTTCATCTTCCTGTAGGGCGAAGggaggaccgtgagagacgcagctgaacaacgcgccccagacacggcacttt                        >  1:509584‑M1/1‑25 (MQ=255)
                                                 gACTTCATCTTCCTGTAGGGCGAAGggaggaccgtgagagacgcagctgaacaacgcgccccagacacggcaattt                        <  2:790555‑M1/76‑52 (MQ=255)
                                                  aCTTCATCTTCCTGTAGGGCGAAGggaggaccgtgagagacgcagctgaacaacgcgccccagacacggcaatttt                       >  2:1155165‑M1/1‑24 (MQ=255)
                                                   cTTCATCTTCCTGTAGGGCGAAGggaggaccgtgagagacgcagctgaacaacgcgccccagacacggcaattttt                      <  1:37797‑M1/76‑54 (MQ=255)
                                                    ttCATCTTCCTGTAGGGCGAAGggaggaccgtgagagacgcagctgaacaacgcgccccagacacggcaattttta                     >  1:346772‑M1/1‑22 (MQ=255)
                                                       aTCTTCCTGTAGGGCGAAGggaggaccgtgagagacgcagctgaacaacgcgccccagacacggcaatttttatct                  <  2:887602‑M1/76‑58 (MQ=255)
                                                         cTTCCTGTAGGGCGAAGggaggaccgtgagagacgcagctgaacaacgcgccccagacacggcaatttttatcttc                >  2:1035429‑M1/1‑17 (MQ=255)
                                                          ttCCTGTAGGGCGAAGggaggaccgtgagagacgcagctgaacaacgcgccccagacacggcaatttttatcttcg               <  2:1058808‑M1/76‑61 (MQ=255)
                                                                 aGGGCGAAGggaggaccgtgagagacgcagctgaacaacgcgccccagacacggcaatttttatcttcgaaataga        <  1:1492212‑M1/76‑68 (MQ=255)
                                                                   ggCGAAGggaggaccgtgagagacgcagctgaacaacgcgccccagacacggcaatttttatcttcgaaatagaac      >  1:1391876‑M1/1‑7 (MQ=255)
                                                                       aaGggaggaccgtgagagacgcagctgaacaacgcgccccagacacggcaatttttatcttcgaaatagaactgcc  >  2:393450‑M1/1‑3 (MQ=255)
                                                                       aaGggaggaccgtgagagacgcagctgaacaacgcgccccagacacggcaatttttatcttcgaaatagaactgcc  >  2:1035010‑M1/1‑3 (MQ=255)
                                                                       aaGggaggaccgtgagagacgcagctgaacaacgcgccccagacacggcaatttttatcttcgaaatagaactgcc  >  2:585495‑M1/1‑3 (MQ=255)
                                                                       aaGggaggaccgtgagagacgcagctgaacaacgcgccccagacacggcaatttttatcttcgaaatagaactgcc  >  2:600989‑M1/1‑3 (MQ=255)
                                                                       aaGggagaaccgtgagagacgcagctgaacaacgcgccccagacacggcaatttttatcttcgaaatagaactgcc  >  1:49293‑M1/1‑3 (MQ=255)
                                                                         |                                                                         
CATCGCAGCGTGCGGTGATTTCTTCCGGCGTCAGCCAGCAAACTTCACTGACTTCATCTTCCTGTAGGGCGAAGGGACCGTGAGAGACGCAGCTGAACAACGCGCCCCAGACACGGCAATTTTTATCTTCGAAATAGAACTGCCCGT  >  REL606/2362315‑2362461

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 20 ≤ ATCG/ATCG < 28 ≤ ATCG/ATCG < 33 ≤ 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.