Predicted mutation
evidence seq id position mutation freq annotation gene description
JC JC REL606 1,821,525 IS150 (–) +3 bp 95.0% 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)41 (1.000) 28/64 NT 89.1% noncoding (1977/1977 nt) IS150 repeat region
?REL606 = 1821527 5 (0.110)coding (573/1041 nt) ynjI predicted inner membrane protein
* ? REL606 1821525 =0 (0.000)55 (1.340) 33/64 NT 100% coding (575/1041 nt) ynjI predicted inner membrane protein
?REL606 3651091 = NA (NA)noncoding (1/1443 nt) IS150 repeat region

CGTGTTTAACTGTCCAACTTTTTGGGGTCAGTACACC...................................  >  REL606/590437‑590471
..................................ACCATGGAGCCAATACCTTTAAAGCAAACATGGAGGAA  <  REL606/1821527‑1821491
                                                                        
   GTTTAACTGTCCAACTTTTTGGGGTCAGTACACCAT                                   >  1:5480529/1‑36 (MQ=1)
   GTTTAACTGTCCAACTTTTTGGGGTCAGTACACCAT                                   <  1:5345524/36‑1 (MQ=1)
   GTTTAACTGTCCAACTTTTTGGGGTCAGTACACCAT                                   >  1:660915/1‑36 (MQ=1)
    TTTAACTGTCCAACTTTTTGGGGTCAGTACCCCATG                                  <  1:5251465/36‑1 (MQ=0)
    TTTAACTGTCCAACTTTTTGGGGTCAGTACACCATG                                  <  1:1599675/36‑1 (MQ=1)
    TTTAACTGTCCAACTTTTTGGGGTCAGTACACCATG                                  <  1:3908495/36‑1 (MQ=1)
    TTTAACTGTCCAACTTTTTGGGGTCAGTACACCATG                                  >  1:3574340/1‑36 (MQ=1)
      TAACTGTCCAACTTTTTGGGGTCAGTACACCATGGA                                >  1:5073458/1‑36 (MQ=1)
      TAACTGTCCAACTTTTTGGGGTCAGTACACCATGGA                                <  1:310714/36‑1 (MQ=1)
       AACTGTCCAACTTTTTGGGGTCAGTACACCATGGAG                               <  1:3146714/36‑1 (MQ=1)
        ACTGTCCAACTTTTTGGGGTCAGTACACCATGGAGC                              >  1:4599144/1‑36 (MQ=1)
        ACTGTCCAACTTTTTGGGGTCAGTACACCATGGAGC                              <  1:3980702/36‑1 (MQ=1)
        ACTGTCCAACTTTTTGGGGTCAGTACACCATGGAGC                              <  1:968183/36‑1 (MQ=1)
          TGTCCAACTTTTTGGGGTCAGTACACCATGGAGCCA                            >  1:1332250/1‑36 (MQ=1)
           GTCCAACTTTTTGGGGTCAGTACACCATGGAGCCAA                           <  1:5928682/36‑1 (MQ=1)
             CCAACTTTTTGGGGTCAGTACACCATGGAGCCAATA                         <  1:5888107/36‑1 (MQ=1)
             CCAACTTTTTGGGGTCAGTACACCATGGAGCCAATA                         <  1:4370302/36‑1 (MQ=1)
             CCAACTTTTTGGGGTAAGTACACCATGGAGCCAATA                         <  1:4552603/36‑1 (MQ=0)
                 CTTTTTGGGGTCAGTACACCATGGAGCCAATACCTT                     <  1:244038/36‑1 (MQ=38)
                  TTTTTGGGGTCAGTACACCATGGAGCCAATACCTTT                    <  1:3116255/36‑1 (MQ=38)
                   TTTTGGGGTCAGTACACCATGGAGCCAATACCTTTA                   >  1:3329205/1‑36 (MQ=38)
                   TTTTGGGGTCAGTACACCATGGAGCCAATACCTTTA                   <  1:3203603/36‑1 (MQ=38)
                   TTTTGGGGTCAGTACACCATGGAGCCAATACCTTTA                   <  1:403225/36‑1 (MQ=38)
                     TTGGGGTCAGTACACCATGGAGCCAATACCTTTAAA                 <  1:813051/36‑1 (MQ=40)
                     TTGGGGTCAGTACACCATGGAGCCAATACCTTTAAA                 >  1:5524429/1‑36 (MQ=40)
                       GGGGTCAGTACACCATGGAGCCAATACCTTTAAAGC               >  1:1455654/1‑36 (MQ=35)
                       GGGGTCAGTACACCATGGAGCCAATACCTTTAAAGC               <  1:855712/36‑1 (MQ=35)
                       GGGGTCAGTACACCATGGAGCCAATACCTTTAAAGC               >  1:662170/1‑36 (MQ=35)
                        GGGTCAGTACACCATGGAGCCAATACCTTTAAAGCA              >  1:1207454/1‑36 (MQ=32)
                        GGGTCAGTACACCATGGAGCCAATACCTTTAAAGCA              >  1:5744827/1‑36 (MQ=32)
                         GGTCAGTACACCATGGAGCCAATACCTTTAAAGCAA             >  1:4824470/1‑36 (MQ=1)
                         GGTCAGTACACCATGGAGCCAATACCTTTAAAGCAA             <  1:825250/36‑1 (MQ=1)
                          gtaAGTACACCATGGAGCCAATACCTTTAAAGCAAA            <  1:3887684/33‑1 (MQ=0)
                          GTCAGTACACCATGGAGCCAATACCTTTAAAGCAAA            <  1:4495235/36‑1 (MQ=1)
                          GTCAGTACACAATGGAGCCAATACCTTTAAAGCAAA            <  1:3763532/36‑1 (MQ=0)
                            CAGTACACCATGGAGCCAATACCTTTAAAGCAAACA          <  1:306759/36‑1 (MQ=1)
                              GTACACCATGGAGCCAATACCTTTAAAGCAAACATG        >  1:5644853/1‑36 (MQ=1)
                              GTACACCATGGAGCCAATACCTTTAAAGCAAACATG        >  1:5103398/1‑36 (MQ=1)
                                ACACCATGGAGCCAATACCTTTAAAGCAAACATGGA      <  1:2859617/36‑1 (MQ=1)
                                ACACCATGGAGCCAATACCTTTAAAGCAAACATGGA      <  1:1160114/36‑1 (MQ=1)
                                 CACCATGGAGCCAATACCTTTAAAGCAAACATGGAG     >  1:4906228/1‑36 (MQ=1)
                                                                        
CGTGTTTAACTGTCCAACTTTTTGGGGTCAGTACACC...................................  >  REL606/590437‑590471
..................................ACCATGGAGCCAATACCTTTAAAGCAAACATGGAGGAA  <  REL606/1821527‑1821491

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