Predicted mutation
evidence seq id position mutation freq annotation gene description
JC JC REL606 2,448,493 IS186 (+) +6 bp 66.5% intergenic (+24/+21) nupC → / ← yfeA nucleoside (except guanosine) transporter/predicted diguanylate cyclase

New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? REL606 2448493 =30 (0.470)59 (0.660) 32/42 NT 58.4% intergenic (+24/+26) nupC/yfeA nucleoside (except guanosine) transporter/predicted diguanylate cyclase
?REL606 = 2774196 NA (NA)noncoding (1343/1343 nt) IS186 repeat region
* ? REL606 = 244849832 (0.500)64 (0.710) 34/42 NT 58.8% intergenic (+29/+21) nupC/yfeA nucleoside (except guanosine) transporter/predicted diguanylate cyclase
?REL606 2772854 = NA (NA)noncoding (1/1343 nt) IS186 repeat region

Only 100 of 300 total aligned reads displayed.

CTGTAAGACCATACATAAAAAAGCCGGGGATAATTCCCCGGTTTTTTTACGCCTGTTAATCAGCTAATG  >  REL606/2448464‑2448532
                                  |                                  
ggttcaacccttaagttagctcttatgggATAAtt                                    <  1:1328599‑M1/6‑1 (MQ=255)
ggttcaacccttaagttagcgcttatgggATAAtt                                    <  1:1350345‑M1/6‑1 (MQ=255)
cTGTAAGACCATACATAAAAAAGCCGGGGATAAtt                                    <  1:3938081/35‑1 (MQ=255)
cTGTAAGACCATACATAAAAAAGCCGGGGATAAtt                                    >  1:7714225/1‑35 (MQ=255)
cTGTAAGACCATACATAAAAAAGCCGGGGATAAtt                                    >  1:8948967/1‑35 (MQ=255)
cTGTAAGACCATACATAAAAAAGCCGGGGATAAtt                                    <  1:19032853/35‑1 (MQ=255)
  ttcaacccttaagttagcgcttatgggATAATTcc                                  >  1:15301037‑M1/28‑35 (MQ=255)
    caacccttaagttagcgcttatgggATAATTcccc                                >  1:19825005‑M1/26‑35 (MQ=255)
    aaGCCCATACATAAAAAAGCCGGGGATAATTcccc                                <  1:19702743/35‑1 (MQ=255)
    aaGACCATACATAAAAAAGCCGGGGATAATTcccc                                <  1:2828101/35‑1 (MQ=255)
    aaGACCATACATAAAAAAGCCGGGGATAATTcccc                                <  1:17004888/35‑1 (MQ=255)
     aacccttaagttagcgcttatgggATAATTCCCCg                               >  1:18362053‑M1/25‑35 (MQ=255)
     aacccttaagttagcgcttatgggATAATTCCCCg                               >  1:17130580‑M1/25‑35 (MQ=255)
     aGACCATACATAAAAAAGCCGGGGATAATTCCCCg                               >  1:19880050/1‑35 (MQ=255)
     aGACCATACATAAAAAAGCCGGGGATAATTCCCCg                               <  1:3727615/35‑1 (MQ=255)
     aGACCATACATAAAAAAGCCGGGGATAATTCCCCg                               <  1:14381800/35‑1 (MQ=255)
      acccttaagttagcgcttatgggATAATTCCCCgg                              >  1:18139901‑M1/24‑35 (MQ=255)
       cccttaagttagcgcttatgggATAATTCCCCGGt                             <  1:11425113‑M1/13‑1 (MQ=255)
        ccttaagttagcgcttatgggATAATTCCCCGGtt                            <  1:12909675‑M1/14‑1 (MQ=255)
         cttaagttagcgcttatgggATAATTCCCCGGttt                           >  1:1623639‑M1/21‑35 (MQ=255)
         cttaagttagcgcttatgggATAATTCCCCGGttt                           <  1:2267930‑M1/15‑1 (MQ=255)
         catacataAAAAAGCCGGGGATAATTCCCCGGttt                           <  1:5135638/35‑1 (MQ=255)
         catacataAAAAAGCCGGGGATAATTCCCCGGttt                           >  1:12429540/1‑35 (MQ=255)
          ttaagttagcgcttatgggATAATTCCCCGGtttt                          >  1:9217002‑M1/20‑35 (MQ=255)
          atacataAAAAAGCCGGGGATAATTcccataagcg                          <  1:9595250‑M1/35‑11 (MQ=255)
          atacataAAAAAGCCGGGGATAATTCCCCGGtttt                          <  1:7419672/35‑1 (MQ=255)
          atacataAAAAAGCCGGGGATAATTCCCCGGtttt                          >  1:19431837/1‑35 (MQ=255)
           tacataAAAAAGCCGGGGATAATTcccataagcgc                         <  1:17704988‑M1/35‑12 (MQ=255)
           tacataAAAAAGCCGGGGATAATTCCCCGGttttt                         >  1:13410668/1‑35 (MQ=255)
           taagttagcgcttatgggATAATTCCCCGGttttt                         <  1:15814655‑M1/17‑1 (MQ=255)
           taagttagcgcttatgggATAATTCCCCGGttttt                         >  1:17223589‑M1/19‑35 (MQ=255)
            acataAAAAAGCCGGGGATAATTcccataagcgct                        <  1:14786978‑M1/35‑13 (MQ=255)
            acataAAAAAGCCGGGGATAATTCCCCGGtttttt                        <  1:4535044/35‑1 (MQ=255)
            aagttagcgcttatgggATAATTCCCCGGtttttt                        >  1:20345336‑M1/18‑35 (MQ=255)
             cataAAAAAGCCGGGGATGATTcccataagcgcta                       >  1:139947‑M1/1‑22 (MQ=255)
             cataAAAAAGCCGGGGATGATTcccataagcgcta                       <  1:19409856‑M1/35‑14 (MQ=255)
             cataAAAAAGCCGGGGATAATTcccataagcgcta                       >  1:18646502‑M1/1‑22 (MQ=255)
             agttagcgcttatgggATAATTCCCCGGttttttt                       <  1:234345‑M1/19‑1 (MQ=255)
             agttagcgcttatgggATAATTCCCCGGttttttt                       <  1:204309‑M1/19‑1 (MQ=255)
              gttagcgcttatgggATAATTCCCCGGTTTTTTTa                      <  1:8657160‑M1/20‑1 (MQ=255)
               ttagcgcttatgggATAATTCCCCGGTTTTTTTAc                     <  1:7855353‑M1/21‑1 (MQ=255)
               ttagcgcttatgggATAATTCCCCGGTTTTTTTAc                     <  1:2587811‑M1/21‑1 (MQ=255)
               ttagcgcttatgggATAATTCCCCGGTTTTTTTAc                     >  1:6366113‑M1/15‑35 (MQ=255)
               ttagcgcttatgggATAATTCCCCGGTTTTTTTAc                     <  1:17432747‑M1/21‑1 (MQ=255)
                tagcgcttatgggATAATTCCCCGGTTTTTTTACg                    <  1:3648455‑M1/22‑1 (MQ=255)
                aaaaaaGCCGGGGATAATTcccataagcgctaact                    >  1:3109192‑M1/1‑19 (MQ=255)
                aaaaaaGCCGGGGATAATTCCCCGGTTTTTTTACg                    >  1:6068266/1‑35 (MQ=255)
                aaaaaaGCCGGGGATAATTCCCCGGTTTTTTTACg                    <  1:2743146/35‑1 (MQ=255)
                 agcgcttatgggATAATTCCCCGGTTTTTTTACGc                   >  1:13550861‑M1/13‑35 (MQ=255)
                 agcgcttatgggATAATTCCCCGGTTTTTTTACGc                   >  1:11810127‑M1/13‑35 (MQ=255)
                 agcgcttatgggATAATTCCCCGGTTTTTTTACGc                   <  1:15743302‑M1/23‑1 (MQ=255)
                 aaaaaGCCGGGGATAATTcccataagcgctaactt                   <  1:6507339‑M1/35‑18 (MQ=255)
                  gcgcttatgggATAATTCCCCGGTTTTTTTACGcc                  >  1:18932778‑M1/12‑35 (MQ=255)
                   aaaGCCGGGGATAATTCCCCGGTTTTTTTACGCCt                 <  1:15137944/35‑1 (MQ=255)
                    gcttatgggATAATTCCCCGGTTTTTTTACGCCTg                >  1:3428254‑M1/10‑35 (MQ=255)
                    aaGCCGGGGATAATTcccataagcgctaacttaag                >  1:9139407‑M1/1‑15 (MQ=255)
                    aaGCCGGGGATAATTCCCCGGTTTTTTTACGCCTg                <  1:11594496/35‑1 (MQ=255)
                     cttatgggATAATTCCCCGGTTTTTTTACGCCTGt               >  1:17902106‑M1/9‑35 (MQ=255)
                     cttatgggATAATTCCCCGGTTTTTTTACGCCTGt               >  1:9395986‑M1/9‑35 (MQ=255)
                      gCCGGGGATAATTCCCCGGTTTTTTTACGCCTGtt              >  1:15402996/1‑35 (MQ=255)
                      gCCGGGGATAATTCCCCGGTTTTTTTACGCCTGtt              >  1:16367244/1‑35 (MQ=255)
                       tatgggATAATTCCCCGGTTTTTTTACGCCTGTTa             <  1:1101830‑M1/29‑1 (MQ=255)
                        cGGTGATAATTcccataagcgctaacttaagggtt            >  1:8531378‑M1/1‑11 (MQ=255)
                        cGGGGATAATTcccataagcgctaacttaagggtt            >  1:2061989‑M1/1‑11 (MQ=255)
                        cGGGGATAATTcccataagcgctaacttaagggtt            <  1:20117976‑M1/35‑25 (MQ=255)
                        cGGGGATAATTCCCCGGTTTTTTTACGCCTGTTaa            >  1:16906840/1‑35 (MQ=255)
                        atgggATAATTCCCCTGTTTTTTTACGCCTGTTaa            >  1:6222865‑M1/6‑35 (MQ=255)
                         tgggATAATTCCCCGGTTTTTTTACGCCTGTTAAt           <  1:18444975‑M1/31‑1 (MQ=255)
                         tgggATAATTCCCCGGTTTTTTTACGCCTGTTAAt           <  1:14295695‑M1/31‑1 (MQ=255)
                         tgggATAATTCCCCGGTTTTTTTACGCCTGTTAAt           >  1:19316554‑M1/5‑35 (MQ=255)
                         tgggATAATTCCCCGGTTTTTTTACGCCTGTTAAt           >  1:6573409‑M1/5‑35 (MQ=255)
                         ggggATAATTcccataagcgctaacttaagggttg           <  1:7014372‑M1/35‑26 (MQ=255)
                          gggATAATTcccataagcgctaacttaagggttgt          >  1:11732822‑M1/1‑9 (MQ=255)
                           ggATAATTcccataagcgctaacttaagggttgtg         >  1:7481101‑M1/1‑8 (MQ=255)
                           ggATAATTcccataagcgctaacttaagggttgtg         <  1:12675952‑M1/35‑28 (MQ=255)
                           ggATAATTCCCCGGTTTTTTTACGCCTGTTAATCa         <  1:4144605‑M1/33‑1 (MQ=255)
                           ggATAATTCCCCGGTTTTTTTACGCCTGTTAATCa         >  1:2869436‑M1/3‑35 (MQ=255)
                           ggATAATTCCCCGGTTTTTTTACGCCTGTTAATCa         <  1:15685924‑M1/33‑1 (MQ=255)
                            gATAATTcccataagcgctaacttaagggttgtgg        <  1:7316999‑M1/35‑29 (MQ=255)
                            gATAATTcccataagcgctaacttaagggttgtgg        <  1:16547122‑M1/35‑29 (MQ=255)
                            gATAATTCCCCGGTTTTTTTACGCCTGTTAATCAg        <  1:8166078‑M1/34‑1 (MQ=255)
                            gATAATTCCCCGGTTTTTTTACGCCTGTTAATCAg        <  1:10145800‑M1/34‑1 (MQ=255)
                             aTAATTCCCCGGTTTTTTTACGCCTGTTAATCAGc       >  1:13777010/1‑35 (MQ=255)
                             aTAATTCCCCGGTTTTTTTACGCCTGTTAATCAGc       >  1:4695640/1‑35 (MQ=255)
                             aTAATTCCCCGGTTTTTTTACGCCTGTTAATCAGc       >  1:19175267/1‑35 (MQ=255)
                              tAATTcccataagcgctaacttaagggttgtggta      <  1:9365536‑M1/35‑31 (MQ=255)
                              tAATTcccataagcgctaacttaagggttgtggta      <  1:379599‑M1/35‑31 (MQ=255)
                              tAATTCCCCGGTTTTTTTACGCCTGTTAATCAGCt      <  1:12229113/35‑1 (MQ=255)
                               aaTTcccataagcgctaacttaagggttgtggtat     <  1:9811270‑M1/35‑32 (MQ=255)
                               aaTTCCCCGGTTTTTTTACGCCTGTTAATCAGCTa     <  1:14549923/35‑1 (MQ=255)
                               aaTTCCCCGGTTTTTTTACGCCTGTTAATCAGCTa     >  1:7067413/1‑35 (MQ=255)
                               aaTTCCCCGGTTTTTTTACGCCTGTTAATCAGCTa     <  1:5575053/35‑1 (MQ=255)
                                  tcccataagcgctaacttaagggttgtggtattac  <  1:17343877‑M1/35‑35 (MQ=255)
                                  tcccataagcgctaacttaagggttgtggtattac  >  1:13196095‑M1/1‑1 (MQ=255)
                                  tcccataagcgctaacttaagggttgtggtattac  <  1:5235844‑M1/35‑35 (MQ=255)
                                  tcccataagcgctaacttaagggttgtggtattac  <  1:10619169‑M1/35‑35 (MQ=255)
                                  tcccataagcgctaacttaagggttgtggtattac  <  1:5414698‑M1/35‑35 (MQ=255)
                                  tccaataagcgctaacttaagggttgtggtattac  <  1:2935583‑M1/35‑35 (MQ=255)
                                  tccaataagcgctaacttaagggttgtggtattac  <  1:3872060‑M1/35‑35 (MQ=255)
                                  tCCCCGGTTTTTTTACGCCTGTTAATCAGCTAATg  <  1:4942343/35‑1 (MQ=255)
                                  |                                  
CTGTAAGACCATACATAAAAAAGCCGGGGATAATTCCCCGGTTTTTTTACGCCTGTTAATCAGCTAATG  >  REL606/2448464‑2448532

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 32 ≤ 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.