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

Only 100 of 116 total aligned reads displayed.

ATCTCTCGCTGCAAAAGCTGGAGGAAGCCTCACTTTCTGCCCAGTCATCCACCAAAGAAGAAGCCAGCCTG  >  REL606/3015739‑3015809
                                   |                                   
aTCTCTCGCTGCAAAAGCTGGAGGAAGCCTCACttt                                     >  2:2378829/1‑36 (MQ=255)
aTCTCTCGCTGCAAAAGCTGGAGGAAGCCTCACttt                                     <  1:188746/36‑1 (MQ=255)
aTCTCTCGCTGCAAAAGCTGGAGGAAGCCTCACttt                                     >  1:3297775/1‑36 (MQ=255)
  ttcatcgtccagcaaaatgggtgcagtacaCTTTCt                                   <  1:3007280‑M2/6‑1 (MQ=255)
  ctctcGCTGCAAAAGCTGGAGGAAGCCTCACTTTtg                                   >  2:328517‑M1/1‑34 (MQ=255)
   tctcGCTGCAAAAGCTGGAGGAAGCCTCACTTTtgt                                  >  2:1469917‑M1/1‑33 (MQ=255)
   tctcGCTGCAAAAGCTGGAGGAAGCCTCACTTTtgt                                  <  1:1140254‑M1/36‑4 (MQ=255)
   tctcGCTGCAAAAGCTGGAGGAAGCCTCACTTTtgt                                  >  2:1914102‑M1/1‑33 (MQ=255)
   tctcGCTGCAAAAGCTGGAGGAAGCCTCACTTTtgt                                  >  2:372091‑M1/1‑33 (MQ=255)
   tctcGCTGCAAAAGCTGGAGGAAGCCTCACTTTtgt                                  <  2:1300651‑M1/36‑4 (MQ=255)
   tcatcgtccgacaaaatgggtgcagtacaCTTTCTg                                  >  1:1294855‑M2/30‑36 (MQ=255)
   tcatcgtccaacaaaatgggtgcagtacaCTTTCTg                                  <  2:1530587‑M2/7‑1 (MQ=255)
    ctcGCTGCAAAAGCTGGAGGAAGCCTCACTTTtgta                                 <  2:2493685‑M1/36‑5 (MQ=255)
    catcgtccaacaaaatgggtgcagtacaCTTTCTGc                                 <  1:3212360‑M2/8‑1 (MQ=255)
    catcgtccaacaaaatgggtgcagtacaCTTTCTGc                                 >  1:791289‑M2/29‑36 (MQ=255)
     tcGCTGCAAAAGCTGGAGGAAGCCTCACTTTtgtac                                >  1:1176597‑M1/1‑31 (MQ=255)
     tcGCTGCAAAAGCTGGAGGAAGCCTCACTTTtgtac                                >  1:1278072‑M1/1‑31 (MQ=255)
     atcgtccaacaaaatgggtgcagtacaCTTTCTGcc                                >  2:3663076‑M2/28‑36 (MQ=255)
      tcgtccaacaaaatgggtgcagtacaCTTTCTGccc                               >  1:2109187‑M2/27‑36 (MQ=255)
      cGCTGCAAAAGCTGGAGGAAGCCTCACTTTtgtact                               <  1:3464218‑M1/36‑7 (MQ=255)
       gCTGCAAAAGCTGGAGGAAGCCTCACTTTtgtactg                              >  2:127834‑M1/1‑29 (MQ=255)
        gtccaacaaaatgggtgccgtacaCTTTCTGCCCAg                             >  2:2007225‑M2/25‑36 (MQ=255)
        cTGCAAAAGCTGGAGGAAGCCTCACTTTtgtactga                             >  1:2794001‑M1/1‑28 (MQ=255)
         tccaacaaaatgggtgcagtacaCTTTCTGCCCAGt                            >  1:3464077‑M2/24‑36 (MQ=255)
         tccaacaaaatgggtgcagtacaCTTTCTGCCCAGt                            <  1:312984‑M2/13‑1 (MQ=255)
          gCAAAAGCTGGAGGAAGCCTCACTTTtgtactgacc                           <  1:676952‑M1/36‑11 (MQ=255)
          ccaacaaaatgggtgcagtacaCTTTCTGCCCAGTc                           <  2:1109873‑M2/14‑1 (MQ=255)
           caacaaaatgggtgcagtacaCTTTCTGCCCAGTCa                          >  2:2309774‑M2/22‑36 (MQ=255)
            aagagctGGAGGAAGCCTCACTTTtgtactgacccc                         <  1:652336‑M1/33‑13 (MQ=255)
            aacaaaatgggtgcagtacaCTTTCTGCCCAGTCAt                         <  1:1538542‑M2/16‑1 (MQ=255)
            aaaaGCTGGAGGAAGCCTCACTTTtgtactgacccc                         <  1:2299269‑M1/36‑13 (MQ=255)
             acaaaattggtgcagtacaCTTTCTGCCCAGTCATc                        <  2:173679‑M2/17‑1 (MQ=255)
             aaaGCTGGAGGAAGCCTCACTTTtgtactgacccca                        >  2:1959424‑M1/1‑23 (MQ=255)
             aaaGCTGGAGGAAGCCTCACTTTtgtactgacccca                        <  1:2612123‑M1/36‑14 (MQ=255)
              aaGCTGGAGGAAGCCTCACTTTtgtactgaccccaa                       >  2:1567932‑M1/1‑22 (MQ=255)
               aaaatgggtgcagtacaCTTTCTGCCCAGTTATcca                      >  1:3118564‑M2/18‑36 (MQ=255)
               aaaatgggtgcagtacaCTTTCTGCCCAGTCATcca                      >  1:3012721‑M2/18‑36 (MQ=255)
               aaaatgggtgcagtacaCTTTCTGCCCAGTCATcca                      >  1:2688366‑M2/18‑36 (MQ=255)
               aaaatgggtgcagtacaCTTTCTGCCCAGTCATcca                      >  1:26606‑M2/18‑36 (MQ=255)
               aGCTGGAGGAAGCCTCACTTTagtactgaccccaaa                      >  1:3570011‑M1/1‑21 (MQ=255)
                gCTGGAGGAAGCCTCACTTTtgtactgaccccaaaa                     >  2:2477106‑M1/1‑20 (MQ=255)
                gCTGGAGGAAGCCTCACTTTtgtactgaccccaaaa                     >  1:2365250‑M1/1‑20 (MQ=255)
                aaatgggtgcagtacaCTTTCTGCCCAGTCATccac                     <  1:1787012‑M2/20‑1 (MQ=255)
                aaaggggtgcagtacaCTTTCTGCCCAGTCATccac                     <  2:2072391‑M2/20‑1 (MQ=255)
                 aatgggcgcagtacaCTTTCTGCCCAGTCATccacc                    >  2:2628651‑M2/16‑36 (MQ=255)
                  atgggtgcagtacaCTTTCTGCCCAGTCatccacta                   >  1:2764810‑M2/15‑34 (MQ=255)
                   tgggtgcagtaccCTTTCTGCCCAGTCATCCACCaa                  >  2:1784357‑M2/14‑36 (MQ=255)
                   tgggtgcagtacaCTTTCTGCCCAGTCATCCACCaa                  <  2:430658‑M2/23‑1 (MQ=255)
                   ggaggaAGCCTCACTTTtgtactgaccccaaaaagt                  >  1:3520078‑M1/1‑17 (MQ=255)
                   ggaggaAGCCTCACTTTtgtactgaccccaaaaagt                  >  2:3431683‑M1/1‑17 (MQ=255)
                    gggtgcagtacaCTTTCTGCCCAGTCATCCACCaaa                 <  1:892269‑M2/24‑1 (MQ=255)
                     ggtgcagtacaCTTTCTGCCCAGTCATCCAccaacg                >  1:455081‑M2/12‑34 (MQ=255)
                      gtgcagtacaCTTTCTGCCCAGTCATCCACCAaaga               >  2:3676929‑M2/11‑36 (MQ=255)
                      gtgcagtacaCTTTCTGCCCAGTCATCCACCAaaga               >  1:764288‑M2/11‑36 (MQ=255)
                      gtgcagtacaCTTTCTGCCCAGTCATCCACCAaaga               >  1:1122133‑M2/11‑36 (MQ=255)
                      gtgcagtacaCTTTCTGCCCAGTCATCCACCAaaga               >  1:289823‑M2/11‑36 (MQ=255)
                      ggaAGCCTCACTTTtgtactgaccccaaaaagttgg               >  1:2976673‑M1/1‑14 (MQ=255)
                      ggaAGCCTCACTTTtgtactgaccccaaaaagttgg               <  1:3223186‑M1/36‑23 (MQ=255)
                       gaAGCCTCACTTTtgtactgaccccaaaacgttgga              >  2:756664‑M1/1‑13 (MQ=255)
                       gaAGCCTCACTTTtgtactgaccccaaaaagttgga              <  2:2550722‑M1/36‑24 (MQ=255)
                        gcagtacaCTTTCTGCCCAGTCATCCACCAaagaag             <  2:2675430‑M2/28‑1 (MQ=255)
                        aaGCCTCACTTTtgtactgaccccaaaaggttggac             <  2:1616879‑M1/36‑25 (MQ=255)
                         cagtacaCTTTCTGCCCAGTCATCCACCAaagaaga            <  1:91716‑M2/29‑1 (MQ=255)
                         cagtacaCTTTCTGCCCAGTCATCCACCAaagaaga            >  2:411315‑M2/8‑36 (MQ=255)
                         aGCCTCACTTTtgtactgaccccaacaagttggaca            >  1:848799‑M1/1‑11 (MQ=255)
                          gCCTCACTTTtgtactgaccccaaaaagttggacag           >  2:3385812‑M1/1‑10 (MQ=255)
                          gCCTCACTTTtgtactgaccccaaaaagttggacag           >  1:544242‑M1/1‑10 (MQ=255)
                          agtacaCTTTCTGCCCAGTCATCCACCAaagaagaa           <  1:3279116‑M2/30‑1 (MQ=255)
                            tacaCTTTCTGCCCAGTCATCCACCAAAGAAGAAGc         >  1:1779211‑M2/5‑36 (MQ=255)
                            tacaCTTTCTGCCCAGTCATCCACCAAAGAAGAAGc         >  1:3111387‑M2/5‑36 (MQ=255)
                            tacaCTTGCTGCCCAGTCATCCACCAAAGAAGAAGc         <  2:3561696‑M2/32‑1 (MQ=255)
                            cTCACTTTtgtactgaccccaaaaagttggacagtt         >  1:2861377‑M1/1‑8 (MQ=255)
                            cTCACTTTtgtactgaccccaaaaagttggacagtt         <  2:231592‑M1/36‑29 (MQ=255)
                            cTCACTTTtgtactgaccccaaaaagttggacagtt         <  2:1397468‑M1/36‑29 (MQ=255)
                             tCACTTTtgtactgaccccaaaaagttggacagtta        >  1:2435280‑M1/1‑7 (MQ=255)
                             acaCTTTCTGCCCAGTCATCCACCAAAGAAGAAGcc        <  1:2628850‑M2/33‑1 (MQ=255)
                             acaCTTTCTGCCCAGTCATCCACCAAAGAAGAAGcc        >  1:3060729‑M2/4‑36 (MQ=255)
                             acaCTTTCTGCCCAGTCATCCACCAAAGAAGAAGcc        <  1:1645864‑M2/33‑1 (MQ=255)
                             acACTATCAGCCCAGACATCCACCAAAGAAGAagcc        <  1:3209652/35‑1 (MQ=25)
                              caCTTTCTGCCCAGTCATCCACCAAAGACGAAGCCa       >  2:796419‑M2/3‑36 (MQ=255)
                              cACTTTtgtactgaccccaaaaagttggacagttaa       >  1:577402‑M1/1‑6 (MQ=255)
                               aCTTTtgtactgaccccaaaaagttggacagttaaa      >  2:1268307‑M1/1‑5 (MQ=255)
                               aCTTTtgtactgaccccaaaaagttggacagttaaa      <  2:289993‑M1/36‑32 (MQ=255)
                               aCTTTCTGCCCAGTCATCCACCAAAGAAGAAGCCAg      <  1:77937‑M2/35‑1 (MQ=255)
                               aCTTTCTGCCCAGTCATCCACCAAAGAAGAAGCCAg      >  2:2734471‑M2/2‑36 (MQ=255)
                               aCTTTCTGCCCAGTCATCCACCAAAGAAGAAGCCAg      >  1:1928816‑M2/2‑36 (MQ=255)
                                cTTTtgtactgaccccaaaaagttggacagttaaac     <  1:2692004‑M1/36‑33 (MQ=255)
                                cTTTCTGCCCAGTCATCCACCAAAGAAGAagccagc     >  2:1262244/1‑36 (MQ=255)
                                cTTTCTGCCCAGTCATCCACCAAAGAAGAagccagc     >  1:2504498/1‑36 (MQ=255)
                                cTTTCTGCCCAGTCATCCACCAAAGAAGAagccagc     <  1:2673107/36‑1 (MQ=255)
                                 ttttgtactgaccccaaaaagttggacagttaaaca    >  2:3061769‑M1/1‑3 (MQ=255)
                                  tttgtactgaccccaaaaagttggacagttaaacac   <  1:2969682‑M1/36‑35 (MQ=255)
                                  ttCTGCCCAGTCATCCACCAAAGAAGAAGCCAGCCt   <  2:2743274/36‑1 (MQ=255)
                                  ttCTGCCCAGTCATCCACCAAAGAAGAAGCCAGCCt   <  1:825625/36‑1 (MQ=255)
                                   ttgtactgaccccaaaaagttggacagttaaacacg  >  2:1480648‑M1/1‑1 (MQ=255)
                                   ttgtactgaccccaaaaagttggacagttaaacacg  <  1:490564‑M1/36‑36 (MQ=255)
                                   ttgtactgaccccaaaaagttggacagttaaacacg  >  1:1916519‑M1/1‑1 (MQ=255)
                                   ttgtactgaccccaaaaagttggacagttaaacacg  >  2:1255031‑M1/1‑1 (MQ=255)
                                   ttgtactgaccccaaaaagttggacagttaaacacg  >  2:526166‑M1/1‑1 (MQ=255)
                                   tCTGCCCAGTCATCCACCAAAGAAGAAGCCAGCCTg  <  2:2942323/36‑1 (MQ=255)
                                   |                                   
ATCTCTCGCTGCAAAAGCTGGAGGAAGCCTCACTTTCTGCCCAGTCATCCACCAAAGAAGAAGCCAGCCTG  >  REL606/3015739‑3015809

Base quality scores:  ATCG  < 3 ≤  ATCG  < 12 ≤  ATCG  < 20 ≤  ATCG  < 29 ≤  ATCG  < 32 ≤  ATCG