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

ATTCCATTTCATCGTCCAACAAAATGGGTGCAGTACA...................................  <  REL606/588531‑588495
...................................CACTTTCTGCCCAGTCATCCACCAAAGAAGAAGCCAG  >  REL606/3015771‑3015805
                                                                        
  TCCATTTCATCGTCCAACAAAATGGGTGCAGTACAC                                    >  1:32994/1‑36 (MQ=1)
  TCCATTTCATCGTCCAACAAAATGGGTGCAGTACAC                                    >  1:5575559/1‑36 (MQ=1)
  TCCATTTCATCGTCCAACAAAATGGGTGCAGTACAC                                    >  1:5302910/1‑36 (MQ=1)
    CATTTCATCGTCCAACAAAATGGGTGCAGTACACcg                                  <  1:5030201/36‑3 (MQ=1)
    CATTTCATCGTCCAACAAAATGGGTGCAGTACACcg                                  <  1:4089607/36‑3 (MQ=1)
    CATTTCATCGTCCAACAAAATGGGTGCAGTACACTT                                  <  1:5001079/36‑1 (MQ=1)
     ATTTCATCGTCCAACAAAATGGTTGCAGTACACTTT                                 <  1:2633218/36‑1 (MQ=0)
     ATTTCATCGTCCAACAAAATGGGTGCAGTACACTTT                                 >  1:505086/1‑36 (MQ=1)
     ATTTCATCGTCCAACAAAATGGGTGCAGTACACTTT                                 <  1:3674647/36‑1 (MQ=1)
      TTTCATCGTCCAACAAAATGGGTGCAGTACACTTTC                                <  1:5659277/36‑1 (MQ=1)
      TTTCATCGTCCAACAAAATGGGTGCAGTACACTTTC                                <  1:4291068/36‑1 (MQ=1)
       TTCATCGTCCAACAAAATGGGTGCAGTACACTTTCT                               >  1:5668941/1‑36 (MQ=1)
       TTCATCGTCCAACAAAATGGGTGCAGTACACTTTCT                               <  1:3384874/36‑1 (MQ=1)
       TTCATCGTCCAACAAAATGGGTGCAGTACACTTTCT                               <  1:4939874/36‑1 (MQ=1)
        TCATCGTCCAACAAAATGGGTGCAGTCCACTTTCTG                              <  1:4489903/36‑1 (MQ=0)
        TCATCGTCCAACAAAATGGGTGCAGTACACTTTCTG                              <  1:3326452/36‑1 (MQ=1)
        TCATCGTCCAACAAAATGGGTGCAGTACACTTTCTG                              <  1:2286918/36‑1 (MQ=1)
        TCATCGTCCAACAAAATGGGTGCAGTACACTTTCTG                              >  1:4962357/1‑36 (MQ=1)
         CATCGTCCAACAAAATGGGTGCAGTACACTTTCTGC                             >  1:3720376/1‑36 (MQ=1)
         CATCGTCCAACAAAATGGGTGCAGTACACTTTCTGC                             <  1:1513001/36‑1 (MQ=1)
          ATCGTCCAACAAAATGGGTGCAGTACACTTTCTGCC                            >  1:81350/1‑36 (MQ=1)
           TCGTCCAACAAAATGGGTGCAGTCCACTTTCTGCCC                           <  1:3234553/36‑1 (MQ=12)
           TCGTCCAACAAAATGGGTGCAGTACACTTTCTGCCC                           <  1:5408728/36‑1 (MQ=32)
            CGTCCAACAAAATGGGTGCAGTACACTTTCTGCCCA                          >  1:4324547/1‑36 (MQ=35)
            CGTCCAACAAAATGGGTGCAGTACACTTTCTGCCCA                          <  1:538007/36‑1 (MQ=35)
            CGTCCAACAAAATGGGTGCAGTACACTTTCTGCCCA                          <  1:361038/36‑1 (MQ=35)
              TCCAACAAAATGGGTGCAGTACACTTTCTGCCCAGT                        >  1:1969056/1‑36 (MQ=38)
               CCAACAAAATTGGTGCAGTACACTTTCTGCCCAGTC                       <  1:4802819/36‑1 (MQ=255)
               CCAACAAAATGGGTGCAGTACGCTTTCTGCCCAGTC                       <  1:5227793/36‑1 (MQ=255)
               CCAACAAAATGGGTGCAGTACACTTTCTGCCCAGTC                       <  1:2431871/36‑1 (MQ=38)
                 AACAAAATGGGTGCAGTACACTTTCTGCCCAGTCAT                     <  1:1480623/36‑1 (MQ=255)
                 AACAAAATGGGTGCAGTACACTTTCTGCCCAGTCAT                     <  1:443466/36‑1 (MQ=255)
                 AACAAAATGGGTGCAGTACACTTTCTGCCCAGTCAT                     <  1:5116024/36‑1 (MQ=255)
                 AACAAAATGGGTGCAGTACACTTTCTGCCCAGTCAT                     <  1:4639783/36‑1 (MQ=255)
                  ACAAAATGGGTGCAGTACACTTTCTGCCCAGTCATC                    >  1:605637/1‑36 (MQ=255)
                  ACAAAATGGGTGCAGTACACTTTCTGCCCAGTCATC                    >  1:1266429/1‑36 (MQ=255)
                   CAAAATGGGTGCAGTACACTTTCTGCCCAGTCATCC                   <  1:5858404/36‑1 (MQ=255)
                   CAAAATGGGTGCAGTACACTTTCTGCCCAGTCATCC                   >  1:2865298/1‑36 (MQ=255)
                   CAAAATGGGTGCAGTACACTTTCTGCCCAGTCATCC                   >  1:3341192/1‑36 (MQ=255)
                   CAAAATGAGTGCAGTACACTTTCTGCCCAGTCATCC                   <  1:126203/36‑1 (MQ=255)
                    AAAATGGGTGCAGTACACTTTCTGCCCAGTCATCCA                  >  1:2209821/1‑36 (MQ=255)
                    AAAATGGGTGCAGTACACTTTCTGCCCAGTCATCCA                  <  1:4604750/36‑1 (MQ=255)
                     AAATGGGTGCAGTACACTTTCTGCCCAGTCATCCAC                 <  1:3030726/36‑1 (MQ=255)
                     AAATGGGTGCAGTACACTTTCTGCCCAGTCATCCAC                 <  1:3980988/36‑1 (MQ=255)
                       ATGGGTGCAGTACACTTTCTGCCCAGTCATCCACCA               >  1:1935402/1‑36 (MQ=255)
                       ATGGGTGCAGTACACTTTCTGCCCAGTCATCCACCA               >  1:629254/1‑36 (MQ=255)
                       ATGGGTGCAGTACACTTTCTGCCCAGTCATCCACCA               >  1:3745824/1‑36 (MQ=255)
                        TGGGTGCAGTACACTTTCTGCCCAGTCATCCACCAA              <  1:1777862/36‑1 (MQ=255)
                        TGGGTGCAGTACACTTTCTGCCCAGTCATCCACCAA              <  1:5558911/36‑1 (MQ=255)
                         GGGTGCAGTACACTTTCTGCCCAGTCATCCACCAAA             >  1:5236612/1‑36 (MQ=255)
                          GGTGCAGTACACTTTCTGCCCAGTCATCCACCAAAG            <  1:5090981/36‑1 (MQ=255)
                          GGTGCAGTACACTTTCTGCCCAGTCATCCACCAAAG            <  1:2999480/36‑1 (MQ=255)
                           GTGCAGTACACTTTCTGCCCAGTCATCCACCAAAGA           <  1:3118425/36‑1 (MQ=255)
                           GTGCAGTACACTTTCTGCCCAGCCATCCACCAAAGA           >  1:46851/1‑36 (MQ=255)
                            TGCAGTACACTTTCTGCCCAGTCATCCACCAAAGAA          <  1:5304408/36‑1 (MQ=255)
                             GCAGTACACTTTCTGCCCAGGCATCCACCAAAGAAG         >  1:283614/1‑36 (MQ=255)
                              CAGTACACTTTCTGCCCAGTCATCCACCAAAGAAag        >  1:5359568/1‑34 (MQ=255)
                               AGTACACTTTCTGCCCAGTCATCCACCAAAGAAGAA       >  1:3762210/1‑36 (MQ=255)
                                GTACACTTTCTGCCCAGTCATCCACCAAAGAAGAAG      >  1:2890612/1‑36 (MQ=255)
                                  ACACTTTCTGCCCATTCATCCACCAAAGAAGAAGCC    >  1:1767284/1‑36 (MQ=255)
                                  ACACTTTCTGCCCAGTCATCCACCAAAGAAGAAGCC    <  1:3912605/36‑1 (MQ=255)
                                  ACACTTTCTGCCCAGTCATCCACCAAAGAAGAAGCC    >  1:1935244/1‑36 (MQ=255)
                                                                        
ATTCCATTTCATCGTCCAACAAAATGGGTGCAGTACA...................................  <  REL606/588531‑588495
...................................CACTTTCTGCCCAGTCATCCACCAAAGAAGAAGCCAG  >  REL606/3015771‑3015805

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