Predicted mutation
evidence seq id position mutation freq annotation gene description
JC JC REL606 1,764,698 IS1 (–) +9 bp 52.5% 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 48.6% noncoding (1/768 nt) IS1 repeat region
?REL606 1764698 = 19 (0.440)intergenic (‑318/‑15) ppsA/ydiA phosphoenolpyruvate synthase/hypothetical protein
* ? REL606 = 176470619 (0.470)24 (0.420) 22/66 NT 55.8% intergenic (‑326/‑7) ppsA/ydiA phosphoenolpyruvate synthase/hypothetical protein
?REL606 = 1974039 NA (NA)noncoding (768/768 nt) IS1 repeat region

TTTTAATTTATGCTTTCATAGAATTATGTCTGCATCACGGGAAGAACAAAATGGATAATGCTGTTGATCG  >  REL606/1764663‑1764732
                                   |                                  
ttttAATTTATGCTTTCATAGAATTATGTCTGCATc                                    >  1:5599018/1‑36 (MQ=255)
 tttAATTTATGCTTTCATAGAATTATGTCTGCATCa                                   <  1:1321385/36‑1 (MQ=255)
  ttAATTTATGCTTTCATAGAATTATGTCTGCATCAc                                  >  1:3536816/1‑36 (MQ=255)
   tAATTTATGCTTTCATAGAATTATGTCTGCATCACg                                 >  1:3120296/1‑36 (MQ=255)
   tAATTTATGCTTTCATAGAATTATGTCTGCATCACg                                 >  1:6257725/1‑36 (MQ=255)
    aaTTTATGCTTTCATAGAATTATGTCTGCATCACgg                                <  1:4212477/36‑1 (MQ=255)
    aaTTTATGCTTTCATAGAATTATGTCTGCATCACgg                                <  1:715239/36‑1 (MQ=255)
     aTTTATGCTTTCATAGAATTATGTCTGCATCACggg                               >  1:5972968/1‑36 (MQ=255)
      tttATGCTTTCATAGAATTATGTCTGCATCACGGga                              >  1:1489919/1‑36 (MQ=255)
      cactaaatcagtaagttggcagcatcaccCACGGga                              >  1:6112959‑M2/30‑36 (MQ=255)
       ttATGCTTTCATAGAATTATGTCTGCATCACGGGAA                             >  1:43580‑M1/1‑36 (MQ=255)
       actaaatcagtaagttggcagcatcaccCACGGgaa                             <  1:4784487‑M2/8‑1 (MQ=255)
        tATGCTTTCATAGAATTATGTCTGCATCACGGGAAG                            >  1:1469307‑M1/1‑36 (MQ=255)
        ctaaatcagtaagttggcagcatcaccCACGGgaag                            >  1:640076‑M2/28‑36 (MQ=255)
        ctaaatcagtaagttggcagcatcaccCACGGgaag                            <  1:790364‑M2/9‑1 (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:2062680/36‑1 (MQ=255)
                       ttATGTCTGCATCACGGGAAGAACAAAATGGATAAt             <  1:890041/36‑1 (MQ=255)
                       ggcagcatcaccCACGTGAAGAACAAAATGGATAAt             <  1:1920264‑M2/24‑1 (MQ=255)
                        tATGTCTGCATCACGGGAAGCACAAAATGGATAATg            >  1:5851087/1‑36 (MQ=255)
                        tATGTCTGCATCACGGGAAGAACAAAATGGATAATg            <  1:878644/36‑1 (MQ=255)
                        tATGTCTGCATCACGGGAAGAACAAAATGGATAATg            >  1:3622721/1‑36 (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:1276842‑M1/36‑21 (MQ=255)
                            tCTGCATCACGGGAAGggtaatgactccaacttatt        <  1:6777160‑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:7906835‑M1/36‑23 (MQ=255)
                              tGCATCACGGGAAGggtaatgactccaacttattga      <  1:2099107‑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:4043803‑M1/1‑10 (MQ=255)
                                  tCACGGGAAGggtaatgactccaacttattgatagt  <  1:1975244‑M1/36‑27 (MQ=255)
                                  cCACGGGAAGAACAAAATGGATAATGCTGTTGATCg  <  1:5046047‑M2/35‑1 (MQ=255)
                                   |                                  
TTTTAATTTATGCTTTCATAGAATTATGTCTGCATCACGGGAAGAACAAAATGGATAATGCTGTTGATCG  >  REL606/1764663‑1764732

Base quality scores:  ATCG  < 3 ≤  ATCG  < 20 ≤  ATCG  < 27 ≤  ATCG  < 32 ≤  ATCG  < 37 ≤  ATCG