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

CGTGTTTAACTGTCCAACTTTTTGGGGTCAGTACAAACTTTGATTTATAGTCAGGTGGGGCTTTTCTCTG  >  REL606/3652499‑3652568
                                  |                                   
cttgtttAACTGTCCAACTTTTTGGGGTCAGTACAa                                    <  1:2807421‑M2/34‑2 (MQ=255)
 gtgtTTAACTGTCCAACTTTTTGGGGTCAGTACAaa                                   <  1:4826587‑M2/36‑3 (MQ=255)
 gtgtTTAACTGTCCAACTTTTTGGGGTCAGTACAaa                                   >  1:2987721‑M2/1‑34 (MQ=255)
  tgtTTAACTGTCCAACTTTTTGGGGTCAGTACAAAc                                  <  1:2129808/36‑1 (MQ=35)
   gtTTAACTGTCCAACTTTTTGGGGTCAGTACAaaag                                 <  1:4279700‑M2/36‑5 (MQ=255)
   gtTTAACTGTCCAACTTTTTGGGGTCAGTACAaaag                                 <  1:1389932‑M2/36‑5 (MQ=255)
   gtTTAACTGTCCAACTTTTTGGGGTCAGTACAaaag                                 >  1:2465218‑M2/1‑32 (MQ=255)
   gtTTAACTGTCCAACTTTTTGGGGTCAGTACAAACt                                 <  1:5315288/36‑1 (MQ=38)
    tttAACTGTCCAACTTTTTGGGGTCAGTACAaaagt                                <  1:410310‑M2/36‑6 (MQ=255)
    tttAACTGTCCAACTTTTTGGGGTCAGTACAaaagt                                <  1:1847956‑M2/36‑6 (MQ=255)
     ttAACTGTCCAACTTTTTTGGGTCAGTACAAACttt                               <  1:1570240/36‑1 (MQ=255)
     ttAACTGTCCAACTTTTTGGGGTCAGTACAAACttt                               >  1:628629/1‑36 (MQ=255)
      tAACTGTCCAACTTTTTGGGGTCAGTACAAACTTTg                              >  1:4075555/1‑36 (MQ=255)
       aaCTGTCCAACTTTTTGGGGTCAGTACAaaagtgag                             >  1:2476383‑M2/1‑28 (MQ=255)
       aaCTGTCCAACTTTTTGGGGTCAGTACAaaactgag                             >  1:3654675‑M2/1‑28 (MQ=255)
       aaCTGTCCAACTTTTTGGGGTCAGTACAAACTTTGa                             >  1:2191119/1‑36 (MQ=255)
       aaCTGTACAACTTTTTGGGGTCAGTACAaaagtgag                             <  1:2595381‑M2/36‑9 (MQ=255)
        aCTGTCCAACTTTTTGGGGTCAGTACAaaagtgagg                            <  1:4060139‑M2/36‑10 (MQ=255)
        aCTGTCCAACTTTTTGGGGTCAGTACAAACTTTGAt                            <  1:876687/36‑1 (MQ=255)
          tGTCCAACTTTTTGGGGTCAGTACAaaagtgaggct                          <  1:4970273‑M2/36‑12 (MQ=255)
          tGTCCAACTTTTTGGGGTCAGTACAaaagtgaggct                          >  1:4551445‑M2/1‑25 (MQ=255)
          tGTCCAACTTTTTGGGGTCAGTACAAACTTTGAttt                          <  1:5082237/36‑1 (MQ=255)
          tGTCCAACTTTTTGGGGTCAGTACAAACTTTGAttt                          <  1:63592/36‑1 (MQ=255)
          tGTCCAACTTTTTGGGGTCAGTACAAACTTTGAttt                          >  1:1786235/1‑36 (MQ=255)
           gTCCAACTTTTTGGGGTCAGTACAaaagtgaggctt                         >  1:3382003‑M2/1‑24 (MQ=255)
            tCCAACTTTTTGGGGTCAGTACAAACTTTGATTtat                        <  1:5108031/36‑1 (MQ=255)
            tCCAACTTCTTGGGGTCAGTACAaaagtgaggcttc                        <  1:2876506‑M2/36‑14 (MQ=255)
             ccAACTTTTTGGGGTCAGTACAaaagtgaggcttcc                       <  1:5473186‑M2/36‑15 (MQ=255)
             ccAACTTTTTGGGGTCAGTACAACCTTTGATTtata                       <  1:5245846/36‑1 (MQ=255)
             ccAACTTTTTGGGGTCAGTACAAACTTTGATTtata                       <  1:1458900/36‑1 (MQ=255)
             ccAACTTTTTGGGGTCAGTACAAACTTTGATTtata                       <  1:5341831/36‑1 (MQ=255)
               aaCTTTTTGGGGTCAGTACAaaagtgaggcttcctc                     >  1:1309324‑M2/1‑20 (MQ=255)
               aaCTTTTTGGGGTCAGTACAaaagtgaggcttcctc                     <  1:4892485‑M2/36‑17 (MQ=255)
               aaCTTTTTGGGGTCAGTACAaaagtgaggcttcccc                     >  1:4188176‑M2/1‑20 (MQ=255)
               aaCTTTTTGGGGTCAGTACAAACTTTGATTTATAGt                     <  1:5208186/36‑1 (MQ=255)
               aaCTTTTTGGGGTCAGTACAAACTTTGATTTATAGt                     <  1:5159090/36‑1 (MQ=255)
               aaATTTTTGGGGTCAGTACAaaagtgaggcttcctc                     >  1:775796‑M2/1‑20 (MQ=255)
                aCTTTTTGGGGTCAGTACAaaagtgaggcttcctcc                    >  1:4773693‑M2/1‑19 (MQ=255)
                aCTTTTTGGGGTCAGTACAaaagtgaggcttcctcc                    >  1:681979‑M2/1‑19 (MQ=255)
                aCTTTTTGGGGTCAGTACAAACTTTGATTTATAGTc                    <  1:3084600/36‑1 (MQ=255)
                  tttttGGGGTCAGTACAaaagtgaggcttcctccag                  <  1:2474881‑M2/36‑20 (MQ=255)
                   ttttGGGGTCAGTACAaaagtgaggcttcctccagc                 <  1:2150569‑M2/36‑21 (MQ=255)
                   ttttGGGGTCAGTACAaaagtgaggcttcctccagc                 <  1:3514421‑M2/36‑21 (MQ=255)
                   ttttGGGGTCAGTACAaaagtgaggcttcctccagc                 >  1:1498344‑M2/1‑16 (MQ=255)
                   ttttGGGGTCAGTACAAACTTTGATTTATAGTCAgg                 >  1:3049913/1‑36 (MQ=255)
                   ttttGGGGTCAGTACAAACTTTGATTTATAGTCAgg                 >  1:2853017/1‑36 (MQ=255)
                   ttttGGGGCCAGTACAAACTTTGATTTATAGTCAgg                 >  1:1724702/1‑36 (MQ=255)
                    tttGGGGTCAGTACAaaagtgaggcttcctccagct                <  1:5243616‑M2/36‑22 (MQ=255)
                    tttGGGGTCAGTACAaaagtgaggcttcctccagct                >  1:555464‑M2/1‑15 (MQ=255)
                    tttGGGGTCAGTACAAACTTTGATTTATAGTCAGGt                <  1:1509926/36‑1 (MQ=255)
                     ttGGGGTCAGTACAaaagtgaggcttcctccagctt               >  1:557101‑M2/1‑14 (MQ=255)
                     ttGGGGTCAGTACAaaagtgaggcttcctccagctt               >  1:971262‑M2/1‑14 (MQ=255)
                     ttGGGGTCAGTACAaaagtgaggcttcctccagctt               >  1:1272556‑M2/1‑14 (MQ=255)
                     taggggtcAGTACAaaagtgaggcttcctccagctt               >  1:893765‑M2/3‑14 (MQ=255)
                      tGGGGTCAGTACAaaagtgaggcttcctccagcttt              <  1:3515120‑M2/36‑24 (MQ=255)
                      tGGGGTCAGTACAAACTTTGATTTATAGTCAGGTgg              <  1:3714777/36‑1 (MQ=255)
                       ggggTCAGTACAaaagtgaggcttcctccagctttt             >  1:2733800‑M2/1‑12 (MQ=255)
                       ggggTCAGTACAaaagtgaggcttcctccagctttt             <  1:3938734‑M2/36‑25 (MQ=255)
                       ggggTCAGTACAaaagtgaggcttcctccagctttt             <  1:4636560‑M2/36‑25 (MQ=255)
                       ggggTCAGTACAaaagtgaggcttcctccagctttt             >  1:3575707‑M2/1‑12 (MQ=255)
                       ggggTCAGTACAAACTTTGATTTATAGTCAGGTggg             <  1:4025114/36‑1 (MQ=255)
                       ggggTCAGTACAAACTTTGATTTATAGTCAGGTggg             >  1:3883202/1‑36 (MQ=255)
                       ggggTAAGTACAAACTTTGATTTATAGTCCGGTggg             <  1:538050/36‑1 (MQ=38)
                         ggTCAGTACAaaagtgaggcttcctccagcttttgc           >  1:2313750‑M2/1‑10 (MQ=255)
                          gTCAGTACAaaagtgaggcttcctccagcttttgca          <  1:4481969‑M2/36‑28 (MQ=255)
                          gTCAGTACAaaagtgaggcttcctccagcttttgca          <  1:5210242‑M2/36‑28 (MQ=255)
                          gTCAGTACAACCTTTGATTTATAGTCAGGTGGGGCt          <  1:3426586/36‑1 (MQ=255)
                          gTCAGTACAAACTTTGATTTATAGTCAGGTGGGGCt          >  1:4179216/1‑36 (MQ=255)
                          gTCAGTACAAACTTTGATTTATAGTCAGGTGGGGCt          >  1:3034226/1‑36 (MQ=255)
                           tCAGTACAaaagtgaggcttcctccagcttttgcag         <  1:802282‑M2/36‑29 (MQ=255)
                           tCAGTACAaaagtgaggcttcctccagcttttgcag         >  1:473191‑M2/1‑8 (MQ=255)
                            cAGTACAaaagtgaggcttcctccagcttttgcagc        <  1:5087770‑M2/36‑30 (MQ=255)
                            cAGTACAaaagtgaggcttcctccagcttttgcagc        >  1:3915431‑M2/1‑7 (MQ=255)
                            cAGTACAaaagtgaggcttcctccagcttttgcagc        >  1:775927‑M2/1‑7 (MQ=255)
                            cAGTACAAACTTTGATTTATAGTCAGGTGGGGCttt        <  1:4132656/36‑1 (MQ=255)
                             aGTACAaaagtgaggcttcctccagcttttgcagcg       <  1:4478402‑M2/36‑31 (MQ=255)
                             aGTACAAACTTTGATTTATAGTCAGGTGGGGCtttt       <  1:3681883/36‑1 (MQ=255)
                             aGTACAAACTTTGATTTATAGTCAGGTGGGGCtttt       <  1:3927858/36‑1 (MQ=255)
                             aGTACAAACTTTGATTTATAGTCAGGTGGGGCtttt       <  1:5437268/36‑1 (MQ=255)
                              gTACAaaagtgaggcttcctccagcttttgcagcga      <  1:1626819‑M2/36‑32 (MQ=255)
                              gTACAAACTTTGATTTATAGTCAGGTGGGGCTTTtc      <  1:2071123/36‑1 (MQ=255)
                              gTACAAACTTTGATTTATAGTCAGGTGGGGCTTTtc      <  1:5588852/36‑1 (MQ=255)
                              gTACAAACTTTGATTTATAGTCAGGTGGGGCTTTtc      >  1:239842/1‑36 (MQ=255)
                               tACAaaagtgaggcttcctccagcttttgcagcgag     >  1:4009279‑M2/1‑4 (MQ=255)
                               tACAaaagtgaggcttcctccagcttttgcagagag     <  1:5893108‑M2/36‑33 (MQ=255)
                               tACAAACTTTGATTTATAGTCAGGTGGGGCTTTtct     >  1:4250862/1‑36 (MQ=255)
                                aCAaaagtgaggcttcctccagcttttgcagcgaga    <  1:1322040‑M2/36‑34 (MQ=255)
                                aCAaaagtgaggcttcctccagcttttgcagcgaga    >  1:4847499‑M2/1‑3 (MQ=255)
                                aCAAACTTTGATTTATAGTCAGGTGGGGCTTTtctc    <  1:2620850/36‑1 (MQ=255)
                                 cAaaagtgaggcttcctccagcttttgcagcgagag   >  1:5783680‑M2/1‑2 (MQ=255)
                                 cAaaagtgaggcttcctccagcttttgcagcgagag   >  1:2228618‑M2/1‑2 (MQ=255)
                                 cAAACTTTGATTTATAGTCAGGTGGGGCTTTTctct   <  1:843177/36‑1 (MQ=255)
                                 cAAACTTTGATTTATAGTCAGGTGGGGCTTTTctct   >  1:361552/1‑36 (MQ=255)
                                  aaaagtgaggcttcctccagcttttgcagcgagaga  <  1:5524381‑M2/36‑36 (MQ=255)
                                  aaaagtgaggcttcctccagcttttgcagcgagaga  <  1:1689590‑M2/36‑36 (MQ=255)
                                  aaaagtgaggcttcctccagcttttgcagcgagaga  <  1:3248500‑M2/36‑36 (MQ=255)
                                  aaaCTTTGATTTATAGACAGATGGGGATTTTatctt  >  1:4789098/1‑35 (MQ=12)
                                  |                                   
CGTGTTTAACTGTCCAACTTTTTGGGGTCAGTACAAACTTTGATTTATAGTCAGGTGGGGCTTTTCTCTG  >  REL606/3652499‑3652568

Base quality scores:  ATCG  < 3 ≤  ATCG  < 16 ≤  ATCG  < 23 ≤  ATCG  < 32 ≤  ATCG  < 37 ≤  ATCG