Predicted mutation
evidence seq id position mutation freq annotation gene description
JC NC_000913 3,992,314 92 bp x 2 62.1% duplication cyaA → adenylate cyclase

New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913 3992314 =36 (0.480)74 (1.020) 54/248 NT 62.1% coding (1162/2547 nt) cyaA adenylate cyclase
?NC_000913 = 3992405 56 (0.780)coding (1253/2547 nt) cyaA adenylate cyclase

Only 100 of 170 total aligned reads displayed.

GATAACCGCGCTAACTGGAAGATTGATCAGGTGCGTGAGGCGCACAACGAGTTGCTCGACGCGATGATGCAGAGCTACCGTAATCTGATCCGCTTTGCGCGTCGCAATAACCTTAGCGTCTCCGCCAGTCCGCAGGATATCGGCGTGCTGACGCGTAAGCTGTATGCCGCGTTTGAAGCATTACCAGGTAAAGTGACGCTGGTAAACCCGCAGATTTCACCCGATCTCTCGGAACCGAATCTGACCTTTATTTATGTGCCGCCGGGCCGGGCTAACCGTTCAGGTTGGTATCTGTAT  >  NC_000913/3992257‑3992553
                                                                                                                                                    |                                                                                                                                                    
gATAACCGCGCTAACTGGAAGATTGATCAGGTGCGTGAGGCGCACAACGAGTTGCTCGACGCGATGATGCAGAGCTACCGTAATCTGATCCGCTTTGCGCGTCGCAATAACCTTAGCGTCTCCGCCAGTCCGCAGGATATCGGCGTGCt                                                                                                                                                      >  4:355962/1‑149 (MQ=255)
ctttgccgctcgcaataacctatgcgtctgcgccagtccgcagaaaatcggtgtggtGACGCGATGATGCAGAGCTACCGTAATCTGATCCGCTTTGCGCGTCGCAATAACCTTAGCGTCTCCGCCAGTCCGCAGGATATCGGCGTGCt                                                                                                                                                      >  2:5436‑M1/58‑149 (MQ=255)
      cgcgCTAACTGGAAGATTGATCAGGTGCGTGAGGCGCACAACGAGTTGCTCGACGCGATGATGCAGAGCTACCGTAATCTGATCCGCTTTGCGCGTCGCAATAACCTTAGCGTCTCCGCCAGTCCGCAGGATATCGGCGTGCTgacgcg                                                                                                                                                <  1:112998‑M2/149‑7 (MQ=255)
       cgtcgcaataaccttagcgtctccgccagtccgcaggatatcggcgtgctGACGCGATGATGCAGAGCTACCGTAATCTGATCCGCTTTGCGCGTCGCAATAACCTTAGCGTCTCCGCCAGTCCGCAGGATATCGGCGTGCTGACGCGt                                                                                                                                               >  4:214778‑M1/51‑149 (MQ=255)
         gcTAACTGGAAGATTGATCAGGTGCGTGAGGCGCACAACGAGTTGCTCGACGCGATGATGCAGAGCTACCGTAATCTGATCCGCTTTGCGCGTCGCAATAACCTTAGCGTCTCCGCCAGTCCGCAGGATATCGGCGTGCTgacgcgatg                                                                                                                                             <  2:72239‑M2/149‑10 (MQ=255)
         gcTAACTGGAAGATTGATCAGGTGCGTGAGGCGCACAACGAGTTGCTCGACGCGATGATGCAGAGCTACCGTAATCTGATCCGCTTTGCGCGTCGCAATAACCTTAGCGTCTCCGCCAGTCCGCAGGATATCGGCGTGCTGACGCGTaa                                                                                                                                             <  2:55114/149‑1 (MQ=255)
          cTAACTGGAAGATTGATCAGGTGCGTGAGGCGCACAACGAGTTGCTCGACGCGATGATGCAGAGCTACCGTAATCTGATCCGCTTTGCGCGTCGCAATAACCTTAGCGTCTCCGCCAGTCCGCAGGATATCGGCGTGCTgacgcgatga                                                                                                                                            <  3:150056‑M2/149‑11 (MQ=255)
            caataaccttagcgtctccgccagtccgcaggatatcggcgtgctGACGCGATGATGCAGAGCTACCGTAATCTGATCCGCTTTGCGCGTCGCAATAACCTTAGCGTCTCCGCCAGTCCGCAGGATATCGGCGTGCTGACGCGTAAGCt                                                                                                                                          >  4:503817‑M1/46‑149 (MQ=255)
             aataaccttagcgtctccgccagtccgcaggatatcggcgtgctGACGCGATGATGCAGAGCTACCGTAATCTGATCCGCTTTGCGCGTCGCAATAACCTTAGCGTCTCCGCCAGTCCGCAGGATATCGGCGTGCTGACGCGTAAGCTg                                                                                                                                         <  4:134072‑M1/105‑1 (MQ=255)
              ataaccttagcgtctccgccagtccgcaggatatcggcgtgctGACGCGATGATGCAGAGCTACCGTAATCTGATCCGCTTTGCGCGTCGCAATAACCTTAGCGTCTCCGCCAGTCCGCAGGATATCGGCGTGCTGACGCGTAAGCTGt                                                                                                                                        <  3:359779‑M1/106‑1 (MQ=255)
               taaccttagcgtctccgccagtccgcaggatatcggcgtgctGACGCGATGATGCAGAGCTACCGTAATCTGATCCGCTTTGCGCGTCGCAATAACCTTAGCGTCTCCGCCAGTCCGCAGGATATCGGCGTGCTGACGCGTAAGCTGTa                                                                                                                                       <  3:387760‑M1/107‑1 (MQ=255)
                ggAAGATTGATCAGGTGCGTGAGGCGCACAACGACTTGCTCGACGCGATGATGCAGCGCTACCGTAATCTGATCCGCTTTGCGCGTCGCAATAACCTTAGCGTCTGCGCCAGTCCGCAGGATATCGGCGTGCTgacgcgatgatgcaga                                                                                                                                      >  4:426686‑M2/1‑133 (MQ=255)
                 gAAGATTGATCAGGTGCGTGAGGCGCACAACGAGTTGCTCGACGCGATGATGCAGAGCTACCGTAATCTGATCCGCTTTGCGCGTCGCAATAACCTTAGCGTCTCCGCCAGTCCGCAGGATATCGGCGTGCTgacgcgatgatgcagag                                                                                                                                     <  3:621139‑M2/149‑18 (MQ=255)
                        gATCAGGTGCGTGAGGCGCACAACGAGTTGCTCGACGCGATGATGCAGAGCTACCGTAATCTGATCCGCTTTGCGCGTCGCAATAACCTTAGCGTCTCCGCCACTCCGCAGGATATCGGCGTGCTgacgcgatgatgcagatctaccgt                                                                                                                              <  3:594973‑M2/149‑25 (MQ=255)
                          tCAGGTGCGTGAGGCGCACAACGAGTTGCTCGACGCGATGATGCAGAGCTACCGTAATCTGATCCGCTTTGCGCGTCGCAATAACCTTAGCGTCTCCGCCAGTCCGCAGGATATCGGCGTGCTgacgcgatgatgcagagctaccgtaa                                                                                                                            <  2:276383‑M2/149‑27 (MQ=255)
                           cAGGTGCGTGAGGCGCACAACGAGTTGCTCGACGCGATGATGCAGAGCTACCGTAATCTGATCCGCTTTGCGCGTCGCAATAACCTTAGCGTCTCCGCCAGTCCGCAGGATATCGGCGTGCTgacgcgatgatgcagagctaccgtaat                                                                                                                           >  2:450709‑M2/1‑122 (MQ=255)
                             ccgccagtccgcaggatatcggcgtgctGACGCGATGATGCAGAGCTACCGTAATCTGATCCGCTTTGCGCGTCGCAATAACCTTAGCGTCTCCGCCAGTCCGCAGGATATCGGCGTGCTGACGCGTAAGCTGTATGCCGCGTTTGAAg                                                                                                                         >  4:514805‑M1/29‑149 (MQ=255)
                               tGCGTGAGGCGCACAACGAGTTGCTCGACGCGATGATGCAGAGCTACCGTAATCTGATCCGCTTTGCGCGTCGCAATAACCTTAGCGTCTCCGCCAGTCCGCAGGATATCGGCGTGCTgacgcgatgatgcagagctaccgtaatctga                                                                                                                       >  1:483938‑M2/1‑118 (MQ=255)
                               gccagtccgcaggatatcggcgtgctGACGCGATGATGCAGAGCTACCGTAATCTGATCCGCTTTGCGCGTCGCAATAACCTTAGCGTCTCCGCCAGTCCGCAGGATATCGGCGTGCTGACGCGTAAGCTGTAt                                                                                                                                      <  1:191096‑M1/108‑1 (MQ=255)
                                gCGTGAGGCGCACAACGAGTTGCTCGACGCGATGATGCAGAGCTACCGTAATCTGATCCGCTTTGCGCGTCGCAATAACCTTAGCGTCTCCGCCAGTCCGCAGGATATCGGCGTGCTgacgcgatgatgcagagctaccgtaatctgat                                                                                                                      >  1:72239‑M2/1‑117 (MQ=255)
                                    tccgcaggatatcggcgtgctGACGCGATGATGCAGAGCTACCGTAATCTGATCCGCTTTGCGCGTCGCAATAACCTTAGCGTCTCCGCCAGTCCGCAGGATATCGGCGTGCTGACGCGTAAGCTGTATGCCGCGTTTGAAGCATTAcc                                                                                                                  >  3:134072‑M1/22‑149 (MQ=255)
                                    gAGGCGCACAACGAGTTGCTCGACGCGATGATGCAGAGCTACCGTAATCTGATCCGCTTTGCGCGTCGCAATAACCTTAGCGTCTCCGCCAGTCCGCAGGATATCGGCGTGCTgacgcgatgatgcagagctaccgtaatctgatccgc                                                                                                                  <  4:445478‑M2/149‑37 (MQ=255)
                                     ccgcaggatatcggcgtgctGACGCGATGATGCAGAGCTACCGTAATCTGATCCGCTTTGCGCGTCGCAATAACCTTAGCGTCTCCGCCAGTCCGCAGGATATCGGCGTGCTGACGCGTAAGCTGTATGCCGCGTTTGAAGCATTACCa                                                                                                                 >  2:311350‑M1/21‑149 (MQ=255)
                                     aGGCGCACAACGAGTTGCTCGACGCGATGATGCAGAGCTACCGTAATCTGATCCGCTTTGCGCGTCGCAATAACCTTAGCGTCTCCGCCAGTCCGCAGGATATCGGCGTGCTgacgcgatgatgcagagctaccgtaatctgatccgct                                                                                                                 <  4:610849‑M2/149‑38 (MQ=255)
                                      cgcaggatatcggcgtgctGACGCGATGATGCAGAGCTACCGTAATCTGATCCGCTTTGCGCGTCGCAATAACCTTAGCGTCTCCGCCAGTCCGCAGGATATCGGCGTGCTGACGCGTAAGCTGTATGCCGCGTTTGAAGCATTACCAg                                                                                                                >  4:621139‑M1/20‑149 (MQ=255)
                                      cgcaggatatcggcgtgctGACGCGATGATGCAGAGCTACCGTAATCTGATCCGCTTTGCGCGTCGCAATAACCTTAGCGTCTCCGCCAGTCCGCAGGATATCGGCGTGCTGACGCGTAAGCTGTATGCCGCGTTTGAAGCATTACCAg                                                                                                                >  3:544185‑M1/20‑149 (MQ=255)
                                      cgcaggaaatcggcgtcctGACGCGATGATGCAGAGCTACCGTAATCTGATCCGCTTTGCGCGTCGCAATAACCTTAGCGTCTCCGCCAGTCCGCAGGATATCGGCGTGCTGACGCGTAAGCTGTATGCCGCGTTTGAAGCATTACCAg                                                                                                                >  3:233844‑M1/20‑149 (MQ=255)
                                         gcACAACGAGTTGCTCGACGCGATGATGCAGAGCTACCGTAATCTGATCCGCTTTGCGCGTCGCAATAACCTTAGCGTCTCCGCCAGTCCGCAGGATATCGGCGTGCTgacgcgatgatgcagagctaccgtaatctgatccgctttgc                                                                                                             <  4:629976‑M2/149‑42 (MQ=255)
                                            caACGAGTTGCTCGACGCGATGATGCAGAGCTACCGTAATCTGATCCGCTTTGCGCGTCGCAATAACCTTAGCGTCTCCGCCAGTCCGCAGGATATCGGCGTGCTgacgcgatgatgcagagctaccgtaatctgatccgctttgcgcg                                                                                                          <  4:160606‑M2/149‑45 (MQ=255)
                                             aaCGAGTTGCTCGACGCGATGATGCAGAGCTACCGTAATCTGATCCGCTTTGCGCGTCGCAATTACCTTAGCGTCTCCGCCAGTCCGCAGGATATCGGCGTGCTgacgcgatgatgcagagctaccgtaatctgatccgctttgcgcgt                                                                                                         <  4:544185‑M2/149‑46 (MQ=255)
                                                  gcgtgctGACGCGATGATGCAGAGCTACCGTAATCTGATCCGCTTTGCGCGTCGCAATAACCTTAGCGTCTCCGCCAGTCCGCAGGATATCGGCGTGCTGACGCGTAAGCTGTATGCCGCGTTTGAAGCATTACCAGGTAAAGTGACGc                                                                                                    <  4:342517‑M1/142‑1 (MQ=255)
                                                  gTTGCTCGACGCGATGATGCAGAGCTACCGTAATCTGATCCGCTTTGCGCGTCGCAATAACCTTAGCGTCTCCGCCAGTCCGCAGGATATCGGCGTGCTgacgcgatgatgcagagctaccgtaatctgatccgctttgcgcgtcgcaa                                                                                                    <  3:503817‑M2/149‑51 (MQ=255)
                                                  gTTGCTCGACGCGATGATGCAGAGCTACCGTAATCTGATCCGCTTTGCGCGTCGCAATAACCTTAGCGTCTCCGCCAGTCCGCAGGATATCGGCGTGCTgacgcgatgatgcagagctaccgtaatctgatccgctttgcgcgtcgcaa                                                                                                    <  2:324472‑M2/149‑51 (MQ=255)
                                                      cTCGACGCGATGATGCAGAGCTACCGTAATCTGATCCGCTTTGCGCGTCGCAATAACCTTAGCGTCTCCGCCAGTCCGCAGGATATCGGCGTGCTgaagcgatgatgcagagctaccgtaatctgatccgctttgcgcgtcgcaataac                                                                                                <  2:667703‑M2/149‑55 (MQ=255)
                                                       tCGACGCGATGATGCAGAGCTACCGTAATCTGATCCGCTTTGCGCGTCGCAATAACCTTAGCGTCTCCGCCAGTCCGCAGGATATCGGCGTGCTGACGCGTAAGCTGTATGCCGCGTTTGAAGCATTACCAGGTaaa                                                                                                           >  1:568643/1‑137 (MQ=255)
                                                       tCGACGCGATGATGCAGAGCTACCGTAATCTGATCCGCTTTGCGCGTCGCAATAACCTTAGCGTCTCCGCCAGTCCGCAGGATATCGGCGTGCTGACGCGTAAGCTGTATGCCGCGTTTGAAGCATTACCAGGTAAAGTGACGCTGGTa                                                                                               >  2:179754/1‑149 (MQ=255)
                                                         GACGCGATGATGCAGAGCTACCGTAATCTGATCCGCTTTGCGCGTCGCAATAACCTTAGCGTCTCCGCCAGTCCGCAGGATATCGGCGTGCTGACGCGTAAGCTGTATGCCGCGTTTGAAGCATTACCAGGTAAAGTGACGCTGGTaaa                                                                                             <  3:339300‑M1/149‑1 (MQ=255)
                                                          aCGCGATGATGCAGAGCTACCGTAATCTGATCCGCTTTGCGCGTCGCAATAACCTTAGCGTCTCCGCCAGTCCGCAGGATATCGGCGTGCTgacgcgatgatgc                                                                                                                                         >  3:619666‑M2/1‑91 (MQ=255)
                                                          aCGCGATGATGCAGAGCTACCGTAATCTGATCCGCTTTGCGCGTCGCAATAACCTTAGCGTCTCCGCCAGTCCGCAGGATATCGGCGTGCTgacgcgatgatgc                                                                                                                                         >  1:489027‑M2/1‑91 (MQ=255)
                                                          ACGCGATGATGCAGAGCTACCGTAATCTGATCCGCTTTGCGCGTCGCAATAACCTTAGCGTCTCCGCCAGTCCGCAGGATATCGGCGTGCTGACGCGTAAGCTGTATGCCGCGTTTGAAGCATTACCAGGTAAAGTGACGCTGGTAAAc                                                                                            <  1:96391‑M1/149‑1 (MQ=255)
                                                            gcgATGATGCAGAGCTACCGTAATCTGATCCGCTTTGCGCGTCGCAATAACCTTAGCGTCTCCGCCAGTCCGCAGGATATCGGCGTGCTgacgcgatgatgcagagctaccgtaatctgatccgctttgcgcgtcgcaataaccttagc                                                                                          >  1:324472‑M2/1‑89 (MQ=255)
                                                            gcgATGATGCAGAGCTACCGTAATCTGATCCGCTTTGCGCGTCGCAATAACCTTAGCGTCTCCGCCAGTCCGCAGGATATCCGCGTGCTgacgcgatgatgcagagctaccgtaatctgatccgctttgcgcgtcgcaataaccttagc                                                                                          >  4:47959‑M2/1‑89 (MQ=255)
                                                            GCGATGATGCAGAGCTACCGTAATCTGATCCGCTTTGCGCGTCGCAATAACCTTAGCGTCTCCGCCAGTCCGCAGGATATCGGCGTGCTGACGCGTAAGCTGTATGCCGCGTTTGAAGCATTACCAGGTAAAGTGACGCTGGTAAAccc                                                                                          >  2:112998‑M1/1‑149 (MQ=255)
                                                              gatgatGCAGAGCTACCGTAATCTGATCCGCTTTGCGCGTCGCAATAACCTTAGCGTCTCCGCCAGTCCGCAGGATATCGGCGTGCTgacgcgatgatgcagagctaccgtaatctgatccgctttgcgcgtcgcaataaccttagcgt                                                                                        <  3:214778‑M2/149‑63 (MQ=255)
                                                                    gCAGAGCTACCGTAATCTGATCCGCTTTGCGCGTCGCAATAACCTTACCGTCTCCGCCAGTCCGCAGGATATCGGCGTGCTgacgcgatgatgcagagctaccgtaatctgatccgctttgcgcgtcgcaataaccttagcgtctccgc                                                                                  >  3:479652‑M2/1‑81 (MQ=255)
                                                                     cAGAGCTACCGTAATCTGATCCGCTTTGCGCGTCGCAATAACCTTAGCGTCTCCGCCAGTCCGCAGGATATCGGCGTGCTGACGCGTAAGCTGTATGCCGCGTTTGAAg                                                                                                                         <  4:475398/109‑1 (MQ=255)
                                                                     cAGAGCTACCGTAATCTGATCCGCTTTGCGCGTCGCAATAACCTTAGCGTCTCCGCCAGTCCGCAGGATATCGGCGTGCTGACGCGTAAGCTGTATGCCGCGTTTGAAg                                                                                                                         >  3:475398/1‑109 (MQ=255)
                                                                      agagCTACCGTAATCTGATCCGCTTTGCGCGTCGCAATAACCTTAGCGTCTCCGCCAGTCCGCAGGATATCGGCGTGCTgacgcgatgatgcagagctaccgtaatctgatccgctttgcgcgtcgcaataaccttagcgtctccgcca                                                                                >  3:647683‑M2/1‑79 (MQ=255)
                                                                      agagCTACCGTAATCTGATCCGCTTTGCGCGTCGCAATAACCTTAGCGTCTCCGCCAGTCCGCAGGATATCGGCGTGCTGACGCGTAAGCTGTATGCCGCGTTTGAAGCATTACCAGGTAAAGTGACGCTGGTAAACCCGCAGATTTCa                                                                                >  3:650640/1‑149 (MQ=255)
                                                                       gagCTACCGTAATCTGATCCGCTTTGCGCGTCGCAATAACCTTAGCGTCTCCGCCAGTCCGCAGGATATCGGCGTGCTGACGCGTAAGCTGTATGCCGCGTTTGAAGCATTACCAGGTAAAGTGACGCTGGTAAACCCGCAGATTTCAc                                                                               <  2:476967/149‑1 (MQ=255)
                                                                           tACCGTAATCTGATCCGCTTTGCGCGTCGCAATAACCTTAGCGTCTCCGCCAGTCCGCAGGATATCGGCGTGCTGACGCGTAAGCTGTATGCCGCGTTTGAAGCATTACCAGGTAAAGTGACGCTGGTAAACCCGCAGATTTCACCCGa                                                                           >  1:26155/1‑149 (MQ=255)
                                                                             acgTAATCTGATCCGCTTTGCGCGTCGCTATAACCTTAGCGTCACCGCCAGTCCGCAGGATATCGGCGTGCTgacgcgatgatgcagagctaccgtaatctgatccgctttgcgcgtcgcaataaccttagcgtctccgacagtccgca                                                                         >  4:594973‑M2/2‑72 (MQ=255)
                                                                              cGTAATCTGATCCGCTTTGCGCGTCGCAATAACCTTAGCGTCTCCGCCAGTCCGCAGGATATCGGCGTGCTGACGCGTAAGCTGTATGCCGCGTTTGAAGCATTACCAGGTAAAGTGACGCTGGTAAACCCGCAGATTTCACCCGAtct                                                                        >  2:598913/1‑149 (MQ=255)
                                                                              cGTAATCTGATCCGCTTTGCGCGTCGCAATAACCTTAGCGTCTCCGCCAGTCCGCAGGATATCGGCGTGCTGACGCGTAAGCTGTATGCCGCGTTTGAAGCATTACCAGGTAAAGTGACGCTGGTAAACCCGCAGATTTCACCCGAtct                                                                        >  2:424486/1‑149 (MQ=255)
                                                                               gTAATCTGATCCGCTTTGCGCGTCGCAATAACCTTAGCGTCTCCGCCAGTCCGCAGGATATCGGCGTGCTGACGCGTAAGCTGTATGCCGCGTTTGAAGCATTACCAGGTAAAGTGACGCTGGTAAACCCGCAGATTTCACCCGAtctc                                                                       >  2:268137/1‑149 (MQ=255)
                                                                                tAATCTGATCCGCTTTGCGCGTCGCAATAACCTTAGCGTCTCCGCCAGTCCGCAGGATATCGGCGTGCTGACGCGTAAGCTGTATGCCGCGTTTGAAGCATTACCAGGTAAAGTGACGCTGGTAAACCCGCAGATTTCACCCGAtctct                                                                      <  3:342517/149‑1 (MQ=255)
                                                                                   tCTGATCCGCTTTGCGCGTCGCAATAACCTTAGCGTCTCCGCCAGTCCGCAGGATATCGGCGTGCTgacgcgatgatgcagagctaccgtaatctgatccgctttgcgcgtcgcaataaccttagc                                                                                          >  3:79707‑M2/1‑66 (MQ=255)
                                                                                   tCTGATCCGCTTTGCGCGTCGCAATAACCTTAGCGTCTCCGCCAGTCCGCAGGATATCGGCGTGCTGACGCGTAAGCTGTATGCCGCGTTTGAAGCATTACCAGGTAAAGTGACGCTGGTAAACCCGCAGATTTCACCCGATCTCTCgg                                                                   >  2:500301/1‑149 (MQ=255)
                                                                                     tGATCCGCTTTGCGCGTCGCAATAACCTTAGCGTCTCCGCCAGTCCGCAGGATATCGGCGTGCTGACGCGTAAGCTGTATGCCGCGTTTGAAGCATTACCAGGTAAAGTGACGCTGGTAAACCCGCAGATTTCACCCGATCTCTCGGaa                                                                 >  1:54401/1‑149 (MQ=255)
                                                                                     tGATCCGCTTTGCGCGTCGCAATAACCTTAGCGTCTCCGCCAGTCCGCAGGATATCGGCGTGCTGACGCGTAAGCTGTATGCCGCGTTTGAAGCATTACCAGGTAAAGTGACGCTGGTAAACCCGCAGATTTCACCCGATCTCTCGGaa                                                                 >  3:465767/1‑149 (MQ=255)
                                                                                      gATCCGCTTTGCGCGTCGCAATAACCTTAGCGTCTCCGCCAGTCCGCAGGATATCGGCGTGCTGACGCGTAAGCTGTATGCCGCGTTTGAAGCATTACCAGGTAAAGTGACGCTGGTAAACCCGCAGATTTCACCCGATCTCTCGGAAc                                                                <  1:246251/149‑1 (MQ=255)
                                                                                       aTCCGCTTTGCGCGTCGCAATAACCTTAGCGTCTCCGCCAGTCCGCAGGATATCGGCGTGCTgacgcgatgatgcagagctaccgtaatctgatccgctttgcgcgtcgcaataaccttagcgtctccgccagtccgcaggatatcggc                                                               >  3:276692‑M2/1‑62 (MQ=255)
                                                                                        tCCGCTTTGCGCGTCGCAATAACCTTAGCGTCTCCGCCAGTCCGCAGGATATCGGCGTGCTGACGCGTAAGCTGTATGCCGCGTTTGAAGCATTACCAGGTAAAGTGACGCTGGTAAACCCGCAGATTTCACCCGATCTCTCGGAACCg                                                              >  1:646876/1‑149 (MQ=255)
                                                                                           gCTTTGCGCGTCGCAATAACCTTAGCGTCTCCGCCAGTCCGCAGGATATCGGCGTGCTGACGCGTAAGCTGTATGCCGCGTTTGAAGCATTACCAGGTAAAGTGACGCTGGTAAACCCGCAGATTTCACCCGATCTCTCGGAACCGAAt                                                           <  1:397342/149‑1 (MQ=255)
                                                                                           gCTTTGCGCGTCGCAATAACCTTAGCGTCTCCGCCAGTCCGCAGGATATCGGCGTGCTGACGCGTAAGCTGTATGCCGCGTTTGAAGCATTACCAGGTAAAGTGACGCTGGTAAACCCGCAGATTTCACCCGATCTCTCGGAACCGAAt                                                           <  1:281874/149‑1 (MQ=255)
                                                                                            cTTTGCGCGTCGCAATAACCTTAGCGTCTCCGCCAGTCCGCAGGATATCGGCGTGCTGACGCGTAAGCTGTATGCCGCGTTTGAAGCATTACCAGGTAAAGTGACGCTGGTAAACCCGCAGATTTCACCCGATCTCTCGGAACCGAATc                                                          <  4:557826/149‑1 (MQ=255)
                                                                                            cTTTGCCGCTCGCAATAACCTATGCGTCTGCGCCAGTCCGCAGAAAATCGGTGTGGTgacgcgatgatgcagagctaccgtaatctgatccgctttgcgcgtcgcaataaccttagcgtctccgccagtccgcaggatatcggcgtgct                                                          >  2:5436‑M2/1‑57 (MQ=255)
                                                                                                   cgTCGCAATAACCTTAGCGTCTCCGCCAGTCCGCAGGATATCGGCGTGCTgacgcgatgatgcagagctaccgtaatctgatccgctttgcgcgtcgcaataaccttagcgtctccgccagtccgcaggatatcggcgtgctgacgcgt                                                   >  4:214778‑M2/1‑50 (MQ=255)
                                                                                                    gTCGCAATAACCTTAGCGTCTCCGCCAGTCCGCAGGATATCGGCGTGCTGACGCGTAAGCTGTATGCCGCGTTTGAAGCATTACCAGGTAAAGTGACGCTGGTAAACCCGCAGATTTCACCCGATCTCTCGGAACCGAATCTGACCttt                                                  <  3:103644/149‑1 (MQ=255)
                                                                                                      cGCAATAACCTTAGCGTCTCCGCCAGTCCGCAGGATATCGGCGTGCTgacgcgatgatgcagagctaccgtaatctgatccgctttgcgcgtcgcaataaccttagcgtctccgccagtccgcaggatatcggcgtgctgacgcgtaag                                                >  1:128378‑M2/1‑47 (MQ=255)
                                                                                                      cGCAATAACCTTAGCGTCTCCGCCAGTCCGCAGGATATCGGCGTGCTgacgcgatgatgcagagctaccgtaatctgatccgctttgcgcgtcgcaataaccttagcgtctccgccagtccgcaggatatcggcgtgctgacgcgtaag                                                <  2:128378‑M2/149‑103 (MQ=255)
                                                                                                      cGCAATAACCTTAGCGTCTCCGCCAGTCCGCAGGATATCGGCGTGCTgacgcgatgatgcagagctaccgtaatctgatccgctttgcgcgtcgcaataaccttagcgcctccgccagtccgcaggatatcggcgtgctgacgcgtaag                                                <  1:501797‑M2/149‑103 (MQ=255)
                                                                                                      cGCAATAACCTTAGCGTCTCCGCCAGTCCGCAGGATATCGGCGTGCTGACGCGTAAGCTGTATGCCGCGTTTGAAGCATTACCAGGTAAAGTGACGCTGGTAAACCCGCAGATTTCACCCGATCTCTCGGAACCGAATCTGACCtttat                                                <  1:216932/149‑1 (MQ=255)
                                                                                                        cAATAACCTTAGCGTCTCCGCCAGTCCGCAGGATATCGGCGTGCTGACGCGTAAGCTGTATGCCGCGTTTGAAGCATTACCAGGTAAAGTGACGCTGGTAAACCCGCAGATTTCACCCGATCTCTCGGAACCGAATCTGACCtttattt                                              >  1:122419/1‑149 (MQ=255)
                                                                                                         aaTAACCTTAGCGTCTCCGCCAGTCCGCAGGATATCGGCGTGCTgacgcgatgatgcagagctaccgtaatctgatccgctttgcgcgtcgcaataaccttagcgtctccgccagtccgcaggatatcggcgtgctgacgcgtaagctg                                             <  4:134072‑M2/149‑106 (MQ=255)
                                                                                                           tAACCTTAGCGTCTCCGCCAGTCCGCAGGATATCGGCGTGCTgacgcgatgatgcagagctaccgtaatctgatccgctttgcgcgtcgcaataaccttagcgtctccgccagtccgcaggatatcggcgtgctgacgcgtaagctgta                                           <  4:103644‑M2/149‑108 (MQ=255)
                                                                                                           tAACCTTAGCGTCTCCGCCAGTCCGCAGGATATCGGCGTGCTgacgcgatgatgcagagctaccgtaatctgatccgctttgcgcgtcgcaataaccttagcgtctccgccagtccgcaggatatcggcgtgctgacgcgtaagctgta                                           <  3:387760‑M2/149‑108 (MQ=255)
                                                                                                              ccTTAGCGTCTCCGCCAGTCCGCAGGATATCGGCGTGCTgacgcgatgatgcagagctaccgtaatctgatccgctttgcgcgtcgcaataaccttagcgtctccgccagtccgcaggatatcggcgtgctgacgcgtaagctgtatgc                                        <  2:170401‑M2/149‑111 (MQ=255)
                                                                                                               cTTAGCGTCTCCGCCAGTCCGCAGGATATCGGCGTGCTGACGCGTAAGCTGTATGCCGCGTTTGAAGCATTACCAGGTAAAGTGACGCTGGTAAACCCGCAGATTTCACCCGATCTCTCGGAACCGAATCTGACCTTTATTTATGTgcc                                       >  1:474711/1‑149 (MQ=255)
                                                                                                               cTTAGCGTCTCCGCCAGTCCGCAGGATATCGGCGTGCTGACGCGTAAGCTGTATGCCGCGTTTGAAGCATTACCAGGTAAAGTGACGCTGGTAAACCCGCAGATTTCACCCGATCTCTCGGAACCGAATCTGACCTTTATTTATGTgcc                                       >  2:316752/1‑149 (MQ=255)
                                                                                                               cTTAGCGTCTCCGCCAGTCCGCAGGATATCGGCGTGCTGACGCGTAAGCTGTATGCCGCGTTTGAAGCATTACCAGGTAAAGTGACGCTGGTAAACCCGCAGATTTCACCCGATCTCTCGGAACCGAATCTGACCTTTATTTATGTgcc                                       >  2:591503/1‑149 (MQ=255)
                                                                                                                ttAGCGTCTCCGCCAGTCCGCAGGATATCGGCGTGCTGACGCGTAAGCTGTATGCCGCGTTTGAAGCATTACCAGGTAAAGTGACGCTGGTAAACCCGCAGATTTCACCCGATCTCTCGGAACCGAATCTGACCTttatttat                                            <  4:66864/143‑1 (MQ=255)
                                                                                                                ttAGCGTCTCCGCCAGTCCGCAGGATATCGGCGTGCTGACGCGTAAGCTGTATGCCGCGTTTGAAGCATTACCAGGTAAAGTGACGCTGGTAAACCCGCAGATTTCACCCGATCTCTCGGAACCGAATCTGACCTttatttat                                            >  3:66864/1‑143 (MQ=255)
                                                                                                                ttAGCGTCTCCGCCAGTCCGCAGGATATCGGCGTGCTGACGCGTAAGCTGTATGCCGCGTTTGAAGCATTACCAGGTAAAGTGACGCTGGTAAACCCGCAGATTTCACCCGATCTCTCGGAACCGAATCTGACCTTTATTTATGTgccg                                      >  4:248709/1‑149 (MQ=255)
                                                                                                                       ctcCGCCAGTCCGCAGGATATCGGCGTGCTgacgcgatgatgcagagctaccgtaatctgatccgctttgcgcgtcgcaataaccttagcgtctccgccagtccgcaggatatcggcgtgctgacgcgtaagctgtatgccgcgtttga                               >  3:400140‑M2/1‑30 (MQ=255)
                                                                                                                           gCCAGTCCGCAGGATATCGGCGTGCTgacgcgatgatgcagagctaccgtaatctgatccgctttgcgcgtcgcaataaccttagcgtctccgccagtccgcaggatatcggcgtgctgacgcgtaagctgtatgccgcgtttgaagca                           <  4:648940‑M2/149‑124 (MQ=255)
                                                                                                                           gCCAGTCCGCAGGATATCGGCGTGCTgacgcgatgatgcagagctaccgtaatctgatccgctttgcgcgtcgcaataaccttagcgtctccgccagtccgcaggatatcggcgtgctgacgcgtaagctgtat                                          >  2:191096‑M2/1‑26 (MQ=255)
                                                                                                                           gCCAGTCCGCAGGATATCGGCGTGCTgacgcgatgatgcagagctaccgtaatctgatccgctttgcgcgtcgcaataaccttagcgtctccgccagtccgcaggatatcggcgtgctgacgcgtaagctgtat                                          <  1:191096‑M2/134‑109 (MQ=255)
                                                                                                                            ccAGTCCGCAGGATATCGGCGTGCTgacgcgatgatgcagagctaccgtaatctgatccgctttgcgcgtcgcaataaccttagcgtctccgccagtccgcaggatatcggcgtgctgacgcgtaagctgtatgccgcgtttgaagcat                          >  2:334645‑M2/1‑25 (MQ=255)
                                                                                                                              aGTCCGCAGGATATCGGCGTGCTGACGCGTAAGCTGTATGCCGCGTTTGAAGCATTACCAGGTAAAGTGACGCTGGTAAACCCGCAGATTTCACCCGATCTCTCGGAACCGAATCTGACCTTTATTTATGTGCCGCCGGGCCGGGCTaa                        <  4:465767/149‑1 (MQ=255)
                                                                                                                              aGTCCGCAGGATATCGGCGTGCTGACGCGTAAGCTGTATGCCGCGTTTGAAGCATTACCAGGTAAAGTGACGCTGGTAAACCCGCAGATTTCACCCGATCTCTCGGAACCGAATCTGACCTTTATTTATGTGCCGCCGGGCCGGGCTaa                        <  1:268137/149‑1 (MQ=255)
                                                                                                                               gTCCGCAGGATATCGGCGTGCTGACGCGTAAGCTGTATGCCGCGTTTGAAGCATTACCAGGTAAAGTGACGCTGGTAAACCCGCAGATTTCACCCGATCTCTCGGAACCGAATCTGACCTTTATTTATGTGCCGCCGGGCCGGGCTAAc                       <  4:359779/149‑1 (MQ=255)
                                                                                                                                tCCGCAGGATATCGGCGTGCTgacgcgatgatgcagagctaccgtaatctgatccgctttgcgcgtcgcaataaccttagcgtctccgccagtccgcaggatatcggcgtgctgacgcgtaagctgtatgccgcgtttgaagcattacc                      <  4:572613‑M2/149‑129 (MQ=255)
                                                                                                                                  cGCAGGATATCGGCGTGCTgacgcgatgatgcagagctaccgtaatctgatccgctttgcgcgtcgcaataaccttagcgtctccgccagtccgcaggatatcggcgtgctgacgcgtaagctgtatgccgcgtttgaagcattaccag                    >  3:544185‑M2/1‑19 (MQ=255)
                                                                                                                                  cGCAGGATATCGGCGTGCTgacgcgatgatgcagagctaccgtaatctgatccgctttgcgcgtcgcaataaccttagcgtctccgccagtccgcaggatatcggcgtgctgacgcgtaagctgtatgccgcgtttgaagcattaccag                    >  4:387760‑M2/1‑19 (MQ=255)
                                                                                                                                      ggATATCGGCGTGCTGACGCGTAAGCTGTATGCCGCGTTTGAAGCATTACCAGGTAAAGTGACGCTGGTAAACCCGCAGATTTCACCCGATCTCTCGGAACCGAATCTGACCTTTATTTATGTGCCGCCGGGCCGGGCTAACCGTTCAg                <  1:424486/149‑1 (MQ=255)
                                                                                                                                          atCGGCGTCTgacgcgatgatgcagagctaccgtaatctgatccgctttgcgcgtcgcaataaccttagcgtctccgccagtccgcaggatatcggcgtgctgacgcgtaagctgtatgccgcgtttgaagcattaccaggtaaagtga           >  1:333480‑M2/1‑10 (MQ=255)
                                                                                                                                              gCGTGCTgacgcgatgatgcagagctaccgtaatctgatccgctttgcgcgtcgcaataaccttagcgtctccgccagtccgcaggatatcggcgtgc                                                           >  1:23250‑M2/1‑7 (MQ=255)
                                                                                                                                              gCGTGCTgacgcgatgatgcagagctaccgtaatctgatccgctttgcgcgtcgcaataaccttagcgtctccgccagtccgcaggatatcggcgtgc                                                           <  2:23250‑M2/98‑92 (MQ=255)
                                                                                                                                                    tGACGCGTAAGCTGTATGCCGCGTTTGAAGCATTACCAGGTAAAGTGACGCTGGTAAACCCGCAGATTTCACCCGATCTCTCGGAACCGAATCTGACCTTTATTTATGTGCCGCCGGGCCGGGCTAACCGTTCAGGTTGGTATCTGtat  >  4:625503/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
GATAACCGCGCTAACTGGAAGATTGATCAGGTGCGTGAGGCGCACAACGAGTTGCTCGACGCGATGATGCAGAGCTACCGTAATCTGATCCGCTTTGCGCGTCGCAATAACCTTAGCGTCTCCGCCAGTCCGCAGGATATCGGCGTGCTGACGCGTAAGCTGTATGCCGCGTTTGAAGCATTACCAGGTAAAGTGACGCTGGTAAACCCGCAGATTTCACCCGATCTCTCGGAACCGAATCTGACCTTTATTTATGTGCCGCCGGGCCGGGCTAACCGTTCAGGTTGGTATCTGTAT  >  NC_000913/3992257‑3992553

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