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

Only 100 of 121 total aligned reads displayed.

ACATAAAACACTATCAATAAGTTGGAGTCATTACCCTCGACTGGTAACCATCCGATTGAAACGTTGCTGC  >  REL606/1974005‑1974074
                                  |                                   
aCATAAAACACTATCAATAAGTTGGAGTCATTACCc                                    <  1:2458080‑M2/36‑2 (MQ=255)
aCATAAAACACTATCAATAAGTTGGAGTCATTACCc                                    <  1:3008077‑M2/36‑2 (MQ=255)
aCATAAAACACTATCAATAAGTTGGAGTCATTACCc                                    <  1:3410304‑M2/36‑2 (MQ=255)
aCATAAAACACTATCAATAAGTTGGAGTCATTACCc                                    <  1:3734371‑M2/36‑2 (MQ=255)
aCATAAAACACTATCAATAAGTTGGAGTCATTACCc                                    >  1:4761762‑M2/1‑35 (MQ=255)
aCATAAAACACTATCAATAAGTTGGAGTCATTACCc                                    >  1:4982619‑M2/1‑35 (MQ=255)
aCATAAAACACTATCAATAAGTTGGAGTCATTACCc                                    <  1:516215‑M2/36‑2 (MQ=255)
aCATAAAACACTATCAATAAGTTGGAGTCATTACCc                                    >  1:5188539‑M2/1‑35 (MQ=255)
aCATAAAACACTATCAATAAGTTGGAGTCATTACCc                                    <  1:5722688‑M2/36‑2 (MQ=255)
aCATAAAACACTATCAATAAGTTGGAGTCATTACCc                                    >  1:1604643‑M2/1‑35 (MQ=255)
aCATAAAACACTATCAATAAGTTGGAGTCATTACCc                                    <  1:588374‑M2/36‑2 (MQ=255)
aCATAAAACACTATCAATAAGTTGGAGTCATTACCc                                    <  1:593979‑M2/36‑2 (MQ=255)
aCATAAAACACTATCAATAAGTTGGAGTCATTACCc                                    <  1:6223999‑M2/36‑2 (MQ=255)
aCATAAAACACTATCAATAAGTTGGAGTCATTACCc                                    <  1:7859628‑M2/36‑2 (MQ=255)
aCATAAAACACTATCAATAAGTTGGAGTCATTACCc                                    <  1:1271475‑M2/36‑2 (MQ=255)
aCATAAAACACTATCAATAAGTTGGAGTCATTACCc                                    <  1:895780‑M2/36‑2 (MQ=255)
 cATAAAACACTATCAATTAGCTGGAGTCATTTCCCt                                   >  1:3404107/1‑36 (MQ=11)
   tAAAACACTATCAATAAGTTGGAGTCATTACCCTCg                                 >  1:850847/1‑36 (MQ=38)
    aaaaCACTATCAATCAGTTGGAGTCATTACCcttcc                                <  1:2382221‑M2/36‑6 (MQ=255)
    aaaaCACTATCAATAAGTTGGAGTCATTACCCTCGa                                >  1:961562/1‑36 (MQ=40)
    aaaaCAATATCAATAAGTTGGAGACATCACCCTcca                                >  1:5896671/1‑34 (MQ=11)
    aaaaAACTATCAATAAGTTGGAGTCATTACCCTCGa                                >  1:5987234/1‑36 (MQ=255)
     aaaCACTATCAATAAGTTGGAGTCATTACCcttccc                               >  1:86446‑M2/1‑30 (MQ=255)
     aaaCACTATCAATAAGTTGGAGTCATTACCCTCGAc                               >  1:2457135/1‑36 (MQ=255)
     aaaCACTATCAATAAGTTGGAGTCATTACCCTCGAc                               <  1:614480/36‑1 (MQ=255)
      aaCACTATCAATAAGTTGGAGTCATTACCCTCGACt                              <  1:1433345/36‑1 (MQ=255)
       acacTATCAATAAGTTGGAGTCATTACCcttcccgt                             >  1:3577214‑M2/1‑28 (MQ=255)
       acacTATCAATAAGTTGGAGTCATTACCCTCGACTg                             <  1:1648880/36‑1 (MQ=255)
       acacTATCAATAAGTTGGAGTCATTACCCTCGACTg                             >  1:4774118/1‑36 (MQ=255)
       acacTATCAATAAGTTGGAGTCATTACCCTCGACTg                             >  1:6137225/1‑36 (MQ=255)
        cacTATCAATAAGTTGGAGTCATTACCCTCGACTgg                            >  1:60503/1‑36 (MQ=255)
         acTATCAATAAGTTGGAGTCATTACCcttcccgtga                           >  1:4043803‑M2/1‑26 (MQ=255)
         acTATCAATAAGTTGGAGTCATTACCcttcccgtga                           <  1:1975244‑M2/36‑11 (MQ=255)
         acTATCAATAAGTTGGAGTCATTACCCTCGACTGGt                           <  1:3656574/36‑1 (MQ=255)
          cTATCAATAAGTTGGAGTCATTACCCTCGACTGGTa                          <  1:704100/36‑1 (MQ=255)
          cTATCAATAAGTTGGAGTCATTACCCTCGACTGGTa                          >  1:5632391/1‑36 (MQ=255)
          cTAACAATAAGTTGGAGTCAATACCCTCGACTGGTa                          <  1:7829981/36‑1 (MQ=21)
           tATCAATAAGTTGGAGTCATTACCcttcacgtgatg                         <  1:2342855‑M2/36‑13 (MQ=255)
           tATCAATAAGTTGGAGTCATTACCCTCGACTGGTaa                         <  1:7830058/36‑1 (MQ=255)
           tATCAATAAGTTGGAGTCATTACCCTCGACTGGTaa                         <  1:1704051/36‑1 (MQ=255)
           tATCAATAAGTTGGAGTCATTACCCTCGACTGGTaa                         >  1:335273/1‑36 (MQ=255)
            aTCAATAAGTTGGAGTCATTACCcttcccgtgatgc                        >  1:5501470‑M2/1‑23 (MQ=255)
            aTCAATAAGTTGGAGTCATTACCCTCGACTGGTAAc                        <  1:7040052/36‑1 (MQ=255)
            aTCAATAAGTTGGAGTCATTACCCTCGACTGGTAAc                        <  1:2554140/36‑1 (MQ=255)
            aTCAATAAGTTGGAATCATTACCCTCGACTGGTAAc                        <  1:2236947/36‑1 (MQ=255)
             tCAATAAGTTGGAGTCATTACCcttcccgtgatgca                       <  1:2099107‑M2/36‑15 (MQ=255)
             tCAATAAGTTGGAGTCATTACCcttcccgtgatgca                       <  1:7906835‑M2/36‑15 (MQ=255)
              cAATAAGTTGGAGTCATTACCCTCGACTGGTAACCa                      >  1:5300466/1‑36 (MQ=255)
               aaTAATTTGGAGTCATTACCCTCGACTGGTAACCAt                     <  1:7183765/36‑1 (MQ=255)
               aaTAAGTTGGAGTCATTACCcttcccgtgatgcaga                     <  1:1276842‑M2/36‑17 (MQ=255)
               aaTAAGTTGGAGTCATTACCcttcccgtgatgcaga                     <  1:6777160‑M2/36‑17 (MQ=255)
               aaTAAGTTGGAGTCATTACCCTCGACTGGTAACCAt                     >  1:7559446/1‑36 (MQ=255)
                aTAAGTTGGAGTCATTACCcttcccgtgatgcagac                    >  1:5833057‑M2/1‑19 (MQ=255)
                aTAAGTTGGAGTCATTACCCTCGACTGGTAACCATc                    >  1:5989380/1‑36 (MQ=255)
                 tAAGTTGGAGTCATTACCcttcccgtgatgcagaca                   <  1:2267589‑M2/36‑19 (MQ=255)
                 tAAGTTGGAGTCATTACCCTCGACTGGTAACCATcc                   <  1:530814/36‑1 (MQ=255)
                 tAAGTTGGAGTCATTACCCTCGACTGGTAACCATcc                   <  1:7805453/36‑1 (MQ=255)
                  aaGTTGGAGTCATTACCcttcccgtgatgcagacat                  <  1:7172494‑M2/36‑20 (MQ=255)
                  aaGTTGGAGTCATTACCcttcccgtgatgcagacat                  >  1:2230274‑M2/1‑17 (MQ=255)
                  aaGTTGGAGTCATTACCCTCGACTGGTAACCATCCg                  <  1:3306889/36‑1 (MQ=255)
                   aGTTGGAGTCATTACCCTCGACTGGTAACCATCCGa                 >  1:3427696/1‑36 (MQ=255)
                   aGTTGGAGTCATTACCCTCGACTGGTAACCATCCGa                 >  1:1442381/1‑36 (MQ=255)
                   aGTTGGAGTCATTACCCTCGACTGGTAACCATCCGa                 >  1:3111744/1‑36 (MQ=255)
                    gtgggAGTCATTACCCTCGACTGGTAACCATCCGAt                <  1:1056114/33‑1 (MQ=255)
                    gTTTGAGTCATTACCCTCGACTGGTAACCATCCGAt                >  1:2636065/1‑36 (MQ=255)
                    gTTGGAGTCATTACCcttcccgtgatgcagacataa                >  1:2580310‑M2/1‑15 (MQ=255)
                    gTTGGAGTCATTACCCTCGACTGGTAACCATCCGAt                <  1:7633556/36‑1 (MQ=255)
                    gTTGGAGTCATTACCCTCGACTGGTAACCATCCGAt                >  1:6157623/1‑36 (MQ=255)
                    gTTGGAGTCATTACCCTCGACTGGTAACCATCCGAt                >  1:3813737/1‑36 (MQ=255)
                     ttGGAGTCATTACCCTCGACTGGTACCCATCCGAtt               <  1:1056188/36‑1 (MQ=255)
                     ttGGAGTCATTACCCTCGACTGGTAACCATCCGAtt               >  1:413698/1‑36 (MQ=255)
                     ttGGAGTCATTACCCTCGACTGGTAACCATCCGAtt               >  1:484360/1‑36 (MQ=255)
                      tGGAGTCATTACCCTCGACTGGTAACCATCCGATTg              >  1:1478658/1‑36 (MQ=255)
                      tGGAGTCATTACCCTCGACTGGTAACCATCCGATTg              >  1:6438069/1‑36 (MQ=255)
                       ggAGTCATTACCCTCGACTGGTAACCATCCGATTGa             <  1:3783439/36‑1 (MQ=255)
                        gAGTCATTACCcttcccgtgatgcagacataattct            <  1:1618273‑M2/36‑26 (MQ=255)
                        gAGTCATTACCCTCGACTGGTAACCATCCGATTGaa            >  1:6332700/1‑36 (MQ=255)
                         aGTCATTACCcttcccgtgatgcagacataattcta           >  1:3385321‑M2/1‑10 (MQ=255)
                         aGTCATTACCCTCGACTGGTAACCATCCGATTGaaa           <  1:1011139/36‑1 (MQ=255)
                         aGTCATTACCCTCGACTGGTAAACATCAGATTGaaa           >  1:6791684/1‑36 (MQ=37)
                          gTCATTACCcttcccgtgatgcagacataattctat          >  1:3094612‑M2/1‑9 (MQ=255)
                          gTCATTACCcttcccgtgatgcagacataattctat          <  1:5498435‑M2/36‑28 (MQ=255)
                          gTCATTACCCTCGACTGGTACCCATCCGATTGAAAc          <  1:7573862/36‑1 (MQ=255)
                          gTCATTACCCTCGACTGGTAACCATCCGATTGAAAc          <  1:6555701/36‑1 (MQ=255)
                          gTCATTACCCTCGACTGGTAACCATCCGATTGAAAc          >  1:6236573/1‑36 (MQ=255)
                           tCATTACCCTCGACTGGTAACCATCCGATTGAAACg         <  1:5038944/36‑1 (MQ=255)
                           tCATTACCCTCGACTGGTAACCATCCGATTGAAACg         >  1:4540118/1‑36 (MQ=255)
                            cATTACCCTCGACTGGTAACCATCCGATTGAAACGt        <  1:4084524/36‑1 (MQ=255)
                            cATTACCCTCGACTGGTAACCATCCGATTGAAACGt        <  1:1915552/36‑1 (MQ=255)
                              ttACCcttcccgtgatgcagacataattctatgaaa      >  1:2119132‑M2/1‑5 (MQ=255)
                              ttACCCTCGACTGGTAACCATCCGATTGAAACGTtg      >  1:2496224/1‑36 (MQ=255)
                               tACCCTCTACTTGTAACCATCCGATTGAAACGTtgc     <  1:2560058/36‑1 (MQ=38)
                               tACCCTCGACTGGTAACCATCCGATTGAAACGTtgc     <  1:5254580/36‑1 (MQ=255)
                               tACCCTCGACTGGTAACCATCCGATTGAAACGTtgc     >  1:6506675/1‑36 (MQ=255)
                                aCCCTCGACTGGTAACCATCCTATTGAAACGTtgct    <  1:1926303/36‑1 (MQ=255)
                                aCCCTCGACTGGTAACCATCCGATTGAAACGTtgct    <  1:593443/36‑1 (MQ=255)
                                 cccttcccgtgatgcagacataattctatgaaagca   >  1:3623005‑M2/1‑2 (MQ=255)
                                 cccTCGACTGGTAACCATCCGATTGAAACGTtgctg   >  1:7006612/1‑36 (MQ=255)
                                  ccttcccgtgatgcagacataattctatgaaagcat  >  1:4073314‑M2/1‑1 (MQ=255)
                                  ccTCGACTGGTAACCATCCGATTGAAACGTtgctgc  <  1:2953604/36‑1 (MQ=255)
                                  |                                   
ACATAAAACACTATCAATAAGTTGGAGTCATTACCCTCGACTGGTAACCATCCGATTGAAACGTTGCTGC  >  REL606/1974005‑1974074

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