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

CATTTCATCGTCCAACAAAATGGGTGCAGTACAGGTTGATAGAAATATCGCCTGGGGCTTTTGTCCAT  >  REL606/2775845‑2775912
                                |                                   
cATTTCATCGTCCAACAAAATGGTAGCAGTACAGGt                                  >  1:3045433/1‑36 (MQ=11)
      aTCGTGCGACAAAATGGGTGCAGTACAGGTTGATAg                            <  2:833501/36‑1 (MQ=18)
      aTCGTCCAACAAAATGGGGGTAGTACAGGTCGAttg                            >  2:485289/1‑34 (MQ=14)
       tCGTCCAACAAAATGGGTGCAGTACAGGTTGATAGa                           <  1:3100740/36‑1 (MQ=255)
        cGTCCATCATGATGGGTGCCGTGCAGGTTGATAgta                          <  1:1463930/36‑3 (MQ=9)
        cGTCCAACAAAATGGGTGCAGTACAGGTTGATAGaa                          >  2:1264660/1‑36 (MQ=255)
         gTCCAACAAAATGGGTGCAGTACAGGTTGATAGaaa                         <  1:156173/36‑1 (MQ=255)
         gTCCAACAAAATGGGTGCAGTACAGGTTGATAGaaa                         <  2:2079666/36‑1 (MQ=255)
          tCCTACAGAATGGGTGCAGTACAGGTTGATAGAAat                        <  2:3580019/36‑1 (MQ=25)
          tCCAACAAAATGGGTGCAGTACAGGTTGATAGAAat                        <  1:3333768/36‑1 (MQ=255)
            caacaaAATGGGTGCAGTACAGGTTGATAGAAatat                      <  2:1707059/36‑1 (MQ=255)
            caacaaAATGGGTGCAGTACAGGTTGATAGAAatat                      <  2:2629300/36‑1 (MQ=255)
             aacaaAATGGGTGCAGTACAGGTTGATAGAAATATc                     >  1:1843326/1‑36 (MQ=255)
             aacaaAATGGGTGCAGTACAGGTTGATAGAAATATc                     >  1:894172/1‑36 (MQ=255)
             aacaaAATGGGTGCAGTACAGGTTGATAGAAACATc                     >  2:578100/1‑36 (MQ=255)
               caaAATGGGTGCAGTACAGGTTGATAGAAATATCGc                   <  2:1648107/36‑1 (MQ=255)
                aaaaTGGGTGCAGTACAGGTTGATAGAAATATCGcc                  >  1:1833756/1‑36 (MQ=255)
                aaaaTGGGTGCAGTACAGGTTGATAGAAATATCAcc                  >  1:101053/1‑36 (MQ=255)
                  aaTGGGTGCAGTACAGGTTGATAGAAATATCGCCTg                <  1:3359118/36‑1 (MQ=255)
                  aaTGGGTGCAGTACAGGTTGATAGAAATATCGCCTg                >  2:2549125/1‑36 (MQ=255)
                    tGGGTGCAGTACAGGTTGATAGAAATATCGCCTggg              >  1:1873813/1‑36 (MQ=255)
                     gggTGCAGTACAGGTTGATAGAAATATCGCCTgggg             <  2:3509995/36‑1 (MQ=255)
                        tGCAGTACAGGTTGATAGAAATATCGCCTGGGGCtt          <  2:3177660/36‑1 (MQ=255)
                         gCAGTGCTGTTTGATTGAAATATCGCCTGGGGCttt         <  2:3049427/36‑1 (MQ=19)
                         gCAGTACAGGTTGATAGAAATATCGCCTGGGGCttt         <  1:848925/36‑1 (MQ=255)
                         gCAGTACAGGTTGATAGAAATATCGCCTGGGGCttt         <  1:2221371/36‑1 (MQ=255)
                           aGTACAGGTTGATAGAAATATCGCCTGGGGGTTTTg       >  2:258623/1‑36 (MQ=255)
                           aGTACAGGTTGATAGAAATATCGCCTGGGGCTTTTg       >  1:1149831/1‑36 (MQ=255)
                            gTACATGTTGATAGAAATATCGCCTGGGGCTTTTGt      <  1:374849/36‑1 (MQ=255)
                            gTACAGGTTGATAGAAATATCGCCTGGGGCTTTTGt      <  1:3434025/36‑1 (MQ=255)
                               cAGGTTGATAGAAATATCGCCTGGGGCTTTTTTCCa   >  1:2855589/1‑36 (MQ=255)
                               cAGGTTGATAGAAATATCGCCTGGGGCTTTTGTCCa   >  1:1452388/1‑36 (MQ=255)
                                aGGTTGATAGAAATATCGCCTGGGGCTTTTGTCCAt  >  2:1205291/1‑36 (MQ=255)
                                |                                   
CATTTCATCGTCCAACAAAATGGGTGCAGTACAGGTTGATAGAAATATCGCCTGGGGCTTTTGTCCAT  >  REL606/2775845‑2775912

Base quality scores:  ATCG  < 3 ≤  ATCG  < 8 ≤  ATCG  < 18 ≤  ATCG  < 29 ≤  ATCG  < 32 ≤  ATCG