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

CTCTCGCTGCAAAAGCTGGAGGAAGCCTCACTTT...................................  >  REL606/3015741‑3015774
..................................TGTACTGACCCCAAAAAGTTGGACAGTTAAACACG  <  REL606/3894996‑3894962
                                                                     
CTCTCGCTGCAAAAGCTGGAGGAAGCCTCACTTTTG                                   >  2:328517/1‑36 (MQ=1)
 TCTCGCTGCAAAAGCTGGAGGAAGCCTCACTTTTGT                                  <  2:1300651/36‑1 (MQ=32)
 TCTCGCTGCAAAAGCTGGAGGAAGCCTCACTTTTGT                                  >  2:1469917/1‑36 (MQ=32)
 TCTCGCTGCAAAAGCTGGAGGAAGCCTCACTTTTGT                                  >  2:1914102/1‑36 (MQ=32)
 TCTCGCTGCAAAAGCTGGAGGAAGCCTCACTTTTGT                                  >  2:372091/1‑36 (MQ=32)
 TCTCGCTGCAAAAGCTGGAGGAAGCCTCACTTTTGT                                  <  1:1140254/36‑1 (MQ=32)
  CTCGCTGCAAAAGCTGGAGGAAGCCTCACTTTTGTA                                 <  2:2493685/36‑1 (MQ=35)
   TCGCTGCAAAAGCTGGAGGAAGCCTCACTTTTGTAC                                >  1:1278072/1‑36 (MQ=38)
   TCGCTGCAAAAGCTGGAGGAAGCCTCACTTTTGTAC                                >  1:1176597/1‑36 (MQ=38)
    CGCTGCAAAAGCTGGAGGAAGCCTCACTTTTGTACT                               <  1:3464218/36‑1 (MQ=38)
     GCTGCAAAAGCTGGAGGAAGCCTCACTTTTGTACTG                              >  2:127834/1‑36 (MQ=38)
     GCTGCAAAAGCTGGAGGAAGCCTCACTCTTGTACTG                              >  1:3234402/1‑36 (MQ=255)
      CTGCAAAAGCTGGAGGAAGCCTCACTTTTGTACTGA                             >  1:2794001/1‑36 (MQ=38)
      CTGCAAAAGATGGAGGAAGCCTCACTTTTGTACTGA                             >  1:1716309/1‑36 (MQ=255)
        GCAAAAGCTGGAGGAAGCCTCACTTTTGTACTGACC                           <  1:676952/36‑1 (MQ=38)
          aagAGCTGGAGGAAGCCTCACTTTTGTACTGACCCC                         <  1:652336/33‑1 (MQ=255)
          AAAAGCTGGATGAAGCCTCACTTTTGTACTGACCCC                         <  1:681680/36‑1 (MQ=255)
          AAAAGCTGGAGGAAGCCTCACTTTTGTACTGACCCC                         <  1:2299269/36‑1 (MQ=38)
           AAAGCTGGAGGAAGCCTCACTTTTGTACTGACCCCA                        >  2:1959424/1‑36 (MQ=38)
           AAAGCTGGAGGAAGCCTCACTTTTGTACTGACCCCA                        <  1:2612123/36‑1 (MQ=38)
            AAGCTGGAGGAAGCCTCACTTTTGTACTGACCCCAA                       >  2:1567932/1‑36 (MQ=40)
             AGCTGGAGGAAGCCTCACTTTAGTACTGACCCCAAA                      >  1:3570011/1‑36 (MQ=255)
              GCTGGAGGAAGCCTCACTTTTGTACTGACCCCAAAA                     >  1:2365250/1‑36 (MQ=255)
              GCTGGAGGAAGCCTCACTTTTGTACTGACCCCAAAA                     >  2:2477106/1‑36 (MQ=255)
                TGGAGGAAGCCTCACTTTTGTACTGACCCCAAAAAG                   >  1:3094531/1‑36 (MQ=255)
                 GGAGGAAGCCTCACTTTTGTACTGACCCCAAAAAGT                  >  2:3431683/1‑36 (MQ=255)
                 GGAGGAAGCCTCACTTTTGTACTGACCCCAAAAAGT                  >  1:3520078/1‑36 (MQ=255)
                    GGAAGCCTCACTTTTGTACTGACCCCAAAAAGTTGG               <  1:3223186/36‑1 (MQ=255)
                    GGAAGCCTCACTTTTGTACTGACCCCAAAAAGTTGG               >  1:2976673/1‑36 (MQ=255)
                     GAAGCCTCACTTTTGTACTGACCCCAAAACGTTGGA              >  2:756664/1‑36 (MQ=255)
                     GAAGCCTCACTTTTGTACTGACCCCAAAAAGTTGGA              >  2:3699850/1‑36 (MQ=255)
                     GAAGCCTCACTTTTGTACTGACCCCAAAAAGTTGGA              <  2:2550722/36‑1 (MQ=255)
                      AAGCCTCACTTTTGTACTGACCCCAAAAGGTTGGAC             <  2:1616879/36‑1 (MQ=255)
                       AGCCTCACTTTTGTACTGACCCCAACAAGTTGGACA            >  1:848799/1‑36 (MQ=255)
                       AGCCTCACTTTTGTACTGACCCCAAAAAGTTGGACA            >  2:1372317/1‑36 (MQ=255)
                        GCCTCACTTTTGTACTGACCCCAAAAAGTTGGACAG           >  2:3385812/1‑36 (MQ=255)
                        GCCTCACTTTTGTACTGACCCCAAAAAGTTGGACAG           >  1:544242/1‑36 (MQ=255)
                          CTCACTTTTGTACTGACCCCAAAAAGTTGGACAGTT         <  2:231592/36‑1 (MQ=255)
                          CTCACTTTTGTACTGACCCCAAAAAGTTGGACAGTT         >  1:2861377/1‑36 (MQ=255)
                          CTCACTTTTGTACTGACCCCAAAAAGTTGGACAGTT         <  2:1397468/36‑1 (MQ=255)
                           TCACTTTTGTACTGACCCCAAAAAGTTGGACAGTTA        >  1:2435280/1‑36 (MQ=255)
                            CACTTTTGTACTGACCCCAAAAAGTTGGACAGTTAA       >  1:577402/1‑36 (MQ=255)
                             ACTTTTGTACTGACCCCAAAAAGTTGGACAGTTAAA      >  2:1268307/1‑36 (MQ=40)
                             ACTTTTGTACTGACCCCAAAAAGTTGGACAGTTAAA      >  1:396217/1‑36 (MQ=40)
                             ACTTTTGTACTGACCCCAAAAAGTTGGACAGTTAAA      <  2:289993/36‑1 (MQ=40)
                              CTTTTGTACTGACCCCAAAAAGTTGGACAGTTAAAC     <  1:2692004/36‑1 (MQ=38)
                               TTTTGTACTGACCCCAAAAAGTTGGACAGTTAAACA    >  2:3061769/1‑36 (MQ=35)
                                TTTGTACTGACCCCAAAAAGTTGGACAGTTAAACAC   <  1:2969682/36‑1 (MQ=32)
                                 TTGTACTGACCCCAAAAAGTTGGACAGTTAAACACG  <  1:490564/36‑1 (MQ=1)
                                 TTGTACTGACCCCAAAAAGTTGGACAGTTAAACACG  <  1:263649/36‑1 (MQ=1)
                                 TTGTACTGACCCCAAAAAGTTGGACAGTTAAACACG  >  2:1255031/1‑36 (MQ=1)
                                 TTGTACTGACCCCAAAAAGTTGGACAGTTAAACACG  <  2:1101184/36‑1 (MQ=1)
                                 TTGTACTGACCCCAAAAAGTTGGACAGTTAAACACG  >  1:1916519/1‑36 (MQ=1)
                                 TTGTACTGACCCCAAAAAGTTGGACAGTTAAACACG  >  2:526166/1‑36 (MQ=1)
                                 TTGTACTGACCCCAAAAAGTTGGACAGTTAAACACG  >  2:1480648/1‑36 (MQ=1)
                                                                     
CTCTCGCTGCAAAAGCTGGAGGAAGCCTCACTTT...................................  >  REL606/3015741‑3015774
..................................TGTACTGACCCCAAAAAGTTGGACAGTTAAACACG  <  REL606/3894996‑3894962

Base quality scores:  ATCG  < 3 ≤  ATCG  < 11 ≤  ATCG  < 19 ≤  ATCG  < 30 ≤  ATCG  < 32 ≤  ATCG