New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? minE 957613 =51 (0.980)3 (0.070) 3/90 NT 7.0% coding (2693/2796 nt) pqqL predicted peptidase
?minE 957637 = 35 (0.770)coding (2669/2796 nt) pqqL predicted peptidase

TGTTAACATTCTCTACCGTCATTTGTTTCAACAATTGCTCCTGCTCAGTCCATGCTGCAGGATCGTCATATTGAATAAGACTATTTACAATAGTGTTCGCTAATTGTTGAACGCTACGCTGTTGGATATCGAG  >  minE/957575‑957707
                                                              |                                                                      
tgtTAACATTCTCTACCGTCATTTGTTTCAACAATTGCTCCTGCTCAGTCCATGCTGCAGGATCGTCatat                                                                >  1:1935926/1‑71 (MQ=255)
    aaCATTCTCTACCGTCATTTGTTTCAACAATTGCTCCTGCTCAGTCCATGCTGCAGGATCGTCATATTGa                                                             <  1:52451/70‑1 (MQ=255)
     aCATTCTCTACCGTCATTTGTTTCAACAATTGCTCCTGCTCAGTCCATGCTGCAGGATCGTCATATTGAAt                                                           <  1:1979973/71‑1 (MQ=255)
      cATTCTCTACCGTCATTTGTTTCAACAATTGCTCCTGCTCAGTCCATGCTGCAGGATCGTCATATTGAAt                                                           >  1:3374844/1‑70 (MQ=255)
       aTTCTCTACCGTCATTTGTTTCAACAATTGCTCCTGCTCAGTCCATGCTGCAGGATCGTCATATTGAATaa                                                         <  1:2962759/71‑1 (MQ=255)
       aTTCTCTACCGTCATTTGTTTCAACAATTGCTCCTGCTCAGTCCATGCTGCAGGATCGTCATATTGAATaa                                                         <  1:3735992/71‑1 (MQ=255)
             tACCGTCATTTGTTTCAACAATTGCTCCTGCTCAGTCCATGCTGCAGGATCGTCATATTGAATAAGACt                                                     <  1:401965/69‑1 (MQ=255)
               ccGTCATTTGTTTCAACAATTGCTCCTGCTCAGTCCATGCTGCAGGATCGTCATATTGAATAAGACTAtt                                                  <  1:1157115/70‑1 (MQ=255)
               ccGTCATTTGTTTCAACAATTGCTCCTGCTCAGTCCATGCTGCAGGATCGTCATATTGAATAAGACTAtt                                                  <  1:681145/70‑1 (MQ=255)
                 gTCATTTGTTTCAACAATTGCTCCTGCTCAGTCCATGCTGCAGGATCGTCATATTGAATAAGACTATTTAc                                               <  1:3272556/71‑1 (MQ=255)
                   cATTTGTTTCAACAATTGCTCCTGCTCAGTCCATGCTGCAGGATCGTCATATTGAATAAGACTATTTACaa                                             >  1:2580715/1‑71 (MQ=255)
                   cATTTGTTTCAACAATTGCTCCTGCTCAGTCCATGCTGCAGGATCGTCATATTGAATAAGACTATTTACaa                                             >  1:3135054/1‑71 (MQ=255)
                    aTTTGTTTCAACAATTGCTCCTGCTCAGTCCATGCTGCAGGATCGTCATATTGAATAAGACt                                                     <  1:1872733/62‑1 (MQ=255)
                    aTTTGTTTCAACAATTGCTCCTGCTCAGTCCATGCTGCAGGATCGTCATATTGAATAAGACTATTTACaa                                             <  1:1100099/70‑1 (MQ=255)
                    aTTTGTTTCAACAATTGCTCCTGCTCAGTCCATGCTGCAGGATCGTCATATTGAATAAGACTATTTACaa                                             <  1:2178685/70‑1 (MQ=255)
                    aTTTGTTTCAACAATTGCTCCTGCTCAGTCCATGCTGCAGGATCGTCATATTGAATAAGACTATTTACaa                                             <  1:717242/70‑1 (MQ=255)
                    aTTTGTTTCAACAATTGCTCCTGCTCAGTCCATGCTGCAGGATCGTCATATTGAATAAGACTATTTACaa                                             <  1:1441221/70‑1 (MQ=255)
                    aTTTGTTTCAACAATTGCTCCTGCTCAGTCCATGCTGCAGGATCGTCATATTGAATAAGACTATTTACaa                                             <  1:385746/70‑1 (MQ=255)
                    aTTTGTTTCAACAATTGCTCCTGCTCAGTCCATGCTGCAGGATCGTCATATTGAATAAGACTATTTACaa                                             <  1:1384710/70‑1 (MQ=255)
                        gTTTCAACAATTGCTCCTGCTCAGTCCATGCTGCAGGATCGTCATATTGAATAAGACTATTTACAATAgtg                                        >  1:2038115/1‑71 (MQ=255)
                        gTTTCAACAATTGCTCCTGCTCAGTCCATGCTGCAGGATCGTCATATTGAATAAGACTATTTACAATAgt                                         <  1:1669390/70‑1 (MQ=255)
                         tttCAACAATTGCTCCTGCTCAGTCCATGCTGCAGGATCGTCATATTGAATAAGACTATTTACAATAgtgt                                       >  1:653566/1‑71 (MQ=255)
                         tttCAACAATTGCTCCTGCTCAGTCCATGCTGCAGGATCGTCATATTGAATAAGACTATTTACAATAgtgt                                       <  1:2522405/71‑1 (MQ=255)
                          ttCAACAATTGCTCCTGCTCAGTCCATGCTGCAGGATCGTCATATTGAAt                                                           <  1:3538756/50‑1 (MQ=255)
                           tCAACAATTGCTCCTGCTCAGTCCATGCTGCAGGATCGTCATATTGAATAAGACTATTTACAATAg                                          <  1:3076477/66‑1 (MQ=255)
                            caacaaTTGCTCCTGCTCAGTCCATGCTGCAGGATCGTCATATTGAATAAGACTATTTACAATAGTGTTc                                     <  1:3952366/70‑1 (MQ=255)
                             aacaaTTGCTCCTGCTCAGTCCATGCTGCAGGATCGTCATATTGAATaa                                                         <  1:1133419/49‑1 (MQ=255)
                             aacaaTTGCTCCTGCTCAGTCCATGCTGCAGGATCGTCATATTGAATAAGACTATTTACAATAGTGTTCGc                                   >  1:3116328/1‑71 (MQ=255)
                              acaaTTGCTCCTGCTCAGTCCATGCTGCAGGATCGTCATATTGAATAAGACTATTTACAATAGTGTTCGc                                   >  1:1012953/1‑70 (MQ=255)
                                    gCTCCTGCTCAGTCCATGCTGCAGGATCGTCATATTGAATAAGACTATTTACAATa                                           >  1:890985/1‑56 (MQ=255)
                                      tcctgcagcatggactgagcaggaTCGTCATATTGaa                                                            <  1:3694331‑M2/13‑1 (MQ=255)
                                       CCTGCTCAGTCCATGCTGCAGGATCGTCATATTGAATAAGACTATTTACAATAGTGTTCGCTAAttgttg                          <  1:2357699‑M1/70‑1 (MQ=255)
                                        ctgcagcatggactgagcaggaTCGTCATATTGAATa                                                          >  1:2604336‑M2/23‑37 (MQ=255)
                                        CTGCTCAGTCCATGCTGCAGGATCGTCATATTGAATAAGACTATTTAc                                               >  1:3190479‑M1/1‑48 (MQ=255)
                                        CTGCTCAGTCCATGCTGCAGGATCGTCATATTGAATAAGACTATTTACAATAGTGTTCGCTAATTGTTGaa                        <  1:2949720‑M1/71‑1 (MQ=255)
                                        CTGCTCAGTCCATGCTGCAGGATCGTCATATTGAATAAGACTATTTACAATAGTGTTCGCTAATTGTTGaa                        >  1:610267‑M1/1‑71 (MQ=255)
                                         TGCTCAGTCCATGCTGCAGGATCGTCATATTGAATAAGACTATTTACAATAGTGTTCGCTAATTGTTGAAc                       <  1:3051049‑M1/71‑1 (MQ=255)
                                          GCTCAGTCCATGCTGCAGGATCGTCATATTGAATAAGACTATTTACAATAGTGTTCGCTAATTGTTGAACg                      <  1:2001824‑M1/71‑1 (MQ=255)
                                           cagcatggactgagcaggaTCGTCATATTGAATAAGACt                                                     >  1:445808‑M2/20‑39 (MQ=255)
                                           CTCAGTCCATGCTGCAGGATCGTCATATTGAATAAGACTa                                                    >  1:1990194‑M1/1‑40 (MQ=255)
                                           CTCAGTCCATGCTGCAGGATCGTCATATTGAATAAGACTATTTAc                                               <  1:1252651‑M1/45‑1 (MQ=255)
                                            tCAGTCCATGCTGCAGGATCGTCATATTGAATAAGACTATTTACAATAGTGTTCGCTa                                 <  1:2118237/58‑1 (MQ=255)
                                            tCAGTCCATGCTGCAGGATCGTCATATTGAATAAGACTATTTACAATAGTGTTCGCTa                                 <  1:3614863/58‑1 (MQ=255)
                                             cAGTCCATGCTGCAGGATCGTCATATTGAATAAGACTATTTACAATAGTGTTc                                     <  1:1083119/53‑1 (MQ=255)
                                             cAGTCCATGCTGCAGGATCGTCATATTGAATAAGACTATTTACAATAGTGTTCGCTAATTGTTGAacgcta                   >  1:2969582/1‑71 (MQ=255)
                                                  cATGCTGCAGGATCGTCATATTGAATAAGACTATTTACAATAGTGTTCGCTAATTGTTGAACGCTACGCTg              <  1:3245032/71‑1 (MQ=255)
                                                   aTGCTGCAGGATCGTCATATTGAATAAGACTATTTACAATAgtg                                        <  1:3263566/44‑1 (MQ=255)
                                                   aTGCTGCAGGATCGTCATATTGAATAAGACTATTTACAATAgtg                                        <  1:1651637/44‑1 (MQ=255)
                                                    tgctgcAGGATCGTCATATTGAATAAGACTATTTACAATAGTGTTCGCTAATTGTTGAACGCTACGCTGtt            >  1:3513179/1‑71 (MQ=255)
                                                    tgctgcAGGATCGTCATATTGAATAAGACTATTTACAATAGTGTTCGCTAATTGTTGAACGCTACGCTGtt            >  1:2123492/1‑71 (MQ=255)
                                                           ggaTCGTCATATTGAATAAGACTATTTACAATAGTGTTCGCTAATTGTTGAACGCTACGCTGTTGGATATc     <  1:468527‑M2/68‑1 (MQ=255)
                                                            gaTCGTCATATTGAATAAGACTATTTACAATAGTGTTCGCTAAttgttg                          <  1:33271‑M2/47‑1 (MQ=255)
                                                            gaTCGTCATATTGAATAAGACTATTTACAATAGTGTTCGCTAAttgttg                          <  1:767254‑M2/47‑1 (MQ=255)
                                                            gaTCGTCATATTGAATAAGACTATTTACAATAGTGTTCGCTAATTGTTGAACGCTACGCTGt             <  1:1813647‑M2/60‑1 (MQ=255)
                                                            gaTCGTCATATTGAATAAGACTATTTACAATAGTGTTCGCTAATTGTTGAACGCTACGCTGTTGGATATc     >  1:1396762‑M2/3‑70 (MQ=255)
                                                              tCGTCATATTGAATAAGACTATTTACAATAGTGTTCGCTAATTGTTGAACGCTACGCTGTTGGATATCGAg  <  1:1885048/71‑1 (MQ=255)
                                                              |                                                                      
TGTTAACATTCTCTACCGTCATTTGTTTCAACAATTGCTCCTGCTCAGTCCATGCTGCAGGATCGTCATATTGAATAAGACTATTTACAATAGTGTTCGCTAATTGTTGAACGCTACGCTGTTGGATATCGAG  >  minE/957575‑957707

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