Predicted mutation
evidence seq id position mutation annotation gene description
JC REL606 3,571,131 +GTACA coding (332/435 nt) uspA → universal stress global response regulator

New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? REL606 = 35711292 (0.040)48 (0.880)
+CAGTA
39/136 0.3 96.2% coding (330/435 nt) uspA universal stress global response regulator
?REL606 3571130 = 2 (0.030)coding (331/435 nt) uspA universal stress global response regulator

GAGCGGCAGCGGCGACCTGGGCCAAGTTCTGGTCGATGCAATCAAGAAATACGATATGGATCTGGTGGTTTGTGGTCACCACCAGGACTTCTGGAGCAAACTGATGTCTTCCGCACGTCAGCTGATCAACACCGTTC  >  REL606/3571054‑3571190
                                                                           |                                                             
gAGCGGCAGCGGCGACCTGGGCCAAGTTCTGGTCGATGCAATCAAGAAATACGATATGGATCTGGTGGTTTGTGGt                                                               <  2:1710191/76‑1 (MQ=255)
  gcggcagcggcGACCTGGGCCAAGTTCTGGTCGATGCAATCAAGAAATACGATATGGATCTGGTGGTTTGTGGTca                                                             <  1:806606‑M1/76‑3 (MQ=255)
  gcggcagcggcGACCTGGGCCAAGTTCTGGTCGATGCAATCAAGAAATACGATATGGATCTGGTGGTTTGTGGTca                                                             <  1:806584‑M1/76‑3 (MQ=255)
   cggcagcggcGACCTGGGCCAAGTTCTGGTCGATGCAATCAAGAAATACGATATGGATCTGGTGGTTTGTGGTcag                                                            <  1:982285‑M1/76‑4 (MQ=255)
       agcggcGACCTGGGCCAAGTTCTGGTCGATGCAATCAAGAAATACGATATGGATCTGGTGGTTTGTGGTcagtaca                                                        >  1:854604‑M1/1‑69 (MQ=255)
          ggcgACCTGGGCCAAGTTCTGGTCGATGCAATCAAGAAATACGATATGGATCTGGTGGTTTGTGGTcagtacacca                                                     <  1:1806764‑M1/76‑11 (MQ=255)
           gcgACCTGGGCCAAGTTCTGGTCGATGCAATCAAGAAATACGATATGGATCTGGTGGTTTGTGGTcagtacaccac                                                    <  1:1083323‑M1/76‑12 (MQ=255)
            cgACCTGGGCCAAGTTCTGGTCGATGCAATCAAGAAATACGATATGGATCTGGTGGTTTGTGGTcagtacaccacc                                                   >  1:1142578‑M1/1‑64 (MQ=255)
             gACCTGGGCCAAGTTCTGGTCGATGCAATCAAGAAATACGATATGGATCTGGTGGTTTGTGGTcagtacaccacca                                                  <  2:301899‑M1/76‑14 (MQ=255)
              aCCTGGGCCAAGTTCTGGTCGATGCAATCAAGAAATACGATATGGATCTGGTGGTTTGTGGTcagtacaccaccag                                                 <  1:43785‑M1/76‑15 (MQ=255)
               ccTGGGCCAAGTTCTGGTCGATGCAATCAAGAAATACGATATGGATTTGGTGGTTTGTGGTcagtacaccaccagg                                                <  1:1438507‑M1/76‑16 (MQ=255)
                cTGGGCCAAGTTCTGGTCGATGCAATCAAGAAATACGATATGGATCTGGTGGTTTGTGGTcagtacaccaccagga                                               >  1:1651649‑M1/1‑60 (MQ=255)
                 tGGGCCAAGTTCTGGTCGATGCAATCAAGAAATACGATATGGATCTGGTGGTTTGTGGTcagtacaccaccaggac                                              <  2:1346584‑M1/76‑18 (MQ=255)
                  gggCCAAGTTCTGGTCGATGCAATCAAGAAATACGATATGGATCTGGTGGTTTGTGGTcagtacaccaccaggact                                             <  2:1036295‑M1/76‑19 (MQ=255)
                   ggCCAAGTTCTGGTCGATGCAATCAAGAAATACGATATGGATCTGGTGGTTTGTGGTcagtacaccaccaggactt                                            <  2:1190717‑M1/76‑20 (MQ=255)
                     ccAAGTTCTGGTCGATGCAATCAAGAAATACGATATGGATCTGGTGGTTTGTGGTcagtacaccaccaggacttct                                          <  2:1362008‑M1/76‑22 (MQ=255)
                       aaGTTCTGGTCGATGCAATCAAGAAATACGATATGGATCTGGTGGTTTGTGGTcagtacaccaccaggacttctgg                                        <  2:1625909‑M1/76‑24 (MQ=255)
                        aGTTCTGGTCGATGCAATCAAGAAATACGATATGGATCTGGTGGTTTGTGGTcagtacaccaccaggacttctgga                                       >  1:1738394‑M1/1‑52 (MQ=255)
                        aGTTCTGGTCGATGCAATCAAGAAATACGATATGGATCTGGTGGTTTGTGGTcagtacaccaccaggacttctgga                                       <  1:93334‑M1/76‑25 (MQ=255)
                         gTTCTGGTCGATGCAATCAAGAAATACGATATGGATCTGGTGGTTTGTGGTcagtacaccaccaggacttctggag                                      <  2:675115‑M1/76‑26 (MQ=255)
                             tGGTCGATGCAATCAAGAAATACGATATGGATCTGGTGGTTTGTGGTcagtacaccaccaggacttctggagcaaa                                  >  1:1653620‑M1/1‑47 (MQ=255)
                              ggTCGATGCAATCAAGAAATACGATATGGATCTGGTGGTTTGTGGTcagtacaccaccaggacttctggagcaaac                                 >  2:22497‑M1/1‑46 (MQ=255)
                                   aTGCAATCAAGAAATACGATATGGATCTGGTGGTTTGTGGTcagtacaccaccaggacttctggagcaaactgatg                            >  1:435401‑M1/1‑41 (MQ=255)
                                    tGCAATCAAGAAATACGATATGGATCTGGTGGTTTGTGGTcagtacaccaccaggacttctggagcaaactgatgt                           <  2:1140570‑M1/76‑37 (MQ=255)
                                    tGCAATCAAGAAATACGATATGGATCTGGTGGTTTGTGGTcagtacaccaccaggacttctggagcaaactgatgt                           <  1:1474011‑M1/76‑37 (MQ=255)
                                     gCAATCAAGAAATACGATATGGATCTGGTGGTTTGTGGTcagtacaccaccaggacttctggagcaaactgatgtc                          >  2:1244897‑M1/1‑39 (MQ=255)
                                     gCAATCAAGAAATACGATATGGATCTGGTGGTTTGTGGTcagtacaccaccaggacttctggagcaaactgatgtc                          <  2:10773‑M1/76‑38 (MQ=255)
                                        aTCAAGAAATACGATATGGATCTGGTGGTTTGTGGTCACCACCAGGACTTCTGGAGCAAACTGATGTCTTCCGCAc                       >  1:529438/1‑76 (MQ=255)
                                         tCAAGAAATACGATATGGATCTGGTGGTTTGTGGTcagtacaccaccaggacttctggagcaaactgatgtcttcc                      <  1:298618‑M1/76‑42 (MQ=255)
                                          cAAGAAATACGATATGGATCTGGTGGTTTGTGGTcagtacaccaccaggacttctggagcaaactgatgtcttccg                     <  2:945667‑M1/76‑43 (MQ=255)
                                           aaGAAATACGATATGGATCTGGTGGTTTGTGGTgagtacaccaccaggacttctggagcaaactgatgtcttccgc                    <  2:473189‑M1/76‑44 (MQ=255)
                                           aaGAAATACGATATGGATCTGGTGGTTTGTGGTcagtacaccaccaggacttctggagcaaactgatgtcttccgc                    >  1:940276‑M1/1‑33 (MQ=255)
                                           aaGAAATACGATATGGATCTGGTGGTTTGTGGTcagtacaccaccaggacttctggagcaaactgatgtcttccgc                    >  2:1355104‑M1/1‑33 (MQ=255)
                                             gAAATACGATATGGATCTGGTGGTTTGTGGTcagtacaccaccaggacttctggagcaaactgatgtcttccgcac                  <  1:443803‑M1/76‑46 (MQ=255)
                                             gAAATACGATATGGATCTGGTGGTTTGTGGTcagtacaccaccaggacttctggagcaaactgatgtcttccgcac                  >  2:784905‑M1/1‑31 (MQ=255)
                                              aaaTACGATATGGATCTGGTGGTTTGTGGTcagtacaccaccaggacttctggagcaaactgatgtcttccgcacg                 >  1:267665‑M1/1‑30 (MQ=255)
                                                  aCGATATGGATCTGGTGGTTTGTGGTcagtacaccaccaggacttctggagcaaactgatgtcttccgcacgtcag             >  1:834474‑M1/1‑26 (MQ=255)
                                                  aCGATATGGATCTGGTGGTTTGTGGTcagtacaccaccaggacttctggagcaaactgatgtcttccgcacgtcag             <  1:1212095‑M1/76‑51 (MQ=255)
                                                   cGATATGGATCTGGTGGTTTGTGGTcagtacaccaccaggacttctggagcaaactgatgtcttccgcacgtcagc            <  1:414057‑M1/76‑52 (MQ=255)
                                                   cGATATGGATCTGGTGGTTTGTGGTcagtacaccaccaggacttctggagcaaactgatgtcttccgcacgtcagc            <  2:1599441‑M1/76‑52 (MQ=255)
                                                     atatGGATCTGGTGGTTTGTGGTcagtacaccaccaggacttctggagcaaactgatgtcttccgcacgtcagctg          >  1:304818‑M1/1‑23 (MQ=255)
                                                     atatGGATCTGGTGGTTTGTGGTcagtacaccaccaggacttctggagcaaactgatgtcttccgcacgtcagctg          >  2:1033508‑M1/1‑23 (MQ=255)
                                                      tatGGATCTGGTGGTTTGTGGTcagtacaccaccaggacttctggagcaaactgatgtcttccgcacgtcagctga         >  1:257187‑M1/1‑22 (MQ=255)
                                                      tatGGATCTGGTGGTTTGTGGTcagtacaccaccaggacttctggagcaaactgatgtcttccgcacgtcagctga         >  1:1562972‑M1/1‑22 (MQ=255)
                                                      tatGGATCTGGTGGTTTGTGGTcagtacaccaccaggacttctggagcaaactgatgtcttccgcacgtcagctga         >  2:1492320‑M1/1‑22 (MQ=255)
                                                      tatGGATCTGGTGGTTTGTGGTcagtacaccaccaggacttctggagcaaactgatgtcttccgcacgtcagctga         >  1:350490‑M1/1‑22 (MQ=255)
                                                       atGGATCTGGTGGTTTGTGGTcagtacaccaccaggacttctggagcaaactgatgtcttccgcacgtcagctgat        <  2:1399018‑M1/76‑56 (MQ=255)
                                                         ggATCTGGTGGTTTGTGGTcagtacaccaccaggacttctggagcaaactgatgtcttccgcacgtcagctgatca      >  2:658463‑M1/1‑19 (MQ=255)
                                                         ggATCTGGTGGTTTGTGGTcagtacaccaccaggacttctggagcaaactgatgtcttccgcacgtcagctgatca      >  1:1028149‑M1/1‑19 (MQ=255)
                                                          gATCTGGTGGTTTGTGGTcagtacaccaccaggacttctggagcaaactgatgtcttccgcacgtcagctgatcaa     >  1:1793745‑M1/1‑18 (MQ=255)
                                                          gATCTGGTGGTTTGTGGTCACCACCAGGACTTCTGGAGCAAACTGATGTCTTCCGCACGTCAGCTGATCAACACCg     >  2:171054/1‑76 (MQ=255)
                                                           aTCTGGTGGTTTGTGGTcagtacaccaccaggacttctggagcaaactgatgtcttccgcacgtcagctgatcaac    >  2:254203‑M1/1‑17 (MQ=255)
                                                           aTCTGGTGGTTTGTGGTcagtacaccaccaggacttctggagcaaactgatgtcttccgcacgtcagctgatcaac    >  2:821261‑M1/1‑17 (MQ=255)
                                                             cTGGTGGTTTGTGGTcagtacaccaccaggacttctggagcaaactgatgtcttccgcacgtcagctgatcaacac  >  2:580604‑M1/1‑15 (MQ=255)
                                                                           |                                                             
GAGCGGCAGCGGCGACCTGGGCCAAGTTCTGGTCGATGCAATCAAGAAATACGATATGGATCTGGTGGTTTGTGGTCACCACCAGGACTTCTGGAGCAAACTGATGTCTTCCGCACGTCAGCTGATCAACACCGTTC  >  REL606/3571054‑3571190

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