Predicted mutation
evidence seq id position mutation freq annotation gene description
JC JC REL606 3,015,771 IS150 (–) +4 bp 98.3% coding (1192‑1195/1677 nt) ECB_02816 → KpsD protein

New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? REL606 588495 =NA (NA)62 (1.460) 36/66 NT 98.4% noncoding (1/1977 nt) IS150 repeat region
?REL606 3015771 = 1 (0.090)coding (1192/1677 nt) ECB_02816 KpsD protein
* ? REL606 = 30157741 (0.050)54 (1.280) 37/66 NT 98.2% coding (1195/1677 nt) ECB_02816 KpsD protein
?REL606 = 3652533 NA (NA)noncoding (1443/1443 nt) IS150 repeat region

Only 100 of 119 total aligned reads displayed.

ATCTCTCGCTGCAAAAGCTGGAGGAAGCCTCACTTTCTGCCCAGTCATCCACCAAAGAAGAAGCCAGCCT  >  REL606/3015739‑3015808
                                   |                                  
atttcatcgtccaacaaaatggttgcagtacaCttt                                    <  1:2633218‑M2/4‑1 (MQ=255)
atttcatcgtccaacaaaatgggtgcagtacaCttt                                    <  1:3674647‑M2/4‑1 (MQ=255)
 tttcatcgtccaacaaaatgggtgcagtacaCTTTc                                   <  1:4291068‑M2/5‑1 (MQ=255)
 tttcatcgtccaacaaaatgggtgcagtacaCTTTc                                   <  1:5659277‑M2/5‑1 (MQ=255)
 tctctcGCTGCAAAAGCTGGAGGAAGCCTCACTTTt                                   <  1:3248500‑M1/36‑2 (MQ=255)
 tctctcGCTGCAAAAGCTGGAGGAAGCCTCACTTTt                                   <  1:1689590‑M1/36‑2 (MQ=255)
  ttcatcgtccaacaaaatgggtgcagtacaCTTTCt                                  <  1:3384874‑M2/6‑1 (MQ=255)
  ttcatcgtccaacaaaatgggtgcagtacaCTTTCt                                  >  1:5668941‑M2/31‑36 (MQ=255)
  ttcatcgtccaacaaaatgggtgcagtacaCTTTCt                                  <  1:4939874‑M2/6‑1 (MQ=255)
  ctctcGCTGCAAAAGCTGGAGGAAGCCTCACTTTtg                                  >  1:2228618‑M1/1‑34 (MQ=255)
   tctcGCTGCAAAAGCTGGAGGAAGCCTCACTTTtgt                                 >  1:4847499‑M1/1‑33 (MQ=255)
   tctcGCTGCAAAAGCTGGAGGAAGCCTCACTTTtgt                                 <  1:1322040‑M1/36‑4 (MQ=255)
   tcatcgtccaacaaaatgggtgcagtccaCTTTCTg                                 <  1:4489903‑M2/7‑1 (MQ=255)
   tcatcgtccaacaaaatgggtgcagtacaCTTTCTg                                 >  1:4962357‑M2/30‑36 (MQ=255)
   tcatcgtccaacaaaatgggtgcagtacaCTTTCTg                                 <  1:3326452‑M2/7‑1 (MQ=255)
   tcatcgtccaacaaaatgggtgcagtacaCTTTCTg                                 <  1:2286918‑M2/7‑1 (MQ=255)
    ctcTCTGCAAAAGCTGGAGGAAGCCTCACTTTtgta                                <  1:5893108‑M1/36‑5 (MQ=255)
    ctcGCTGCAAAAGCTGGAGGAAGCCTCACTTTtgta                                >  1:4009279‑M1/1‑32 (MQ=255)
    catcgtccaacaaaatgggtgcagtacaCTTTCTGc                                <  1:1513001‑M2/8‑1 (MQ=255)
    catcgtccaacaaaatgggtgcagtacaCTTTCTGc                                >  1:3720376‑M2/29‑36 (MQ=255)
     tcGCTGCAAAAGCTGGAGGAAGCCTCACTTTtgtac                               <  1:1626819‑M1/36‑6 (MQ=255)
     atcgtccaacaaaatgggtgcagtacaCTTTCTGcc                               >  1:81350‑M2/28‑36 (MQ=255)
      tcgtccaacaaaatgggtgcagtccaCTTTCTGccc                              <  1:3234553‑M2/10‑1 (MQ=255)
      tcgtccaacaaaatgggtgcagtacaCTTTCTGccc                              <  1:5408728‑M2/10‑1 (MQ=255)
      cGCTGCAAAAGCTGGAGGAAGCCTCACTTTtgtact                              <  1:4478402‑M1/36‑7 (MQ=255)
       gCTGCAAAAGCTGGAGGAAGCCTCACTTTtgtactg                             <  1:5087770‑M1/36‑8 (MQ=255)
       gCTGCAAAAGCTGGAGGAAGCCTCACTTTtgtactg                             >  1:775927‑M1/1‑29 (MQ=255)
       gCTGCAAAAGCTGGAGGAAGCCTCACTTTtgtactg                             >  1:3915431‑M1/1‑29 (MQ=255)
       cgtccaacaaaatgggtgcagtacaCTTTCTGCCCa                             <  1:538007‑M2/11‑1 (MQ=255)
        cTGCAAAAGCTGGAGGAAGCCTCACTTTtgtactga                            >  1:473191‑M1/1‑28 (MQ=255)
        cTGCAAAAGCTGGAGGAAGCCTCACTTTtgtactga                            <  1:802282‑M1/36‑9 (MQ=255)
         tccaacaaaatgggtgcagtacaCTTTCTGCCCAGt                           >  1:1969056‑M2/24‑36 (MQ=255)
         tGCAAAAGCTGGAGGAAGCCTCACTTTtgtactgac                           <  1:5210242‑M1/36‑10 (MQ=255)
         tGCAAAAGCTGGAGGAAGCCTCACTTTtgtactgac                           <  1:4481969‑M1/36‑10 (MQ=255)
          gCAAAAGCTGGAGGAAGCCTCACTTTtgtactgacc                          >  1:2313750‑M1/1‑26 (MQ=255)
          ccaacaaaattggtgcagtacaCTTTCTGCCCAGTc                          <  1:4802819‑M2/14‑1 (MQ=255)
          ccaacaaaatgggtgcagtacgCTTTCTGCCCAGTc                          <  1:5227793‑M2/14‑1 (MQ=255)
          ccaacaaaatgggtgcagtacaCTTTCTGCCCAGTc                          <  1:2431871‑M2/14‑1 (MQ=255)
            aacaaaatgggtgcagtacaCTTTCTGCCCAGTCAt                        <  1:443466‑M2/16‑1 (MQ=255)
            aacaaaatgggtgcagtacaCTTTCTGCCCAGTCAt                        <  1:4639783‑M2/16‑1 (MQ=255)
            aacaaaatgggtgcagtacaCTTTCTGCCCAGTCAt                        <  1:5116024‑M2/16‑1 (MQ=255)
            aaaaGCTGGAGGAAGCCTCACTTTtgtactgacccc                        >  1:3575707‑M1/1‑24 (MQ=255)
            aaaaGCTGGAGGAAGCCTCACTTTtgtactgacccc                        >  1:2733800‑M1/1‑24 (MQ=255)
            aaaaGCTGGAGGAAGCCTCACTTTtgtactgacccc                        <  1:3938734‑M1/36‑13 (MQ=255)
             acaaaatgggtgcagtacaCTTTCTGCCCAGTCATc                       >  1:605637‑M2/20‑36 (MQ=255)
             acaaaatgggtgcagtacaCTTTCTGCCCAGTCATc                       >  1:1266429‑M2/20‑36 (MQ=255)
             aaaGCTGGAGGAAGCCTCACTTTtgtactgacccca                       <  1:3515120‑M1/36‑14 (MQ=255)
              caaaatgggtgcagtacaCTTTCTGCCCAGTCATcc                      >  1:2865298‑M2/19‑36 (MQ=255)
              caaaatgggtgcagtacaCTTTCTGCCCAGTCATcc                      <  1:5858404‑M2/18‑1 (MQ=255)
              caaaatgagtgcagtacaCTTTCTGCCCAGTCATcc                      <  1:126203‑M2/18‑1 (MQ=255)
              aaGCTGGAGGAAGCCTCACTTTtgtactgaccccta                      >  1:893765‑M1/1‑22 (MQ=255)
              aaGCTGGAGGAAGCCTCACTTTtgtactgaccccaa                      >  1:1272556‑M1/1‑22 (MQ=255)
              aaGCTGGAGGAAGCCTCACTTTtgtactgaccccaa                      >  1:557101‑M1/1‑22 (MQ=255)
               aaaatgggtgcagtacaCTTTCTGCCCAGTCATcca                     >  1:2209821‑M2/18‑36 (MQ=255)
               aaaatgggtgcagtacaCTTTCTGCCCAGTCATcca                     <  1:4604750‑M2/19‑1 (MQ=255)
               aGCTGGAGGAAGCCTCACTTTtgtactgaccccaaa                     >  1:555464‑M1/1‑21 (MQ=255)
               aGCTGGAGGAAGCCTCACTTTtgtactgaccccaaa                     <  1:5243616‑M1/36‑16 (MQ=255)
                gCTGGAGGAAGCCTCACTTTtgtactgaccccaaaa                    >  1:1498344‑M1/1‑20 (MQ=255)
                gCTGGAGGAAGCCTCACTTTtgtactgaccccaaaa                    <  1:3514421‑M1/36‑17 (MQ=255)
                aaatgggtgcagtacaCTTTCTGCCCAGTCATccac                    <  1:3030726‑M2/20‑1 (MQ=255)
                aaatgggtgcagtacaCTTTCTGCCCAGTCATccac                    <  1:3980988‑M2/20‑1 (MQ=255)
                  atgggtgcagtacaCTTTCTGCCCAGTCATccacca                  >  1:1935402‑M2/15‑36 (MQ=255)
                  atgggtgcagtacaCTTTCTGCCCAGTCATccacca                  >  1:629254‑M2/15‑36 (MQ=255)
                  atgggtgcagtacaCTTTCTGCCCAGTCATccacca                  >  1:3745824‑M2/15‑36 (MQ=255)
                   tgggtgcagtacaCTTTCTGCCCAGTCATCCACCaa                 <  1:1777862‑M2/23‑1 (MQ=255)
                   tgggtgcagtacaCTTTCTGCCCAGTCATCCACCaa                 <  1:5558911‑M2/23‑1 (MQ=255)
                   ggaggaAGCCTCACTTTtgtactgaccccaaaaagt                 >  1:681979‑M1/1‑17 (MQ=255)
                   ggaggaAGCCTCACTTTtgtactgaccccaaaaagt                 >  1:4773693‑M1/1‑17 (MQ=255)
                    ggggaagCCTCACTTTtgtactgaccccaaaaagtt                >  1:4188176‑M1/3‑16 (MQ=255)
                    gaggaAGCCTCACTTTtgtactgaccccaaaaattt                >  1:775796‑M1/1‑16 (MQ=255)
                    gaggaAGCCTCACTTTtgtactgaccccaaaaagtt                <  1:4892485‑M1/36‑21 (MQ=255)
                    gaggaAGCCTCACTTTtgtactgaccccaaaaagtt                >  1:1309324‑M1/1‑16 (MQ=255)
                     ggtgcagtacaCTTTCTGCCCAGTCATCCACCAaag               <  1:5090981‑M2/25‑1 (MQ=255)
                      gtgcagtacaCTTTCTGCCCAGTCATCCACCAaaga              <  1:3118425‑M2/26‑1 (MQ=255)
                      gtgcagtacaCTTTCTGCCCAGCCATCCACCAaaga              >  1:46851‑M2/11‑36 (MQ=255)
                      ggaAGCCTCACTTTtgtactgaccccaaaaagttgg              <  1:5473186‑M1/36‑23 (MQ=255)
                       tgcagtacaCTTTCTGCCCAGTCATCCACCAaagaa             <  1:5304408‑M2/27‑1 (MQ=255)
                       gaAGCCTCACTTTtgtactgaccccaagaagttgga             <  1:2876506‑M1/36‑24 (MQ=255)
                        gcagtacaCTTTCTGCCCAGGCATCCACCAaagaag            >  1:283614‑M2/9‑36 (MQ=255)
                        aaGCCTCACTTTtgtactgaccccaaaaagttggac            >  1:3382003‑M1/1‑12 (MQ=255)
                         cagtacaCTTTCTGCCCAGTCATCCACcaaagaaag           >  1:5359568‑M2/8‑34 (MQ=255)
                         aGCCTCACTTTtgtactgaccccaaaaagttggaca           >  1:4551445‑M1/1‑11 (MQ=255)
                         aGCCTCACTTTtgtactgaccccaaaaagttggaca           <  1:4970273‑M1/36‑26 (MQ=255)
                          agtacaCTTTCTGCCCAGTCATCCACCAaagaagaa          >  1:3762210‑M2/7‑36 (MQ=255)
                           gtacaCTTTCTGCCCAGTCATCCACCAaagaagaag         >  1:2890612‑M2/6‑36 (MQ=255)
                           ccTCACTTTtgtactgaccccaaaaagttggacagt         <  1:4060139‑M1/36‑28 (MQ=255)
                            cTCAGTTTtgtactgaccccaaaaagttggacagtt        >  1:3654675‑M1/1‑8 (MQ=255)
                            cTCACTTTtgtactgaccccaaaaagttgtacagtt        <  1:2595381‑M1/36‑29 (MQ=255)
                            cTCACTTTtgtactgaccccaaaaagttggacagtt        >  1:2476383‑M1/1‑8 (MQ=255)
                             acaCTTTCTGCCCAGTCATCCACCAAAGAAGAAGcc       >  1:1935244‑M2/4‑36 (MQ=255)
                              caCTTTCTGCCCAGTCATCCACCAAAGAAGAAGCCa      >  1:1487903‑M2/3‑36 (MQ=255)
                              caCTTTCTGCCCAGTCATCCACAAAAGAAGAAGCCa      >  1:4269836‑M2/3‑36 (MQ=255)
                              cACTTTCTGCAAAGTAATCCACCAAAGAAGAagcca      <  1:359731/36‑1 (MQ=25)
                               aCTTTtgtactgaccccaaaaagttggacagttaaa     <  1:410310‑M1/36‑32 (MQ=255)
                               aCTTTtgtactgaccccaaaaagttggacagttaaa     <  1:1847956‑M1/36‑32 (MQ=255)
                               aCTTTCTGCCCAGTCATCAACCAAAGAAGAAGCCAg     >  1:2688501‑M2/2‑36 (MQ=255)
                                cTTTtgtactgaccccaaaaagttggacagttaaac    >  1:2465218‑M1/1‑4 (MQ=255)
                                cTTTtgtactgaccccaaaaagttggacagttaaac    <  1:4279700‑M1/36‑33 (MQ=255)
                                cTTTtgtactgaccccaaaaagttggacagttaaac    <  1:1389932‑M1/36‑33 (MQ=255)
                                  tttgtactgaccccaaaaagttggacagttaaacac  >  1:2987721‑M1/1‑2 (MQ=255)
                                   |                                  
ATCTCTCGCTGCAAAAGCTGGAGGAAGCCTCACTTTCTGCCCAGTCATCCACCAAAGAAGAAGCCAGCCT  >  REL606/3015739‑3015808

Base quality scores:  ATCG  < 3 ≤  ATCG  < 14 ≤  ATCG  < 22 ≤  ATCG  < 30 ≤  ATCG  < 37 ≤  ATCG