Predicted mutation
evidence seq id position mutation freq annotation gene description
JC JC REL606 3,015,771 IS150 (–) +4 bp 96.7% 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 96.9% noncoding (1/1977 nt) IS150 repeat region
?REL606 3015771 = 2 (0.040)coding (1192/1677 nt) ECB_02816 KpsD protein
* ? REL606 = 30157742 (0.050)54 (1.280) 37/66 NT 96.4% coding (1195/1677 nt) ECB_02816 KpsD protein
?REL606 = 3652533 NA (NA)noncoding (1443/1443 nt) IS150 repeat region

Only 100 of 126 total aligned reads displayed.

TGAATCTCTCGCTGCAAAAGCTGGAGGAAGCCTCACTTTCTGCCCAGTCATCCACCAAAGAAGAAGCCAGC  >  REL606/3015736‑3015806
                                   |                                   
tccatttcatcgtccaacaaaatgggtgcagtacac                                     >  1:32994‑M2/36‑36 (MQ=255)
tccatttcatcgtccaacaaaatgggtgcagtacac                                     >  1:5302910‑M2/36‑36 (MQ=255)
tccatttcatcgtccaacaaaatgggtgcagtacac                                     >  1:5575559‑M2/36‑36 (MQ=255)
 gAATCTCTCGCTGCAAAAGCTGGAGGAAGCCTCACt                                    <  1:4216820/36‑1 (MQ=255)
 gAATCTCTCGCTGCAAAAGCTGGAGGAAGCCTCACt                                    >  1:1635166/1‑36 (MQ=255)
  catttcatcgtccaacaaaatgggtgcagtacaccg                                   <  1:4089607‑M2/3‑3 (MQ=255)
  catttcatcgtccaacaaaatgggtgcagtacaccg                                   <  1:5030201‑M2/3‑3 (MQ=255)
  catttcatcgtccaacaaaatgggtgcagtacaCtt                                   <  1:5001079‑M2/3‑1 (MQ=255)
  aaTCTCTCGCTGTACAAGCTGGAGGAAGCCTCACtt                                   <  1:360969/36‑1 (MQ=37)
  aaTCTCTCGCTGCAAAAGCTGGAGGAAGCCTCACTT                                   <  1:5329493‑M1/36‑1 (MQ=255)
   atttcatcgtccaacaaaatgggtgcagtacaCttt                                  >  1:505086‑M2/33‑36 (MQ=255)
    tttcatcgtccaacaaaatgggtgcagtacaCTTTc                                 <  1:4291068‑M2/5‑1 (MQ=255)
    tctctcGCTGCAAAAGCTGGAGGAAGCCTCACTTTt                                 <  1:1689590‑M1/36‑2 (MQ=255)
    tctctcGCTGCAAAAGCTGGAGGAAGCCTCACTTTt                                 <  1:5524381‑M1/36‑2 (MQ=255)
    tctctcGCTGCAAAAGCTGGAGGAAGCCTCACTTTt                                 <  1:3248500‑M1/36‑2 (MQ=255)
     ttcatcgtccaacaaaatgggtgcagtacaCTTTCt                                <  1:4939874‑M2/6‑1 (MQ=255)
     ttcatcgtccaacaaaatgggtgcagtacaCTTTCt                                <  1:3384874‑M2/6‑1 (MQ=255)
     ttcatcgtccaacaaaatgggtgcagtacaCTTTCt                                >  1:5668941‑M2/31‑36 (MQ=255)
     ctctcGCTGCAAAAGCTGGAGGAAGCCTCACTTTtg                                >  1:2228618‑M1/1‑34 (MQ=255)
     ctctcGCTGCAAAAGCTGGAGGAAGCCTCACTTTtg                                >  1:5783680‑M1/1‑34 (MQ=255)
      tctcGCTGCAAAAGCTGGAGGAAGCCTCACTTTtgt                               >  1:4847499‑M1/1‑33 (MQ=255)
      tcatcgtccaacaaaatgggtgcagtccaCTTTCTg                               <  1:4489903‑M2/7‑1 (MQ=255)
      tcatcgtccaacaaaatgggtgcagtacaCTTTCTg                               <  1:2286918‑M2/7‑1 (MQ=255)
      tcatcgtccaacaaaatgggtgcagtacaCTTTCTg                               <  1:3326452‑M2/7‑1 (MQ=255)
      tcatcgtccaacaaaatgggtgcagtacaCTTTCTg                               >  1:4962357‑M2/30‑36 (MQ=255)
       ctcTCTGCAAAAGCTGGAGGAAGCCTCACTTTtgta                              <  1:5893108‑M1/36‑5 (MQ=255)
       ctcGCTGCAAAAGCTGGAGGAAGCCTCACTTTtgta                              >  1:4009279‑M1/1‑32 (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)
          gCTGCAAAAGCTGGAGGAAGCCTCACTTTtgtactg                           <  1:5087770‑M1/36‑8 (MQ=255)
          cgtccaacaaaatgggtgcagtacaCTTTCTGCCCa                           <  1:538007‑M2/11‑1 (MQ=255)
          cgtccaacaaaatgggtgcagtacaCTTTCTGCCCa                           >  1:4324547‑M2/26‑36 (MQ=255)
          cgtccaacaaaatgggtgcagtacaCTTTCTGCCCa                           <  1:361038‑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)
             ccaacaaaatgggtgcagtacgCTTTCTGCCCAGTc                        <  1:5227793‑M2/14‑1 (MQ=255)
               aacaaaatgggtgcagtacaCTTTCTGCCCAGTCAt                      <  1:4639783‑M2/16‑1 (MQ=255)
               aacaaaatgggtgcagtacaCTTTCTGCCCAGTCAt                      <  1:5116024‑M2/16‑1 (MQ=255)
               aacaaaatgggtgcagtacaCTTTCTGCCCAGTCAt                      <  1:443466‑M2/16‑1 (MQ=255)
               aacaaaatgggtgcagtacaCTTTCTGCCCAGTCAt                      <  1:1480623‑M2/16‑1 (MQ=255)
               aaaaGCTGGAGGAAGCCTCACTTTtgtactgacccc                      >  1:2733800‑M1/1‑24 (MQ=255)
               aaaaGCTGGAGGAAGCCTCACTTTtgtactgacccc                      <  1:3938734‑M1/36‑13 (MQ=255)
                acaaaatgggtgcagtacaCTTTCTGCCCAGTCATc                     >  1:1266429‑M2/20‑36 (MQ=255)
                aaaGCTGGAGGAAGCCTCACTTTtgtactgacccca                     <  1:3515120‑M1/36‑14 (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:3980988‑M2/20‑1 (MQ=255)
                   aaatgggtgcagtacaCTTTCTGCCCAGTCATccac                  <  1:3030726‑M2/20‑1 (MQ=255)
                    cTGGAGGAAGCCTCACTTTtgtactgaccccaaaaa                 <  1:2474881‑M1/36‑18 (MQ=255)
                     atgggtgcagtacaCTTTCTGCCCAGTCATccacca                >  1:629254‑M2/15‑36 (MQ=255)
                     atgggtgcagtacaCTTTCTGCCCAGTCATccacca                >  1:3745824‑M2/15‑36 (MQ=255)
                     atgggtgcagtacaCTTTCTGCCCAGTCATccacca                >  1:1935402‑M2/15‑36 (MQ=255)
                      tgggtgcagtacaCTTTCTGCCCAGTCATCCACCaa               <  1:5558911‑M2/23‑1 (MQ=255)
                      ggaggaAGCCTCACTTTtgtactgaccccaaaaagt               >  1:4773693‑M1/1‑17 (MQ=255)
                      ggaggaAGCCTCACTTTtgtactgaccccaaaaagt               >  1:681979‑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:2999480‑M2/25‑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)
                            aGCCTCACTTTtgtactgaccccaaaaagttggaca         >  1:4551445‑M1/1‑11 (MQ=255)
                             agtacaCTTTCTGCCCAGTCATCCACCAaagaagaa        >  1:3762210‑M2/7‑36 (MQ=255)
                              gtacaCTTTCTGCCCAGTCATCCACCAaagaagaag       >  1:2890612‑M2/6‑36 (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)
                                acaCTTTCTGCCCATTCATCCACCAAAGAAGAAGcc     >  1:1767284‑M2/4‑36 (MQ=255)
                                acaCTTTCTGCCCAGTCATCCACCAAAGAAGAAGcc     >  1:1935244‑M2/4‑36 (MQ=255)
                                acaCTTTCTGCCCAGTCATCCACCAAAGAAGAAGcc     <  1:3912605‑M2/33‑1 (MQ=255)
                                 caCTTTCTGCCCAGTCATCCACCAAAGAAGAAGCCa    >  1:1487903‑M2/3‑36 (MQ=255)
                                 caCTTTCTGCCCAGTCATCCACCAAAGAAGAAGCCa    <  1:699593‑M2/34‑1 (MQ=255)
                                 cACTTTCTGCAAAGTAATCCACCAAAGAAGAagcca    <  1:359731/36‑1 (MQ=25)
                                  aCTTTtgtactgaccccaaaaagttggacagttaaa   <  1:1847956‑M1/36‑32 (MQ=255)
                                  aCTTTtgtactgaccccaaaaagttggacagttaaa   <  1:410310‑M1/36‑32 (MQ=255)
                                  aCTTTCTGCCCAGTCATCAACCAAAGAAGAAGCCAg   >  1:2688501‑M2/2‑36 (MQ=255)
                                   cTTTtgtactgaccccaaaaagttggacagttaaac  <  1:1389932‑M1/36‑33 (MQ=255)
                                   cTTTtgtactgaccccaaaaagttggacagttaaac  >  1:2465218‑M1/1‑4 (MQ=255)
                                   cTTTtgtactgaccccaaaaagttggacagttaaac  <  1:4279700‑M1/36‑33 (MQ=255)
                                   |                                   
TGAATCTCTCGCTGCAAAAGCTGGAGGAAGCCTCACTTTCTGCCCAGTCATCCACCAAAGAAGAAGCCAGC  >  REL606/3015736‑3015806

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