Predicted mutation
evidence seq id position mutation annotation gene description
JC JC REL606 1,916,224 IS1 (–) +9 bp [pykA] [pykA]

New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? REL606 1916224 =1 (0.020)62 (0.960) 50/200 0.3 98.3% intergenic (+127/‑1) yebK/pykA predicted DNA‑binding transcriptional regulator/pyruvate kinase
?REL606 1973272 = NA (NA)noncoding (1/768 nt) IS1 repeat region
* ? REL606 = 19162321 (0.020)68 (1.060) 53/200 0.2 98.4% coding (8/1443 nt) pykA pyruvate kinase
?REL606 2034326 = NA (NA)noncoding (768/768 nt) IS1 repeat region

CAATGCCCGGCATACTCATGTCCATTAGCACCACGTCAACGGCATTTGCCCGACACCACTTAACGGCGTCTTCACCGCACGATGCCTCACCGACGACTTTGGTGATGCTGCCAACTTACTGATTTAGTGTATGATGGTGTTTTTGAGGTGCTCCAGTGGCTTCTGTTTCTATCAGCTGTCCCTCCTGTTCAGCTACT  >  REL606/1973172‑1973368
                                                                                                    |                                                                                                
cAATGCCCGGCATACTCATGTCCATTAGCACCACGTCAACGGCATTTGCCCGACACCACTTAACGGCGTCTTCACCGCACGATGCCTCACCGACGACTTTg                                                                                                  >  2:1356502/1‑101 (MQ=255)
 aaTTCCCGGCATACTCATGTCCATTAGCACCACGTCAACGGCATTTGCCCGACACCACTTAACGGCGTCTTCACCGCACGATGCCTCACCGACGACTTTgg                                                                                                 <  2:853198/101‑1 (MQ=255)
     cccGGCATACTCATGTCCATTAGCACCACGTCAACGGCATTTGCCCGACACCACTTAACGGCGTCTTCACCGCACGATGCCTCACCGACGACTTTGGTGAt                                                                                             <  1:451068/101‑1 (MQ=255)
      ccGGCATGCTGATTGCCATTAGCACCACGTCAACGGCATTTGCCTGACGCCACTTAACGGTGTCTTCCGCGCTCGATTCCTTGCCGACGACTTTGGTGAtg                                                                                            <  2:824270/101‑1 (MQ=255)
      ccGGCATACTCATGTCCATTAGCACCACGTCAACGGCATTTGCCCGACACCACTTAACGGCGTCTTCACCGCACGATGCCTCACCGACGACTTTGGTGAtg                                                                                            <  2:208532/101‑1 (MQ=255)
       cGGCATACTCATGTCCATTAGCACCACGTCAACGGCATTTGCCCGACACCACTTAACGGCGTCTTCACCGCACGATGCCTCACCGACGACTTTGGTGAtgc                                                                                           <  1:900260/101‑1 (MQ=255)
         gCATACTCATGTCCATTAGCACCACGTCAACGGCATTTGCCCGACACCACTTAACGGCGTCTTCACCGCACGATGCCTCACCGACGACTTTGGTGAtgctg                                                                                         <  1:929186/101‑1 (MQ=255)
           aTACTCATGTCCATTAGCACCACGTCAACGGCATTTGCCCGACACCACTTAACGGCGTCTTCACCGCACGATGCCTCACCGACGACTTTGGTGATGCTGcc                                                                                       >  2:360737/1‑101 (MQ=255)
            tACTCATGTCCATTAGCACCACGTCAACGGCATTTGCCCGACACCACTTAACGGCGTCTTCACCGCACGATGCCTCACCGACGACTTTGGTGATGCTGCCa                                                                                      >  1:1343764/1‑101 (MQ=255)
              cTCATGTCCATTAGCACCACGTCAACGGCATTTGCCCGACACCACTTACCGGCGTCTTCACCGCACGATGCCTCACCGACGACTTTGGTGATGCTGCCAAc                                                                                    <  2:505730/101‑1 (MQ=255)
              cTCATGTCCATTAGCACCACGTCAACGGCATTTGCCCGACACCACTTAACGGCGTCTTCACCGCACGATGCCTCACCGACGACTTTGGTGATGCTGCCAAc                                                                                    <  1:1052917/101‑1 (MQ=255)
                cATGTCCATTAGCACCACGTCAACGGCATTTGCCCGACACCACTTAACGGCGTCTTCACCGCACGATGCCTCACCGACGACTTTGGTGATGCTGCCAACtt                                                                                  <  2:237266/101‑1 (MQ=255)
                 aTGTCCATTAGCACCACGTCAACGGCATTTGCCCGACACCACTTAACGGCGTCTTCACCGCACGATGCCTCACCGACGACTTTGGTGATGCTGCCAACTTa                                                                                 >  2:502101/1‑101 (MQ=255)
                   gTCCATTAGCACCACGTCAACGGCATTTGCCCGACACCACTTAACGGCGTCTTCACCGCACGATGCCTCACCGACGACTTTGGTGATGCTGCCAACTTACt                                                                               <  1:1266360/101‑1 (MQ=255)
                    tCCATTATCGCCACGTCAACGGCATTTGCCCGACACCTCTGAACGGCGTCTTCACCGCACGATTCCTCACCGACGACTTTGGTGATGCTGCCAACTTGCTg                                                                              <  2:1536206/101‑1 (MQ=255)
                     ccATTAGCACCACGTCAACGGCATTTGCCCGACACCACTTAACGGCGTCTTCACCGCACGATGCCTCACCGACGACTTTGGTGATGCTGCCAACTTACTGa                                                                             >  1:338483/1‑101 (MQ=255)
                      cATTAGCACCACGTCAACGGCATTTGCCCGACACCACTTAACGGCGTCTTCACCGCACGATGCCTCACCGACGACTTTGGTGATGCTGCCAACTTACTGAt                                                                            <  1:60752/101‑1 (MQ=255)
                      cATTAGCACCACGTCAACGGCATTTGCCCGACACCACTTAACGGCGTCTTCACCGCACGATGCCTCACCGACGACTTTGGTGATGCTGCCAACTTACTGAt                                                                            >  2:1038034/1‑101 (MQ=255)
                        ttAGCACCACGTCAACGGCATTTGCCCGACACCACTTAACGGCATCTTCACCGCACGATGCCTCACCGACGACTTTGGTGATGCTGCCAACTTACTGAttt                                                                          <  2:425475/101‑1 (MQ=255)
                         tAGCACCACGTCAACGGCATTTGCCCGACACCACTTAACGGCGTCTTCACCGCACGATGCCTCACCGACGACTTTGGTGATGCTGCCAACTTACTGATTTa                                                                         <  1:884963/101‑1 (MQ=255)
                          aGCACCACGTCAACGGCATTTGCCCGACACCACTTAACGGCGTCTTCACCGCACGATGCCTCACCGACGACTTTGGTGATGCTGCCAACTTACTGATTTAg                                                                        <  1:457110/101‑1 (MQ=255)
                           gCACCACGTCAACGGCATTTGCCCGACACCACTTAACGGCGTCTTCACCGCACGATGCCTCACCGACGACTTTGGTGATGCTGCCAACTTACTGATTTAgt                                                                       <  2:341926/101‑1 (MQ=255)
                           gCACCACGTCAACGGCATTTGCCCGACACCACTTAACGGCGTCTTCACCGCACGATGCCTCACCGACGACTTTGGTGATGCTGCCAACTTACTGATTTAgt                                                                       <  2:1404381/101‑1 (MQ=255)
                            caccacGTCAACGGCATTTGCCCGACACCACTTAACGGCGTCTTCACCGCACGATGCCTCACCGACGACTTTGGTGATGCTGCCAACTTACTGATTTAgtg                                                                      >  1:1003163/1‑101 (MQ=255)
                             accacGTCAACGGCATTTGCCCGACACCACTTAACGGCGTCTTCACCGCACGATGCCTCACCGACGACTTTGGTGATGCTGCCAACTTACTGATTTAgtgt                                                                     >  1:312659/1‑101 (MQ=255)
                               cacGTCAACGGCATTTGCCCGACACCACTTAACGGCGTCTTCACCGCACGATGCCTCACCGACGACTTTGGTGATGCTGCCAACTTACTGATTTAGTGTat                                                                   >  2:229170/1‑101 (MQ=255)
                                acGTCAACGGCATTTGCCCGACACCACTTAACGGCGTCTTCACCGCACGATGCCTCACCGACGACTTTGGTGATGCTGCCAACTTACTGATTTAGTGTatg                                                                  >  1:1025547/1‑101 (MQ=255)
                                  gTCAACGGCATTTGCCCGACACCACTTAACGGCGTCTTCACCGCACGATGCCTCACCGACGACTTTGGTGATGCTGCCAACTTACTGATTTAGTGTatgat                                                                >  2:1368106/1‑101 (MQ=255)
                                  gTCAACGGCATTTGCCCGACACCACTTAACGGCGTCTTCACCGCACGATGCCTCACCGACGACTTTGGTGATGCTGCCAACTTACTGATTTAGTGTatgat                                                                <  1:617754/101‑1 (MQ=255)
                                    cAACGGCATTTGCCCGACACCACTTAACGGCGTCTTCACCGCACGATGCCTCACCGACGACTTTGGTGATGCTGCCAACTTACTGATTTAGTGTATGATgg                                                              <  2:633815/101‑1 (MQ=255)
                                      aCGGCATTTGCCCGACACCACTTAACGGCGTCTTCACCGCACGATGCCTCACCGACGACTTTGGTGATGCTGCCAACTTACTGATTTAGTGTATGATGgtg                                                            <  2:1592786/101‑1 (MQ=255)
                                           aTTTGCCCGACACCACTTAACGGCGTCTTCACCGCACGATGCCTCACCGACGACTTTGGTGATGCTGCCAACTTACTGATTTAGTGTATGATGGTGttttt                                                       >  1:712136/1‑101 (MQ=255)
                                           aTTTGCCCGACACCACTTAACGGCGTCTTCACCGCACGATGCCTCACCGACGACTTTGGTGATGCTGCCAACTTACTGATTTAGTGTATGATGGTGttttt                                                       <  1:252580/101‑1 (MQ=255)
                                           aTTTGCCCGACACCACTTAACGGCGTCTTCACCGCACGATGCCTCACCGACGACTTTGGTGATGCTGCCAACTTACTGATTTAGTGTATGATGGGGttttt                                                       >  1:357776/1‑101 (MQ=255)
                                           aTTTGCCCGACACCACTTAACGGCGTCTTCACCGCAAGATGCCTCACCGACGACTTTGGTGATGCTGCCAACTTACTGATTTAGTGTATGATGGTGttttt                                                       >  1:1697601/1‑101 (MQ=255)
                                            tttGCCCGACACCACTTAACGGCGTCTTCACCGCACGATGCCTCGCCGACGACTTTGGTGATGCTGCCAACTTACTGATTTAGTGTATGATGGTGTTTTTg                                                      <  2:1626125/101‑1 (MQ=255)
                                                 ccGACACCACTTAACGGCGTCTTCACCGCACGATGCCTCACCGACGACTTTGGTGATGCTGCCAACTTACTGATTTAGTGTATGATGGTGTTTTTGAGGTg                                                 <  1:1690759/101‑1 (MQ=255)
                                                   gACACCACTTAACGGCGTCTTCACCGCACGATGCCTCACCGACGACTTTGGTGATGCTGCCAACTTACTGATTTAGTGTATGATGGTGTTTTTGAGGTGCt                                               <  1:419967/101‑1 (MQ=255)
                                                    acacCACTTAACGGCGTCTTCACCGCACGATGCCTCACCGACGACTTTGGTGATGCTGCCAACTTACTGATTTAGTGTATGATGGTGTTTTTGAGGTGCTc                                              >  2:265021/1‑101 (MQ=255)
                                                      accacTTAACGGCGTCTTCACCGCACGATGCCTCACCGACGACTTTGGTGATGCTGCCAACTTACTGATTTAGTGTATGATGGTGTTTTTGAGGTGCTCCa                                            <  2:712131/101‑1 (MQ=255)
                                                      accacTTAACGGCGTCTTCACCGCACGATGCCTCACCGACGACTTTGGTGATGCTGCCAACTTACTGATTTAGTGTATGATGGTGTTTTTGAGGTGCTCCa                                            >  1:471541/1‑101 (MQ=255)
                                                          cTTAACGGCGTCTTCCCCGCACGATGCCTCACCGACGACTTTGGTGATGCTGCCAACTTACTGATTTAGTGTATGATGGTGTTTTTGAGGTGCTCCAGTgg                                        <  1:360737/101‑1 (MQ=255)
                                                            tAACGGCGTCTTCACCGCACGATGCCTCACCGACGACTTTGGTGATGCTGCCAACTTACTGATTTAGTGTATGATGGTGTTTTTGAGGTGCTCCAGTGGCt                                      <  1:846561/101‑1 (MQ=255)
                                                             aaCGGCGTCTTCACCGCACGATGCCTCACCGACGACTTTGGTGATGCTGCCAACTTACTGATTTAGTGTATGATGGTGTTTTTGAGGTGCTCCAGTGGCtt                                     >  1:1202203/1‑101 (MQ=255)
                                                                  cGTCTTCACCGCACGATGCCTCACCGACGACTTTGGTGATGCTGCCAACTTACTGATTTAGTGTATGATGGTGTTTTTGAGGTGCTCCAGTGGCTTCTGtt                                >  2:898253/1‑101 (MQ=255)
                                                                   gTCTTCACCGCACGATGCCTCACCGACGACTTTGGTGATGCTGCCAACTTACTGATTTAGTGTATGATGGTGTTTTTGAGGTGCTCCAGTGGCTTCTGttt                               <  1:959890/101‑1 (MQ=255)
                                                                   gTCTTCACCGCACGATGCCTCACCGACGACTTTGGTGATGCTGCCAACTTACTGATTTAGTGTATGATGGTGTTTTTGAGGTGCTCCAGTGGCTTCTGttt                               >  1:281695/1‑101 (MQ=255)
                                                                         aCCGCACGATGCCTCACCGACGACTTTGGTGATGCTGCCAACTTACTGATTTAGTGTATGATGGTGTTTTTGAGGTGCTCCAGTGGCTTCTGTTTCTATCa                         <  1:426757/101‑1 (MQ=255)
                                                                         aCCGCACGATGCCTCACCGACGACTTTGGTGATGCTGCCAACTTACTGATTTAGTGTATGATGGTGTTTTTGAGGTGCTCCAGTGGCTTCTGTTTCTATCa                         <  1:883872/101‑1 (MQ=255)
                                                                              aCGATGCCTCACCGACGACTTTGGTGATGCTGCCAACTTACTGATTTAGTGTATGATGGTGTTTTTGAGGTGCTCCAGTGGCTTCTGTTTCTATCAGCTGt                    <  2:309918/101‑1 (MQ=255)
                                                                                  tGCCTCACCGACGACTTTGGTGATGCTGCCAACTTACTGATTTAGTGTATGATGGTGTTTTTGAGGTTCTCCAGTGGCTTCTGTTTCTATCAGCTGTCcct                >  2:1509167/1‑101 (MQ=255)
                                                                                  tGCCTCACCGACGACTTTGGTGATGCTGCCAACTTACTGATTTAGTGTATGATGGTGTTTTTGAGGTGCTCCAGTGGCTTCTGTTTCTATCAGCTGTCcct                >  1:524847/1‑101 (MQ=255)
                                                                                  tGCCTCACCGACGACTTTGGTGATGCTGCCAACTTACTGATTTAGTGTATGATGGTGTTTTTGAGGTGCTCCAGTGGCTTCTGTTTCTATCAGCTGTCcct                >  1:140381/1‑101 (MQ=255)
                                                                                   gCCTCACCGACGACTTTGGTGATGCTGCCAACTTACTGATTTAGTGTATGATGGTGTTTTTGAGGTGCTCCAGTGGCTTCTGTTTCTATCAGCTGTCcctc               <  2:785599/101‑1 (MQ=255)
                                                                                     cTCACCGACGACTTTGGTGATGCTGCCAACTTACTGATTTAGTGTATGATGGTGTTTTTGAGGTGCTCCAGTGGCTTCTGTTTCTATCAGCTGTCcctcct             >  1:520454/1‑101 (MQ=255)
                                                                                     cTCACCGACGACTTTGGTGATGCTGCCAACTTACTGATTTAGTGTATGATGGTGTTTTTGAGGTGCTCCAGTGGCTTCTGTTTCTATCAGCTGTCcctcct             >  2:785454/1‑101 (MQ=255)
                                                                                     cTCACCGACGACTTTGGTGATGCTGCCAACTTACTGATTTAGTGTATGATGGTGTTTTTGAGGTGCTCCAGTGGCTTCTGTTTCTATCAGCTGTCcctcct             >  2:177365/1‑101 (MQ=255)
                                                                                     cTCACCGACGACTTTGGTGATGCTGCCAACTTACTGATTTAGTGTATGATGGTGTTTTTGAGGTGCTCCAGTGGCTTCTGTTTCTATCAGATGTCcctccc             >  2:14820/1‑100 (MQ=255)
                                                                                            acgacTTTGGTGATGCTGCCAACTTACTGATTTAGTGTATGATGGTGTTTTTGAGGTGCTCCAGTGGCTTCTGTTTCTATCAGCTGTCCCTCCTGTTCAGc      >  2:1184394/1‑101 (MQ=38)
                                                                                               acTTTGGTGATGCTGCCAACTTACTGATTTAGTGTATGATGGTGTTTTTGAGGTGCTCCAGTGGCTTCTGTTTCTATCAGCTGTCCCTCCTGTTCAGCTAc   >  2:49262/1‑101 (MQ=33)
                                                                                                cTTTGGTGATGCTGCCAACTTACTGATTTAGTGTATGATGGTGTTTTTGAGGTGCTCCAGTGGCTTCTGTTTCTATCAGCTGTCCCTCCTGTTCAGCTACt  <  2:312659/101‑1 (MQ=32)
                                                                                                    |                                                                                                
CAATGCCCGGCATACTCATGTCCATTAGCACCACGTCAACGGCATTTGCCCGACACCACTTAACGGCGTCTTCACCGCACGATGCCTCACCGACGACTTTGGTGATGCTGCCAACTTACTGATTTAGTGTATGATGGTGTTTTTGAGGTGCTCCAGTGGCTTCTGTTTCTATCAGCTGTCCCTCCTGTTCAGCTACT  >  REL606/1973172‑1973368

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