New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913 = 451418382 (1.060)3 (0.040) 3/262 NT 3.6% coding (126/903 nt) fecB ferric citrate ABC transporter periplasmic binding protein
?NC_000913 = 4514201 80 (1.040)coding (108/903 nt) fecB ferric citrate ABC transporter periplasmic binding protein
Rejected: Frequency below/above cutoff threshold.

GACGGACTGCCACGGTTTCAGGTGCGCACGCACTTCGGGCAGGATGCGTTTTGCATCGTTATCGTCGGCAATACCGATCGGGATGACGTCCACGGCGGCCAGCGCATCGGCGAACGAGAGTTCCAGCACCACAATCCGTTGTGGCGTTTTTTCGAGTGTAAACGTGCCGTGTTCGTCCTGAACCGTGGCGGCAAAGGCGTGGCTGATCACCAGCAGCAGGCCTGCAAAAAGAAAACGGATAAATGCCAACATAATCACATTCCAGCTAAAAGCCCGGCAAGCCGGGCGTTAAC  >  NC_000913/4514057‑4514349
                                                                                                                                                |                                                                                                                                                    
gACGGACTGCCACGGTTTCAGGTGCGCACGCACTTCGGGCAGGATGCGTTTTGCATCGTTATCGTCGGCAATACCGATCGGGATGACGTCCACGGCGGCCAGCGCATCGGCGAACGAGAGTTCCAGCACCACAATCCGTTGTGGCGttt                                                                                                                                                  >  2:103371/1‑149 (MQ=255)
 aCGGACTGCCACGGTTTCAGGTGCGCACGCACTTCGGGCAGGATGCGTTTTGCATCGTTATCGTCGGCAATACCGATCGGGATGACGTCCACGGCGGCCAGCGCATCGGCGAACGAGAGTTCCAGCACCACAATCCGTTGTGGCGtttt                                                                                                                                                 >  4:496420/1‑149 (MQ=255)
  cGGACTGCCACGGTTTCAGGTGCGCACGCACTTCGGGCAGGATGCGTTTTGCATCGTTATCGTCGGCAATACCGATCGGGATGACGTCCACGGCGGCCAGCGCATCGGCGAACGAGAGTTCCAGCACCACGATCCGTTGTGGCGttttt                                                                                                                                                <  2:182934/149‑1 (MQ=255)
   ggACTGCCACGGTTTCAGGTGCGCACGCACTTCGGGCAGGATGCGTTTTGCATCGTTATCGTCGGCAATACCGATCGGGATGACGTCCACGGCGGCCAGCGCATCGGCGAACGAGAGTTCCAGCACCACAATCCGTTGTGGCGtttttt                                                                                                                                               >  3:521570/1‑149 (MQ=255)
       tGCCACGGTTTCAGGTGCGCACGCACTTCGGGCAGGATGCGTTTTGCATCGTTATCGTCGGCAATACCGATCGGGATGACGTCCACGGCGGCCAGCGCATCGGCGAACGAGAGTTCCAGCACCACAATCCGTTGTGGCGTTTTTTCGAg                                                                                                                                           >  1:325550/1‑149 (MQ=255)
         ccACGGTTTCAGGTGCGCACGCACTTCGGGCAGGATGCGTTTTGCATCGTTATCGTCGGCAATACCGATCGGGATGACGTCCACGGCGGCCAGCGCATCGGCGAACGAGAGTTCCAGCACCACAATCCGTTGTGGCGTTTTTTCGAgtg                                                                                                                                         >  1:237235/1‑149 (MQ=255)
          cACGGTTTCAGGTGCGCACGCACTTCGGGCAGGATGCGTTTTGCATCGTTATCGTCGGCAATACCGATCGGGATGACGTCCACGGCGGCCAGCGCATCGGCGAACGAGAGTTCCAGCACCACAATCCGTTGTGGCGTTTTTTCGAgtgt                                                                                                                                        >  1:123126/1‑149 (MQ=255)
           aCGGTTTCAGGTGCGCACGCACTTCGGGCAGGATGCGTTTTGCATCGTTATCGTCGGCAATACCGATCGGGATGACGTCCACGGCGGCCAGCGCATCGGCGAACGAGAGTTCCAGCACCACAATCCGTTGTGGCGTTTTTTCGAGTGTa                                                                                                                                       >  4:6761/1‑149 (MQ=255)
              gTTTCAGGTGCGCACGCACTTCGGGCAGGATGCGTTTTGCATCGTTATCGTCGGCAATACCGATCGGGATGACGTCCACGGCGGCCAGCGCATCGGCGAACGAGAGTTCCAGCACCACAATCCGTTGTGGCGTTTTTTCGAGTGTAAAc                                                                                                                                    >  4:366169/1‑149 (MQ=255)
                ttCAGGTGCGCACGCACTTCGGGCAGGATGCGTTTTGCATCGTTATCGTCGGCAATACCGATCGGGATGACGTCCACGGCGGCCAGCGCATCGGCGAACGAGAGTTCCAGCACCACAATCCGTTGTGGCGTTTTTTCGAGTGTAAACGt                                                                                                                                  >  3:422462/1‑149 (MQ=255)
                     gTGCGCACGCACTTCGGGAAGGATGCGTTTTGCATCGTTATCGTCGGCAATACCGATCGGGATGACGTCCACGGCGGCCAGCGCATCGGCGAACGAGAGTTCCAGCACCACAATCCGTTGTGGCGTTTTTTCGAGTGTAAACGTGCCgt                                                                                                                             <  2:701143/149‑1 (MQ=255)
                            cgcacTTCGGGCAGGATGCGTTTTGCATCGTTATCGTCGGCAATACCGATCGGGATGACGTCCACGGCGGCCAGCGCATCGGCGAACGAGAGTTCCAGCACCACAATCCGTTGTGGCGTTTTTTCGAGTGTAAACGTGCCGTGTTCGTc                                                                                                                      <  3:303259/149‑1 (MQ=255)
                            cgcacTTCGGGCAGGATGCGTTTTGCATCGTTATCGTCGGCAATACCGATCGGGATGACGTCCACGGCGGCCAGCGCATCGGCGAACGAGAGTTCCAGCACCACAATCCGTTGTGGCGTTTTTTCGAGTGTAAACGTGCCGTGTTCGTc                                                                                                                      <  2:237235/149‑1 (MQ=255)
                            cgcacTTCGGGCAGGATGCGTTTTGCATCGTTATCGTCGGCAATACCGATCGGGATGACGTCCACGGCGGCCAGCGCATCGGCGAACGAGAGTTCCAGCACCACAATCCGTTGTGGCGTTTTTTCGAGTGTAAACGTGCCGTGTTCGTc                                                                                                                      <  1:369657/149‑1 (MQ=255)
                             gcacTTCGGGCAGGATGCGTTTTGCATCGTTATCGTCGGCAATACCGATCGGGATGACGTCCACGGCGGCCAGCGCATCGGCGAACGAGAGTTCCAGCACCACAATCCGTTGTGGCGTTTTTTCGAGTGTAAACGTGCCGTGTTCGTcc                                                                                                                     >  2:247605/1‑149 (MQ=255)
                               acTTCGGGCAGGATGCGTTTTGCATCGTTATCGTCGGCAATACCGATCGGGATGACGTCCACGGCGGCCAGCGCATCGGCGAACGAGAGTTCCAGCACCACAATCCGTTGTGGCGTTTTTTCGAGTGTAAACGTGCCGTGTTCGTCCTg                                                                                                                   >  1:1756/1‑149 (MQ=255)
                                 ttCGGGCAGGATGCGTTTTGCATCGTTATCGTCGGCAATACCGATCGGGATGACGTCCACGGCGGCCAGCGCATCGGCGAGCGAGAGTTCCAGCACCACAATCCGTTGTGGCGTTTTTTCGAGTGTAAACGTGCCGTGTTCGTCCTGaa                                                                                                                 >  4:65336/1‑149 (MQ=255)
                                     ggCAGGATGCGTTTTGCATCGTTATCGTCGGCAATACCGATCGGGATGACGTCCACGGCGGCCAGCGCATCGGCGAACGAGAGTTCCAGCACCACAATCCGTTGTGGCGTTTTTTCGAGTGTAAACGTGCCGTGTTCGTCCTGAACCGt                                                                                                             >  4:377961/1‑149 (MQ=255)
                                      gCAGGATGCGTTTTGCATCGTTATCGTCGGCAATACCGATCGGGATGACGTCCACGGCGGCCAGCGCATCGGCGAACGAGAGTTCCAGCACCACAATCCGTTGTGGCGTTTTTTCGAGTGTAAACGTGCCGTGTTCGTCCTGAACCGTg                                                                                                            <  4:119338/149‑1 (MQ=255)
                                       cAGGATGCGTTTTGCATCGTTATCGTCGGCAATACCGATCGGGATGACGTCCACGGCGGCCAGCGCATCGGCGAACGAGAGTTCCAGCACCACAATCCGTTGTGGCGTTTTTTCGAGTGTAAACGTGCCGTGTTCGTCCTGAACCGTgg                                                                                                           <  3:480863/149‑1 (MQ=255)
                                         ggATGCGTTTTGCATCGTTATCGTCGGCAATACCGATCGGGATGACGTCCACGGCGGCCAGCGCATCGGCGAACGAGAGTTCCAGCACCACAATCCGTTGTGGCGTTTTTTCGAGTGTAAACGTGCCGTGTTCGTCCTGAACCGTggcg                                                                                                         <  3:163390/149‑1 (MQ=255)
                                         ggATGCGTTTTGCATCGTTATCGTCGGCAATACCGATCGGGATGACGTCCACGGCGGCCAGCGCATCGGCGAACGAGAGTTCCAGCACCACAATCCGTTGTGGCGTTTTTTCGAGTGTAAACGTGCCGTGTTCGTCCTGAACCGTggcg                                                                                                         <  1:168002/149‑1 (MQ=255)
                                          gATGCGTTTTGCATCGTTATCGTCGGCAATACCGATCGGGATGACGTCCACGGCGGCCAGCGCATCGGCGAACGAGAGTTCCAGCACCACAATCCGTTGTGGCGTTTTTTCGAGTGTAAACGTGCCGTGTTCGTCCTGAACCGTggcgg                                                                                                        >  2:118212/1‑149 (MQ=255)
                                          gATGCGTTTTGCATCGTTATCGTCGGCAATACCGATCGGGATGACGTCCACGGCGGCCAGCGCATCGGCGAACGAGAGTTCCAGCACCACAATCCGTTGTGGCGTTTTTTCGAGTGTAAACGTGCCGTGTTCGTCCTGAACCGTggcgg                                                                                                        <  2:21196/149‑1 (MQ=255)
                                              cGTTTTGCATCGTTATCGTCGGCAATACCGATCGGGATGACGTCCACGGCGGCCAGCGCATCGGCGAACGAGAGTTCCAGCACCACAATCCGTTGTGGCGTTTTTTCGAGTGTAAACGTGCCGTGTTCGTCCTGAACCGTGGCGGCaaa                                                                                                    <  4:521570/149‑1 (MQ=255)
                                               gTTTTGCATCGTTATCGTCGGCAATACCGATCGGGATGACGTCCACGGCGGCCAGCGCATCGGCGAACGAGAGTTCCAGCACCACAATCCGTTGTGGCGTTTTTTCGAGTGTAGACGTGCCGTGTTCGTCCTGAACCGTGGCGGCAAAg                                                                                                   >  3:466311/1‑149 (MQ=255)
                                               gTTTTGCATCGTTATCGTCGGCAATACCGATCGGGATGACGTCCACGGCGGCCAGCGCATCGGCGAACGAGAGTTCCAGCACCACAATCCGTTGTGGCGTTTTTTCGAGTGTAAACGTGCCGTGTTCGTCCTGAACCGTGGCGGCAAAg                                                                                                   >  3:498810/1‑149 (MQ=255)
                                               gTTTTGCATCGTTATCGTCGGCAATACCGATCGGGATGACGTCCACGGCGGCCAGCGCATCGGCGAACGAGAGTTCCAGCACCACAATCCGTTGTGGCGTTTTTTCGAGTGTAAACGTGCCGTGTTCGTCCTGAACCGTGGCGGCAAAg                                                                                                   >  3:465935/1‑149 (MQ=255)
                                                ttttGCATCGTTATCGTCGGCAATACCGATCGGGATGACGTCCACGGCGGCCAGCGCATCGGCGAACGAGAGTTCCAGCACCACAATCCGTTGTGGCGTTTTTTCGAGTGTAAACGTGCCGTGTTCGTCCTGAACCGTGGCGGCAAAgg                                                                                                  >  3:521987/1‑149 (MQ=255)
                                                      aTCGTTATCGTCGGCAATACCGATCGGGATGACGTCCACGGCGGCCAGCGCATCGGCGAACGAGAGTTCCAGCACCACAATCCGTTGTGGCGTTTTTTCGAGTGTAAACGTGCCGTGTTCGTCCTGAACCGTGGCGGCAAAGGCGTGGc                                                                                            >  1:594490/1‑149 (MQ=255)
                                                      aTCGTTATCGTCGGCAATACCGATCGGGATGACGTCCACGGCGGCCAGCGCATCGGCGAACGAGAGTTCCAGCACCACAATCCGTTGTGGCGTTTTTTCGAGTGTAAACGTGCCGTGTTCGTCCTGAACCATGGCGGCAAAGGCGTGGc                                                                                            >  1:653177/1‑149 (MQ=255)
                                                       tCGTTATCGTCGGCAATACCGATCGGGATGACGTCCACGGCGGCCAGCGCATCGGCGAACGAGAGTTCCAGCACCACAATCCGTTGTGGCGTTTTTTCGAGTGTAAACGTGCCGTGTTCGTCCTGAACCGTGGCGGCAAAGGCGTGGCt                                                                                           >  3:650630/1‑149 (MQ=255)
                                                          ttATCGTCGGCAATTCCGATCGGGATGACGTCCACGGCGGCCAGCGCATCGGCGAACGAGAGTTCCAGCACCACAATCCGTTGTGGCGTTTTTTCGAGTGTAAACGTGCCGTGTTCGTCCTGAACCGTGGCGGCAAAGGCGTGGCTGAt                                                                                        <  4:249717/149‑1 (MQ=255)
                                                             tcgtcgGCAATACCGATCGGGATGACGTCCACGGCGGCCAGCGCATCGGCGAACGAGAGTTCCAGCACCACAATCCGTTGTGGCGTTTTTTCGAGTGTAAACGTGCCGTGTTCGTCCTGAACCGTGGCGGCAAAGGCGTGGCTGATCAc                                                                                     >  4:256564/1‑149 (MQ=255)
                                                              cgtcgGCAATACCGATCGGGATGACGTCCACGGCGGCCAGCGCATCGGCGAACGAGAGTTCCAGCACCACAATCCGTTGTGGCGTTTTTTCGAGTGTAAACGTGCCGTGTTCGTCCTGAACCGTGGCGGCAAAGGCGTGGCTGATCAcc                                                                                    <  2:487167/149‑1 (MQ=255)
                                                                  ggCAATACCGATCGGGATGACGTCCACGGCGGCCAGCGCATCGGCGAACGAGAGTTCCAGCACCACAATCCGTTGTGGCGTTTTTTCGAGTGTAAACGTGCCGTGTTCGTCCTGAACCGTGGCGGCAAAGGCGTGGCTGATCACcagca                                                                                >  4:500036/1‑149 (MQ=255)
                                                                   gCAATACCGATCGGGATGACGTCCACGGCGGCCAGCGCATCGGCGAAGGAGAGTTCCAGCACCACAATCCGTTGTGGCGTTTTTTCGAGTGTAAACGTGCCGTGTTCGTCCTGTACCGTGGCGGCAAAGGCGTGGCTGATCACcagcag                                                                               >  3:465978/1‑149 (MQ=255)
                                                                   gCAATACCGATCGGGATGACGTCCACGGCGGCCAGCGCATCGGCGAACGAGAGTTCCAGCACCACAATCCGTTGTGGCGTTTTTTCGAGTGTAAACGTGCCGTGTTCGTCCTGAACCGTGGCGGCAAAg                                                                                                   >  4:417458/1‑129 (MQ=255)
                                                                   gCAATACCGATCGGGATGACGTCCACGGCGGCCAGCGCATCGGCGAACGAGAGTTCCAGCACCACAATCCGTTGTGGCGTTTTTTCGAGTGTAAACGTGCCGTGTTCGTCCTGAACCGTGGCGGCAAAg                                                                                                   <  3:417458/129‑1 (MQ=255)
                                                                   gCAATACCGATCGGGATGACGTCCACGGCGGCCAGCGCATCGGCGAACGAGAGTTCCAGCACCACAATCCGTTGTGGCGTTTTTTCGAGTGTAAACGTGCCGTGTTCGTCCTGAACCGTGGCGGCAAAGGCGTGGCTGATCACcagcag                                                                               >  4:81913/1‑149 (MQ=255)
                                                                     aaTACCGATCGGGATGACGTCCACGGCGGCCAGCGCATCGGCGAACGAGAGTTCCAGCACCACAATCCGTTGTGGCGTTTTTTCGAGTGTAAACGTGCCGTGTTCGTCCTGAACCGTGGCGGCAAAGGCGTGGCTGATCACcagcagca                                                                             >  3:332215/1‑149 (MQ=255)
                                                                          cgctCGGGATGACGTCCACGGCGGCCAGCGCATCGGCGAACGAGAGTTCCAGCACCACAATCCGTTGTGGCGTTTTTTCGAGTGTAAACGTGCCGTGTTCGTCCTGAACCGTGGCGGCAAAGGCGTGGCTGATCACCAGCAGCAGGCCt                                                                        <  4:30679/146‑1 (MQ=255)
                                                                          cGATCGGGATGACGTCCACGGCGGCCAGCGCATCGGCGAACGAGAGTTCCAGCACCACAATCCGTTGTGGCGTTTTTTCGAGTGTAAACGTGCCGTGTTCGTCCTGAACCGTGGCGGCAAAGGCGTGGCTGATCACCAGCAGCAGGCCt                                                                        >  1:493941/1‑149 (MQ=255)
                                                                           gATCGGGATGACGTCCACGGCGGCCAGCGCATCGGCGAACGAGAGTTCCAGCACCACAATCCGTTGTGGCGTTTTTTCGAGTGTAAACGTGCCGTGTTCGTCCTGAACCGTGGCGGCAAAGGCGTGGCTGATCACCAGCAGCAGGCCTg                                                                       <  2:623047/149‑1 (MQ=255)
                                                                                 gATGACGTCCACGGCGGCCAGCGCATCGGCGAACGAGAGTTCCAGCACCACAATCCGTTGTGGCGTTTTTTCGAGTGTAAACGTGCCGTGTTCGTCCTGAACCGTGGCGGCAAAGGCGTGGCTGACCACCAGCAGCAGGCCTGCaaaaa                                                                 >  3:82877/1‑149 (MQ=255)
                                                                                 gATGACGTCCACGGCGGCCAGCGCATCGGCGAACGAGAGTTCCACCACCACAATCCGTTGTGGCGTTTTTTCGAGTGTAAACGTGCCGTGTTCGTCCTGAACCGTGGCGGCAAAGGCGTGGCTGATCACCAGCAGCAGGCCTGCaaaaa                                                                 <  3:583006/149‑1 (MQ=255)
                                                                                   tGACGTCCACGGCGGCCAGCGCATCGGCGAACGAGAGTTCCAGCACCACAATCCGTTGTGGCGTTTTTTCGAGTGTAAACGTGCCGTGTTCGTCCTGAACCGTg                                                                                                            <  3:4237/104‑1 (MQ=255)
                                                                                   tGACGTCCACGGCGGCCAGCGCATCGGCGAACGAGAGTTCCAGCACCACAATCCGTTGTGGCGTTTTTTCGAGTGTAAACGTGCCGTGTTCGTCCTGAACCGTg                                                                                                            >  4:4237/1‑104 (MQ=255)
                                                                                   tGACGTCCACGGCGGCCAGCGCATCGGCGAACGAGAGTTCCAGCACCACAATCCGTTGTGGCGTTTTTTCGAGTGTAAACGTGCCGTGTTCGTCCTGAACCGTGGCGGCAAAGGCGTGGCTGATCACCAGCAGCAGGCCTGCAAAAAGa                                                               >  2:228806/1‑149 (MQ=255)
                                                                                      cGTCCACGGCGGCCAGCGCATCGGCGAACGAGAGTTCCAGCACCACAATCCGTTGTGGCGTTTTTTCGAGTGTAAACGTGCCGTGTTCGTCCTGAACCGTGGCGGCAAAGGCGTGGCTGATCACCAGCAGCAGGCCTGCAAAAAGaaaa                                                            <  3:334252/149‑1 (MQ=255)
                                                                                      cGTCCACGGCGGCCAGCGCATCGGCGAACGAGAGTTCCAGCACCACAATCCGTTGTGGCGTTTTTTCGAGTGTAAACGTGCCGTGTTCGTCCTGAACCGTGGCGGCAAAGGCGTGGCTGATCACCAGCAGCAGGCCTGCAAAAAGaaaa                                                            <  4:465935/149‑1 (MQ=255)
                                                                                      cGTCCACGGCGGCCAGCGCATCGGCGAACGAGAGTTCCAGCACCACAATCCGTTGTGGCGTTTTTTCGAGTGTAAACGTGCCGTGTTCGTCCTGAACCGTGGCGGCAAAGGCGTGGCTGATCACCAGCAGCAGGCCTGCAAAAAGaaaa                                                            <  3:81913/149‑1 (MQ=255)
                                                                                       gTCCACGGCGGCCAGCGCATCGGCGAACGAGAGTTCCAGCACCACAATCCGTTGTGGCGTTTTTTCGAGTGTAAACGTGCCGTGTTCGTCCTGAACCGTGGCGGCAAAGGCGTGGCTGATCACCAGCAGCAGGCCTGCAAAAAGAAAAc                                                           >  3:428652/1‑149 (MQ=255)
                                                                                            cggcggCCAGCGCATCGGCGAACGAGAGTTCCAGCACCACAATCCGTTGTGGCGTTTTTTCGAGTGTAAACGTGCCGTGTTCGTCCTGAACCGTGGCGGCAAAGGCGTGGCTGATCACCAGCAGCAGGCCTGCAAAAAGAAAACGGATa                                                      >  4:292295/1‑149 (MQ=255)
                                                                                             ggcggcCAGCGCATCGGCGAACGAGAGTTCCAGCACCACAATCCGTTGTGGCGTTTTTTCGAGTGTAAACGTGCCGTGTTCGTCCTGAACCGTGGCGGCAAAGGCGTGGCTGATCACCAGCAGCAGGCCTGCAAAAAGAAAACGGATaa                                                     >  2:583526/1‑149 (MQ=255)
                                                                                                  ccAGCGCATCGGCGAACGAGAGTTCCAGCACCACAATCCGTTGTGGCGTTTTTTCGAGTGTAAACGTGCCGTGTTCGTCCTGAACCGTGGCGGCAAAGGCGTGGCTGATCACCAGCAGCAGGCCTGCAAAAAGAAAACGGATAAATGcc                                                <  1:103371/149‑1 (MQ=255)
                                                                                                   cAGCGCATCGGCGAACGAGAGTTCCAGCACCACAATCCGTTGTGGCGTTTTTTCGAGTGTAAACGTGCCGTGTTCGTCCTGAACCGTGGCGGCAAAGGCGTGGCTGATCACCAGCAGCAGGCCTGCAAAAAGAAAACGGATAAATGCCa                                               >  4:392510/1‑149 (MQ=255)
                                                                                                         aTCGGCGAACGAGAGTTCCAGCACCACAATCCGTTGTGGCGTTTTTTCGAGTGTAAACGTGCCGTGTTCGTCCTGAACCGTGGCGGCAAAGGCGTGGCTGATCACCAGCAGCAGGCCTGCAAAAAGAAAACGGATAAATGCCAACATaa                                         >  2:565626/1‑149 (MQ=255)
                                                                                                         aTCGGCGAACGAGAGTTCCAGCACCACAATCCGTTGTGGCGTTTTTTCGAGTGTAAACGTGCCGTGTTCGTCCTGAACCGTGGCGGCAAAGGCGTGGCTGATCACCAGCAGCAGGCCTGCAAAAAGAAAACGGATAAATGCCAACATaa                                         >  1:556700/1‑149 (MQ=255)
                                                                                                          tCGGTGAACGAGAGTTCCAGCACCACAATCCGTTGTGGCATTTTTTCGAGTGTAAACGTGCCGTGTTCGTCCTGAACCGTGGCGGCAAACGCGTGGCTGATCACCAGCAGCAGGCCTGCAAAAAGATAACGGATAAATGCCAACATAAt                                        <  4:574015/149‑1 (MQ=255)
                                                                                                          tCGGCGAACGAGAGTTCCAGCACCTCGATCCGTTGTGGCGTTTTTTCGAGTGTAAACGTGCCGTGTTCGTCCTGAACCGTGGCGGCAAAGGCGTGGCTGATCACCAGCAGCAGGCCTGCAAAAAGAAAACGGATAAATGCCAACATAAt                                        <  4:401742/149‑1 (MQ=255)
                                                                                                            ggCGAACGAGAGTTCCAGCACCACAATCCGTTGTGGCGTTTTTTCGAGTGTAAACGTGCCGTGTTCGTCCTGAACCGTGGCGGCAAAGGCGTGGCTGATCACCAGCAGCAGGCCTGCAAAAAGAAAACGGATAAATGCCAACATAATca                                      <  3:292295/149‑1 (MQ=255)
                                                                                                              cGAACGAGAGTTCCAGCACCACAATCCGTTGTGGCGTTTTTTCGAGTGTAAACGTGCCGTGTTCGTCCTGAACCGTGGCGGCAAAGGCGTGGCTGATCACCAGCAGCAGGCCTGCAAAAAGAAAACGGATAAATGCCAACATAATcaca                                    <  1:702625/149‑1 (MQ=255)
                                                                                                              cGAACGAGAGTTCCAGCACCACAATCCGTTGTGGCGTTTTTTCGAGTGTAAACGTGCCGTGTTCGTCCTGAACCGTGGCGGCAAAGGCGTGGCTGATCACCAGCAGCAGGCCTGCAAAAAGAAAACGGATAAATGCCAACATAATcaca                                    <  4:498810/149‑1 (MQ=255)
                                                                                                              cGAACGAGAGTTCCAGCACCACAATCCGTTGTGGCGTTTTTTCGAGTGTAAACGTGCCGTGTTCGTCCTGAACCGTGGCGGCAAAGGCGTGGCTGATCACCAGCAGCAGGCCTGCAAAAAGAAAACGGATAAATGCCAACATAATcaca                                    >  2:468448/1‑149 (MQ=255)
                                                                                                              cGAACGAGAGTTCCAGCACCACAATCCGTTGTGGCGTTTTTTCGAGTGTAAACGTGCCGTGTTCGTCCTGAACCGTGGCGGCAAAGGCGTGGCTGATCACCAGCAGCAGGCCTGCAAAAAGAAAACGGATAAATGCCAACATAATcaca                                    <  3:366169/149‑1 (MQ=255)
                                                                                                               gAACGAGAGTTCCAGCACCACAATCCGTTGTGGCGTTTTTTCGAGTGTAAACGTGCCGTGTTCGTCCTGAACCGTGGCGGCAAAGGCGTGGCTGATCACCAGCAGCAGGCCTGCAAAAAGAAAACGGATAAATGCCAACATAATCACAt                                   >  4:265916/1‑149 (MQ=255)
                                                                                                                 aCGAGAGTTCCAGCACCACAATCCGTTGTGGCGTTTTTTCGAGTGTAAACGTGCCGTGTTCGTCCTGAACCGTGGCGGCAAAGGCGTGGCTGATCACCAGCAGCAGGCCTGCAAAAAGAAAACGGATAAATGCCAACATAATCACATTc                                 >  4:217465/1‑149 (MQ=255)
                                                                                                                  cGAGAGTTCCAGCACCACAATCCGTTGTGGCGTTTTTTCGAGTGTAAACGTGCCGTGTTCGTCCTGAACCGTGGCGGCAAAGGCGTGGCTGATCACCAGCAGCAGGCCTGCAAAAAGAAAACGGATAAATGCCAACATAATCACATTcc                                <  4:680920/149‑1 (MQ=255)
                                                                                                                  cGAGAGTTCCAGCACCACAATCCGTTGTGGCGTTTTTTCGAGTGTAAACGTGCCGTGTTCGTCCTGAACCGTGGCGGCAAAGGCGTGGCTGATCACCAGCAGCAGGCCTGCAAAAAGAAAACGGATAAATGCCAACATAATCACATTcc                                <  1:124624/149‑1 (MQ=255)
                                                                                                                   gagaGTTCCAGCACCACAATCCGTTGTGGCGTTTTTTCGAGTGTAAACGTGCCGTGTTCGTCCTGAACCGTGGCGGCAAAGGCGTGGCTGATCACCAGCAGCAGGCCTGCAAAAAGAAAACGGATAAATGCCAACATAATCACATTCCa                               >  2:287360/1‑149 (MQ=255)
                                                                                                                   gagaGTTCCAGCACCACAATCCGTTGTGGCGTTTTTTCGAGTGTAAACGTGCCGTGTTCGTCCTGAACCGTGGCGGCAAAGGCGTGGCTGATCACCAGCAGCAGGCCTGCAAAAAGAAAACGGATAAATGCCAACATAATCACATTCCa                               >  2:391978/1‑149 (MQ=255)
                                                                                                                     gagTTCCAGCACCACAATCCGTTGTGGCGTTTTTTCGAGTGTAAACGTGCCGTGTTCGTCCTGAACCGTGGCGGCAAAGGCGTGGCTGATCACCAGCAGCAGGCCTGCAAAAAGAAAACGGATAAATGCCAACATAATCACATTCCAGc                             <  2:50142/149‑1 (MQ=255)
                                                                                                                        ttCCAGCACCACAATCCGTTGTGGCGTTTTTTCGATTGTAAACGTGCCGTGTTCGTCCTGAACCGTGGCGGCAAAGGCGTGGCTGATCACCAGCAGCAGGCCTGCAAAAAGAAAACGGATAAATGCCAACATAATCACATTCCAGCTaa                          >  4:630262/1‑149 (MQ=255)
                                                                                                                        ttCCAGCACCACAATCCGTTGTGGCGTTTTTTCGAGTGTAAACGTGCCGTGTTCGTCCTGAACCGTGGCGGCAAAGGCGTGGCTGATCACCAGCAGCAGGCCTGCAAAAAGAAAACGGATAAATGCCAACATAATCACATTCCAGCTaa                          >  4:387581/1‑149 (MQ=255)
                                                                                                                           cAGCACCACAATCCGTTGTGGCgctggaactctcgttcgccgatgcgctggccgccgt                                                                                                                  <  4:102497‑M2/58‑37 (MQ=255)
                                                                                                                           cAGCACCACAATCCGTTGTGGCgctggaactctcgttcgccgatgcgctggccgccgt                                                                                                                  >  3:102497‑M2/1‑22 (MQ=255)
                                                                                                                            aGCACCACAATCCGTTGTGGTgctggaactctcgttcgccgatgcgctggccgccgtggacgtcatcccgatcggtattgccgacgataacgatgcaaaacgcatcctgcccgaagtgcgtgcgcacctgaaaccgtggcagtccgtcg                      <  4:384735‑M2/149‑129 (MQ=255)
                                                                                                                            aGCACCACAATCCGTTGTGGCGTTTTTTCGAGTGTAAACGTGCCGTGTTCGTCCTGAACCGTGGCGGCAAAGGCGTGGCTGATCACCAGCAGCAGGCCTGCAAAAAGAAAACGGATAAATGCCAACATAATCACATTCCAGCTAAAAGc                      >  4:374654/1‑149 (MQ=255)
                                                                                                                                 cacaATCCGTTGTGGCGTTTTTTCGAGTGTAAACGTGCCGTGTTCGTCCTGAACCGTGGCGGCAAAGGCGTGGCTGATCACCAGCAGCAGGCCTGCAAAAAGAAAACGGATAAATGCCAACATAATCACATTCCAGCTAAAAGCCCGGc                 >  3:326727/1‑149 (MQ=255)
                                                                                                                                 cacaATCCGTTGTGGCGTTTTTTCGAGTGTAAACGTGCCGTGTTCGTCCTGAACCGTGGCGGCAAAGGCGTGGCTGATCACCAGCAGCAGGCCTGCAAAAAGAAAACGGATAAATGCCAACATAATCACATTCCAGCTAAAAGCCCGGc                 >  4:380322/1‑149 (MQ=255)
                                                                                                                                 cacaATCCGGTGTTGCGTTTTTGCGAGTGAAAACGTGCGGTGTTCGCCCTGAACCGTGGCGGCAAAGAACTGCCAGATAACCAATAGCAGGCCTGCAAACAGAAAACGGATAAATGCTAACAGAATCACATTCCATGAAAAAGACCGGc                 >  4:274963/1‑149 (MQ=255)
                                                                                                                                  acaATCCGTTGTGGCGTTTTTTCGAGTGTAAACGTGCCGTGTTCGTCCTGAACCGTGGCGGCAAAGGCGTGGCTGATCACCAGCAGCAGGCCTGCAAAAAGAAAACGGATAAATGCCAACATAATCACATTCCAGCTAAAAGCCCGGCa                >  3:396389/1‑149 (MQ=255)
                                                                                                                                     aTCCGTTGTGGCGTTTTTTCGAGTGTAAACGTGCCGTGTTCGTCCTGAACCGTGGCGGCAAAGGCGTGGCTGATCACCAGCAGCAGGCCTGCAAAAAGAAAACGGATAAATGCCAACATAATCACATTCCAGCTAAAAGCCCGGCAAGc             >  1:381232/1‑149 (MQ=255)
                                                                                                                                        cGTTGTGGCGTTTTTTCGAGTGTAAACGTGCCGTGTTCGTCCTGAACCGTGGCGGCAAAGGCGTGGCTGATCACCAGCAGCAGGCCTGCAAAAAGAAAACGGATAAATGCCAACATAATCACATTCCAGCTAAAAGCCCGGCAAGCCgg          <  1:468448/149‑1 (MQ=255)
                                                                                                                                        cGTTGTGGCGTTTTTTCGAGTGTAAACGTGCCGTGTTCGTCCTGAACCGTGGCGGCAAAGGCGTGGCTGATCACCAGCAGCAGGCCTGCAAAAAGAAAACGGATAAATGCCAACATAATCACATTCCAGCTAAAAGCCCGGCAAGCCgg          <  1:565626/149‑1 (MQ=255)
                                                                                                                                             tgGCGTTTTTTCGAGTGTAAACGTGCCGTGTTCGTCCTGAACCGTGGCGGCAAAGGCGTGGCTGATCACCAGCAGCAGGCCTGCAAAAAGAAAACGGATAAATGCCAACATAATCACATTCCAGCTAAAAGCCCGGCAAGCCGGGCGtt     <  2:594490/149‑1 (MQ=255)
                                                                                                                                               gCGTTTTTTCGAGTGTAAACGTGCCGTGTTCGTCCTGAACCGTGGCGGCAAAGGCGTGGCTGATCACCAGCAGCAGGCCTGCAAAAAGAAAACGGATAAATGCCAACATAATCACATTCCAGCTAAAAGCCCGGCAAGCCGGGCGTTaa   >  4:501091/1‑149 (MQ=255)
                                                                                                                                                cGTTTTTTCGAGTGTAAACGTGCCGTGTTCGTCCTGAACCGTGGCGGCAAAGGCGTGGCTGATCACCAGCAGCAGGCCTGCAAAAAGAAAACGGATAAATGCCAACATAATCACATTCCAGCTAAAAGCCCGGCAAGCCGGGCGTTAac  <  3:451649/149‑1 (MQ=255)
                                                                                                                                                |                                                                                                                                                    
GACGGACTGCCACGGTTTCAGGTGCGCACGCACTTCGGGCAGGATGCGTTTTGCATCGTTATCGTCGGCAATACCGATCGGGATGACGTCCACGGCGGCCAGCGCATCGGCGAACGAGAGTTCCAGCACCACAATCCGTTGTGGCGTTTTTTCGAGTGTAAACGTGCCGTGTTCGTCCTGAACCGTGGCGGCAAAGGCGTGGCTGATCACCAGCAGCAGGCCTGCAAAAAGAAAACGGATAAATGCCAACATAATCACATTCCAGCTAAAAGCCCGGCAAGCCGGGCGTTAAC  >  NC_000913/4514057‑4514349

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.