Predicted mutation
evidence seq id position mutation freq annotation gene description
JC JC REL606 969,836 IS150 (+) +3 bp 100% coding (810‑812/2283 nt) pflB ← pyruvate formate lyase I

New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? REL606 = 590471NA (NA)78 (0.820) 44/80 NT 98.6% noncoding (1977/1977 nt) IS150 repeat region
?REL606 969836 = 1 (0.010)coding (812/2283 nt) pflB pyruvate formate lyase I
* ? REL606 664688 =NA (NA)78 (0.870) 44/76 NT 98.6% noncoding (1/1443 nt) IS150 repeat region
?REL606 = 969838 1 (0.010)coding (810/2283 nt) pflB pyruvate formate lyase I

Only 100 of 259 total aligned reads displayed.

GACCTATATGCCTCGTGTTTAACTGTCCAACTTTTTGGGGTCAGTACACCCTGGGGCTGCCACTCCAAACCCGAACAATTTGGATGGTAGTCCCTTCT  >  REL606/590424‑590521
                                               |                                                  
  ccTATATGCCTCGTGTTTAACTGTCCAACTTTTTGGGGTCAGTACACCCTg                                               >  1:10156163/1‑51 (MQ=38)
   cTATATGCCTCGTGTTTAACTGTCCAACTTTTTGGGGTCAGTACACattgg                                              >  1:2124492‑M2/1‑46 (MQ=255)
   cTATATGCCTCGTGTTTAACTGTCCAACTTTTTGGGGTCAGTACACattgg                                              <  1:8426630‑M2/51‑6 (MQ=255)
   cTATATGCCTCGTGTTTAACTGTCCAACTTTTTGGGGTCAGTACACCCTgg                                              <  1:3565644/51‑1 (MQ=38)
   cTATATGCCTCGTGTTTAACTGTCCAACTTTTTGGGGTCAGTACACCCTgg                                              >  1:521504/1‑51 (MQ=40)
   cTATATGCCTCGTGTTTAACTGTCCAACTTTTTGGGGTCAGTACACCCTgg                                              <  1:6853844/51‑1 (MQ=38)
    tataTGCCTCGTGTTTAACTGTCCAACTTTTTGGGGTCAGTACACCCTggg                                             >  1:3114948/1‑51 (MQ=40)
    tataTGCCTCGTGTTTAACTGTCCAACTTTTTGGGGTCAGTACACCCTggg                                             >  1:9871884/1‑51 (MQ=255)
     atatGCCTCGTGTTTAACTGTCCAACTTTTTGGGGTCAGTACACattgggg                                            <  1:9149450‑M2/51‑8 (MQ=255)
      tatGCCTCGTGTTTAACTGTCCAACTTTTTGGGGTCAGTACACattggggc                                           >  1:10616848‑M2/1‑43 (MQ=255)
      tatGCCTCGTGTTTAACTGTCCAACTTTTTGGGGTCAGTACACCCTGGGGc                                           <  1:5231565/51‑1 (MQ=255)
       atGCCTCGTGTTTAACTGTCCAACTTTTTGGGGTCAGTACACCCTGGGGCt                                          <  1:1183745/51‑1 (MQ=255)
       atGCCTCGTGTTTAACTGTCCAACTTTTTGGGGTCAGTACACCCTGGGGCt                                          <  1:3705195/51‑1 (MQ=255)
         gCCTCGTGTTTAACTGTCCAACTTTTTGGGGTCAGTACACattggggctta                                        >  1:6398612‑M2/1‑40 (MQ=255)
         gCCTCGTGTTTAACTGTCCAACTTTTTGGGGTCAGTACACattggggctta                                        <  1:9641294‑M2/51‑12 (MQ=255)
         gCCTCGTGTTTAACTGTCCAACTTTTTGGGGTCAGTACACattggggctta                                        >  1:2096949‑M2/1‑40 (MQ=255)
         gCCTCGTGTTTAACTGTCCAACTTTTTGGGGTCAGTACACattggggctta                                        <  1:9790238‑M2/51‑12 (MQ=255)
         gCCTCGTGTTTAACTGTCCAACTTTTTGGGGTCAGTACACCCTGGGGCTGc                                        >  1:10334395/1‑51 (MQ=255)
         gCCTCGTGTTTAACTGTCCAACTTTTTGGGGTCAGTACACCCTGGGGCTGc                                        >  1:1944371/1‑51 (MQ=255)
          ccTCGTGTTTAACTGTCCAACTTTTTGGGGTCAGTACACattggggcttaa                                       <  1:3672636‑M2/51‑13 (MQ=255)
          ccTCGTGTTTAACTGTCCAACTTTTTGGGGTCAGTACACattggggcttaa                                       >  1:5401252‑M2/1‑39 (MQ=255)
          ccTCGTGTTTAACTGTCCAACTTTTTGGGGTCAGTACACCCTGGGGCTGcc                                       >  1:2992054/1‑51 (MQ=255)
           cTCGTGTTTAACTGTCCAACTTTTTGGGGTCAGTACACCCTGGGGCTGCCa                                      <  1:3328180/51‑1 (MQ=255)
           cTCGTGTTTAACTGTCCAACTTTTTGGGGTCAGTACACCCTGGGGCTGCCa                                      <  1:6136895/51‑1 (MQ=255)
            tCGTGTTTAACTGTCCAACTTTTTGGGGTCAGTACACCCTGGGGCTGCCAc                                     >  1:4980311/1‑51 (MQ=255)
              gtgtTTAACTGTCCAACTTTTTGGGGTCAGTACACattggggcttaagcgc                                   >  1:8880578‑M2/1‑35 (MQ=255)
              gtgtTTAACTGTCCAACTTTTTGGGGTCAGTACACCCTGGGGCTGCCACTc                                   <  1:1253842/51‑1 (MQ=255)
              gtgtTTAACTGTCCAACTTTTTGGGGTCAGTACACCCTGGGGCTGCCACTc                                   >  1:8568193/1‑51 (MQ=255)
               tgtTTAACTGTCCAACTTTTTGGGGTCAGTACACattggggcttaagcgcg                                  <  1:8080060‑M2/51‑18 (MQ=255)
               tgtTTAACTGTCCAACTTTTTGGGGTCAGTACACCCTGGGGCTGCCACTcc                                  >  1:171351/1‑51 (MQ=255)
               tgtTTAACTGTCCAACTTTTTGGGGTCAGTACACCCTGGGGCTGCCACTcc                                  <  1:4445727/51‑1 (MQ=255)
                gtTTAACTGTCCAACTTTTTGGGGTCAGTACACattggggcttaagcgcgc                                 >  1:1003660‑M2/1‑33 (MQ=255)
                gtTTAACTGTCCAACTTTTTGGGGTCAGTACACattggggcttaagcgcgc                                 >  1:6692189‑M2/1‑33 (MQ=255)
                 tttAACTGTCCAACTTTTTGGGGTCAGTACACCCTGGGGCTGCCACTCCaa                                >  1:6001076/1‑51 (MQ=255)
                 tttAACTGTCCAACTTTTTGGGGTCAGTACACCCTGGGGCTGCCACTCCaa                                >  1:10883901/1‑51 (MQ=255)
                 tttAACTGTCCAACTTTTTGGGGTCAGTACACCCTGGGGCTGCCACTCCaa                                >  1:2287266/1‑51 (MQ=255)
                  ttAACTGTCCAACTTTTTGGGGTCAGTACACattggggcttaagcgcgcgc                               <  1:3896983‑M2/51‑21 (MQ=255)
                  ttAACTGTCCAACTTTTTGGGGTCAGTACACattggggcttaagcg                                    >  1:4895184‑M2/1‑31 (MQ=255)
                  ttAACTGTCCAACTTTTTGGGGTCAGTACACattggggcttaagcg                                    >  1:8189043‑M2/1‑31 (MQ=255)
                  ttAACTGTCCAACTTTTTGGGGTCAGTACACattggggcttaagcg                                    >  1:4495882‑M2/1‑31 (MQ=255)
                  ttAACTGTCCAACTTTTTGGGGTCAGTACACCCTGGGGCTGCCACTCCaaa                               <  1:6607591/51‑1 (MQ=255)
                   tAACTGTCCAACTTTTTGGGGTCAGTACACattggggcttaagcgcgcgct                              <  1:2027046‑M2/51‑22 (MQ=255)
                   tAACTGTCCAACTTTTTGGGGTCAGTACACCCTGGGGCTGCCACTCCAAAc                              <  1:2824112/51‑1 (MQ=255)
                   tAACTGTCCAACTTTTTGGGGTCAGTACACCCTGGGGCTGCCACTCCAAAc                              >  1:2575498/1‑51 (MQ=255)
                    aaCTGTCCAACTTTTTGGGGTCAGTACACCCTGGGGCTGCCACTCCAAAcc                             <  1:2756453/51‑1 (MQ=255)
                     aCTGTCCAACTTTTTGGGGTCAGTACACattggggcttaagcgcgcgcttc                            >  1:9251749‑M2/1‑28 (MQ=255)
                     aCTGTCCAACTTTTTGGGGTCAGTACACattggggcttaagcgcgcgcttc                            >  1:6364666‑M2/1‑28 (MQ=255)
                     aCTGTCCAACTTTTTGGGGTCAGTACACattggggcttaagcgcgcgcttc                            >  1:1442235‑M2/1‑28 (MQ=255)
                       tGTCCAACTTTTTGGGGTCAGTACACattggggcttaagcgcgcgcttcaa                          <  1:8356440‑M2/51‑26 (MQ=255)
                       tGTCCAACTTTTTGGGGTCAGTACACCCTGGGGCTGCCACTCCAAACCCGa                          >  1:8776001/1‑51 (MQ=255)
                        gTCCAACTTTTTGGGGTCAGTACACattggggcttaagcgcgcgcttcaac                         <  1:1895223‑M2/51‑27 (MQ=255)
                        gTCCAACTTTTTGGGGTCAGTACACCCTGGGGCTGCCACTCCAAACCCGaa                         >  1:10437437/1‑51 (MQ=255)
                         tCCAACTTTTTGGGGTCAGTACACCCTGGGGCTGCCACTCCAAACCCGAAc                        <  1:2870504/51‑1 (MQ=255)
                           cAACTTTTTGGGGTCAGTACACCCTGGGGCTGCCACTCCAAACCCGAACaa                      <  1:4678331/51‑1 (MQ=255)
                           cAACTTTTTGGGGTCAGTACACCCTGGGGCTGCCACTCCAAACCCGAACaa                      <  1:7412440/51‑1 (MQ=255)
                            aaCTTTTTGGGGTCAGTACACattggggcttaagcgcgcgcttcaactcac                     <  1:6938500‑M2/51‑31 (MQ=255)
                            aaCTTTTTGGGGTCAGTACACattggggcttaagcgcgcgcttcaactcac                     <  1:6496282‑M2/51‑31 (MQ=255)
                            aaCTTTTTGGGGTCAGTACACattggggcttaagcgcgcgcttcaactcac                     <  1:5822169‑M2/51‑31 (MQ=255)
                            aaCTTTTTGGGGTCAGTACACattggggcttaagcgcgcgcttcaactcac                     <  1:4275553‑M2/51‑31 (MQ=255)
                            aaCTTTTTGGGGTCAGTACACCCTGGGGCTGCCACTCCAAACCCGAACAAt                     >  1:5728620/1‑51 (MQ=255)
                            aaCTTTTTGGGGTCAGTACACCCTGGGGCTGCCACTCCAAACCCGAACAAt                     <  1:7959003/51‑1 (MQ=255)
                            aaCTTTTTGGGGTCAGTACACCCTGGGGCTGCCACTCCAAACCCGAACAAt                     >  1:9504791/1‑51 (MQ=255)
                            aaCTTTTTGGGGTCAGTACACCCTGGGGCTGCCACTCCAAACCCGAACAAt                     >  1:3753573/1‑51 (MQ=255)
                             aCTTTTTGGGGTCAGTACACattggggcttaagcgcgcgcttcaactcacc                    <  1:9550063‑M2/51‑32 (MQ=255)
                             aCTTTTTGGGGTCAGTACACCCTGGGGCTGCCACTCCAAACCCGAACAAtt                    <  1:5684127/51‑1 (MQ=255)
                              cTTTTTGGGGTCAGTACACattggggcttaagcgcgcgcttcaactcaccc                   <  1:4758722‑M2/51‑33 (MQ=255)
                               tttttGGGGTCAGTACACattggggcttaagcgcgcgcttcaactcacccc                  <  1:426337‑M2/51‑34 (MQ=255)
                               tttttGGGGTCAGTACACattggggcttaagcgcgcgcttcaactcacccc                  <  1:6836304‑M2/51‑34 (MQ=255)
                                ttttGGGGTCAGTACACCCTGGGGCTGCCACTCCAAACCCGAACAAt                     <  1:4037911/47‑1 (MQ=255)
                                ttttGGGGTCAGTACACCCTGGGGCTGCCACTCCAAACCCGAACAATTTgg                 >  1:2660387/1‑51 (MQ=255)
                                 tttGGGGTCAGTACACattggggcttaagcgcgcgcttcaactcacccctg                <  1:5640748‑M2/51‑36 (MQ=255)
                                 tttGGGGTCAGTACACattggggcttaagcgcgcgcttcaactcacccctg                <  1:10537789‑M2/51‑36 (MQ=255)
                                 tttGGGGTCAGTACACCCTGGGGCTGCCACTCCAAACCCGAACAATTTGGa                >  1:5023078/1‑51 (MQ=255)
                                 tttGGGGTCAGTACACCCTGGGGCTGCCACTCCAAACCCGAACAATTTGGa                <  1:8862699/51‑1 (MQ=255)
                                  ttGGGGTCAGTACACattggggcttaagcgcgcgcttcaactcacccctgc               <  1:8998031‑M2/51‑37 (MQ=255)
                                  ttGGGGTCAGTACACattggggcttaagcgcgcgcttcaactcacccctgc               >  1:5788684‑M2/1‑15 (MQ=255)
                                   tGGGGTCAGTACACattggggcttaagcgcgcgcttcaactcacccctgca              >  1:404386‑M2/1‑14 (MQ=255)
                                   tGGGGTCAGTACACCCTGGGGCTGCCACTCCAAACCCGAACAATTTGGATg              <  1:7029537/51‑1 (MQ=255)
                                    ggggTCAGTACACattggggcttaagcgcgcgcttcaactcacccctgcaa             >  1:4335306‑M2/1‑13 (MQ=255)
                                    ggggTCAGTACACCCTGGGGCTGCCACTCCAAACCCGAACAATTTGGATgg             <  1:9094839/51‑1 (MQ=255)
                                     gggTCAGTACACattggggcttaagcgcgcgcttcaactcacccctgcaat            >  1:1872016‑M2/1‑12 (MQ=255)
                                     gggTCAGTACACattggggcttaagcgcgcgcttcaactcacccctgcaat            <  1:7160641‑M2/51‑40 (MQ=255)
                                      ggTCAGTACACCCTGGGGCTGCCACTCCAAACCCGAACAATTTGGATGGTa           <  1:1638856/51‑1 (MQ=255)
                                       gTCAGTACACattggggcttaagcgcgcgcttcaactcacccctgcaatat          >  1:2355436‑M2/1‑10 (MQ=255)
                                       gTCAGTACACattggggcttaagcgcgcgcttcaactcacccctgcaatat          >  1:4807175‑M2/1‑10 (MQ=255)
                                       gTCAGTACACCCTGGGGCTGCCACTCCAAACCCGAACAATTTGGATGGTAg          <  1:8023653/51‑1 (MQ=255)
                                       gTCAGTACACCCTGGGGCTGCCACTCCAAACCCGAACAATTTGGATGGTAg          <  1:7593254/51‑1 (MQ=255)
                                        tCAGTACACattggggcttaagcgcgcgcttcaactcacccctgcaatatt         <  1:3106259‑M2/51‑43 (MQ=255)
                                         cAGTACACattggggcttaagcgcgcgcttcaactcacccctgcaatattt        <  1:1988575‑M2/51‑44 (MQ=255)
                                          aGTACACattggggcttaagcgcgcgcttcaactcacccctgcaatatttt       <  1:7638693‑M2/51‑45 (MQ=255)
                                          aGTACACCCTGGGGCTGCCACTCCAAACCCGAACAATTTGGATGGTAGTcc       >  1:3829854/1‑51 (MQ=255)
                                           gTACACCCTGGGGCTGCCACTCCAAACCCGAACAATTTGGATGGTAg          >  1:7266687/1‑47 (MQ=255)
                                           gTACACCCTGGGGCTGCCACTCCAAACCCGAACAATTTGGATGGTAGTccc      <  1:9123388/51‑1 (MQ=255)
                                            tACACattggggcttaagcgcgcgcttcaactcacccctgcaatattttgt     <  1:10836098‑M2/51‑47 (MQ=255)
                                            tACACCCTGGGGCTGCCACTCCAAACCCGAACAATTTGGATGGTAGTCCct     >  1:86421/1‑51 (MQ=255)
                                              cacattggggcttaagcgcgcgcttcaactcacccctgcaatattttgtcc   >  1:4006105‑M2/1‑3 (MQ=255)
                                               acCCTGGGGCTGCCACTCCAAACCCGAACAATTTGGATGGTAGTCCcttct  >  1:916885/1‑51 (MQ=255)
                                               acCCTGGGGCTGCCACTCCAAACCCGAACAATTTGGATGGTAGTCCcttct  >  1:1281307/1‑51 (MQ=255)
                                               acCCTGGGGCTGCCACTCCAAACCCGAACAATTTGGATGGTAGTCCcttct  >  1:7757644/1‑51 (MQ=255)
                                               acCCTGGGGCTGCCACTCCAAACCCGAACAATTTGGATGGTAGTCCcttct  >  1:4572930/1‑51 (MQ=255)
                                               |                                                  
GACCTATATGCCTCGTGTTTAACTGTCCAACTTTTTGGGGTCAGTACACCCTGGGGCTGCCACTCCAAACCCGAACAATTTGGATGGTAGTCCCTTCT  >  REL606/590424‑590521

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