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

AATCTCTCGCTGCAAAAGCTGGAGGAAGCCTCACTTT..................................  >  REL606/3015738‑3015774
...................................TTTGTACTGACCCCAAAAAGTTGGACAGTTAAACAC  <  REL606/3652533‑3652500
                                                                       
  TCTCTCGCTGCAAAAGCTGGAGGAAGCCTCACTTTT                                   <  1:3248500/36‑1 (MQ=255)
  TCTCTCGCTGCAAAAGCTGGAGGAAGCCTCACTTTT                                   <  1:5524381/36‑1 (MQ=255)
  TCTCTCGCTGCAAAAGCTGGAGGAAGCCTCACTTTT                                   <  1:1689590/36‑1 (MQ=255)
   CTCTCGCTGCAAAAGCTGGAGGAAGCCTCACTTTTG                                  >  1:2228618/1‑36 (MQ=255)
   CTCTCGCTGCAAAAGCTGGAGGAAGCCTCACTTTTG                                  >  1:5783680/1‑36 (MQ=255)
    TCTCGCTGCAAAAGCTGGAGGAAGCCTCACTTTTGT                                 >  1:4847499/1‑36 (MQ=255)
    TCTCGCTGCAAAAGCTGGAGGAAGCCTCACTTTTGT                                 <  1:1322040/36‑1 (MQ=255)
     CTCTCTGCAAAAGCTGGAGGAAGCCTCACTTTTGTA                                <  1:5893108/36‑1 (MQ=255)
     CTCGCTGCAAAAGCTGGAGGAAGCCTCACTTTTGTA                                >  1:4009279/1‑36 (MQ=255)
      TCGCTGCAAAAGCTGGAGGAAGCCTCACTTTTGTAC                               <  1:1626819/36‑1 (MQ=255)
       CGCTGCAAAAGCTGGAGGAAGCCTCACTTTTGTACT                              <  1:4478402/36‑1 (MQ=255)
        GCTGCAAAAGCTGGAGGAAGCCTCACTTTTGTACTG                             >  1:775927/1‑36 (MQ=255)
        GCTGCAAAAGCTGGAGGAAGCCTCACTTTTGTACTG                             >  1:3915431/1‑36 (MQ=255)
        GCTGCAAAAGCTGGAGGAAGCCTCACTTTTGTACTG                             <  1:5087770/36‑1 (MQ=255)
         CTGCAAAAGCTGGAGGAAGCCTCACTTTTGTACTGA                            <  1:802282/36‑1 (MQ=255)
         CTGCAAAAGCTGGAGGAAGCCTCACTTTTGTACTGA                            >  1:473191/1‑36 (MQ=255)
          TGCAAAAGCTGGAGGAAGCCTCACTTTTGTACTGAC                           <  1:4481969/36‑1 (MQ=255)
          TGCAAAAGCTGGAGGAAGCCTCACTTTTGTACTGAC                           <  1:5210242/36‑1 (MQ=255)
           GCAAAAGCTGGAGGAAGCCTCACTTTTGTACTGACC                          >  1:2313750/1‑36 (MQ=255)
             AAAAGCTGGAGGAAGCCTCACTTTTGTACTGACCCC                        >  1:2733800/1‑36 (MQ=38)
             AAAAGCTGGAGGAAGCCTCACTTTTGTACTGACCCC                        <  1:4636560/36‑1 (MQ=38)
             AAAAGCTGGAGGAAGCCTCACTTTTGTACTGACCCC                        >  1:3575707/1‑36 (MQ=38)
             AAAAGCTGGAGGAAGCCTCACTTTTGTACTGACCCC                        <  1:3938734/36‑1 (MQ=38)
              AAAGCTGGAGGAAGCCTCACTTTTGTACTGACCCCA                       <  1:3515120/36‑1 (MQ=38)
               AAGCTGGAGGAAGCCTCACTTTTGTACTGACCCCta                      >  1:893765/1‑34 (MQ=25)
               AAGCTGGAGGAAGCCTCACTTTTGTACTGACCCCAA                      >  1:557101/1‑36 (MQ=35)
               AAGCTGGAGGAAGCCTCACTTTTGTACTGACCCCAA                      >  1:971262/1‑36 (MQ=35)
               AAGCTGGAGGAAGCCTCACTTTTGTACTGACCCCAA                      >  1:1272556/1‑36 (MQ=35)
                AGCTGGAGGAAGCCTCACTTTTGTACTGACCCCAAA                     >  1:555464/1‑36 (MQ=32)
                AGCTGGAGGAAGCCTCACTTTTGTACTGACCCCAAA                     <  1:5243616/36‑1 (MQ=32)
                 GCTGGAGGAAGCCTCACTTTTGTACTGACCCCAAAA                    <  1:2150569/36‑1 (MQ=1)
                 GCTGGAGGAAGCCTCACTTTTGTACTGACCCCAAAA                    <  1:3514421/36‑1 (MQ=1)
                 GCTGGAGGAAGCCTCACTTTTGTACTGACCCCAAAA                    >  1:1498344/1‑36 (MQ=1)
                  CTGGAGGAAGCCTCACTTTTGTACTGACCCCAAAAA                   <  1:2474881/36‑1 (MQ=1)
                    GGAGGAAGCCTCACTTTTGTACTGACCCCAAAAAGT                 >  1:681979/1‑36 (MQ=1)
                    GGAGGAAGCCTCACTTTTGTACTGACCCCAAAAAGT                 >  1:4773693/1‑36 (MQ=1)
                     ggGGAAGCCTCACTTTTGTACTGACCCCAAAAAGTT                >  1:4188176/3‑36 (MQ=1)
                     GAGGAAGCCTCACTTTTGTACTGACCCCAAAAATTT                >  1:775796/1‑36 (MQ=0)
                     GAGGAAGCCTCACTTTTGTACTGACCCCAAAAAGTT                <  1:4892485/36‑1 (MQ=1)
                     GAGGAAGCCTCACTTTTGTACTGACCCCAAAAAGTT                >  1:1309324/1‑36 (MQ=1)
                       GGAAGCCTCACTTTTGTACTGACCCCAAAAAGTTGG              <  1:5473186/36‑1 (MQ=1)
                        GAAGCCTCACTTTTGTACTGACCCCAAGAAGTTGGA             <  1:2876506/36‑1 (MQ=0)
                         AAGCCTCACTTTTGTACTGACCCCAAAAAGTTGGAC            >  1:3382003/1‑36 (MQ=1)
                          AGCCTCACTTTTGTACTGACCCCAAAAAGTTGGACA           <  1:4970273/36‑1 (MQ=1)
                          AGCCTCACTTTTGTACTGACCCCAAAAAGTTGGACA           >  1:4551445/1‑36 (MQ=1)
                            CCTCACTTTTGTACTGACCCCAAAAAGTTGGACAGT         <  1:4060139/36‑1 (MQ=1)
                             CTCAGTTTTGTACTGACCCCAAAAAGTTGGACAGTT        >  1:3654675/1‑36 (MQ=0)
                             CTCACTTTTGTACTGACCCCAAAAAGTTGTACAGTT        <  1:2595381/36‑1 (MQ=0)
                             CTCACTTTTGTACTGACCCCAAAAAGTTGGACAGTT        >  1:2476383/1‑36 (MQ=1)
                                ACTTTTGTACTGACCCCAAAAAGTTGGACAGTTAAA     <  1:1847956/36‑1 (MQ=1)
                                ACTTTTGTACTGACCCCAAAAAGTTGGACAGTTAAA     <  1:410310/36‑1 (MQ=1)
                                 CTTTTGTACTGACCCCAAAAAGTTGGACAGTTAAAC    >  1:2465218/1‑36 (MQ=1)
                                 CTTTTGTACTGACCCCAAAAAGTTGGACAGTTAAAC    <  1:1389932/36‑1 (MQ=1)
                                 CTTTTGTACTGACCCCAAAAAGTTGGACAGTTAAAC    <  1:4279700/36‑1 (MQ=1)
                                                                       
AATCTCTCGCTGCAAAAGCTGGAGGAAGCCTCACTTT..................................  >  REL606/3015738‑3015774
...................................TTTGTACTGACCCCAAAAAGTTGGACAGTTAAACAC  <  REL606/3652533‑3652500

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