Predicted mutation
evidence seq id position mutation annotation gene description
JC JC REL606 1,916,224 IS1 (–) +9 bp [pykA] [pykA]

New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? REL606 1916224 =1 (0.020)62 (0.960) 50/200 0.3 98.3% intergenic (+127/‑1) yebK/pykA predicted DNA‑binding transcriptional regulator/pyruvate kinase
?REL606 1973272 = NA (NA)noncoding (1/768 nt) IS1 repeat region
* ? REL606 = 19162321 (0.020)68 (1.060) 53/200 0.2 98.4% coding (8/1443 nt) pykA pyruvate kinase
?REL606 2034326 = NA (NA)noncoding (768/768 nt) IS1 repeat region

AATATAACCGGAATAAGAGGTGAGCTCATAAAAATCTCCTGATAACTTTTTTATATATATAGCCACGAAAACTAGTGGGAAGAAGTCTAACAAATATCCGGTAATGACTCCAACTTATTGATAGTGTTTTATGTTCAGATAATGCCCGATGACTTTGTCATGCAGCTCCACCGATTTTGAGAACGACAGCGACTTCCGT  >  REL606/2034227‑2034425
                                                                                                   |                                                                                                   
aaTATAACCGGAATAAGAGGTGAGCTCATAAAAATCTCCTGATAACTTTTTTATATATATAGCCACGAAAACTAGTGGGAAGAAGTCTAACAAATATCCgg                                                                                                    <  1:209637/101‑1 (MQ=255)
aaTATAACCGGAATAAGAGGTGAGCTCATAAAAATCTCCTGATAACTTTTTTATATATATAGCCACGAAAACTAGTGGGAAGAAGTCTAACAAATATCCgg                                                                                                    >  2:855565/1‑101 (MQ=255)
aaTATAACCGGAATAAGAGGTGAGCTCATAAAAATCTCCTGATAACTTTTTTATATATATAGCCACGAAAACTAGTGGGAAGAAGTCTAACAAATATCCgg                                                                                                    >  2:1445598/1‑101 (MQ=255)
 atatAACCGGAATAAGAGGTGAGCTCATAAAAATCTCCTGATAACTTTTTTATATATATAGCCACGAAAACTAGTGGGAAGAAGTCTAACAAATATCCGGt                                                                                                   >  2:19177/1‑101 (MQ=255)
 atatAACCGGAATAAGAGGTGAGATCATAAAAATCTCCTGATAACTTTTTTATATATATAGCCACGAAAACTAGTGGTAAGAAGTCTAACAAATATCCGGt                                                                                                   <  2:775882/101‑1 (MQ=255)
  tataACCGGAATAAGAGGTGAGCTCATAAAAATCTCCTGATAACTTTTTTATATATATAGCCACGAAAACTAGTGGGAAGAAGTCTAACAAATATCCGGTa                                                                                                  >  2:48632/1‑101 (MQ=255)
   ataACCGGAATAAGAGGTGAGCTCATAAAAATCTCCTGATAACTTTTTTATATATATAGCCACGAAAACTAGTGGGAAGAAGTCTAACAAATATCCGGTaa                                                                                                 >  1:871383/1‑101 (MQ=255)
     aaCCGGAATAAGAGGTGAGCTCATAAAAATCTCCTGATAACTTTTTTATATATATAGCCACGAAAACTAGTGGGAAGAAGTCTAACAAATATCCGGTAATg                                                                                               >  2:703788/1‑101 (MQ=255)
     aaCCGGAATAAGAGGTGAGCTCATAAAAATCTCCTGATAACTTTTTTATATATATAGCCACGAAAACTAGTGGGAAGAAGTCTAACAAATATCCGGTAATg                                                                                               >  2:1254335/1‑101 (MQ=255)
       ccGGAATAAGAGGTGAGCTCATAAAAATCTCCTGATAACTTTTTTATATATATAGCCACGAAAACTAGTGGGAAGAAGTCTAACAAATATCCGGTAATGAc                                                                                             <  1:1303462/101‑1 (MQ=255)
         ggAATAAGAGGTGAGCTCATAAAAATCTCCTGATAACTTTTTTATATATATAGCCACGAAAACTAGTGGGAAGAAGTCTAACAAATATCCGGTAATGACTc                                                                                           >  1:814971/1‑101 (MQ=255)
          gAATAAGAGGTGAGCTCATAAAAATCTCCTGATAACTTTTTTATATATATAGCCACGAAAACTAGTGGGAAGAAGTCTAACAAATATCCGGTAATGACTcc                                                                                          >  2:713341/1‑101 (MQ=255)
          gAATAAGAGGTGAGCTCATAAAAATCTCCTGATAACTTTTTTATATATATAGCCACGAAAACTAGTGGGAAGAAGTCTAACAAATATCCGGTAATGACTcc                                                                                          >  1:1477581/1‑101 (MQ=255)
            attagagGTGAGCTCATAAAAATCTCCTGATAACTTTTTTATATATATAGCCACGAAAACTAGTGGGAAGAAGTCTAACAAATATCCGGTAATGACTCCaa                                                                                        >  1:1426465/4‑101 (MQ=255)
            aTAAGAGGTGAGCTCATAAAAATCTCCTGATAACTTTTTTATATATATAGCCACGAAAACTAGTGGGAAGAAGTCTAACAAATATCCGGTAATGACTCCaa                                                                                        >  2:678038/1‑101 (MQ=255)
            aTAAGAGGTGAGCTCATAAAAATCTCCTGATAACTTTTTTATATATATAGCCACGAAAACTAGTGGGAAGAAGTCTAACAAATATCCGGTAATGACTCCaa                                                                                        <  2:66736/101‑1 (MQ=255)
            aTAAGAGGTGAGCTCATAAAAATCTCCTGATAACTTTTTTATATATATAGCCACGAAAACTAGTGGGAAGAAGTCTAACAAATATCCGGTAATGACTCCaa                                                                                        >  2:715212/1‑101 (MQ=255)
              aaGAGGTGAGCTCATAAAAATCTCCTGATAACTTTTTTATATATATAGCCACGAAAACTAGTGGGAAGAAGTCTAACAAATATCCGGTAATGACTCCAACt                                                                                      >  1:1503833/1‑101 (MQ=255)
               agagGTGAGCTCATAAAAATCTCCTGATAACTTTTTTATATATATAGCCACGAAAACTAGTGGGAAGAAGTCTAACAAATATCCGGTAATGACTCCAACtt                                                                                     <  2:948579/101‑1 (MQ=255)
                  ggcgAGCTCATAAAAATCTCCTGATAACTTTTTTATATATATAGCCACGAAAACTAGTGGGAAGAAGTCTAACAAATATCCGGTAATGACTCCAACTTAtt                                                                                  >  2:852560/4‑101 (MQ=255)
                   gTGAGCTCATAAAAATCTCCTGATAACTTTTTTATATATATAGCCACGAAAACTAGTGGGAAGAAGTCTAACAAATATCCGGTAATGACTCCAACTTATTg                                                                                 <  1:495173/101‑1 (MQ=255)
                     gAGCTCATAAAAATCTCCTGATAACTTTTTTATATATATAGCCACGAAAACTAGTGGGAAGAAGTCTAACAAATATCCGGTAATGACTCCAACTTATTGAt                                                                               <  1:763464/101‑1 (MQ=255)
                     gAGCTCATAAAAATCTCCTGATAACTTTTTTATATATATAGCCACGAAAACTAGTGGGAAGAAGTCTAACAAATATCCGGTAATGACTCCAACTTATTGAt                                                                               <  2:1156807/101‑1 (MQ=255)
                          cATAAAAATCTCCTGATAACTTTTTTATATATATAGCCACGAAAACTAGTGGGAAGAAGTCTAACAAATATCCGGTAATGACTCCAACTTATTGATAgtgt                                                                          >  2:586402/1‑101 (MQ=255)
                           aTAAAAATCTCCTGATAACTTTTTTATATATATAGCCACGAAAACTAGTGGGAAGAAGTCTGACAAATATCCGGTAATGACTCCAACTTATTGATAGTGtt                                                                         >  2:422972/1‑101 (MQ=255)
                           aTAAAAATCTCCTGATAACTTTTTTATATATATAGCCACGAAAACTAGTGGGAAGAAGTCTAACAAATATCCGGTAATGACTCCAACTTATTGATAGTGtt                                                                         >  1:425442/1‑101 (MQ=255)
                           aTAAAAATCTCCTGATAACTTTTTTATATATATAGCCACGAAAACTAGTGGGAAGAAGTCTAACAAATATCCGGTAATGACTCCAACTTATTGATAGTGtt                                                                         >  1:1180477/1‑101 (MQ=255)
                           aTAAAAATCTCCTGATAACTTTTTTATATATATAGCCACGAAAACTAGTGGGAAGAAGTCTAACAAATATCCGGTAATGACTCCAACTTATTGATAGTGtt                                                                         >  1:1487438/1‑101 (MQ=255)
                           aTAAAAATCTCCTGATAACTTTTTTATATATATAGCCACGAAAACTAGTGGGAAGAAGTCTAACAAATATCCGGTAATGACTCCAACTTATTGATAGTGtt                                                                         >  2:102792/1‑101 (MQ=255)
                           aTAAAAATCTCCTGATAACTTTTTTATATATATAGCCACGAAAACTAGTGGGAAGAAGTCTAACAAATATCCGGTAATGACTCCAACTTATTGATAGTGtt                                                                         >  1:1134571/1‑101 (MQ=255)
                           aTAAAAATCTCCTGATAACTTTTTTATATATATAGCCACGAAAACTAGTGGGAAGAAGTCTAACAAATATCCGGTAATGACTCCAACTTATTGATAGTGtt                                                                         >  1:161146/1‑101 (MQ=255)
                             aaaaaTCTCCTGATAACTTTTTTATATATATAGCCACGAAAACTAGTGGGAAGAAGTCTAACAAATATCCGGTAATGACTCCAACTTATTGATAGTGtttt                                                                       >  1:1706023/1‑101 (MQ=255)
                               aaaTCTCCTGATAACTTTTTTATATATATAGCCACGAAAAGTAGTGGGAAGAAGTCTAACAAATATCCGGTAATGACTCCGACTTATTGATAGTGTTTTAg                                                                     <  2:409819‑M1/101‑2 (MQ=255)
                                aaTCTCCTGATAACTTTTTTATATATATAGCCACGAAAACTAGTGGGAAGAAGTCTAACAAATATCCGGTAATGACTCCAACTTATTGATAGTGTTTTAgt                                                                    >  1:953309‑M1/1‑99 (MQ=255)
                                aaTCTCCTGATAACTTTTTTATATATATAGCCACGAAAACTAGTGGGAAGAAGTCTAACAAATATCCGGTAATGACTCCAACTTATTGATAGTGTTTTAgt                                                                    <  1:914870‑M1/101‑3 (MQ=255)
                                 aTCTCCTGATAACTTTTTTATATATATAGCCACGAAAACTAGTGGGAAGAAGTCTAACAAATATCCGGTAATGACTCCAACTTATTGATAGTGTTTTAgtt                                                                   >  1:1029915‑M1/1‑98 (MQ=255)
                                     ccTGATAACTTTTTTATATATATAGCCACGAAAACTAGTGGGAAGAAGTCTAACAAATATCCGGTAATGACTCCAACTTATTGATAGTGTTTTAgtttgag                                                               <  1:224285‑M1/101‑8 (MQ=255)
                                      cTGATAACTTTTTTATATATATAGCCACGAAAACTAGTGGGAAGAAGTCTAACAAATATCCGGTAATGACTCCAACTTATTGATAGTGTTTTAgtttgagc                                                              <  1:715217‑M1/101‑9 (MQ=255)
                                        gATAACTTTTTTATATATATAGCCACGAAAACTAGTGGGAAGAAGTCTAACAAATATCCGGTAATGACTCCAACTTATTGATAGTGTTTTAgtttgagcag                                                            >  2:799059‑M1/1‑91 (MQ=255)
                                        gATAACTTTTTTATATATATAGCCACGAAAACTAGTGGGAAGAAGTCTAACAAATATCCGGTAATGACTCCAACTTATTGATAGTGTTTTAgtttgagcag                                                            >  2:105695‑M1/1‑91 (MQ=255)
                                         aTAACTTTTTTATATATATAGCCACGAAAACTAGTGGGAAGAAGTCTAACAAATATCCGGTAATGACTCCAACTTATTGATAGTGTTTTAgtttgagcagg                                                           <  1:1510856‑M1/101‑12 (MQ=255)
                                         aTAACTTTTTTATATATATAGCCACGAAAACTAGTGGGAAGAAGTCTAACAAATATCCGGTAATGACTCCAACTTATTGATAGTGTTTTAgtttgagcagg                                                           <  2:1173601‑M1/101‑12 (MQ=255)
                                         aTAACTTTTTTATATATATAGCCACGAAAACTAGTGGGAAGAAGTCTAACAAATATCCGGTAATGACTCCAACTTATTGATAGTGTTTTAgtttgagcagg                                                           <  2:793598‑M1/101‑12 (MQ=255)
                                         aTAACTTTTTTATATATATAGCCACGAAAACTAGTGGGAAGAAGTCTAACAAATATCCGGTAATGACTCCAACTTATTGATAGTGTTTTAgtttgagcagg                                                           <  1:1445604‑M1/101‑12 (MQ=255)
                                           aaCTTTTTTATATATATAGCCACGAAAACTAGTGGGAAGAAGTCTAACAAATATCCGGTAATGACTCCAACTTATTGATAGTGTTTTAgtttgagcaggcg                                                         >  2:570464‑M1/1‑88 (MQ=255)
                                             cTTTTTTATATATATAGCCACGAAAACTAGTGGGAAGAAGTCTAACAAATATCCGGTAATGACTCCAACTTATTGATAGTGTTTTAgtttgagcaggcggt                                                       >  1:910410‑M1/1‑86 (MQ=255)
                                             cTTTTTTATATATATAGCCACGAAAACTAGTGGGAAGAAGTCTAACAAATATCCGGTAATGACTCCAACTTATTGATAGTGCTTTAgtttgagcaggcggt                                                       >  1:147679‑M1/1‑86 (MQ=255)
                                                 tttATATATATAGCCACGAAAACTAGTGGGAAGAAGTCTAACAAATATCCGGTAATGACTCCAACTTATTGATAGTGTTTTAgtttgagcaggcggttaac                                                   <  1:325727‑M1/101‑20 (MQ=255)
                                                  ttATATATATAGCCACGAAAACTAGTGGGAAGAAGTCTAACAAATATCCGGTAATGACTCCAACTTATTGATAGTGTTTTAgtttgagcaggcggttaacc                                                  >  1:959199‑M1/1‑81 (MQ=255)
                                                   tatatatataGCCACGAAAACTAGTGGGAAGAAGTCTAACAAATATCCGGTAATGACTCCAACTTATTGATAGTGTTTTAgtttgagcaggcggttaaccc                                                 >  1:852611‑M1/1‑80 (MQ=255)
                                                        atataGCCACGAAAACTAGTGGGAAGAAGTCTAACAAATATCCGGTAATGACTCCAACTTATTGATAGTGTTTTAgtttgagcaggcggttaaccccggaa                                            >  2:991595‑M1/1‑75 (MQ=255)
                                                        atataGCCACGAAAACTAGTGGGAAGAAGTCTAACAAATATCCGGTAATGACTCCAACTTATTGATAGTGTTTTAgtttgagcaggcaggtaactccgcaa                                            >  2:1334585‑M1/1‑75 (MQ=255)
                                                          ataGCCACGAAAACTAGTGGGAAGAAGTCTAACAAATATCCGGTAATGACTCCAACTTATTGATAGTGTTTTAgtttgagcaggcggttaaccccggaaag                                          <  2:1182486‑M1/101‑29 (MQ=255)
                                                           taGCCACGAAAACTAGTGGGAAGAAGTCTAACAAATATCCGGTAATGACTCCAACTTATTGATAGTGTTTTAgtttgagcaggcggttaaccccggaaagc                                         <  1:1326990‑M1/101‑30 (MQ=255)
                                                             gCCACGAAAACTAGTGGGAAGAAGTCTAACAAATATCCGGTAATGACTCCCGCTTATTGATAGTGTTTTAgtttgagcaggcgatgaaccccggaaagcgg                                       <  2:1078383‑M1/101‑32 (MQ=255)
                                                             gCCACGAAAACTAGTGGGAAGAAGTCTAACAAATATCCGGTAATGACTCCAACTTATTGATAGTGTTTTAgtttgagcaggcggttaaccccggaaagcgg                                       <  2:1180472‑M1/101‑32 (MQ=255)
                                                              ccACGAAAACTAGTGGGAAGAAGTCTAACAAATATCCGGTAATGACTCCAACTTATTGATAGTGTTTTAgtttgagcaggcggttaaccccggaaagcggc                                      >  1:902094‑M1/1‑69 (MQ=255)
                                                              ccACGAAAACTAGTGGGAAGAAGTCTAACAAATATCCGGTAATGACTCCAACTTATTGATAGTGTTTTAgtttgagcaggcggttaaccccggaaagcggc                                      <  2:814966‑M1/101‑33 (MQ=255)
                                                              ccACGAAAACTAGTGGGAAGAAGTCTAACAAATATCCGGTAATGACTCCAACTTATTGATAGTGTTTTAgtttgagcaggcggttaaccccggaaagcggc                                      >  2:285918‑M1/1‑69 (MQ=255)
                                                                 cGAAAACTAGTGGGAAGAAGTCTAACAAATATCCGGTAATGACTCCAACTTATTGATAGTGTTTTAgtttgagcaggcggttaaccccggaaagcggcgtg                                   <  1:997890‑M1/101‑36 (MQ=255)
                                                                  gAAAACTAGTGGGAAGAAGTCTAACAAATATCCGGTAATGACTCCAACTTATTGATAGTGTTTTAgtttgagcaggcggttaaccccggaaagcggcgtgt                                  >  1:1189818‑M1/1‑65 (MQ=255)
                                                                   aaaaCTAGTGGGAAGAAGTCTAACAAATATCCGGTAATGACTCCAACTTATTGATAGTGTTTTAgtttgagcaggcggttaaccccggaaagcggcgtgtc                                 >  2:871948‑M1/1‑64 (MQ=255)
                                                                     aaCTAGTGGGAAGAAGTCTAACAAATATCCGGTAATGACTCCAACTTATTGATAGTGTTTTAgtttgagcaggcggttaaccccggaaagcggcgtgtcat                               >  2:1245982‑M1/1‑62 (MQ=255)
                                                                     aaCTAGTGGGAAGAAGTCTAACAAATATCCGGTAATGACTCCAACTTATTGATAGTGTTTTAgtttgagcaggcggttaaccccggaaagcggcgtgtcat                               >  1:379150‑M1/1‑62 (MQ=255)
                                                                       cTAGTGGGAAGAAGTCTAACAAATATCCGGTAATGACTCCAACTTATTGATAGTGTTTTAgtttgagcaggcggttaaccccggaaagcggcgtgtcatac                             <  2:1000293‑M1/101‑42 (MQ=255)
                                                                        tAGTGGGAAGAAGTCTAACAAATATCCGGTAATGACTCCAACTTATTGATAGTGTTTTAgtttgagcaggcggttaaccccggaaagcggcgtgtcataca                            <  1:19177‑M1/101‑43 (MQ=255)
                                                                         aGTGGGAAGAAGTCTAACAAATATCCGGTAATGACTCCAACTTATTGATAGTGTTTTAgtttgagcaggcggttaaccccggaaagcggcgtgtcatacag                           <  2:425437‑M1/101‑44 (MQ=255)
                                                                         aGTGGGAAGAAGTCTAACAAATATCCGGTAATGACTCCAACTTATTGATAGTGTTTTAgtttgagcaggcggttaaccccggaaagcggcgtgtcatacag                           >  2:1705256‑M1/1‑58 (MQ=255)
                                                                           tGGGAAGAAGTCTAACAAATATCCGGTAATGACTCCAACTTATTGATAGTGTTTTAgtttgagcaggcggttaaccccggaaagcggcgtgtcatacagcg                         >  1:327500‑M1/1‑56 (MQ=255)
                                                                             ggAAGAAGTCTAACAAATATCCGGTAATGACTCCAACTTATTGATAGTGTTTTAgtttgagcaggcggttaaccccggaaagcggcgtgtcatacagcggc                       >  2:331996‑M1/1‑54 (MQ=255)
                                                                             ggAAGAAGTCTAACAAATATCCGGTAATGACTCCAACTTATTGATAGTGTTTTAgtttgagcaggcggttaaccccggaaagcggcgtgtcatacagcggc                       >  2:1573444‑M1/1‑54 (MQ=255)
                                                                               aagaagTCTAACAAATATCCGGTAATGACTCCAACTTATTGATAGTGTTTTAgtttgagcaggcggttaaccccggaaagcggcgtgtcatacagcggcac                     >  1:756440‑M1/1‑52 (MQ=255)
                                                                               aagaagTCTAACAAATATCCGGTAATGACTCCAACTTATTGATAGTGTTTTAgtttgagcaggcggttaaccccggaaagcggcgtgtcatacagcggcac                     >  2:1111518‑M1/1‑52 (MQ=255)
                                                                                 gaagTCTAACAAATATCCGGTAATGACTCCAACTTATTGATAGTGTTTTAgtttgagcaggcggttaaccccggaaagcggcgtgtcatacagcggcacca                   <  1:1191289‑M1/101‑52 (MQ=255)
                                                                                   agTCTAACAAATATCCGGTAATGACTCCAACTTATTGATAGTGTTTTAgtttgagcaggcggttaaccccggaaagcggcgtgtcatacagcggcaccacc                 >  2:784032‑M1/1‑48 (MQ=255)
                                                                                      cTAACAAATATCCGGTAATGACTCCAACTTATTGATAGTGTTTTAgtttgagcaggcggttaaccccggaaagcggcgtgtcatacagcggcaccaccggt              >  2:1388440‑M1/1‑45 (MQ=255)
                                                                                      cTAACAAATATCCGGTAATGACTCCAACTTATTGATAGTGTTTTAgtttgagcaggcggttaaccccggaaagcggcgtgccatacagcggcaccaccggt              >  1:943218‑M1/1‑45 (MQ=255)
                                                                                        aaCAAATATCCGGTAATGACTCCAACTTATTGATAGTGTTTTAgtttgagcaggcggttaaccccggaaagcggcgtgtcatacagcggcaccaccggtac            <  2:953304‑M1/101‑59 (MQ=255)
                                                                                        aaCAAATATCCGGTAATGACTCCAACTTATTGATAGTGTTTTAgtttgagcaggcggttaaccccggaaagcggcgtgtcatacagcggcaccaccggtac            <  1:1362025‑M1/101‑59 (MQ=255)
                                                                                        aaCAAATATCCGGTAATGACTCCAACTTATTGATAGTGTTTTAgtttgagcaggcggttaaccccggaaagcggcgtgtcatacagcggcaccaccggtac            <  2:1134566‑M1/101‑59 (MQ=255)
                                                                                        aaCAAATATCCGGTAATGACTCCAACTTATTGATAGTGTTTTAgtttgagcaggcggttaaccccggaaagcggcgtgtcatacagcggcaccaccggtac            >  2:87647‑M1/1‑43 (MQ=255)
                                                                                        aaCAAATATCCGGTAATGACTCCAACTTATTGATAGTGTTTTAgtttgagcaggcggttaaccccggaaagcggcgtgtcatacagcggcaccaccggtac            <  2:902089‑M1/101‑59 (MQ=255)
                                                                                         aCAAATATCCGGTAATGACTCCAACTTATTGATAGTGTTTTAgtttgagcaggcggttaaccccggaaagcggcgtgtcatacagcggcaccaccggtacg           <  1:1549899‑M1/101‑60 (MQ=255)
                                                                                               atCCGGTAATGACTCCAACTTATTGATAGTGTTTTAgtttgagcaggcggttaaccccggaaagcggcgtgtcatacagcggcaccaccggtacgccattc     >  1:1350848‑M1/1‑36 (MQ=255)
                                                                                                  cGGTAATGACTCCAACTTATTGATAGTGTTTTAgtttgagcaggcggttaaccccggaaagcggcgtgtcatacagcggcaccaccggtacgccattcatc  >  2:1049011‑M1/1‑33 (MQ=255)
                                                                                                   |                                                                                                   
AATATAACCGGAATAAGAGGTGAGCTCATAAAAATCTCCTGATAACTTTTTTATATATATAGCCACGAAAACTAGTGGGAAGAAGTCTAACAAATATCCGGTAATGACTCCAACTTATTGATAGTGTTTTATGTTCAGATAATGCCCGATGACTTTGTCATGCAGCTCCACCGATTTTGAGAACGACAGCGACTTCCGT  >  REL606/2034227‑2034425

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