Predicted mutation
evidence seq id position mutation freq annotation gene description
JC JC REL606 1,821,525 IS150 (–) +3 bp 96.3% coding (573‑575/1041 nt) ynjI ← predicted inner membrane protein

New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? REL606 = 590471NA (NA)51 (0.900) 31/64 NT 96.2% noncoding (1977/1977 nt) IS150 repeat region
?REL606 = 1821527 2 (0.100)coding (573/1041 nt) ynjI predicted inner membrane protein
* ? REL606 1821525 =2 (0.040)54 (0.960) 37/64 NT 96.4% coding (575/1041 nt) ynjI predicted inner membrane protein
?REL606 3651091 = NA (NA)noncoding (1/1443 nt) IS150 repeat region

TTGTGCCTTAATGGAGCATTTTCTTTCTGGCGACCATG...................................  <  REL606/1821560‑1821525
...................................ATGTACTGCACCCATTTTGTTGGACGATGAAATGGAAT  >  REL606/3651091‑3651127
                                                                         
   TGCCTTAATGGAGCATTTTCTTTCTGGCGACCATGT                                    >  1:5320673/1‑36 (MQ=1)
   TGCCTTAATGGAGCATTTTCTTTCTGGCGACCATGT                                    <  1:8315221/36‑1 (MQ=1)
    GCCTTAATGGAGCATTTTCTTTCTGGCGACCATGTA                                   >  1:940919/1‑36 (MQ=1)
     CCTTAATGGAGCATTTTCTTTCTGGCGACCATGTAC                                  <  1:705822/36‑1 (MQ=1)
     CCTTAATGGAGCATTTTCTTTCTGGCGACCATGTAC                                  >  1:6625882/1‑36 (MQ=1)
     CCTTAATGGAGCATTTTCTTTCTGGCGACCATGTAC                                  <  1:5382006/36‑1 (MQ=1)
      CTTAATGGAGCATTTTCTTTCTGGCGACCATGTACT                                 <  1:6894526/36‑1 (MQ=1)
      CTTAATGGAGCATTTTCTTTCTGGCGACCATGTACT                                 >  1:3013390/1‑36 (MQ=1)
        TAATGGAGCATTTTCTTTCTGGCGACCATGTACTGC                               <  1:6629796/36‑1 (MQ=1)
        TAATGGAGCATTTTCTTTCTGGCGACCATGTACTGC                               >  1:5224545/1‑36 (MQ=1)
         AATGGAGCATTTTCTTTATGGCGACCATGTACTGCA                              >  1:2778848/1‑36 (MQ=0)
          ATGGAGCATTTTCTTTCTGGCGACCATGTACTGCAC                             >  1:5179441/1‑36 (MQ=1)
            GGAGCATTTTCTTTCTGGCGACCATGTACTGCACCC                           >  1:7004431/1‑36 (MQ=1)
             GAGCATTTTCTTTCTGGCGACCATGTACTGCACCCA                          <  1:5656971/36‑1 (MQ=1)
               GCATTTTCTTTCTGGCGACCATGTACTGCACCCATT                        <  1:7902749/36‑1 (MQ=35)
               GCATTTTCTTTCTGGCGACCATGTACTGCACCCATT                        <  1:6576740/36‑1 (MQ=35)
                CATTTTCTTTCTGGCGACCATGTACTGCACCCATTT                       <  1:5026240/36‑1 (MQ=38)
                CATTTTCTTTCTGGCGACCATGTACTGCACCCATTT                       <  1:5201900/36‑1 (MQ=38)
                 ATTTTCTTTCTGGCGACCATGTACTGCACCCATTTT                      <  1:993496/36‑1 (MQ=40)
                 ATTTTCTTTCTGGCGACCATGTACTGCACCCATTTT                      <  1:5787654/36‑1 (MQ=40)
                 ATTTTCTTTCTGGCGACCATGTACTGCACCCATTTT                      <  1:3888706/36‑1 (MQ=40)
                 ATTTTCTTTCTGGCGACCATGTACTGCACCCATTTT                      >  1:4868053/1‑36 (MQ=40)
                  TTTTCTTTCTGGCGACCATGTACTGCACCCATTTTG                     >  1:5054637/1‑36 (MQ=38)
                   TTTCTTTCTGGCGACCATGTACTGCACCCATTTTGT                    <  1:8317570/36‑1 (MQ=35)
                   TTTCTTTCTGGCGACCATGTACTGCACCCATTTTGT                    <  1:1072753/36‑1 (MQ=35)
                   TTTCTTTCTGGCGACCATGTACTGCACCCATTTTGT                    <  1:6811842/36‑1 (MQ=35)
                   TTTCTTTCTGGCGACCATGTACTGCACCCATTTTGT                    <  1:7409687/36‑1 (MQ=35)
                    TTCTTTCTGGCGACCATGTGCTGCACCCATTTTGTT                   <  1:127176/36‑1 (MQ=255)
                    TTCTTTCTGGCGACCATGTACTGCACCCATTTTGTT                   <  1:4579496/36‑1 (MQ=32)
                    TTCTTTCTGGCGACCATGTACTGCACCCATTTTGTT                   <  1:7900082/36‑1 (MQ=32)
                    TTCTTTCTGGCGACCATGTACTGCACCCATTTTGTT                   <  1:5104438/36‑1 (MQ=32)
                     TCTTTCTGGCGACCATGTACTGCACCCATTTTGTTG                  <  1:6701061/36‑1 (MQ=1)
                     TCTTTCTGGCGACCATGTACTGCACCCATTTTGTTG                  >  1:387179/1‑36 (MQ=1)
                      CTTTCTGGCGACCATGTACTGCGCCCATTTTGTTGG                 <  1:964202/36‑1 (MQ=0)
                      CTTTCTGGCGACCATGTACTGCACCCATTTTGTTGG                 <  1:1899667/36‑1 (MQ=1)
                       TTTCTGGCGACCATGTACTGCACCCATTTTGTTGGA                >  1:6437998/1‑36 (MQ=1)
                       TTTCTGGCGACCATGTACTGCACCCATTTTGTTGGA                <  1:7783581/36‑1 (MQ=1)
                       TTTCTGGCGACCATGTACTGCACCCATTTTGTTGGA                >  1:2674739/1‑36 (MQ=1)
                        TTCTGGCGACCATGTACTGCACCCATTTTGTTGGAC               >  1:7391684/1‑36 (MQ=1)
                        TTCTGGCGACCATGTACTGCACCCATTTTGTTGGAC               <  1:5893589/36‑1 (MQ=1)
                         TCTGGCGACCATGTACTGCACCCATTTTGTTGGACG              <  1:4008765/36‑1 (MQ=1)
                         TCTGGCGACCATGTACTGCACCCATTTTGTTGGACG              >  1:2246725/1‑36 (MQ=1)
                             GCGACCATGTACTGCGCCCATTTTGTTGGACGATGA          <  1:8456442/36‑1 (MQ=0)
                             GCGACCATGTACTGCACCCATTTTGTTGGACGATGA          >  1:8192928/1‑36 (MQ=1)
                              acACCATGTACTGCACCCATTTTGTTGGACGATGAA         >  1:5603750/3‑36 (MQ=1)
                              CGACCATGTACTGCACCCATTTTGTTGGACGATGAA         <  1:6347382/36‑1 (MQ=1)
                              CGACCATGTACTGCACCCATTTTGTTGGACGATGAA         >  1:811096/1‑36 (MQ=1)
                              CGACCATGTACTGCACCCATTTTGTTGGACGATGAA         <  1:1036401/36‑1 (MQ=1)
                               GACCATGTACTGCACCCATTTTGTTGGACGATGAAA        >  1:7932173/1‑36 (MQ=1)
                               GACCATGTACTGCACCCATTTTGTTGGACGATGAAA        <  1:8033425/36‑1 (MQ=1)
                               GACCATGTACTGCACCCATTTTGTTGGACGATGAAA        <  1:294064/36‑1 (MQ=1)
                                ACCATGTACTGCACCCATTTTGTTGGACGATGAAAT       >  1:4412053/1‑36 (MQ=1)
                                ACCATGTACTGCACCCATTTTGTTGGACGATGAAAT       >  1:1495376/1‑36 (MQ=1)
                                  CATGTACTGCACCCATTTTGTTGGACGATGAAATGG     >  1:7538639/1‑36 (MQ=1)
                                                                         
TTGTGCCTTAATGGAGCATTTTCTTTCTGGCGACCATG...................................  <  REL606/1821560‑1821525
...................................ATGTACTGCACCCATTTTGTTGGACGATGAAATGGAAT  >  REL606/3651091‑3651127

Base quality scores:  ATCG  < 3 ≤  ATCG  < 17 ≤  ATCG  < 25 ≤  ATCG  < 32 ≤  ATCG  < 37 ≤  ATCG