Predicted mutation
evidence seq id position mutation annotation gene description
MC JC REL606 3,894,997 Δ6,934 bp IS150‑mediated rbsD[yieO] rbsD, rbsA, rbsC, rbsB, rbsK, rbsR, [yieO]

Missing coverage evidence...
   seq id start end size ←reads reads→ gene description
* * ÷ REL606 3893578–3894995 3901930 6936–8353 16 [14] [1] 40 insJ‑5–[yieO] insJ‑5, insK‑5, rbsD, rbsA, rbsC, rbsB, rbsK, rbsR, [yieO]

New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? REL606 = 3652533NA (NA)34 (0.650) 23/66 0.9 97.1% noncoding (1443/1443 nt) IS150 repeat region
?REL606 3901931 = 1 (0.020)coding (489/1428 nt) yieO predicted multidrug or homocysteine efflux system

TGTTTAACTGTCCAACTTTTTGGGGTCAGTACAAACTTTGATTTATAGTCAGGTGGGGCTTTTCTCTG  >  REL606/3652501‑3652568
                                |                                   
tgtTTAACTGTCCAACTTTTTGGGGTCAGTACAAAc                                  >  2:3264113/1‑36 (MQ=35)
 gtTTAACTGTCCAACTTTTTGGGGTCAGTACAAACt                                 >  1:1061626/1‑36 (MQ=38)
  tttAACTGTCCAACTTTTTGGGGTCAGTACAAACtt                                >  2:708158/1‑36 (MQ=40)
  tttAACTGTCCAACTTTTTGGGGTCAGTACAAACtt                                <  1:25821/36‑1 (MQ=40)
   ttAACTGTCCAACTTTTTGGGGTCAGTACAAACttt                               <  1:3608504/36‑1 (MQ=255)
    tAACTGTCCAACTTTTTGGGGTCAGTACAAACTTTg                              <  1:1872431/36‑1 (MQ=255)
    tAACTGTCCAACTTTTTGGGGTCAGTACAAACTTTg                              <  1:2695321/36‑1 (MQ=255)
    tAACTGTCCAACTTTTTGGGGTCAGTACAAACTTTg                              >  1:532705/1‑36 (MQ=255)
     aaCTGTCCAACTTTTTGGGGTCAGTACAAACTTTGa                             >  2:206615/1‑36 (MQ=255)
      aCTGTCCAACTTTTTGGGGTCAGTACAAACTTTGAt                            <  2:4323/36‑1 (MQ=255)
       cTGTCCAACTTTTTGGGGTCAGTACAAACTTTGAtt                           <  1:2221190/36‑1 (MQ=255)
        tGTCCAACTTTTTGGGGTCAGTACAAACTTTGAttt                          <  1:3350917/36‑1 (MQ=255)
        tGTCCAACTTTTTGGGGTCAGTACAAACTTTGAttt                          <  2:2545935/36‑1 (MQ=255)
          tcgaaCTTTTTGGGGTCAGTACAAACTTTGATTtat                        <  1:3240468/33‑1 (MQ=255)
           ccAACTTTTTGGGGTCAGTACAAACTTTGATTtata                       >  2:152383/1‑36 (MQ=255)
            cAACTTTTTGGGGTCAGTACAAACTTTGATTTATAg                      <  2:417568/36‑1 (MQ=255)
            cAACTTTTTGGGGTCAGTACAAACTTTGATTTATAg                      >  2:2821518/1‑36 (MQ=255)
            cAACTTTTTGGGGTCAGTACAAACTTTGACTTATAg                      >  2:2061094/1‑36 (MQ=255)
            cAACTGGTTGGGGTCAGTACAAACTTTGATTTATAg                      <  2:958606/36‑1 (MQ=37)
             aaCTTTTTGGGGTCAGTACAAACTTTGATTTATAGt                     <  1:469920/36‑1 (MQ=255)
             aaCTTTTTGGGGTCAGTACAAACTTTGATTTATAGt                     >  2:2168476/1‑36 (MQ=255)
             aaCTTTTTGGGGTCAGTACAAACTTTGATTTATAGt                     <  1:1600294/36‑1 (MQ=255)
             aaCTTTTTGGGGTCAGTACAAACTTTGATTTATAGt                     <  2:979637/36‑1 (MQ=255)
              aCTTTTTGGGGTCAGTACAAACTTTGATTTATAGTc                    <  1:2659049/36‑1 (MQ=255)
              aCTTTTTGGGGTCAGTACAAACTTTGATTTATAGTc                    <  2:2801615/36‑1 (MQ=255)
              aCTTTTTGGGGTCAGTACAAACTTTGATTTATAGTc                    >  1:2224177/1‑36 (MQ=255)
              aCTTTTTGGGGTCAGTACAAACTATGATTTATAGTc                    >  1:880788/1‑36 (MQ=255)
                tttttGGGGTCAGTACAAACTTTGATTTATAGTCAg                  >  2:3329835/1‑36 (MQ=255)
                tttttGGGGTCAGTACAAACTTTGATTTATAGTCAg                  <  2:1688664/36‑1 (MQ=255)
                tttttGGGGTCAGTACAAACTTTGATTTATAGTCAg                  >  2:1132247/1‑36 (MQ=255)
                 ttttGGGGTCAGTACAAACTTTGATTTATAGTCAgg                 >  1:2208984/1‑36 (MQ=255)
                  tttGGGGTCAGTACAAACTTTGATTTATAGTCAGGt                <  1:1789773/36‑1 (MQ=255)
                  tttGGGGTCAGTACAAACTTTGATTTATAGTCAGGt                <  2:1319288/36‑1 (MQ=255)
                  tttGGGGTCAGTACAAACTTTGATTTATAGTCAGGt                <  1:3732144/36‑1 (MQ=255)
                  tttGGGGTCAGTACAAACTTTGATTTATAGTCAGGt                <  2:3085013/36‑1 (MQ=255)
                     ggggTCAGTACAAACTTTGATTTATAGTCAGGTggg             <  2:3567845/36‑1 (MQ=255)
                       ggTCAGTACAAACTTTGATTTATAGTCAGGTGGGGc           <  1:1028952/36‑1 (MQ=255)
                         tCAGTACAAACTTTGATTTATAGCCAGGTGGGGCtt         >  2:2083782/1‑36 (MQ=255)
                          cAGTACAAACTTTGATTTATAGTCAGGTGGGGCttt        <  2:1858197/36‑1 (MQ=255)
                          cAGTACAAACTTTGATTTATAGTCAGGTGGGGCttt        <  1:152372/36‑1 (MQ=255)
                           aGTACAAACTTTGATTTATAGTCAGGTGGGGCtttt       <  2:2972071/36‑1 (MQ=255)
                            gTACAAACTTTGATTTATAGTCAGGTGGGGCTTTtc      <  2:2052653/36‑1 (MQ=255)
                             tACAAACTTTGATTTATAGTCAGGTGGGGCTTTtct     <  1:2600592/36‑1 (MQ=255)
                              aCAAACTTTGATTTATAGTCAGGTGGGGCTTTtctc    <  1:3009665/36‑1 (MQ=255)
                                aaaCTTTGATTTATAGTCAGGTGGGGCTTTTctctt  >  1:2599701/1‑35 (MQ=255)
                                aaaCTTTGATTTATAGTCAGGTGGGGATTTTCtttg  >  1:2916888/1‑36 (MQ=37)
                                |                                   
TGTTTAACTGTCCAACTTTTTGGGGTCAGTACAAACTTTGATTTATAGTCAGGTGGGGCTTTTCTCTG  >  REL606/3652501‑3652568

Base quality scores:  ATCG  < 3 ≤  ATCG  < 13 ≤  ATCG  < 24 ≤  ATCG  < 30 ≤  ATCG  < 32 ≤  ATCG