Predicted mutation
evidence seq id position mutation freq annotation gene description
JC JC REL606 1,764,698 IS1 (–) +9 bp 44.7% intergenic (‑318/‑7) ppsA ← / → ydiA phosphoenolpyruvate synthase/hypothetical protein

New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? REL606 241257 =NA (NA)18 (0.310) 18/66 NT 40.9% noncoding (1/768 nt) IS1 repeat region
?REL606 1764698 = 26 (0.430)intergenic (‑318/‑15) ppsA/ydiA phosphoenolpyruvate synthase/hypothetical protein
* ? REL606 = 176470626 (0.450)24 (0.420) 22/66 NT 48.0% intergenic (‑326/‑7) ppsA/ydiA phosphoenolpyruvate synthase/hypothetical protein
?REL606 = 1974039 NA (NA)noncoding (768/768 nt) IS1 repeat region

TATGCTTTCATAGAATTATGTCTGCATCACGGGAAGAACAAAATGGATAATGCTGTTGATCGCCACGTTTT  >  REL606/1764671‑1764741
                                   |                                   
tATGCTTTCATAGAATTATGTCTGCATCACGGGAAG                                     >  1:1469307‑M1/1‑36 (MQ=255)
ctaaatcagtaagttggcagcatcaccCACGGgaag                                     <  1:790364‑M2/9‑1 (MQ=255)
ctaaatcagtaagttggcagcatcaccCACGGgaag                                     >  1:640076‑M2/28‑36 (MQ=255)
 taaatcagtaagttggcagcatcaccCACGGgaaga                                    <  1:4218586‑M2/10‑1 (MQ=255)
 aTGCTTTCATAGAATTATGTCTGCATCACGGGAAGg                                    >  1:4073314‑M1/1‑35 (MQ=255)
  tGCTTTCATAGAATTATGTCTGCATCACGGgaagaa                                   <  1:2101834/36‑1 (MQ=255)
  tGCTTTCATAGAATTATGTCTGCATCACGGgaagaa                                   >  1:1884559/1‑36 (MQ=255)
  tGCTTTCATAGAATTATGTCTGCATCACGGGAAGgg                                   >  1:3623005‑M1/1‑34 (MQ=255)
     tttCATAGAATTATGTCTGCATCACGGGAAGggtaa                                >  1:2119132‑M1/1‑31 (MQ=255)
     tcagtaagttggcagcatcaccCACGGGAAGAACaa                                <  1:2978246‑M2/14‑1 (MQ=255)
      ttCATAGAATTATGTCTGCATCACGGGAAGAACaaa                               <  1:6905712/36‑1 (MQ=255)
      ttCATAGAATTATGTCTGCATCACGGGAAGAACaaa                               >  1:5197567/1‑36 (MQ=255)
      cagtaagttggcagcatcaccCACGGGAAGAACaaa                               <  1:7909033‑M2/15‑1 (MQ=255)
         aTAGAATTATGTCTGCATCACGGGAAGggtaatgac                            <  1:5498435‑M1/36‑10 (MQ=255)
         aTAGAATTATGTCTGCATCACGGGAAGggtaatgac                            >  1:3094612‑M1/1‑27 (MQ=255)
         aTAGAATTATGTCTGCATCACGGGAAGAACAAAATg                            <  1:872350/36‑1 (MQ=255)
         aTAGAATTATGTCTGCATCACGGGAAGAACAAAATg                            >  1:7170421/1‑36 (MQ=255)
          tAGAATTATGTCTGCATCACGGGAAGggtaatgact                           >  1:3385321‑M1/1‑26 (MQ=255)
          tAGAATTATGTCTGCATCACGGGAAGAACAAAATgg                           <  1:6509469/36‑1 (MQ=255)
          aagttggcagcatcaccCACGGGAAGAACAAAATgg                           <  1:3474185‑M2/19‑1 (MQ=255)
           aGAATTATGTCTGCATCACGGGAAGggtaatgactc                          <  1:1618273‑M1/36‑12 (MQ=255)
            gttggcagcatcaccCACGGGAAGAACAAAATGGAt                         >  1:4523025‑M2/16‑36 (MQ=255)
              aTTATGTCTGCATCACGGGAGGAACAAAATGGATaa                       >  1:5818857/1‑36 (MQ=255)
               ttATGTCTGCATCACGGGAAGggtaatgactccaac                      >  1:2580310‑M1/1‑21 (MQ=255)
               ttATGTCTGCATCACGGGAAGAACAAAATGGATAAt                      <  1:890041/36‑1 (MQ=255)
               ttATGTCTGCATCACGGGAAGAACAAAATGGATAAt                      <  1:2062680/36‑1 (MQ=255)
               ggcagcatcaccCACGTGAAGAACAAAATGGATAAt                      <  1:1920264‑M2/24‑1 (MQ=255)
                tATGTCTGCATCACGGGAAGCACAAAATGGATAATg                     >  1:5851087/1‑36 (MQ=255)
                tATGTCTGCATCACGGGAAGAACAAAATGGATAATg                     >  1:3622721/1‑36 (MQ=255)
                tATGTCTGCATCACGGGAAGAACAAAATGGATAATg                     <  1:878644/36‑1 (MQ=255)
                gcagcatcaccCACGGGAAGAACAAAATGGATAATg                     >  1:6900385‑M2/12‑36 (MQ=255)
                 aTGTCTGCATCACGGGAAGggtaatgactccaactt                    >  1:2230274‑M1/1‑19 (MQ=255)
                 aTGTCTGCATCACGGGAAGggtaatgactccaactt                    <  1:7172494‑M1/36‑18 (MQ=255)
                  tGTCTGCATCACGGGAAGggtaatgactccaactta                   <  1:2267589‑M1/36‑19 (MQ=255)
                  tGTCTGCATCACGGGAAGAACAAAATGGATAATGCt                   <  1:5326800/36‑1 (MQ=255)
                  agcatcaccCACGGGAAGAACAAAATGGATAATGCt                   >  1:7169104‑M2/10‑36 (MQ=255)
                   gcatcaccCACGGGAAGAACAAAATGGATAATGCTg                  <  1:949364‑M2/28‑1 (MQ=255)
                   gcatcaccCACGGGAAGAACAAAATGGATAATGCTg                  >  1:6176105‑M2/9‑36 (MQ=255)
                   gTCTGCATCACGGGAAGggtaatgactccaacttat                  >  1:5833057‑M1/1‑17 (MQ=255)
                    tCTGCATCACGTGAAGAACAAAATGGATAATGCTGt                 <  1:364866/36‑1 (MQ=255)
                    tCTGCATCACGGGAAGggtaatgactccaacttatt                 <  1:6777160‑M1/36‑21 (MQ=255)
                    tCTGCATCACGGGAAGggtaatgactccaacttatt                 <  1:1276842‑M1/36‑21 (MQ=255)
                    catcaccCACGGGAAGAACAAAATGGATGATGCTGt                 <  1:7012737‑M2/29‑1 (MQ=255)
                     cTGCATCACGGGAAGAACAAAATGGATAATGCTGtt                >  1:4725504/1‑36 (MQ=255)
                     atcaccCACGGGAAGAACAAAATGGATAATGCTGtt                <  1:2073157‑M2/30‑1 (MQ=255)
                      tGCATCACGGGAAGggtaatgactccaacttattga               <  1:2099107‑M1/36‑23 (MQ=255)
                      tGCATCACGGGAAGggtaatgactccaacttattga               <  1:7906835‑M1/36‑23 (MQ=255)
                      tGCATCACGGGAAGAACAAAATGGATAATGCTGTTg               >  1:406901/1‑36 (MQ=255)
                      tGCATCACGGGAAGAACAAAATGGATAATGCTGTTg               <  1:559108/36‑1 (MQ=255)
                       gCATCACGGGAAGggtaatgactccaacttattgat              >  1:5501470‑M1/1‑13 (MQ=255)
                        cATCACGTGAAGggtaatgactccaacttattgata             <  1:2342855‑M1/36‑25 (MQ=255)
                        cATCACGGGAAGggtaatgactccaacttattgata             >  1:5297406‑M1/1‑12 (MQ=255)
                        accCACGGGAAGAACAAAATGGATAATGCTGTTGAt             <  1:7325812‑M2/33‑1 (MQ=255)
                         aTCACGGGAAGAACAAAATGGATAATGCTGTTGATc            <  1:5671749/36‑1 (MQ=255)
                          tCACGGGAAGggtaatgactccaacttattgatagt           <  1:1975244‑M1/36‑27 (MQ=255)
                          tCACGGGAAGggtaatgactccaacttattgatagt           >  1:4043803‑M1/1‑10 (MQ=255)
                          cCACGGGAAGAACAAAATGGATAATGCTGTTGATCg           <  1:5046047‑M2/35‑1 (MQ=255)
                            aCGGGAAGggtaatgactccaacttattgatagtgt         >  1:3577214‑M1/1‑8 (MQ=255)
                            ACGGGAAGAACAAAATGGATAATGCTGTTGATCGcc         >  1:3751856‑M2/1‑36 (MQ=255)
                             cGGGAAGAACAAAATGGATAATGCTGTTGATCGCCa        <  1:634756/36‑1 (MQ=255)
                              gggAAGggtaatgactccaacttattgatagtgttt       >  1:86446‑M1/1‑6 (MQ=255)
                              gggAAGAACAAAATGGATAATGCTGTTGATCGCCAc       >  1:1688009/1‑36 (MQ=255)
                              gggAAGAAAAAAAGGGAAAATGCTGTTGATCGCCAc       <  1:2561030/36‑1 (MQ=25)
                               ggAAGggtaatgactccaactgattgatagtgtttt      <  1:2382221‑M1/36‑32 (MQ=255)
                                gaagaaCAAAATGGATAATGCTGTTGATCGCCACGt     <  1:1992351/36‑1 (MQ=255)
                                gaagaaCAAAATGGATAATGCTGTTGATCGCCACGt     <  1:1604030/36‑1 (MQ=255)
                                gaagaaCAAAATGGATAATGCTGTTGATCGCCACGt     <  1:7399957/36‑1 (MQ=255)
                                  agaaCAAAATGGATAATGCTGTTGATCGCCACGttt   >  1:3522835/1‑36 (MQ=255)
                                  agaaCAAAATGGATAATGCTGTTGATCGCCACGttt   <  1:7783390/36‑1 (MQ=255)
                                   gggtaatgactccaacttattgatagtgttttatgt  >  1:1604643‑M1/1‑1 (MQ=255)
                                   gggtaatgactccaacttattgatagtgttttatgt  <  1:2458080‑M1/36‑36 (MQ=255)
                                   gggtaatgactccaacttattgatagtgttttatgt  <  1:7859628‑M1/36‑36 (MQ=255)
                                   gggtaatgactccaacttattgatagtgttttatgt  <  1:895780‑M1/36‑36 (MQ=255)
                                   gggtaatgactccaacttattgatagtgttttatgt  <  1:7787183‑M1/36‑36 (MQ=255)
                                   gggtaatgactccaacttattgatagtgttttatgt  <  1:1271475‑M1/36‑36 (MQ=255)
                                   gggtaatgactccaacttattgatagtgttttatgt  <  1:3008077‑M1/36‑36 (MQ=255)
                                   gggtaatgactccaacttattgatagtgttttatgt  <  1:3410304‑M1/36‑36 (MQ=255)
                                   gggtaatgactccaacttattgatagtgttttatgt  <  1:3734371‑M1/36‑36 (MQ=255)
                                   gggtaatgactccaacttattgatagtgttttatgt  >  1:4761762‑M1/1‑1 (MQ=255)
                                   gggtaatgactccaacttattgatagtgttttatgt  <  1:6619007‑M1/36‑36 (MQ=255)
                                   gggtaatgactccaacttattgatagtgttttatgt  >  1:4982619‑M1/1‑1 (MQ=255)
                                   gggtaatgactccaacttattgatagtgttttatgt  <  1:6223999‑M1/36‑36 (MQ=255)
                                   gggtaatgactccaacttattgatagtgttttatgt  <  1:516215‑M1/36‑36 (MQ=255)
                                   gggtaatgactccaacttattgatagtgttttatgt  <  1:593979‑M1/36‑36 (MQ=255)
                                   gggtaatgactccaacttattgatagtgttttatgt  <  1:588374‑M1/36‑36 (MQ=255)
                                   gggtaatgactccaacttattgatagtgttttatgt  >  1:5188539‑M1/1‑1 (MQ=255)
                                   gggtaatgactccaacttattgatagtgttttatgt  <  1:5722688‑M1/36‑36 (MQ=255)
                                   gaaCAAAATGGATAATGCTGTTGATCGCCACGtttt  >  1:7036106/1‑36 (MQ=255)
                                   |                                   
TATGCTTTCATAGAATTATGTCTGCATCACGGGAAGAACAAAATGGATAATGCTGTTGATCGCCACGTTTT  >  REL606/1764671‑1764741

Base quality scores:  ATCG  < 3 ≤  ATCG  < 18 ≤  ATCG  < 26 ≤  ATCG  < 32 ≤  ATCG  < 37 ≤  ATCG