Predicted mutation
evidence seq id position mutation annotation gene description
JC REL606 2,563,623 +GTAG coding (41/1014 nt) yfgA ← hypothetical protein

New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? REL606 = 25636181 (0.020)56 (0.900)
+GGTA
41/132 0.3 98.3% coding (46/1014 nt) yfgA hypothetical protein
?REL606 2563619 = 1 (0.020)coding (45/1014 nt) yfgA hypothetical protein

TTCAGGCAAAGTCGCTCGGCAACGGCCTGCTGACTAAGTCCTAGTTGTTCGCGAGCATTACGCAGGCGAGCGCCGGTAGTAAGTGCTTCATTTTGGTCGTGCGTGGCTTCAGTATTCATTCGCTGCAGGTACGTTTAAATGAGAATTG  >  REL606/2563545‑2563692
                                                                         |                                                                          
ttCAGGCAAAGTCGCTCGGCAACGGCCTGCTGACTAAGTCCTAGTTGTTCGCGAGCATTACGCAGGCGAGCGCCgg                                                                          >  1:968586‑M1/1‑74 (MQ=255)
 tCAGGCAAAGTCGCTCGGCAACGGCCTGCTGACTAAGTCCTAGTTGTTCGCGAGCATTACGCAGGCGAGCGCCggt                                                                         <  2:1086989‑M1/76‑4 (MQ=255)
 tCAGGCAAAGTCGCTCGGCAACGGCCTGCTGACTAAGTCCTAGTTGTTCGCGAGCATTACGCAGGCGAGCGCCggt                                                                         <  2:1299965‑M1/76‑4 (MQ=255)
 tCAGGCAAAGCCGCGCGGCAACGGCCTGCTGACTAAGTCCTAGTTGTTCGCGAGCATTACGCAGGCGAGCGCCggt                                                                         <  2:548441‑M1/76‑4 (MQ=255)
  cAGGCAAAGTCGCTCGGCAACGGCCTGCTGACTAAGTCCTAGTTGTTCGCGAGCATTACGCAGGCGAGCGCCggta                                                                        <  1:1860198‑M1/76‑5 (MQ=255)
   aGGCAAAGTCGCTCGGCAACGGCCTGCTGACTAAGTCCTAGTTGTTCGCGAGCATTACGCAGGCGAGCGCCggtag                                                                       <  2:500996‑M1/76‑6 (MQ=255)
   aGGCAAAGTCGCTCGGCAACGGCCTGCTGACTAAGTCCTAGTTGTTCGCGAGCATTACGCAGGCGAGCGCCggtag                                                                       >  2:1172512‑M1/1‑71 (MQ=255)
      cAAAGTCGCTCGGCAACGGCCTGCTGACTAAGTCCTAGTTGTTCGCGAGCATTACGCAGGCGAGCGCCggtaggta                                                                    >  2:1206069‑M1/1‑68 (MQ=255)
         aGTCGCTCGGCAACGGCCTGCTGACTAAGTCCTAGTTGTTCGCGAGCATTACGCAGGCGAGCGCCggtaggtagta                                                                 >  2:313762‑M1/1‑65 (MQ=255)
         aGTCGCTCGGCAACGGCCTGCTGACTAAGTCCTAGTTGTTCGCGAGCATTACGCAGGCGAGCGCCggtaggtagta                                                                 >  2:1337935‑M1/1‑65 (MQ=255)
            cGCTCGGCAACGGTCTGCTGACTAAGTCCTAGTTGTTCGCGAGCATTACGCAGGCGAGCGCCggtagggagtaagt                                                              >  2:571503‑M1/1‑62 (MQ=255)
            cGCTCGGCAACGGCCTGCTGACTAAGTCCTAGTTGTTCGCGAGCATTACGCAGGCGAGCGCCggtaggtagtaagt                                                              >  2:2046940‑M1/1‑62 (MQ=255)
             gCTCGGCAACGGCCTGCTGACTAAGTCCTAGTTGTTCGCGAGCATTACGCAGGCGAGCGCCggtaggtagtaagtg                                                             <  2:1571116‑M1/76‑16 (MQ=255)
                cGGCAACGGCCTGCTGACTAAGTCCTAGTTGTTCGCGAGCATTACGCAGGCGAGCGCCggtaggtagtaagtgctt                                                          >  2:440962‑M1/1‑58 (MQ=255)
                cGGCAACGGCCTGCTGACTAAGTCCTAGTTGTTCGCGAGCATTACGCAGGCGAGCGCCggtaggtagtaagtgctt                                                          >  1:1008505‑M1/1‑58 (MQ=255)
                  gCAACGGCCTGCTGACTAAGTCCTAGTTGTTCGCGAGCATTACGCAGGCGAGCGCCagtaggtagtaagtgcttca                                                        <  2:1078815‑M1/76‑21 (MQ=255)
                    aaCGGCCTGCTGACTAAGTCCTAGTTGTTCGCGAGCATTACGCAGGCGAGCGCCggtaggtagtaagtgcttcatt                                                      >  2:38307‑M1/1‑54 (MQ=255)
                     aCGGCCTGCTGACTAAGTCCTAGTTGTTCGCGAGCATTACGCAGGCGAGCGCCggtaggtagtaagtgcttcattt                                                     >  1:1506405‑M1/1‑53 (MQ=255)
                     aCGCCCTGCTGACTAAGTCCTAGTTGTTCGCGAGCATTACGCAGGCGAGCGCCggtaggtagtaagtgcttcattt                                                     <  1:1155006‑M1/76‑24 (MQ=255)
                      cGGCCTGCTGACTAAGTCCTAGTTGTTCGCGAGCATTACGCAGGCGAGCGCCggtaggtagtaagtgcttcatttt                                                    >  2:356481‑M1/1‑52 (MQ=255)
                         ccTGCTGACTAAGTCCTAGTTGTTCGCGAGCATTACGCAGGCGAGCGCCGGTAGTAAGTGCTTCATTTTGGTcgtg                                                 >  2:1862157/1‑76 (MQ=255)
                          ctgctgACTAAGTCCTAGTTGTTCGCGAGCATTACGCAGGCGAGCGCCggtaggtagtaagtgcttcattttggtc                                                >  1:699250‑M1/1‑48 (MQ=255)
                              tgAGTAAGTCCTAGTTGTTCGCGAGCATTACGCAGGCGAGCGCCggtaggtagtaagtgcttcattttggtcgtgc                                            >  1:317188‑M1/1‑44 (MQ=255)
                              tgACTAAGTCCTAGTTGTTCGCGAGCATTACGCAGGCGAGCGCCggtaggtagtaagtgcttcattttgatcgtgc                                            >  1:1292885‑M1/1‑44 (MQ=255)
                               gACTAAGTCCTAGTTGTTCGCGAGCATTACGCAGGCGAGCGCCggtaggtagtaagtgcttcattttggtcgtgcg                                           >  2:963363‑M1/1‑43 (MQ=255)
                               gACGAAGTCCTAGTGGTTCGCGAGCATTACGCGGGCGGGCGCCggtaggtagtaagtgcttcattttggtcgtgcg                                           <  1:475856‑M1/76‑34 (MQ=255)
                                aCTAAGTCCTAGTTGTTCGCGAGCATTACGCAGGCGAGCGCCggtaggtagtaagtgcttcattttggtcgtgcgt                                          >  2:1361659‑M1/1‑42 (MQ=255)
                                  tAAGTCCTAGTTGTTCGCGAGCATTACGCAGGCGAGCGCCggtaggtagtaagtgcttcattttggtcgtgcgtgg                                        <  1:1462833‑M1/76‑37 (MQ=255)
                                  tAAGTCCTAGTTGTTCGCGAGCATTACGCAGGCGAGCGCCggtaggtagtaagtgcttcattttggtcgtgcgtgg                                        >  1:1537366‑M1/1‑40 (MQ=255)
                                  tAAGGCCTAGTTGTTCGCGAGAATTACGCAGGCGAGCGCCggtaggtagtaagtgcttcattttggtcgtgcgtgg                                        <  1:1462844‑M1/76‑37 (MQ=255)
                                   aaGTCCTAGTTGTTCGCGAGCATTACGCAGGCGAGCGCCggtaggtagtaagtgcttcattttggtcgtgcgtggc                                       >  1:1160389‑M1/1‑39 (MQ=255)
                                   aaGTCCTAGTTGTTCGCGAGCATTACGCAGGCGAGCGCCggtaggtagtaagtgcttcattttggtcgtgcgtggc                                       >  2:1402764‑M1/1‑39 (MQ=255)
                                        cTAGTTGTTCGCGAGCATTACGCAGGCGAGCGCCggtaggtagtaagtgcttcattttggtcgtgcgtggcttcag                                  >  2:915321‑M1/1‑34 (MQ=255)
                                         tAGTTGTTCGCGAGCATTACGCAGGCGAGCGCCggtaggtagtaagtgcttcattttggtcgtgcgtggcttcagt                                 <  2:1828786‑M1/76‑44 (MQ=255)
                                          aGTTGTTCGCGAGCATTACGCAGGCGAGCGCCggtaggtagtaagtgcttcattttggtcgtgcgtggcttcagta                                >  1:2002582‑M1/1‑32 (MQ=255)
                                           gttgttCGCGAGCATTACGCAGGCGAGCGCCggtaggtagtaagtgcttcattttggtcgtgcgtggcttcagtat                               >  2:683727‑M1/1‑31 (MQ=255)
                                            ttgttCGCGAGCATTACGCAGGCGAGCGCCggtaggtagtaagtgcttcattttggtcgtgcgtggcttcagtatt                              >  2:1665219‑M1/1‑30 (MQ=255)
                                            ttgttCGCGAGCATTACGCAGGCGAGCGCCggtaggtagtaagtgcttcattttggtcgtgcgtggcttcagtatt                              <  2:473948‑M1/76‑47 (MQ=255)
                                            ttgttCGCGAGCATTACGCAGGCGAGCGCCggtaggtagtaagtgcttcattttggtcgtgcgtggcttcagtatt                              >  1:733928‑M1/1‑30 (MQ=255)
                                              gttCGCGAGCATTACGCAGGCGAGCGCCggtaggtagtaagtgcttcattttggtcgtgcgtggcttcagtattca                            >  1:1416800‑M1/1‑28 (MQ=255)
                                              gttCGCGAGATTTACGCAGGCGAGCGCCggtaggtagtaagtgcttcattttggtcgtgcgtggcttcagtattca                            <  1:1264586‑M1/76‑49 (MQ=255)
                                               ttCGCGAGCATTACGCAGGCGAGCGCCggtaggtagtaagtgcttcattttggtcgtgcgtggcttcagtattcat                           >  2:1411191‑M1/1‑27 (MQ=255)
                                                    gAGCATTACGCAGGCGAGCGCCggtaggtagtaagtgcttcattttggtcgtgcgtggcttcagtattcattcgct                      >  2:182295‑M1/1‑22 (MQ=255)
                                                     aGCATTACGCAGGCGAGCGCCggtaggtagtaagtgcttcattttggtcgtgcgtggcttcagtattcattcgctg                     <  2:1431652‑M1/76‑56 (MQ=255)
                                                      gCATTACGCGGGCGAGCGCCggtaggtagtaagtgcttcattttggtcgtgcgtggcttcagtattcattcgctgc                    <  1:1635490‑M1/76‑57 (MQ=255)
                                                       cATTACGCAGGCGAGCGCCggtaggtagtaagtgcttcattttggtcgggcggggcttcagtattcattcgctgcc                   >  1:1157382‑M1/1‑19 (MQ=255)
                                                         ttACGCGGGCGAGCGCCggtaggtagtaagtgcttcattttggtcgtgcgtggcttcagtattcattcgctgcagg                 <  1:574896‑M1/76‑60 (MQ=255)
                                                         ttACGCAGGCGAGCGCCggtaggtagtaagtgcttcattttggtcgtgcgtggcttcagtattcattcgctgcagg                 >  2:192372‑M1/1‑17 (MQ=255)
                                                          tACGCAGGCGAGCGCCggtaggtagtaagtgcttcattttggtcgtgcgtggcttcagtattcattcgctgcaggt                >  1:778289‑M1/1‑16 (MQ=255)
                                                            cgcgGGCGAGCGCCggtaggtagtaagtgcttcattttggtcgtgcttggcttcagtattcattcgctgcaggtac              <  2:54698‑M1/76‑63 (MQ=255)
                                                            cgcgGGCGAGCGCCggtaggtagtaagtgcttcattttggtcgtgcgtggcttcagtattcattcgctgcaggtac              <  2:69807‑M1/76‑63 (MQ=255)
                                                            cGCAGGCGAGCGCCggtaggtagtaagtgcttcattttggtcgtgcgtggcttcagtattcattcgctgcaggtac              <  2:614167‑M1/76‑63 (MQ=255)
                                                             gCAGGCGAGCGCCggtaggtagtaagtgtttcattttggtcgtgcgtggcttcagtattcattcgctgcaggtacg             >  1:1404722‑M1/1‑13 (MQ=255)
                                                             gCAGGCGAGCGCCggtaggtagtaagtgcttcattttggtcgtgcgtggcttcagtattcattcgctgcaggtacg             >  1:1712752‑M1/1‑13 (MQ=255)
                                                               aGGCGAGCGCCggtaggtagtaagtgcttcattttggtcgtgcgtggcttcagtattcattcgctgcaggtacgtt           >  2:1470879‑M1/1‑11 (MQ=255)
                                                                 gCGAGCGCCggtaggtagtaagtgcttcattttggtcgtgcgtggcttcagtattcattcgctgcaggtacgttta         >  2:1865936‑M1/1‑9 (MQ=255)
                                                                    aGCGCCggtaggtagtaagtgcttcattttggtcgtgcgtggcttcagtattcattcgctgcaggtacgtttaaat      <  2:1691844‑M1/76‑71 (MQ=255)
                                                                     gcgcCggtaggtagtaagtgcttcattttggtcgtgcgtggcttcagtattcattcgctgcaggtacgtttaaatg     <  1:360189‑M1/76‑72 (MQ=255)
                                                                        ccggtaggtagtaagtgcttcattttggtcgtgcgtggcttcagtattcattcgctgcaggtacgtttaaatgaga  >  2:1099464‑M1/1‑2 (MQ=255)
                                                                         |                                                                          
TTCAGGCAAAGTCGCTCGGCAACGGCCTGCTGACTAAGTCCTAGTTGTTCGCGAGCATTACGCAGGCGAGCGCCGGTAGTAAGTGCTTCATTTTGGTCGTGCGTGGCTTCAGTATTCATTCGCTGCAGGTACGTTTAAATGAGAATTG  >  REL606/2563545‑2563692

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