New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? USA300TCH1516_ALE = 2850830262 (0.730)5 (0.020) 4/132 NT 2.2% intergenic (+29/+90) pcp/USA300HOU_RS14605 Pyrrolidone‑carboxylate peptidase/hypothetical protein
?USA300TCH1516_ALE = 2850851 239 (0.810)intergenic (+50/+69) pcp/USA300HOU_RS14605 Pyrrolidone‑carboxylate peptidase/hypothetical protein
Rejected: Frequency below/above cutoff threshold.

Only 100 of 508 total aligned reads displayed.

CAGCCATTGAAGCTATTTCTAATGATGAGGATTTACACCTAGCCCTGGGCACAACGGAATAGTTATAAGTGCACAAACAAAAGCACCCTAATTTGTTAGTTTAAGTCTAACTTTAGGGTGCTTTTGGATGATTGGTATACAGTTTTACTATAATGATGTATTTTAGCTAAACGACATTGTT  >  USA300TCH1516_ALE/2850740‑2850920
                                                                                          |                                                                                          
ggaagttaaaaaaagaagactaaacaatgtcgtttagctaaaatacatcattatagtaaaactgtataccaatcatcCAAAAGCACCCTaa                                                                                            <  1:4485485‑M1/14‑1 (MQ=255)
ggaagttaaaaaaagaagactaaacaatgtcgtttagctaaaatacatcattatagtaaaactgtataccaatcatcCAAAAGCACCCTaa                                                                                            <  2:6200183‑M1/14‑1 (MQ=255)
    cATTGAAGCTATTTCTAATGATGAGGATTTACACCTAGCCCTGGGCACAACGGAATAGTTATAAGTGCACAAACAAAAGCACCCTAATTTg                                                                                        <  2:4244225/91‑1 (MQ=255)
    cATTGAAGCTATTTCTAATGATGAGGATTTACACCTAGCCCTGGGCACAACGGAATAGTTATAAGTGCACAAACAAAAGCACCCTAATTTg                                                                                        <  1:1000929/91‑1 (MQ=255)
    cATTGAAGCTATTTCTAATGATGAGGATTTACACCTAGCCCTGGGCACAACGGAATAGTTATAAGTGCACAAACAAAAGCACCCTAATTTg                                                                                        <  2:354039/91‑1 (MQ=255)
    cATTGAAGCTATTTCTAATGATGAGGATTTACACCTAGCCCTGGGCACAACGGAATAGTTATAAGTGCACAAACAAAAGCACCCTAATTTg                                                                                        <  1:402207/91‑1 (MQ=255)
     aTTGAAGCTATTTCTAATGATGAGGATTTACACCTAGCCCTGGGCACAACGGAATAGTTATAAGTGCACAAACAAAAGCACCCTAATTTGt                                                                                       <  1:1613582/91‑1 (MQ=255)
     aTTGAAGCTATTTCTAATGATGAGGATTTACACCTAGCCCTGGGCACAACGGAATAGTTATAAGTGCACAAACAAAAGCACCCTAATTTGt                                                                                       <  2:6553815/91‑1 (MQ=255)
      ttGAAGCTATTTCTAATGATGAGGATTTACACCTAGCCCTGGGCACAACGGAATAGTTATAAGTGCACAAACAAAAGCACCCTAATTTGtt                                                                                      >  1:1276347/1‑91 (MQ=255)
       tGAAGCTATTTCTAATGATGAGGATTTACACCTAGCCCTGGGCACAACGGAATAGTTATAAGTGCACAAACAAAAGCACCCTAATTTGTTa                                                                                     >  2:643543/1‑91 (MQ=255)
       tGAAGCTATTTCTAATGATGAGGATTTACACCTAGCCCTGGGCACAACGGAATAGTTATAAGTGCACAAACAAAAGCACCCTAATTTGTTa                                                                                     >  2:1262607/1‑91 (MQ=255)
        gAAGCTATTTCTAATGATGAGGATTTACACCTAGCCCTGGGCACAACGGAATAGTTATAAGTGCACAAACAAAAGCACCCTAATTTGTTAg                                                                                    >  2:3778849/1‑91 (MQ=255)
        gAAGCTATTTCTAATGATGAGGATTTACACCTAGCCCTGGGCACAACGGAATAGTTATAAGTGCACAAACAAAAGCACCCTAATTTGTTAg                                                                                    <  2:1678719/91‑1 (MQ=255)
         aaGCTATTTCTAATGATGAGGATTTACACCTAGCCCTGGGCACAACGGAATAGTTATAAGTGCACAAACAAAAGCACCCTAATTTGTTAGt                                                                                   >  1:5609882/1‑91 (MQ=255)
         aaGCTATTTCTAATGATGAGGATTTACACCTAGCCCTGGGCACAACGGAATAGTTATAAGTGCACAAACAAAAGCACCCTAATTTGTTAGt                                                                                   <  2:4551832/91‑1 (MQ=255)
          aGCTATTTCTAATGATGAGGATTTACACCTAGCCCTGGGCACAACGGAATAGTTATAAGTGCACAAACAAAAGCACCCTAATTTGTTAGtt                                                                                  >  2:5754285/1‑91 (MQ=255)
          aGCTATTTCTAATGATGAGGATTTACACCAAGCCCTGGGCACAACGGAATAGTTATAAGTGCACAAACAAAAGCACCCTAATTTCTTatct                                                                                  >  1:5862511/1‑88 (MQ=255)
            cTATTTCTAATGATGAGGATTTACACCTAGCCCTGGGCACAACGGAATAGTTATAAGTGCACAAACAAAAGCACCCTAATTTGTTAGTTTa                                                                                >  2:2814849/1‑91 (MQ=255)
            cTATTTCTAATGATGAGGATTTACACCTAGCCCTGGGCACAACGGAATAGTTATAAGTGCACAAACAAAAGCACCCTAATTTGTTAGTTTa                                                                                <  2:4914946/91‑1 (MQ=255)
              aTTTCTAATGATGAGGATTTACACCTAGCCCTGGGCACAACGGAATAGTTATAAGTGCACAAACAAAAGCGCCCTAATTTGTTAGTTTAAg                                                                              <  2:3177676/91‑1 (MQ=255)
               tttCTAATGATGAGGATTTACACCTAGCCCTGGGCACAACGGAATAGTTATAAGTGCACAAACAAAAGCACCCTAATTTGTTAGTTTAAGt                                                                             <  2:2892987/91‑1 (MQ=255)
                ttCTAATGATGAGGATTTACACCTAGCCCTGGGCACAACGGAATAGTTATAAGTGCACAAACAAAAGCACCCTAATTTGTTAGTTTAAGTc                                                                            <  2:4412021/91‑1 (MQ=255)
                 tCTAATGATGAGGATTTACACCTAGCCCTGGGCACAACGGAATAGTTATAAGTGCACAAACAAAAGCACCCTAATTTGTTAGTTTAAGTCt                                                                           <  2:2173847/91‑1 (MQ=255)
                     atgatgAGGATTTACACCTAGCCCTGGGCACAACGGAATAGTTATAAGTGCACAAACAAAAGCACCCTAATTTGTTAGTTTAAGTCTAACt                                                                       <  2:3052784/91‑1 (MQ=255)
                     atgatgAGGATTTACACCTAGCCCTGGGCACAACGGAATAGTTATAAGTGCACAAACAAAAGCACCCTAATTTGTTAGTTTAAGTCTAACt                                                                       >  1:4717965/1‑91 (MQ=255)
                       gatgaGGATTTACACCTAGCCCTGGGCACAACGGAATAGTTATAAGTGCACAAACAAAAGCACCCTAATTTGTTAGTTTAAGTCTAACTtt                                                                     <  1:3513602‑M2/91‑3 (MQ=255)
                        caatgtcgtttagctaaaatacatcattatagtaaaactgtataccaatcatcCAAAAGCACCCTAATTTGTTAGTTTAAGTCTAACTTTa                                                                    >  1:4202888‑M1/54‑91 (MQ=255)
                         tgaGGATTTACACCTAGCCCTGGGCACAACGGAATAGTTATAAGTGCACAAACAAAAGCACCCTAATTTGTTAGTTTAAGTCTAACTttag                                                                   <  2:467097‑M2/91‑5 (MQ=255)
                           aGGATTTACACCTAGCCCTGGGCACAACGGAATAGTTATAAGTGCACAAACAAAAGCACCCTAATTTGTTAGTTTAAGTCTAACTttaggg                                                                 >  1:3739256‑M2/1‑85 (MQ=255)
                            ggATTTACACCTAGCCCTGGGCACAACGGAATAGTTATAAGTGCACAAACAAAAGCACCCTAATTTGTTAGTTTAAGTCTAACTttagggt                                                                >  1:3633961‑M2/1‑84 (MQ=255)
                              cgtttagctaaaatacatcattatagtaaaactgtataccaatcatcCAAAAGCACCCTAATTTGTTAGTTTAAGTCTAACTTTAGGGTGc                                                              >  1:601162‑M1/48‑91 (MQ=255)
                              aTTTACACCTAGCCCTGGGCACAACGGAATAGTTATAAGTGCACAAACAAAAGCACCCTAATTTGTTAGTTTAAGTCTAACTttagggtgc                                                              <  1:252876‑M2/91‑10 (MQ=255)
                               tttACACCTAGCCCTGGGCACAACGGAATAGTTATAAGTGCACAAACAAAAGCACCCTAATTTGTTAGTTTAAGTCTAACTttagggtgct                                                             >  2:2056594‑M2/1‑81 (MQ=255)
                                 tACACCTAGCCCTGGGCACAACGGAATAGTTATAAGTGCACAAACAAAAGCACCCTAATTTGTTAGTTTAAGTCTAACTttagggtgcttt                                                           <  1:516860‑M2/91‑13 (MQ=255)
                                 tACACCTAGCCCTGGGCACAACGGAATAGTTATAAGCGCACAAACAAAAGCACCCTAATTTGTTAGTTTAAGTCTAACTttagggtgcttt                                                           >  2:2356251‑M2/1‑79 (MQ=255)
                                   cacCTAGCCCTGGGCACAACGGAATAGTTATAAGTGCACAAACAAAAGCACCCTAATTTGTTAGTTTAAGTCTAACTttagggtgcttttg                                                         >  1:4887471‑M2/1‑77 (MQ=255)
                                   cacCTAGCCCTGGGCACAACGGAATAGTTATAAGTGCACAAACAAAAGCACCCTAATTTGTTAGTTTAAGTCTAACTttagggtgcttttg                                                         <  1:4516750‑M2/91‑15 (MQ=255)
                                    acCTAGCCCTGGGCACAACGGAATAGTTATAAGTGCACAAACAAAAGCACCCTAATTTGTTAGTTTAAGTCTAACTTTAGGGTGCTTTTgg                                                        <  2:2276582/91‑1 (MQ=255)
                                     ccTAGCCCTGGGCACAACGGAATAGTTATAAGTGCACAAACAAAAGCACCCTAATTTGTTAGTTTAAGTCTAACTTTAGGGTGCTTTTGga                                                       >  1:5414323/1‑91 (MQ=255)
                                     ccTAGCCCTGGGCACAACGGAATAGTTATAAGTGCACAAACAAAAGCACCCTAATTTGTTAGTTTAAGTCTAACTTTAGGGTGCTTTTGga                                                       >  1:181905/1‑91 (MQ=255)
                                     ccTAGCCCTGGGCACAACGGAATAGTTATAAGTGCACAAACAAAAGCACCCTAATTTGTTAGTTTAAGTCTAACTTTAGGGTGCTTTTGga                                                       <  2:5147132/91‑1 (MQ=255)
                                          cccTGGGCACAACGGAATAGTTATAAGTGCACAAACAAAAGCACCCTAATTTGTTAGTTTAAGTCTAACTTTAGGGTGCTTTTGGATGAtt                                                  <  1:1116441/91‑1 (MQ=255)
                                           ccTGGGCACAACGGAATAGTTATAAGTGCACAAACAAAAGCACCCTAATTTGTTAGTTTAAGTCTAACTTTAGGGTGCTTTTGGATGATTg                                                 >  2:5617925/1‑91 (MQ=255)
                                           ccTGGGCACAACGGAATAGTTATAAGTGCACAAACAAAAGCACCCTAATTTGTTAGTTTAAGTCTAACTTTAGGGTGCTTTTGGATGATTg                                                 <  1:328274/91‑1 (MQ=255)
                                           ccTGGGCACAACGGAATAGTTATAAGTGCACAAACAAAAGCACCCTAATTTGTTAGTTTAAGTCTAACTTTAGGGTGCTTTTGGATGATTg                                                 >  2:6379749/1‑91 (MQ=255)
                                           ccTGGGCACAACGGAATAGTTATAAGTGCACAAACAAAAGCACCCTAATTTGTTAGTTTAAGTCTAACTTTAGGGTGCTTTTGGATGATTg                                                 >  2:1108185/1‑91 (MQ=255)
                                            cTGGGCACAACGGAATAGTTATAAGTGCACAAACAAAAGCACCCTAATTTGTTAGTTTAAGTCTAACTTTAGGGTGCTTTTGGATGATTgg                                                <  1:2237602/91‑1 (MQ=255)
                                               ggCACAACGGAATAGTTATAAGTGCACAAACAAAAGCACCCTAATTTGTTAGTTTAAGTCTAACTTTAGGGTGCTTTTGGATGATTGGtat                                             <  2:1771064/91‑1 (MQ=255)
                                                gCACAACGGAATAGTTATAAGTGCACAAACAAAAGCACCCTAATTTGTTAGTTTAAGTCTAACTTTAGGGTGCTTTTGGATGATTGGtata                                            <  1:6302217/91‑1 (MQ=255)
                                                 cacaACGGAATAGTTATAAGTGCACAAACAAAAGCACCCTAATTTGTTAGTTTAAGTCTAACTTTAGGGTGCTTTTGGATGATTGGTATAc                                           >  2:6725030/1‑91 (MQ=255)
                                                 cacaACGGAATAGTTATAAGTGCACAAACAAAAGCACCCTAATTTGTTAGTTTAAGTCTAACTTTAGGGTGCTTTTGGATGATTGGTATAc                                           <  2:1438395/91‑1 (MQ=255)
                                                  acaACGGAATAGTTATAAGTGCACAAACAAAAGCACCCTAAagttagacttaaactaacaaattagggtgcttttgtttgtgcacttataa                                          <  1:5030718‑M1/91‑51 (MQ=255)
                                                   caACGGAATAGTTATAAGTGCACAAACAAAAGCACCCTAATTTGTTAGTTTAAGTCTAACTTTAGGGTGCTTTTGGATGATTGGTATACAg                                         >  1:3313555/1‑91 (MQ=255)
                                                   caACGGAATAGTTATAAGTGCACAAACAAAAGCACCCTAATTTGTTAGTTTAAGTCTAACTTTAGGGTGCTTTTGGATGATTGGTATACAg                                         >  2:2685851/1‑91 (MQ=255)
                                                    aaCGGAATAGTTATAAGTGCACAAACAAAAGCACCCTAATTTGTTAGTTTAAGTCTAACTTTAGGGTGCTTTTGGATGATTGGTATACAGt                                        <  1:748838/91‑1 (MQ=255)
                                                    aaCGGAATAGTTATAAGTGCACAAACAAAAGCACCCTAATTTGTTAGTTTAAGTCTAACTTTAGGGTGCTTTTGGATGATTGGTATACAGt                                        <  2:1535864/91‑1 (MQ=255)
                                                    aaCGGAATAGTTATAAGTGCACAAACAAAAGCACCCTAATTTGTTAGTTTAAGACTAACTTTAGGGTGCTTTTGGATGATTGGTATACAGt                                        <  1:5350860/91‑1 (MQ=255)
                                                      cGGAATAGTTATAAGTGCACAAACAAAAGCACCCTAATTTGTTAGTTTAAGTCTAACTTTAGGGTGCTTTTGGATGATTGGTATACAGttt                                      <  1:3376907/91‑1 (MQ=255)
                                                         aaTAGTTATAAGTGCACAAACAAAAGCACCCTAATTTGTTAGTTTAAGTCTAACTTTAGGGTGCTTTTGGATGATTGGTATACAGTTTTAc                                   <  2:3624935/91‑1 (MQ=255)
                                                            aGTTATAAGTGCACAAACAAAAGCACCCTAATTTGTTAGTTTAAGTCTAACTTTAGGGTGCTTTTGGATGATTGGTATACAGTTTTACtat                                >  1:6649885/1‑91 (MQ=255)
                                                            aGTTATAAGTGCACAAACAAAAGCACCCTAATTTGTTAGTTTAAGTCTAACTTTAGGGTGCTTTTGGATGATTGGTATACAGTTTTACtat                                >  1:5498611/1‑91 (MQ=255)
                                                             gTTATAAGTGCACAAACAAAAGCACCCTAATTTGTTAGTTTAAGTCTAACTTTAGGGTGCTTTTGGATGATTGGTATACAGTTTTACtata                               >  2:6773013/1‑91 (MQ=255)
                                                              ttATAAGTGCACAAACAAAAGCACCCTAATTTGTTAGTTTAAGTCTAACTTTAGGGTGCTTTTGGATGATTGGTATACAGTTTTACTATaa                              <  1:5159680/91‑1 (MQ=255)
                                                               tataAGTGCACAAACAAAAGCACCCTAATTTGTTAGTTTAAGTCTAACTTTAGGGTGCTTTTGGATGATTGGTATACAGTTTTACTATAat                             <  2:2503927/91‑1 (MQ=255)
                                                               tataAGTGCACAAACAAAAGCACCCTAATTTGTTAGTTTAAGTCTAACTTTAGGGTGCTTTTGGATGATTGGTATACAGTTTTACTATAat                             <  2:5941696/91‑1 (MQ=255)
                                                                  aaGTGCACAAACAAAAGCACCCTAATTTGTTAGTTTAAGTCTAACTTTAGGGTGCTTTTGGATGATTGGTATACAGTTTTACTATAatgat                          <  1:956576/91‑1 (MQ=255)
                                                                  aaGTGCACAAACAAAAGCACCCTAATTTGTTAGTTTAAGTCTAACTTTAGGGTGCTTTTGGATGATTGGTATACAGTTTTACTATAatgat                          >  1:234736/1‑91 (MQ=255)
                                                                   aGTGCACAAACAAAAGCACCCTAATTTGTTAGTTTAAGTCTAACTTTAGGGTGCTTTTGGATGATTGGTATACAGTTTTACTATAatgatg                         >  2:2736937/1‑91 (MQ=255)
                                                                   aGTGCACAAACAAAAGCACCCTAATTTGTTAGTTTAAGTCTAACTTTAGGGTGCTTTTGGATGATTGGTATACAGTTTTACTATAatgatg                         >  1:1726056/1‑91 (MQ=255)
                                                                    gTGCACAAACAAAAGCACCCTAATTTGTTAGTTTAAGTCTAACTTTAGGGTGCTTTTGGATGATTGGTATACAGTTTTACTATAATGATGt                        <  1:2578374/91‑1 (MQ=255)
                                                                       cacaAACAAAAGCACCCTAATTTGTTAGTTTAAGTCTAACTTTAGGGTGCTTTTGGATGATTGGTATACAGTTTTACTATAATGATGTAtt                     <  1:1911658/91‑1 (MQ=255)
                                                                       cacaAACAAAAGCACCCTAATTTGTTAGTTTAAGTCTAACTTTAGGGTGCTTTTGGATGATTGGTATACAGTTTTACTATAATGATGTAtt                     >  2:4629856/1‑91 (MQ=255)
                                                                           aacaaaAGCACCCTTATTTGTTAGTTTAAGTCTAACTTTAGGGTGCTTTTGGATGATTGGTATACAGTTTTACTATAATGATGTATTTTAg                 <  1:6428637/91‑1 (MQ=255)
                                                                             caaaAGCACCCTAAagttagacttaaactaacaaattagggtgcttttgtttgtgcacttataactattccgttgtgcccagggctaggtg               <  1:2821408‑M1/91‑78 (MQ=255)
                                                                              aaaaGCACCCTAAagttagacttaaactaacaaattagggtgcttttgtttgtgcacttataactattccgttgtgcccagggctaggtgt              <  2:5572181‑M1/91‑79 (MQ=255)
                                                                              aaaaGCACCCTAAagttagacttaaactaacaaattagggtgcttttgtttgtgcacttataactattccgttgtgcccagggctaggtgt              >  2:3857299‑M1/1‑13 (MQ=255)
                                                                               aaaGCACCCTAAagttagacttaaactaacaaattagggtgcttttgtttgtgcacttataactattccgttgtgcccagggctaggtgta             >  1:3130791‑M1/1‑12 (MQ=255)
                                                                               aaaGCACCCTAAagttagacttaaactaacaaattagggtgcttttgtttgtgcacttataactattccgttgtgcccagggctaggtgta             <  1:3205867‑M1/91‑80 (MQ=255)
                                                                               aaaGCACCCTAAagttagacttaaactaacaaattagggtgcttttgattgtgcacttataactattccgttgtgcccagggctaggtgta             >  2:4056208‑M1/1‑12 (MQ=255)
                                                                               aaaGCACCCTAAagttagacttaaactaacaaattagggtgattttgtttgtgcacttataactattccgttgtacccagggctaggtgta             >  1:6929897‑M1/1‑12 (MQ=255)
                                                                               AAAGCACCCTAATTTGTTAGTTTAAGTCTAACTTTAGGGTGCTTTTGGATGATTGGTATACAGTTTTACTATAATGATGTATTTTAGCTaa             <  2:1942066‑M1/91‑1 (MQ=255)
                                                                               AAAGCACCCTAATTTGTTAGTTTAAGTCTAACTTTAGGGTGCTTTTGGATGATTGGTATACAGTTTTACTATAATGATGTATTTTAGCTaa             >  1:1394236‑M1/1‑91 (MQ=255)
                                                                                aaGCACCCTAAagttagacttaaactaacaaattagggtgcttttgtttgtgcacttataactattccgttgtgcccagggctaggtgtaa            <  2:6349875‑M1/91‑81 (MQ=255)
                                                                                aaGCACCCTAAagttagacttaaactaacaaattagggtgcttttgtttgtgcacttataactattccgttgtgcccagggctaggtgtaa            <  1:6151061‑M1/91‑81 (MQ=255)
                                                                                AAGCACCCTAATTTGTTAGTTTAAGTCTAACTTTAGGGTGCTTTTGGATGATTGGTATACAGTTTTACTATAATGATGTATTTTAGCTaaa            <  2:807782‑M1/91‑1 (MQ=255)
                                                                                 aGCACCCTAAagttagacttaaactaacaaattagggtgcttttgtttgtgcacttataactattccgttgtgcccagggctaggtgtaaa           >  2:3265698‑M1/1‑10 (MQ=255)
                                                                                  GCACCCTAATTTGTTAGTTTAAGTCTAACTTTAGGGTGCTTTTGGATGATTGGTATACAGTTTTACTATAATGATGTATTTTAGCTAAACg          <  2:4085070‑M1/91‑1 (MQ=255)
                                                                                   cACCCTAAagttagacttaaactaacaaattagggtgcttttgtttgtgcacttataactattccgttgtgcccagggctaggtgtaaatc         <  1:3861358‑M1/91‑84 (MQ=255)
                                                                                   CACCCTAATTTGTTAGTTTAAGTCTAACTTTAGGGTGCTTTTGGATGATTGGTATACAGTTTTACTATAATGATGTATTTTAGCTAAACGa         <  2:2614609‑M1/91‑1 (MQ=255)
                                                                                   CACCCTAATTTGTTAGTTTAAGTCTAACTTTAGGGTGCTTTTGGATGATTGGTATACAGTTTTACTATAATGATGTATTTTAGCTAAACGa         >  2:683498‑M1/1‑91 (MQ=255)
                                                                                    aCCCTAAagttagacttaaactaacaaattagggtgcttttgtttgtgcacttataactattccgttgtgcccagggctaggtgtaaatcc        <  2:6045017‑M1/91‑85 (MQ=255)
                                                                                    aCCCTAAagttagacgtaaactaacaaattagggtgcttttgtttgtgcacttataactattccgttgtgcccagggctaggtgtaaatcc        <  2:5460980‑M1/91‑85 (MQ=255)
                                                                                    ACCCTAATTTGTTAGTTTAAGTCTAACTTTAGGGTGCTTTTGGATGATTGGTATACAGTTTTACTATAATGATGTATTTTAGCTAAACGAc        <  1:6561486‑M1/91‑1 (MQ=255)
                                                                                     cccTAAagttagacttaaactaacaaattagggtgcttttgtttgtgcacttataactattccgttgtgcccagggctaggtgtaaatcct       >  1:3825829‑M1/1‑6 (MQ=255)
                                                                                     CCCTAATTTGTTAGTTTAAGTCTAACTTTAGGGTGCTTTTGGATGATTGGTATACAGTTTTACTATAATGATGTATTTTAGCTAAACGACa       <  2:3021051‑M1/91‑1 (MQ=255)
                                                                                     CCCTAATTTGTTAGTTTAAGTCTAACTTTAGGGTGCTTTTGGATGATTGGTATACAGTTTTACTATAATGATGTATTTTAGCTAAACGACa       <  2:522182‑M1/91‑1 (MQ=255)
                                                                                     CCCTAATTTGTTAGTTTAAGTCTAACTTTAGGGTGCTTTTGGATGATTGGTATACAGTTTTACTATAATGATGTATTTTAGCTAAACGACa       >  1:4347676‑M1/1‑91 (MQ=255)
                                                                                      ccTAAagttagacttaaactaacaaattagggtgcttttgtttgtgcacttataactattccgttgtgcccagggctaggtgtaaatcctc      <  2:335679‑M1/91‑87 (MQ=255)
                                                                                        tAAagttagacttaaactaacaaattagggtgcttttgtttgtgcacttataactattccgttgtgcccagggctaggtgtaaatcctcat    <  1:5272049‑M1/91‑89 (MQ=255)
                                                                                        tAAagttagacttaaactaacaaattagggtgcttttgtttgtgcacttataactattccgttgtgcccagggctaggtgtaaatcctcat    <  2:3480268‑M1/91‑89 (MQ=255)
                                                                                          |                                                                                          
CAGCCATTGAAGCTATTTCTAATGATGAGGATTTACACCTAGCCCTGGGCACAACGGAATAGTTATAAGTGCACAAACAAAAGCACCCTAATTTGTTAGTTTAAGTCTAACTTTAGGGTGCTTTTGGATGATTGGTATACAGTTTTACTATAATGATGTATTTTAGCTAAACGACATTGTT  >  USA300TCH1516_ALE/2850740‑2850920

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