New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913_3_bme_pgi = 120780513 (0.060)167 (0.940) 87/164 0.2 93.3% coding (305/630 nt) stfP e14 prophage; uncharacterized protein
?NC_000913_3_bme_pgi = 1209602 13 (0.070)pseudogene (18/501 nt) stfE pseudogene, e14 prophage; side tail fiber protein fragment family;Phage or Prophage Related

Only 100 of 425 total aligned reads displayed.

CTTCGCGCGCTGCTGTTAAGCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTTATGATTGGTGGTATTGAGCAATCAGTAATCATTCAGTG  >  NC_000913_3_bme_pgi/1207708‑1207904
                                                                                                 |                                                                                                   
catagccaataagcatgttaaatccccctggtaaaagcaggaatgatgccatttgcttgttatttcttcattttcataaacgTTGGTTTGGGAGAAGGa                                                                                                    >  2:1852983‑M1/83‑99 (MQ=255)
catagccaataagcatgttaaatccccctggtaaaagcaggaatgatgccatttgcttgttatttcttcattttcataaacgTTGGTTTGGGAGAAGGAg                                                                                                   <  1:4067232‑M1/18‑1 (MQ=255)
  tagccaataagcatgttaaatccccctggtaaaagcaggaacgatgccatttgcttgttatttcttcattttcataaacgTTGGTTTGGGAGAAGGAGc                                                                                                  <  1:5194449‑M1/19‑1 (MQ=255)
     cgcgcTGCTGTTAAGCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGctctgca                                                                                              <  1:3742161‑M1/100‑8 (MQ=255)
     ccaataagcatgttaaatccccctggtaaaagcaggaatgatgccatttgcttgttatttcttcattttcataaacgTTGGTTTGGGAGAAGGAGCaaaa                                                                                              <  2:3878801‑M1/23‑1 (MQ=255)
        gctgctGTTAAGCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGctctgcatta                                                                                           <  1:3205420‑M1/100‑11 (MQ=255)
           gctgTTAAGCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGctctgcattacct                                                                                        >  2:2652152‑M1/1‑87 (MQ=255)
           agcatgttaaatccccctggtaaaagcaggaatgatgccatttgcttgttatttcttcattttcataaacgTTGGTTTGGGAGAAGGagcaaagct                                                                                            <  2:4763563‑M1/25‑4 (MQ=255)
            ctgTTAAGCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGctctgcattacctg                                                                                       <  2:3506927‑M1/100‑15 (MQ=255)
             tgTTAAGCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAACCGGCTCTCGAAAACCATGGTTAGGGAGCAGGctctgcattacctgt                                                                                      >  2:794048‑M1/1‑85 (MQ=255)
             catgttaaatccccctggtaaaagcaggaatgatgccatttgcttgttatttcttcattttcataaacgTTGGTTTGGGAGAAGGAGCAAAACTCAAt                                                                                        >  2:1427574‑M1/70‑98 (MQ=255)
              atgttaaatccccctggtaaaagcaggaatgatgccatttgcttgttatttcttcattttcataaacgTTGGTTTGGGAGAAGGAGCAAAACTCAATgca                                                                                     >  2:4697110‑M1/69‑100 (MQ=255)
              atgttaaatccccctggtaaaagcaggaatgatgccatttgcttgttatttcttcattttcataaacgTTGGTTTGGGAGAAGGAGCAAAACTCAATg                                                                                       <  1:5210942‑M1/30‑1 (MQ=255)
               ttAAGCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGctctgcattacctgtt                                                                                     <  1:789916‑M1/99‑17 (MQ=255)
                tAAGCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATgcag                                                                                    >  1:2218979/1‑99 (MQ=255)
                gttaaatccccctggtaaaagcaggaatgatgccatttgcttgttatttcttcattttcataaacgTTGGTTTGGGAGAAGGAGCAAAACTCAATgcag                                                                                    <  1:3340555‑M1/33‑1 (MQ=255)
                 ttaaatccccctggtaaaagcaggaatgatgccatttgcttgttatttcttcattttcataaacgTTGGTTTGGGAGAAGGAGCAAAACTCAATgcagc                                                                                   >  1:991654‑M1/66‑99 (MQ=255)
                 aaGCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGctctgcattacctgttgg                                                                                   <  2:325608‑M1/99‑19 (MQ=255)
                  taaatccacctggtaaaagcaggaatgatgccatttgcttgttatttcttcattttcataaacgTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCaa                                                                                 <  2:5224319‑M1/36‑1 (MQ=255)
                   aaatccccctggtaaaagcaggaatgatgccatttgcttgttatttcttcattttcataaacgTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAAc                                                                                <  2:1128340‑M1/37‑1 (MQ=255)
                    ccGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGctctgcattacctgttggtgtc                                                                               <  1:1123057‑M1/100‑23 (MQ=255)
                    ccGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGctctgcattacctgttggtgtc                                                                               <  1:2066365‑M1/100‑23 (MQ=255)
                    aatccccctggtaaaagcaggaatgatgccatttgcttgttatttcttcattttcataaacgTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACg                                                                               >  2:1601098‑M1/63‑100 (MQ=255)
                       cAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGctctgcattacctgttggtgtcccg                                                                            >  2:5054492‑M1/1‑75 (MQ=255)
                       cAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGctctgcattacctgttggtgtcccg                                                                            <  1:1643039‑M1/100‑26 (MQ=255)
                         aGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGctctgcattacctgttggtgtcccgg                                                                           <  1:4543194‑M1/99‑27 (MQ=255)
                          cctggtaaaagcaggaatgatgccatttgcttgttatttcttcattttcataaacgTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACa                                                                         >  2:3207638‑M1/57‑100 (MQ=255)
                          cctggtaaaagcaggaatgatgccatttgcttgttatttcttcattttcataaacgTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACa                                                                         >  1:1338115‑M1/57‑100 (MQ=255)
                           aaTCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGctctgcattacctgttggtgtcccggttc                                                                        >  1:1914985‑M1/1‑71 (MQ=255)
                               cGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGctctgcattacctgttggtgtcccggttccgtg                                                                    <  2:1536052‑M1/100‑34 (MQ=255)
                               cGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGctctgcattacctgttggtgtcccggttccgtg                                                                    <  2:2319539‑M1/100‑34 (MQ=255)
                               cGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGctctgcattacctgttggtgtcccggttccgtg                                                                    >  1:1002507‑M1/1‑67 (MQ=255)
                                gTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGctctgcattacctgttggtgtcccggttccgtgg                                                                   <  2:1717726‑M1/100‑35 (MQ=255)
                                       ggaatgatgccatttgcttgttatttcttcattttcataaacgTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCg                                                            <  1:3236700‑M1/57‑1 (MQ=255)
                                       gATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGctctgcattacctgttggtgtcccggttccgtggccttcag                                                            <  2:2849731‑M1/100‑42 (MQ=255)
                                       gATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGctctgcattacctgttggtgtcccggttccgtggccttcag                                                            <  1:3481454‑M1/100‑42 (MQ=255)
                                        gaatgatgccatttgcttgttatttcttcattttcataaacgTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCgg                                                           >  2:911245‑M1/43‑100 (MQ=255)
                                        gaatgatgccatttgcttgttatttcttcattttcataaacgTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCgg                                                           >  1:2317121‑M1/43‑100 (MQ=255)
                                        agatcaAAGCGGACGGCACTGTGCAAACTGCTCTCGAAAACCTGGGTTTGGGAGAAGGctctgcattacctgttggtgtcccgtttccgtggccttcag                                                            <  2:5187157‑M1/97‑42 (MQ=255)
                                           tgatgccatttgcttgttatttcttcattttcataaacgTTGGTTTGGGAGAAGGAGCAAGACTCAATGCAGCAACGGCTACATTAGGACGCACCGGttt                                                        >  1:4742212‑M1/40‑100 (MQ=255)
                                            cAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGctctgcattacctgttggtgtcccggttccgtggccttcagccact                                                       >  1:1488251‑M1/1‑54 (MQ=255)
                                             atgccatttgcttgttatttcttcattttcataaacgTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTc                                                       >  1:4378109‑M1/38‑99 (MQ=255)
                                              tgccatttgcttgttatttcttcattttcataaacgTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCa                                                      <  1:1551581‑M1/63‑1 (MQ=255)
                                                ccatttgcttgttatttcttcattttcataaacgTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAg                                                   >  2:3990154‑M1/35‑100 (MQ=255)
                                                 catttgcttgttatttcttccttttcataaacgTTGGTTTGGGGGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGc                                                  >  1:3846914‑M1/34‑100 (MQ=255)
                                                  atttgcttgttatttcttcattttcataaacgTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCt                                                 >  1:4928338‑M1/33‑100 (MQ=255)
                                                    aTGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGctctgcattacctgttggtgtcccggttccgtggccttcagccactccgccaac                                               <  1:1436518‑M1/100‑55 (MQ=255)
                                                    aTGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCta                                                >  2:1787356/1‑99 (MQ=255)
                                                     tgcttgttatttcgtcattttcatagacgTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCAt                                                     >  2:2135351‑M1/30‑93 (MQ=255)
                                                     tgcttgttatctcttcattttcataaacgTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCtat                                               >  2:3404782‑M1/30‑99 (MQ=255)
                                                       gCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGctctgcattacctgttggtgtcccggttccgtggccttcagccactccgccaacag                                             <  1:5050607‑M1/99‑57 (MQ=255)
                                                           ttatttcttcattttcataaacgTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGt                                         >  1:4252191‑M1/24‑99 (MQ=255)
                                                             tgCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGctctgcattacctgttggtgtcccggttccgtggccttcagccactccgccaacaggctggct                                      >  1:5264303‑M1/1‑37 (MQ=255)
                                                             tgCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGctctgcattacctgttggtgtcccggttccgtggccttcagccactccgccaacaggctggc                                       <  1:1819239‑M1/99‑63 (MQ=255)
                                                             tgCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGctctgcattacctgttggtgtcccggttccgtggccttcagccactccgccaacaggctggc                                       <  2:3777332‑M1/99‑63 (MQ=255)
                                                             tgCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGctctgcattacctgttggtgtcccggttccgtggccttcagccactccgccaacaggctggc                                       <  2:5275430‑M1/99‑63 (MQ=255)
                                                             tgCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTTa                                       <  2:1994380/99‑1 (MQ=255)
                                                              gCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGctctgcattacctgttggtgtcccggttccgtggccttcagccactccgccaacaggctggct                                      >  1:3682131‑M1/1‑36 (MQ=255)
                                                               cAAACGGCTCTCGAAAACCTTGGTTTGGTATAATGctctgcattacctgttggtttcccggttccgtggccttcagccactccgccaacaggctggct                                      <  1:1723181‑M1/98‑64 (MQ=255)
                                                                    ggCTCTCGAAAACCTTGGTTTGGGAGAAGGctctgcattacctgttggtgtcccggttccgtggccttcagccactccgccaacaggctggctgaaatgc                               <  1:2120481‑M1/100‑71 (MQ=255)
                                                                     gCTCTCGAAAACCTTGGTTTGGGAGAAGGctctgcattacctgttggtgtcccggttccgtggccttcagccactccgccaacaggctggctgacatgc                               >  2:3639506‑M1/1‑29 (MQ=255)
                                                                     attttcataaacgTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTTATGATtggtg                              >  1:2567335‑M1/14‑100 (MQ=255)
                                                                      ctctCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTTATGATtggtgg                             <  2:4656657/100‑1 (MQ=255)
                                                                        ttcataaacgTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTTATGATTggtggt                            >  2:5152844‑M1/11‑99 (MQ=255)
                                                                        ttcataaacgTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTTATGATTGGTGGTa                           <  1:3027062‑M1/90‑1 (MQ=255)
                                                                         tcataaacgTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTTATGATTGGTGGTa                           >  1:4687158‑M1/10‑99 (MQ=255)
                                                                         tcataaacgTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTTATGATTGGTGGTa                           >  2:4278284‑M1/10‑99 (MQ=255)
                                                                         tcGAAAACCTTGGTTTGGGAGAAGGctctgcattacctgttggtgtcccggttccgtggccttcagccactccgccaacaggctggctgaaatgcaatgg                          >  2:2035361‑M1/1‑25 (MQ=255)
                                                                         tcGAAAACCTTGGTTTGGGAGAAGGctctgcattacctgttggtgtcccggttccgtggccttcagccactccgccaacaggctggctgaaatgcaatg                           >  2:3958324‑M1/1‑25 (MQ=255)
                                                                           gAAAACCTTGGTTTGGGAGAAGGctctgcattacctgttggtgtcccggttccgtggccttcagccactccgccaacaggctggctgaaatgcaatggtg                        >  1:891770‑M1/1‑23 (MQ=255)
                                                                           gAAAACCTTGGTTTGGGAGAAGGctctgcattacctgttggtgtcccggttccgtggccttcagccactccgccaacaggctggctgaaatgcaatgg                          <  2:448537‑M1/98‑76 (MQ=255)
                                                                           ataaacgTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTTATGATTGGTGGTAtt                         <  1:4859353‑M1/92‑1 (MQ=255)
                                                                            aaaaCCTTGGTTTGGGAGAAGGctctgcattacctgttggtgtcccggttccgtggccttcagccactccgccaacaggctggctgaaatgcaatggtgc                       <  1:2437122‑M1/100‑79 (MQ=255)
                                                                              agccTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACGTTAGGACGCACCGGTTTCATAGCTATACCGGTTATGATTGGTGGTATTGAGc                     <  2:4858167/98‑1 (MQ=255)
                                                                              aacgTTGGTTTGGGAGAAGGAGCAGAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTTATGATTGGTGGTATTGAg                      >  2:756823‑M1/5‑99 (MQ=255)
                                                                              aacgTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTTATGATTGGTGGTATTGAg                      <  1:2369549‑M1/95‑1 (MQ=255)
                                                                                cgTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTTATGATTGGTGGTATTGAGCaa                   <  1:494965‑M1/98‑1 (MQ=255)
                                                                                ccTTGGTTTGGGAGAAGGctctgcattacctgttggtgtcccggttccgtggccttcagccactccgccaacaggctggctgaaatgcaatggtgccgct                   >  2:3072161‑M1/1‑18 (MQ=255)
                                                                                ccTTGGTTTGGGAGAAGGctctgcattacctgttggtgtcccggttccgtggccttcagccactccgccaacaggctggctgaaatgcaatggtgccgct                   >  1:3250743‑M1/1‑18 (MQ=255)
                                                                                 cTTGGTTTGGGAGAAGGctctgcattacctgttggtgtcccggttccgtggccttcagccactccgccaacaggctggctgaaatgcaatggtgccgctt                  >  1:2837095‑M1/1‑17 (MQ=255)
                                                                                  TTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTTATGATTGGTGGTATTGAGCAATc                 >  1:768779‑M1/1‑100 (MQ=255)
                                                                                   tGGTTTGGGAGAAGGctctgcattacctgttggtgtcccggttccgtggccttcagccactccgccaacaggctggctgaaatgcaatggtgccgctttt                >  1:2616169‑M1/1‑15 (MQ=255)
                                                                                    ggTTTGGGAGAAGGctctgcattacctgttggtgtcccggttccgtggccttcagccactccgccaacaggctggctgaaatgcaatggtgccgcttttt               >  1:3387574‑M1/1‑14 (MQ=255)
                                                                                    GGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTTATGATTGGTGGTATTGAGCAATCAg               <  2:666920‑M1/100‑1 (MQ=255)
                                                                                     gTTTGGGAGAAGGctctgcattacctgttggtgtcccggttccgtggccttcagccactccgccaacaggctggctgaaatgcaatggtgccgctttttc              <  1:272986‑M1/100‑88 (MQ=255)
                                                                                     gTTTGGGAGAAGGctctgcattacctgttagtgtcccggttccgtggccttcagccactccgccaacaggctggctgaaatgcaatggtgccgctttttc              <  1:4104157‑M1/100‑88 (MQ=255)
                                                                                     GTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTTATGATTGGTGGTATTGAGCAATCAGt              >  2:3603742‑M1/1‑100 (MQ=255)
                                                                                      tttGGGAGAAGGctctgcattacctgttggtgtcccggttccgtggccttcagccactccgccaacaggctggctgaaatgcaatggtgccgctttttct             >  2:452763‑M1/1‑12 (MQ=255)
                                                                                      TTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTTATGATTGGTGGTATTGAGCAATCAGt              <  1:362464‑M1/99‑1 (MQ=255)
                                                                                      TTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTTATGATTGGTGGTATTGAGCAATCAGTa             >  1:3193840‑M1/1‑100 (MQ=255)
                                                                                         gggAGAAGGctctgcattacctgttggtgtcccggttccgtggccttcagccactccgccaacaggctggctgaaatgcaatggtgccgctttttctgct          >  1:1861063‑M1/1‑9 (MQ=255)
                                                                                         gggAGAAGGctctgcattacctgttggtgtcccggttccgtggccttcagccactccgccaacaggctggctgaaatgcaatggtgccgctttttctgct          >  1:3978113‑M1/1‑9 (MQ=255)
                                                                                         gggAGAAGGctctgcattacctgttggtgtcccggttccgtggccttcagccactccgccaacaggctggctgaaatgcaatggtgccgctttttctgc           <  1:2926683‑M1/99‑91 (MQ=255)
                                                                                          ggAGAAGGctctgcgttacctgttggtgtcccggttccgtggccttcagccactccgccaacaggctggctgaaatgcaatggtgccgctttttctg            <  2:2489424‑M1/97‑90 (MQ=255)
                                                                                          ggAGAAGGctctgcattacctgttggtgtcccggttccgtggccttcagccactccgccaacaggctggctgaaatgcaatggtgccgctttttctgctg         <  2:2273722‑M1/100‑93 (MQ=255)
                                                                                          ggAGAAGGctctgcattacctgttggtgtcccggttccgtggccttcagccactccgccaacaggctggctgaaatgcaatggtgccgctttttctgct          <  2:1328754‑M1/99‑92 (MQ=255)
                                                                                            agaAGGctctgcattacctgttggtgtcccggttccgtggccttcagccactccgccaacaggctggctgaaatgcaatggtgccgctttttctgctgaa       <  2:1215074‑M1/100‑95 (MQ=255)
                                                                                            agaAGGctctgcattacctgtgtgtgtcccggttccgtggccttcagccactccgccaacaggctggctgaaatgcaatggtgccgctttttctgctgaa       <  1:3096926‑M1/100‑95 (MQ=255)
                                                                                              aaGGctctgcattacctgttggtgtcccggttccgtggccttcagccactccgccaacaggctggctgaaatgcaatggtgccgctttttctgctgaaga     <  1:593170‑M1/100‑97 (MQ=255)
                                                                                                 gctctgcattacctgttggtgtcccggttccgtggccttcagccactccgccaacaggctggctgaaatgcaatggtgccgctttttctgctgaagaata  >  1:1245230‑M1/1‑1 (MQ=255)
                                                                                                 |                                                                                                   
CTTCGCGCGCTGCTGTTAAGCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTTATGATTGGTGGTATTGAGCAATCAGTAATCATTCAGTG  >  NC_000913_3_bme_pgi/1207708‑1207904

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