Predicted mutation
evidence seq id position mutation annotation gene description
JC JC REL606 3,015,771 IS150 (–) +4 bp 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 = 2775877NA (NA)49 (0.940) 33/66 0.2 98.0% noncoding (1/1443 nt) IS150 repeat region
?REL606 3015771 = 1 (0.200)coding (1192/1677 nt) ECB_02816 KpsD protein
* ? REL606 = 30157741 (0.150)55 (1.000) 32/70 0.3 98.2% coding (1195/1677 nt) ECB_02816 KpsD protein
?REL606 = 3894996 NA (NA)noncoding (1443/1443 nt) IS150 repeat region

CATTTCATCGTCCAACAAAATGGGTGCAGTACA...................................  >  REL606/2775845‑2775877
...............................CACTTTCTGCCCAGTCATCCACCAAAGAAGAAGCCAG  >  REL606/3015771‑3015805
                                                                    
CATTTCATCGTCCAACAAAATGGGCGCAGTACACTT                                  >  1:3707352/1‑36 (MQ=0)
 ATTTCATCGTCCAACAAAATGGGTGCAGTACACcgg                                 <  2:2431258/36‑4 (MQ=1)
 ATTTCATCGTCCAACAAAATGGGTGCAGTACACcgg                                 >  2:3225739/1‑33 (MQ=1)
 ATTTCATCGTCCAACAAAATGGGTGCAGTACACaac                                 >  2:1686687/1‑33 (MQ=0)
   TTCATCGTCCAGCAAAATGGGTGCAGTACACTTTCT                               <  1:3007280/36‑1 (MQ=0)
    TCATCGTCCGACAAAATGGGTGCAGTACACTTTCTG                              >  1:1294855/1‑36 (MQ=0)
    TCATCGTCCAACAAAATGGGTGCAGTACACTTTCTG                              <  2:1530587/36‑1 (MQ=1)
     CATCGTCCAACAAAATGGGTGCAGTACACTTTCTGC                             <  1:3212360/36‑1 (MQ=32)
     CATCGTCCAACAAAATGGGTGCAGTACACTTTCTGC                             >  1:791289/1‑36 (MQ=32)
      ATCGTCCAACAAAATGGGTGCAGTACACTTTCTGCC                            >  2:3663076/1‑36 (MQ=35)
       TCGTCCAACAAAATGGGTGCAGTACACTTTCTGCCC                           >  1:2109187/1‑36 (MQ=38)
         GTCCAACAAAATGGGTGCCGTACACTTTCTGCCCAG                         >  2:2007225/1‑36 (MQ=255)
         GTCCAACAAAATGGGTGCAGTACACTTTCTGCCCAG                         <  2:1626024/36‑1 (MQ=38)
          TCCAACAAAATGGGTGCAGTACACTTTCTGCTCAGT                        <  2:2657334/36‑1 (MQ=255)
          TCCAACAAAATGGGTGCAGTACACTTTCTGCCCAGT                        <  1:312984/36‑1 (MQ=38)
          TCCAACAAAATGGGTGCAGTACACTTTCTGCCCAGT                        >  1:3464077/1‑36 (MQ=38)
           CCAACAAAATGGGTGCAGTACACTTTCTGCCCAGTC                       <  2:1109873/36‑1 (MQ=38)
            CAACAAAATGGGTGCAGTACACTTTCTGCCCAGTCA                      >  2:3248431/1‑36 (MQ=40)
            CAACAAAATGGGTGCAGTACACTTTCTGCCCAGTCA                      >  2:2309774/1‑36 (MQ=40)
             AACAAAATGGGTGCAGTACACTTTCTGCCCAGTCAT                     <  1:1538542/36‑1 (MQ=255)
              ACAAAATTGGTGCAGTACACTTTCTGCCCAGTCATC                    <  2:173679/36‑1 (MQ=255)
                AAAATGGGTGCAGTACACTTTCTGCCCAGTTATCCA                  >  1:3118564/1‑36 (MQ=255)
                AAAATGGGTGCAGTACACTTTCTGCCCAGTCATCCA                  >  1:3012721/1‑36 (MQ=255)
                AAAATGGGTGCAGTACACTTTCTGCCCAGTCATCCA                  <  1:2841118/36‑1 (MQ=255)
                AAAATGGGTGCAGTACACTTTCTGCCCAGTCATCCA                  >  1:2688366/1‑36 (MQ=255)
                AAAATGGGTGCAGTACACTTTCTGCCCAGTCATCCA                  >  1:26606/1‑36 (MQ=255)
                 AAATGGGTGCAGTACACTTTCTGCCCAGTCATCCAC                 >  1:2235403/1‑36 (MQ=255)
                 AAATGGGTGCAGTACACTTTCTGCCCAGTCATCCAC                 <  1:1787012/36‑1 (MQ=255)
                 AAAGGGGTGCAGTACACTTTCTGCCCAGTCATCCAC                 <  2:2072391/36‑1 (MQ=255)
                  AATGGGCGCAGTACACTTTCTGCCCAGTCATCCACC                >  2:2628651/1‑36 (MQ=255)
                   ATGGGTGCAGTACACTTTCTGCCCAGTCATCCACta               >  1:2764810/1‑34 (MQ=255)
                    TGGGTGCAGTACCCTTTCTGCCCAGTCATCCACCAA              >  2:1784357/1‑36 (MQ=255)
                    TGGGTGCAGTACACTTTCTGCCCAGTCATCCACCAA              <  2:430658/36‑1 (MQ=255)
                     GGGTGCAGTACACTTTCTGCCCAGTCATCCACCAAA             <  1:892269/36‑1 (MQ=255)
                      GGTGCAGTACACTTTCTGCCCAGTCATCCACCAAcg            >  1:455081/1‑34 (MQ=255)
                       GTGCAGTACACTTTCTGCCCAGTCATCCACCAAAGA           >  1:764288/1‑36 (MQ=255)
                       GTGCAGTACACTTTCTGCCCAGTCATCCACCAAAGA           >  1:289823/1‑36 (MQ=255)
                       GTGCAGTACACTTTCTGCCCAGTCATCCACCAAAGA           >  2:3676929/1‑36 (MQ=255)
                       GTGCAGTACACTTTCTGCCCAGTCATCCACCAAAGA           >  1:1122133/1‑36 (MQ=255)
                         GCAGTACACTTTCTGCCCAGTCATCCACCAAAGAAG         <  2:2675430/36‑1 (MQ=255)
                          CAGTACACTTTCTGCCCAGTCATCCACCAAAGAAGA        <  1:91716/36‑1 (MQ=255)
                          CAGTACACTTTCTGCCCAGTCATCCACCAAAGAAGA        >  2:411315/1‑36 (MQ=255)
                           AGTACACTTTCTGCCCAGTCATCCACCAAAGAAGAA       <  1:3279116/36‑1 (MQ=255)
                             TACACTTTCTGCCCAGTCATCCACCAAAGAAGAAGC     >  1:3111387/1‑36 (MQ=255)
                             TACACTTTCTGCCCAGTCATCCACCAAAGAAGAAGC     >  1:1779211/1‑36 (MQ=255)
                             TACACTTGCTGCCCAGTCATCCACCAAAGAAGAAGC     <  2:3561696/36‑1 (MQ=255)
                              ACACTTTCTGCCCAGTCATCCACCAAAGAAGAAGCC    >  1:3060729/1‑36 (MQ=255)
                              ACACTTTCTGCCCAGTCATCCACCAAAGAAGAAGCC    <  1:2628850/36‑1 (MQ=255)
                              ACACTTTCTGCCCAGTCATCCACCAAAGAAGAAGCC    <  1:1645864/36‑1 (MQ=255)
                                                                    
CATTTCATCGTCCAACAAAATGGGTGCAGTACA...................................  >  REL606/2775845‑2775877
...............................CACTTTCTGCCCAGTCATCCACCAAAGAAGAAGCCAG  >  REL606/3015771‑3015805

Base quality scores:  ATCG  < 3 ≤  ATCG  < 12 ≤  ATCG  < 22 ≤  ATCG  < 30 ≤  ATCG  < 32 ≤  ATCG