New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913_3_bme_pgi = 287771453 (1.580)6 (0.190) 4/266 4.1 11.0% intergenic (+96/+855) iap/ygbF aminopeptidase in alkaline phosphatase isozyme conversion/CRISPR adaptation ssRNA endonuclease
?NC_000913_3_bme_pgi = 2877718 47 (1.470)intergenic (+100/+851) iap/ygbF aminopeptidase in alkaline phosphatase isozyme conversion/CRISPR adaptation ssRNA endonuclease
Rejected: Coverage evenness skew score above cutoff.
Rejected: Frequency below/above cutoff threshold.

TGCGGATAATGCTACCTCTGGTGAAGGAGTTGGCGAAGGCGTCTTGATGGGTTTGAAAATGGGAGCTGGGAGTTCTACCGCAGAGGCGGGGGAACTCCAAGTGATATCCATCATCGCATCCAGTGCGCCCGGTTTATCCCCGCTGATGCGGGGAACACCAGCGTCAGGCGTGAAATCTCACCGTCGTTGCCGGTTTATCCCTGCTGGCGCGGGGAACTCTCGGTTCAGGCGTTGCAAACCTGGCTACCGGGCGGTTTATCCCCGCTAACGCGGGGAACTCGTAGTCCAT  >  NC_000913_3_bme_pgi/2877572‑2877860
                                                                                                                                              |                                                                                                                                                  
tGCGGATAATGCTACCTCTGGTGAAGGAGTTGGCGAAGGCGTCTTGATGGGTTTGAAAATGGGAGCTGGGAGTTCTACCGCAGAGGCGGGGGAACTCCAAGTGATATCCATCATCGCATCCAGTGCGCCCGGTTTATCCCCGCtgctgc                                                                                                                                              >  2:343943/1‑149 (MQ=255)
  cGGATAATGCTACCTCTGGTGAAGGAGTTGGCGAAGGCGTCTTGATGGGTTTGAAAATGGGAGCTTGGAGTTCTACCGCAGAGGCGGGGGAACTCCAAGTGATATCCATCATCGCATCCAGTGCGCCCGGTTTATCCCCGCTGATGCgg                                                                                                                                            <  1:587136/149‑1 (MQ=255)
  cGGATAATGCTACCTCTGGTGAAGGAGTTGGCGAAGGCGTCTTGATGGGTTTGAAAATGGGAGCTGGGAGTTCTACCGCAGAGGCGGGGGAACTCCAAGTGATATCCATCATCGCATCCAGTGCGCCCGGTTTATCCCCGCTGATGCgg                                                                                                                                            <  1:263892/149‑1 (MQ=255)
   ggATAATGTTACCTCTTGTGAAGGTGTTGGGGGAGGCGTCTTGATGGGTTTGAAAATGGGAGCTGGGAGTTCTACCGCAGAGGCGGGGGAACTTCAAGTGATATCCATCATCGCATCCAGTGCGCCCGGTTTATCCCCGCTGATGCggg                                                                                                                                           <  2:51611/149‑1 (MQ=255)
      tAATGCTACCTCTGGTGAAGGAGTTGGCGAAGGCGTCTTGATGGGTTTGAAAATGGGAGCTGGGAGTTCTACCGCAGAGGCGGGGGAACTCCAAGTGATATCCATCATCGCATCCAGTGCGCCCGGTTTATCCCCGCTGATGCGGGGaa                                                                                                                                        <  1:161587/149‑1 (MQ=255)
            tACCTCTGGTGAAGGAGTTGGCGAAGGCGTCTTGATGGGTTTGAAAATGGGAGCTGGGAGTTCTACCGCAGAGGCGGGGGAACTCCAAGTGATATCCATCATCGCATCCAGTGCGCCCGGTTTATCCCCGCTGCTGCCGGGAACACCAg                                                                                                                                  >  2:69277/1‑149 (MQ=255)
               ctctGGTGAAGGAGTTGGCGAAGGCGTCTTGGTGGGTTTGAAAATGGGAGCTGGGAGTTCTACCGCAGAGGCGGGGGAACTCCAAGTGATATCCATCATCGCATCCAGTGCGCCCGGTTTATCCCCGCTGATGCgg                                                                                                                                            >  2:343878/1‑136 (MQ=255)
               ctctGGTGAAGGAGTTGGCGAAGGCGTCTTGGTGGGTTTGAAAATGGGAGCTGGGAGTTCTACCGCAGAGGCGGGGGAACTCCAAGTGATATCCATCATCGCATCCAGTGCGCCCGGTTTATCCCCGCTGATGCgg                                                                                                                                            <  1:343878/136‑1 (MQ=255)
                    gTGAAGGAGTTGGCGAAGGCGTCTTGATGGGTTTGAAAATGGGAGCTGGGAGTTCTACCGCAGAGGCGGGGGAACTCCAAGTGATATCCATCATCGCATCCAGTGCGCCCGGTTTATCCCCGCTGATGCGCGGAACACCAGCGTCAGGc                                                                                                                          >  2:347812/1‑149 (MQ=255)
                       aaGGAGTTGGCGAAGGCGTCTTGATGGGTTTGAAAATGGGAGCTGGGAGTTCTACCGCAGAGGCGGGGGAACTCCAAGTGATATCCATCATCGCATCCAGTGCGCCCGGTTTATCCCCGCTGATGCGGGGAACACCAGCGTCGGGGGTg                                                                                                                       >  2:241764/1‑149 (MQ=255)
                         gtaGTTGGGGAAGGACTCTTGGTGGTTTTGCAGATGGGATCTGGGAGTTCGACCACAGAGACGGAGAAACTCCAAATGACATCCATCATCGCATACAGTGCGCCCGGTTTATCCCCGTTGATGCGCGGAACACCAGCGTCAGGCGTGaa                                                                                                                     <  2:424282/147‑1 (MQ=255)
                             ttGGCGAAGGCGTCTTGATGGGTTTGAAAATGGGAGCTGGGAGTTCTACCGCAGAGGCGGGGGAACTCCAAGTGATATCCATCATCGCATCCAGTGCGCCCGGTTTATCCCCGCTGATGCGGGGAACACCAGCGGCAGGCGTGAAAtct                                                                                                                 >  1:278077/1‑149 (MQ=255)
                                      gCGTCTTGATGGGTTTGAAAATGGGAGCTGGGAGTTCTACCGCAGAGGCGGGGGAACTCCAAGTGATATCCATCATCGCATCCAGTGCGCCCGGTTTATCCCCGCTGATGCGGGGAACACCAGCGTCAGGCGTGAAATCTCACcgtcgt                                                                                                        >  2:579757/1‑149 (MQ=255)
                                        gTCTTGATGGGTTTGAAAATGGGAGCTGGGAGTTCTACCGCAGAGGCGGGGGAACTCCAAGTGATATCCATCATCGCATCCAGTGCGCCCGGTTTATCCCCGCTGATGCGGGGAACACCAGCGTCAGGCATGAAATCTCACCGTCGTTg                                                                                                      >  1:502751/1‑149 (MQ=255)
                                          cTTGATGGGTTTGAAAATGGGAGCTGGGAGTTCTACCGCAGAGGCGGGGGAACTCCAAGTGATATCCATCATCGCATCCAGTGCGCCCGGTTTATCCCCGCTGATGCGGGGAACACCAGCGTCAGGCGTGAAATCTCACCGTCGTTGcc                                                                                                    >  1:23795/1‑149 (MQ=255)
                                           ttGATGGGTTTGAAAATGGGAGCTGGGAGTTCTACCGCCGAGGCGGGGGAACTCCAAGTGATATCCATCATCGCATCCAGTGCGCCCGGTTTATCCCCGCTGATGCGGGGAACACCAGCGTCAGGCGTGAAATCTCACCGTCGTTGCCg                                                                                                   >  1:590788/1‑149 (MQ=255)
                                           ttGATGGGTTTGAAAATGGGAGCTGGGAGTTCTACCGCAGAGGCGGGGGAACTCCAAGTGATATCCATCATCGCATCCAGTGCGCCCGGTTTATCCCCGCTGATGCGGGGAACACCAGCGTCAGGCGTGAAATCTCACCGTCGTTGCCg                                                                                                   <  2:23795/149‑1 (MQ=255)
                                           ttGATGGGTTTGAAAATGGCAGCTGGGAGTTCTACCGCCGTCGCGAGGGAACTCCAAGTGATATCCATCATCGCATCCAGTGCGCCCCGTTTATCCCCGCTGATGCGGGGAACCCCAGCGTCCGGCGTGAAACCCACCCGTCGGTGCCg                                                                                                   >  2:306962/1‑149 (MQ=255)
                                                   tttGAAAATGGGAGCTGGGAGTTCTACCGCAGAGGCGGGGGAACTCCAAGTGATATCCATCATCGCATCCAGTGCGCCCGGTTTATCCCCGCTGATGCGGGGAACACCAGCGTCAGGCGTGAAATCTCACCGTCGTTGCCGGTTTATcc                                                                                           >  2:421009/1‑149 (MQ=255)
                                                   tttGAAAATGGGAGCTGGGAGTTCTACCGCAGAGGCGGGGGAACTCCAAGTGATATCCATCATCGCATCCAGTGCGCCCGGTTTATCCCCGCTGATGCGGGGAACACCAGCGTCAGGCGTGAAATCTCACCGTCGTTGCCGGTTTATcc                                                                                           >  1:506519/1‑149 (MQ=255)
                                                       aaaaTGGGAGCTGGGAGTTCTACCGCAGAGGCGGGGGAACTCCAAGTGATATCCATCATCGCATCCAGTGCGCCCGGTTTATCCCCGCTGATGCGGGGAACACCAGCGTCAGGCGTGAAATCTCAcc                                                                                                             >  1:94459/1‑127 (MQ=255)
                                                       aaaaTGGGAGCTGGGAGTTCTACCGCAGAGGCGGGGGAACTCCAAGTGATATCCATCATCGCATCCAGTGCGCCCGGTTTATCCCCGCTGATGCGGGGAACACCAGCGTCAGGCGTGAAATCTCAcc                                                                                                             <  2:94459/127‑1 (MQ=255)
                                                       aaaaTGGGAGCTGGGAGTTCTACCGCAGAGGCGGGGGAACTCCAAGTGATATCCATCATCGCATCCAGTGCGCCCGGTTTATCCCCGCTGATGCGGGGAACACCAGCGTCAGGCGTGAAATCTCACCGTCGTTGCCGGTTTATCCctgc                                                                                       <  1:174856/149‑1 (MQ=255)
                                                           tGGGAGCTGGGAGTTCTACCGCAGAGGCGGGGGAACTCCAAGTGATATCCATCATCGCATCCAGTGCGCCCGGTTTATCCCCGCTGATGc                                                                                                                                              <  1:570048/90‑1 (MQ=255)
                                                           tGGGAGCTGGGAGTTCTACCGCAGAGGCGGGGGAACTCCAAGTGATATCCATCATCGCATCCAGTGCGCCCGGTTTATCCCCGCTGATGc                                                                                                                                              >  2:570048/1‑90 (MQ=255)
                                                               atcTGGTAGTTCTACCGCAGAGGCGTGGTAACTCCAAGTGATATCCATCATCGCATCCAGTGCGCCCGGTTTATCCCCGCTGAGGCGGGGAACACCAGCGTCAGGCGTGAAATCTCACCGTCGTTGCCGGTTTATCCCTGCTGGCGCgg                                                                               <  1:112696/147‑1 (MQ=255)
                                                                gCTGGGAGTTCTACCGCAGAGGCGGGGGAACTCCAAGTGATATCCATCATCGCATCCAGTGCGCCCGGTTTATCCCCGCTGATGCGGGGAACACCAGCGTCAGGCGTGAAATCTCACCGTCGTTGCCGGTTTATCCCTGCTGGcgcgtg                                                                              >  1:543049/1‑147 (MQ=255)
                                                                      aGTTCTACCGCAGAGGCGGGGGAACTCCAGGGGATATCCATCATCGCATCCAGTGCGCCCGGTTTATCCCCGCTGATGCGGGGAACACCAGCGTCAGGCGTGAAATCTCACCGTCGTTGCCGGTTTATCCCTGCTGGCGCGGGGAActc                                                                        >  1:445087/1‑149 (MQ=255)
                                                                        ttCTACGGCAGAGGCGGGGGAGCTCCAAGTGATATCCATCATCGCATCCAGTGCGCGCGGTTTATCCCCGCTGATGCGGGGAACACCAGCGTCAGGCGTGAAATCTCACCGTCGTTGCCGGTTTATCCCTGCTGGCGCGGGGAACtctc                                                                      <  2:147336/149‑1 (MQ=255)
                                                                         tCTACCGCAGAGGCGGGGGAACTCCAAGTGATATCCATCATCGCATCCAGTGCGCCCGGTTTATCCCCGCTGATGCGGGGAACACCAGCGTCAGGCGTGAAATCTCACCGTCGTTGCCGGTTTATCCCTGCTGGCGCGGGGAACTCTCg                                                                     <  1:579757/149‑1 (MQ=255)
                                                                         tCTACCGCAGAGGCGGGGGAACTCCAAGTGATATCCATCATCGCATCCAGTGCGCCCGGTTTATCCCCGCTGATGCGGGGAACACCAGCGTCAGGCGTGAAATCTCACCGTCGTTGCCGGTTTATCCCTGCTGGCGCGGGGAACTCTCg                                                                     <  1:547032/149‑1 (MQ=255)
                                                                               gCAGAGGCGGTGGAACTCCAAGTGATATCCATCATCGCATCCAGTGCGCCCGGTTTATCCCCGCTGATGCGGGGAACACCAGCGGTAGGCGTGAAATCTCTCCGTCGTTGCCGGTTTATCCCTGCTGGCGCGGGGAACTCTCGGTTCAg                                                               <  1:526427/149‑1 (MQ=255)
                                                                                 agagGCGGGGGTACTCCAAGTGATATCCATCATCGCATCCAGTGCGCCCGGTTTATCCCCGCTGATGCGGGGAACACCAGCGTCAGGCGTGAAATCTCACCGTCGTTGCCGGTTTATCCCTGCTGGCGCGGGGAACTCTCGGTTCAGGc                                                             <  2:27217/149‑1 (MQ=255)
                                                                                    ggCGGGGGAACTCCAAGTGATATCCATCATCGCATCCAGTGCGCCCGGTTTATCCCCGCTGATGCGGGGAACACCAGCGTCAGGCGTGAAATCTCACCGTCGTTGCCGGTTTATCCCTGCTGGCGCGGGGAACTCTCGGTTCAGGCGtt                                                          <  1:421009/149‑1 (MQ=255)
                                                                                     gCGGGGGAACTCCAAGTGATATCCATCATCGCATCCAGTGCGCCCGGTTTATCCCCGCTGATGCGGGGAACACCAGCGTCa                                                                                                                             <  2:485390/81‑1 (MQ=255)
                                                                                     gCGGGGGAACTCCAAGTGATATCCATCATCGCATCCAGTGCGCCCGGTTTATCCCCGCTGATGCGGGGAACACCAGCGTCa                                                                                                                             >  1:485390/1‑81 (MQ=255)
                                                                                                          tCCATCATCGCATCCAGTGCGCCCGGTTTATCCCCGCTGATGCGGGGAACACCAGCGTCAGGCGTGAAATCTCACCGTCGTTGCCGGTTTATCCCTGCTGGCGCGGGGAACTCTCGGTTCCGGCCTTGCAAACCTGGCTACCGCGCGGt                                    >  2:239062/1‑149 (MQ=255)
                                                                                                           ccATCATCGCATCCAGTGCGCCCGGTTTATCCCCGCTGATGCGGGGAACACCAGCGTCAGGCGTGAAATCTCACCGTCGTTGCCGGTTTATCCCTGCTGGCGGGCGGAACTCTCGGTTCAGGCGTTGAAAACCTGGCTACCGGGCGGtt                                   >  1:465519/1‑149 (MQ=255)
                                                                                                             atcatcGCATCCAGTGCGCCCGGTTTATCCCCGCTGATGCGGGGAACACCAGCGTCAGGCGTGAAATCTCACCGTCGTTGCCGGTTTATCCCTGCTGGCGCGGGGGACTCTCGGTTCAGGCGTTGCAAACCTGGCTACCGGGCGGTTTa                                 >  2:11090/1‑149 (MQ=255)
                                                                                                                 tcGCATCCAGTGCGCCCGGTTTATCCCCGCTGATGCGGGGAACACCAGCGTCAGGCGTGAAATCTCACCGTCGTTGCCGGTTTATCCCTGCTGGCGCGGGGAACTCTCGGTTCCGGCGTTGCAAACCCGGCGACCGGGCGGGTCATccc                             >  2:183179/1‑149 (MQ=255)
                                                                                                                  cGCATCCAGTGCGCCCGGTTTATCCCCGCGGATGCGGGGAACACCAGCGTCAGGCGTGAAATCTCACCGTCGTTGCCGGTTTATCCCTGCTGGCGCGGGGAACTCTCGGTTCAGGCGTTGCAAACCTGGCTACCGGGCGGTTTATcccc                            <  2:445087/149‑1 (MQ=255)
                                                                                                                     aTCCAGTGCGCCCGGTTTATCCCCGCTGATGCGGGGAACACCAGCGTCAGGCGTGAAATCTCACCGTCGTTGCCGGTTTATCCCTGCTGGCGCGGGGAACTCTCGGTTCAGGCGTTGCAAACCTGGCTACCGGGCGGTTTATCCCCGCt                         >  1:319555/1‑149 (MQ=255)
                                                                                                                         aGTGCGCCCGGTTTATCCCCGCTGATGCGGGGAACACCAGCGTCAGGCGTGAAATCTCACCGTCGTTGCCGGTTTATCCCTGCTGGCGCggg                                                                              >  2:306100/1‑92 (MQ=255)
                                                                                                                         aGTGCGCCCGGTTTATCCCCGCTGATGCGGGGAACACCAGCGTCAGGCGTGAAATCTCACCGTCGTTGCCGGTTTATCCCTGCTGGCGCggg                                                                              <  1:306100/92‑1 (MQ=255)
                                                                                                                            tcgcCCGCTTTATCCCCGCTGATGCGGGGAACACCAGCGTCAGTCGTGAAATCTCACCGTCGTTGCCGGTTTATCCTGGCTGGCGCGGGGAACTCTCGGTTCAGGCGTTGCAAACCTGGCTACCGGGCGGTTTATCCCCGCTAACGCgg                  <  1:316132/148‑1 (MQ=255)
                                                                                                                             gggccGGTTTGTCCCCGCTGGTGCGGGGAACACCAGCGTCAGGGGTGAAATCTTACCGTCTTTGCCGGTTTATTCCTGATGGCGCGGGGAACTCTCGGTTCAGGCGTTGCAAACCTTGCTACCGGGCGGTTTATCCCCGCTAACGCggg                 <  2:294141/146‑1 (MQ=255)
                                                                                                                             ggcCAGGTTTATCCCCGCTGATGCGGGGAACATCATAGTCAGGTGTGAAATTTCACCGTCGTTGCCGGTTTATCCCTGCTGGAGCGGGGAACTCTCGGTTCAGGCGTTGCAAACCTGGCTACCGGGCGGTTTATCCCCGCTAACGCggg                 <  1:39281/148‑1 (MQ=255)
                                                                                                                             cgcCGGGTTCATCCCCGCTGATGCGGGGGACACGGGCGTCAGGCGTGAAATCTCACCGTCGTTGCCGGTTTATCCCTGCTGGCGCGGGGAACTCTCGGTTCAGGCGTTGCAAACCTGGCTTCCGGGCGGTTTATCCCCGCTAACGCggg                 <  2:10322/149‑1 (MQ=255)
                                                                                                                               cccGGTTTATCCCCGCTGATGCGGGGAACACCAGCGTCAGGCGTGAAATCTCACCGTCGTTGCCGGTTTATCCCTGCTGGCGCGGGGAACTCTCGGTTCCGGCGTTGCAAACCTGGCTACCGGGCGGTTTATCCCCGCTAACGCggggc               >  1:105509/1‑148 (MQ=255)
                                                                                                                                ccGGTTTATCCCCGCTGATGCGGGGAACACCAGCGTCAGGCGTGAAATCTCACCGTCGTTGCCGGTTTATCCCTGCTGGCGCGGGGAACTCTCGGTTCAGGCGTTGCAAACCTGGCTACCGGGCGGTTTATCCCCGCTAACGCGGGGaa              >  1:82691/1‑149 (MQ=255)
                                                                                                                                 cGGTTTATCCCCGCTGATGCGGGGAACACCAGCGTCAGGCGTGAAATCTCACCGTCGTTGCCGGTTTATCCCTGCTGGCGCGGGGAACTCTCGGTTCAGGCGTTGCAAACCTGGCTACCGGACGGTTTATCCCCGCTAACGCCGGGAAc             >  1:7406/1‑149 (MQ=255)
                                                                                                                                          cccGCTGATGCGGGGAACACCAGCGTCAGGCGTAAAATCTTACCGTCGTTGACGGTTTATCCCCGCAGGCGCGGGGATCTCTCGGTTCTGTCGTTGCTAAGCTGACTTTCGGGCGGTTTATCCCCCCTAGTCCGGGGAACTctttatcg    >  2:566795/1‑142 (MQ=255)
                                                                                                                                            cgcgGATGCGGGGACCACCAGCGTCAGGCGTGAAAGTTCACCGTCGTTGCCGGTTTATGCCTGCTGGCGCGGGGAACTCTCGGTTCAGGCGTTGCAAACCTGGCTACCGGGCGGTTTATCCCCGCTAACGCGGGGAACTCGTAGTCcat  <  2:311088/149‑1 (MQ=255)
                                                                                                                                              |                                                                                                                                                  
TGCGGATAATGCTACCTCTGGTGAAGGAGTTGGCGAAGGCGTCTTGATGGGTTTGAAAATGGGAGCTGGGAGTTCTACCGCAGAGGCGGGGGAACTCCAAGTGATATCCATCATCGCATCCAGTGCGCCCGGTTTATCCCCGCTGATGCGGGGAACACCAGCGTCAGGCGTGAAATCTCACCGTCGTTGCCGGTTTATCCCTGCTGGCGCGGGGAACTCTCGGTTCAGGCGTTGCAAACCTGGCTACCGGGCGGTTTATCCCCGCTAACGCGGGGAACTCGTAGTCCAT  >  NC_000913_3_bme_pgi/2877572‑2877860

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