Predicted mutation
evidence seq id position mutation freq annotation gene description
JC JC REL606 3,015,771 IS150 (–) +3 bp 97.6% coding (1192‑1194/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)109 (1.160) 37/44 NT 97.9% noncoding (1/1443 nt) IS150 repeat region
?REL606 3015771 = 2 (0.020)coding (1192/1677 nt) ECB_02816 KpsD protein
* ? REL606 = 30157733 (0.040)97 (0.980) 40/46 NT 96.6% coding (1194/1677 nt) ECB_02816 KpsD protein
?REL606 = 3894996 NA (NA)noncoding (1443/1443 nt) IS150 repeat region

Only 100 of 308 total aligned reads displayed.

GAATCTCTCGCTGCAAAAGCTGGAGGAAGCCTCACTTTCTGCCCAGTCATCCACCAAAGAAGAAGCCAG  >  REL606/3015737‑3015805
                                  |                                  
gAATCTCTCGCTGCAAAAGCTGGAGGAAGCCTCAc                                    <  1:18869849/35‑1 (MQ=255)
gAATCTCTCGCTGCAAAAGCTGGAGGAAGCCTCAc                                    <  1:2246633/35‑1 (MQ=255)
  aTCTCTCGCTGCAAAAGCTGGAGGAAGCCTCACtt                                  <  1:1934452/35‑1 (MQ=255)
  aTCTCTCGCTGCAAAAGCTGGAGGAAGCCTCACtt                                  >  1:512730/1‑35 (MQ=255)
  aTCTCTCGCTGCAAAAGCTGGAGGAAGCCTCACtt                                  <  1:17089884/35‑1 (MQ=255)
    ttcatcgtccaacaaaatgggtgcagtacaCTTTc                                >  1:10585671‑M2/31‑35 (MQ=255)
     tcatcgtccaacaagatgggtgcagtacaCTTTCt                               >  1:8237403‑M2/30‑35 (MQ=255)
     tcatcgtccaacaaaatgggtgcagtacaCTTTCt                               >  1:12178769‑M2/30‑35 (MQ=255)
     tcatcgtccaacaaaatgggtgcagtacaCTTTCt                               <  1:9103824‑M2/6‑1 (MQ=255)
      catcgtccaacaaaatgggtgcagtacaCTTTCTg                              <  1:20125633‑M2/7‑1 (MQ=255)
      catcgtccaacaaaatgggtgcagtacaCTTTCTg                              <  1:19019247‑M2/7‑1 (MQ=255)
      catcgtccaacaaaatgggtgcagtacaCTTTCTg                              >  1:8322173‑M2/29‑35 (MQ=255)
      catcgtccaacaaaatgggtgcagtacaCTTTCTg                              >  1:8960367‑M2/29‑35 (MQ=255)
      catcgtccaacaaaatgggtgcagtacaCTTTCTg                              <  1:12350698‑M2/7‑1 (MQ=255)
      catcgtccaacaaaatgggtgcagtacaCTCTctt                              <  1:10129657‑M2/7‑2 (MQ=255)
       tcGCTGCAAAAGCTGGAGGAAGCCTCACTTtgtac                             >  1:17982405‑M1/1‑30 (MQ=255)
       tcGCTGCAAAAGCTGGAGGAAGCCTCACTTtgtac                             <  1:3346154‑M1/35‑6 (MQ=255)
       atcgtccaacaaaatgggtgcagtacaCTTTCTGc                             >  1:18142849‑M2/28‑35 (MQ=255)
        tcgtccaacaaaatgggtgcagtacaCTTTCTGcc                            <  1:20611137‑M2/9‑1 (MQ=255)
        cGCTGCAAAAGCTGGAGGAAGCCTCACTTtgtact                            <  1:10612558‑M1/35‑7 (MQ=255)
        cGCTGCAAAAGCTGGAGGAAGCCTCACTTtgtact                            <  1:12609396‑M1/35‑7 (MQ=255)
        cGCTGCAAAAGCTGGAGGAAGCCTCACTTtgtact                            <  1:9338198‑M1/35‑7 (MQ=255)
         gCTGCAAAAGCTGGAGGAAGCCTCACTTtgtactg                           <  1:13367379‑M1/35‑8 (MQ=255)
         cgtccaacaaaatgggtgcagtacaCTTTCTGccc                           >  1:14605376‑M2/26‑35 (MQ=255)
           tccaacaaaatgggtgcagtacaCTTTCTGCCCAg                         <  1:11039746‑M2/12‑1 (MQ=255)
           tGCAAAAGCTGGAGGAAGCCTCACTTtgtactgac                         >  1:2810173‑M1/1‑26 (MQ=255)
           tGCAAAAGCTGGAGGAAGCCTCACTTtgtactgac                         <  1:7783829‑M1/35‑10 (MQ=255)
            gCAAAAGCTGGAGGAAGCCTCACTTtgtactgacc                        <  1:15532594‑M1/35‑11 (MQ=255)
            gCAAAAGCTGGAGGAAGCCTCACTTtgtactgacc                        >  1:11668331‑M1/1‑25 (MQ=255)
             caacaaaatgggtgcagtacaCTTTCTGCCCAGTc                       <  1:10430561‑M2/14‑1 (MQ=255)
             cAAAAGCTGGAGGAAGCCTCACTTtgtactgaccc                       >  1:12931173‑M1/1‑24 (MQ=255)
             cAAAAGCTGGAGGAAGCCTCACTTtgtactgaccc                       >  1:5352437‑M1/1‑24 (MQ=255)
             cAAAAGCTGGAGGAAGCCTCACTTtgtactgaccc                       >  1:11189305‑M1/1‑24 (MQ=255)
             cAAAAGCTGGAGGAAGCCTCACTTtgtactgaccc                       >  1:17540841‑M1/1‑24 (MQ=255)
              aacaaaatgggtgcagtacaCTTTCTGCCCAGTCa                      >  1:7545592‑M2/21‑35 (MQ=255)
              aaaaGCTGGAGGAAGCCTCACTTtgtactgacccc                      <  1:12251268‑M1/35‑13 (MQ=255)
              aaaaGCTGGAGGAAGCCTCACTTtgtactgacccc                      <  1:4499658‑M1/35‑13 (MQ=255)
               acaaaatgggtgcagtacaCTTTCTGCCCAGTCAt                     >  1:19286715‑M2/20‑35 (MQ=255)
               aaaGCTGGAGGAAGCCTCACTTtgtactgacccca                     <  1:621475‑M1/35‑14 (MQ=255)
               aaaGCTGGAGGAAGCCTCACTTtgtactgacccca                     <  1:5961968‑M1/35‑14 (MQ=255)
               aaaGCTGGAGGAAGCCTCACTTtgtactgacccca                     >  1:20249177‑M1/1‑22 (MQ=255)
               aaaGCGGGAGGAAGCCTCACTTtgtactgacccca                     <  1:20401881‑M1/35‑14 (MQ=255)
                caaaatgggtgcagtacaCTTTCTGCCCAGTCATc                    <  1:6658111‑M2/17‑1 (MQ=255)
                aaGCTGGAGGAAGCCTCACTTtgtactgaccccaa                    <  1:19489234‑M1/35‑15 (MQ=255)
                 aaaatgtgtgcagtacaCTTTCTGCCCAGTCATcc                   >  1:15094267‑M2/18‑35 (MQ=255)
                 aaaatgggtgcagtacaCTTTCTGCCCAGTCATcc                   <  1:13671114‑M2/18‑1 (MQ=255)
                 aaaatgggtgcagtacaCTTTCTGCCCAGTCATcc                   <  1:5710199‑M2/18‑1 (MQ=255)
                 aaaatgggtgcagtacaCTTTCTGCCCAGTCATcc                   >  1:9695032‑M2/18‑35 (MQ=255)
                 aaaatgggtgcagtacaCTTTCTGCCCAGTCATcc                   >  1:14968577‑M2/18‑35 (MQ=255)
                 aaaatgggtgcagtacaCTTTCTGCACAGTCATcc                   <  1:7277826‑M2/18‑1 (MQ=255)
                 aaaatgggtgcagtacaCTTTCTGCACAGTCATcc                   <  1:383443‑M2/18‑1 (MQ=255)
                 aGCTGGAGGAAGCCTCACTTtgtactgaccccaaa                   <  1:17429807‑M1/35‑16 (MQ=255)
                  gCTGGAGGAAGCCTCACTTtgtactgaccccaaaa                  <  1:12468061‑M1/35‑17 (MQ=255)
                  gCTGGAGGAAGCCTCACTTtgtactgaccccaaaa                  >  1:17809497‑M1/1‑19 (MQ=255)
                  aaatgggtgcagtacaCTTTCTGCCCAGTCATcca                  <  1:19625677‑M2/19‑1 (MQ=255)
                   gatgggtgcagtacaCTTTCTGCCCAGTCATccac                 <  1:163360‑M2/20‑1 (MQ=255)
                   cTGGAGGAAGCCTCACTTtgtactgaccccaaaaa                 <  1:16849819‑M1/35‑18 (MQ=255)
                   cTGGAGGAAGCCTCACTTtgtactgaccccaaaaa                 >  1:19997188‑M1/1‑18 (MQ=255)
                   aattggtgcagtacaCTTTCTGCCCAGTCATccac                 >  1:7176387‑M2/16‑35 (MQ=255)
                   aatgggtgcagtacaCTTTCTGCCCAGTCATccac                 >  1:6693452‑M2/16‑35 (MQ=255)
                    tGGAGGAAGCCTCACTTtgtactgaccccaaaaag                >  1:4751912‑M1/1‑17 (MQ=255)
                     tgggtgcagtacaCTTTCTGCCCAGTCATccacca               >  1:18611485‑M2/14‑35 (MQ=255)
                     tgggtgcagtacaCTTTCTGCCCAGTCATccacca               >  1:11417250‑M2/14‑35 (MQ=255)
                      gggtgcagtacaCTTTCTGCCCAGTCATCCACCaa              >  1:9573095‑M2/13‑35 (MQ=255)
                       ggtgcagtacaCTTTCTGCCCAGTCATCCACCaaa             <  1:7319094‑M2/24‑1 (MQ=255)
                       ggtgcagtacaCTTTCTGCCCAGTCATCCACCaaa             >  1:2351296‑M2/12‑35 (MQ=255)
                       ggtgcagtacaCTTTCTGCCCAGTCATCCACCaaa             >  1:16040514‑M2/12‑35 (MQ=255)
                       ggtgcagtacaCTTTCTGCCCAGTCATCCACCaaa             >  1:10977020‑M2/12‑35 (MQ=255)
                         gaAGCCTCACTTtgtactgaccccaaaaagttgga           >  1:9969419‑M1/1‑12 (MQ=255)
                          gcagtacaCTTTCTGCCCAGTCATCCACCAaagaa          >  1:2393980‑M2/9‑35 (MQ=255)
                          aaGCCTCACTTtgtactgaccccaaaaagttggac          <  1:1284857‑M1/35‑25 (MQ=255)
                          aaGCCTCACTTtgtactgaccccaaaaagttggac          <  1:14048265‑M1/35‑25 (MQ=255)
                           cagtacaCTTTCTGCCCAGTCATCCACCAaagaag         >  1:17722129‑M2/8‑35 (MQ=255)
                           cagtacaCTTTCTGCCCAGTCATCCACCAaagaag         >  1:20451152‑M2/8‑35 (MQ=255)
                           aGCCTCACTTtgtactgaccccaaaaagttggaca         <  1:9754955‑M1/35‑26 (MQ=255)
                           aGCCTCACTTTCTGCCCAGTCATCCACCAaagaag         >  1:4030529/1‑35 (MQ=255)
                            gcatcacTTtgtactgaccccaaaaagttggacag        <  1:16596857‑M1/32‑27 (MQ=255)
                             gtacaCTTTCTGCCCAGTCATCCACCAaagaagaa       >  1:311293‑M2/6‑35 (MQ=255)
                             gtacaCTTTCTGCCCAGTCATCCACCAaagaagaa       <  1:6923071‑M2/30‑1 (MQ=255)
                              tacaCTTTCTGCCCAGTCATCCACCAaagaagaag      >  1:1562547‑M2/5‑35 (MQ=255)
                              tacaCTTTCTGCCCAGTCATCCACCAaagaagaag      >  1:6322809‑M2/5‑35 (MQ=255)
                              nacaCTTTCTGCCCAGTCATCCACCAaagaagaag      >  1:18451743‑M2/5‑35 (MQ=255)
                              cTCACTTtgtactgaccccaaaaagttggacagtt      <  1:17325662‑M1/35‑29 (MQ=255)
                              cTCACTTtgtactgaccccaaaaagttggacagtt      >  1:13160530‑M1/1‑7 (MQ=255)
                               acaCTTTCTGCCCAGTCATCCACCAAAGAAGAAGc     <  1:18973625‑M2/32‑1 (MQ=255)
                               acaCTTTCTGCCCAGTCATCCACCAAAGAAGAAGc     <  1:8743355‑M2/32‑1 (MQ=255)
                                caCTTTCTGCCCAGTCATCCACCAAAGAAGAAGcc    >  1:8836653‑M2/3‑35 (MQ=255)
                                caCTTTCTGCCCAGTCATCCACCAAAGAAGAAGcc    >  1:4622795‑M2/3‑35 (MQ=255)
                                cACTTtgtacttaccccaaaaagttggacagttaa    >  1:13848976‑M1/1‑5 (MQ=255)
                                cACTTtgtactgaccccaaaaagttggacagttaa    <  1:3720998‑M1/35‑31 (MQ=255)
                                cACTTtgtactgaccccaaaaagttggacagttaa    <  1:4202828‑M1/35‑31 (MQ=255)
                                 ccTTTCTGCCCAGTCATCCACAAAAGAAGAagcca   <  1:19849577/34‑1 (MQ=255)
                                 aCTTtgtactgaccccaaaaagttggacagttaaa   >  1:7884569‑M1/1‑4 (MQ=255)
                                 aCTTTCTGCCCAGTCATCCACCAAAGAAGAAGCCa   >  1:15919065‑M2/2‑35 (MQ=255)
                                  cTTtgtactgaccccaaaaagttggacagttaaac  >  1:18550631‑M1/1‑3 (MQ=255)
                                  cTTtgtactgaccccaaaaagttggacagttaaac  <  1:643702‑M1/35‑33 (MQ=255)
                                  cTTtgtactgaccccaaaaagttggacagttaaac  <  1:5562958‑M1/35‑33 (MQ=255)
                                  cTTtgtactgaccccaaaaagttggacagttaaac  >  1:9221751‑M1/1‑3 (MQ=255)
                                  cTTTCTGCCCAGTCATCCACCAAAGAAGAagccag  >  1:4864385/1‑35 (MQ=255)
                                  cTTTCTGCCCAGTCATCCACCAAAGAAGAagccag  <  1:10833767/35‑1 (MQ=255)
                                  |                                  
GAATCTCTCGCTGCAAAAGCTGGAGGAAGCCTCACTTTCTGCCCAGTCATCCACCAAAGAAGAAGCCAG  >  REL606/3015737‑3015805

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 29 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: 
Reads not counted as support for junction
read_name Not counted due to insufficient overlap past the breakpoint.
read_name Not counted due to not crossing MOB target site duplication.