New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913_3_bme_pgi = 394668848 (1.060)7 (0.160) 7/278 4.9 13.2% intergenic (+81/‑12) rrlC/rrfC 23S ribosomal RNA of rrnC operon/5S ribosomal RNA of rrnC operon
?NC_000913_3_bme_pgi = 4171601 NA (NA)intergenic (+57/‑36) rrlB/rrfB 23S ribosomal RNA of rrnB operon/5S ribosomal RNA of rrnB operon
Rejected: Coverage evenness skew score above cutoff.
Rejected: Frequency below/above cutoff threshold.

AGGCCGGGTGTGTAAGCGCAGCGATGCGTTGAGCTAACCGGTACTAATGAACCGTGAGGCTTAACCTTACAACGCCGAAGATGTTTTGGCGGATGAGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGATGAAACAGAATTTGCCTGGCGGCCGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAG‑‑TGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATTAAGCAGTAAG  >  NC_000913_3_bme_pgi/3946540‑3946835
                                                                                                                                                    |                                                                                                                                                     
aGGCCGGGTGTGTAAGCGCAGCGATGCGTTGAGCTAACCGGTACTAATGAACCGTGAGGCTTAACCTTACAACGCCGAATATGTTTTGGCGGATGAGAGAAGATTTTCAGCCTGATACAGATTCAATCAGAACGCAGAAGCGGTCtgat                                                                                                                                                       <  1:609836/149‑1 (MQ=17)
aGGCCGGGTGTGCAAGCGCAGCGATGCGTTTAGCTAACCGGTACTAATGAACCGTGAGGCTTAACCTTACAACGCCGAAGATGTTTTGGCGGATGAGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCtgat                                                                                                                                                       <  2:445848/149‑1 (MQ=17)
    cGGGTGTGTAAGCGCAGCGATGCGTTGAGCTAACCGGTACTAATGAACCGTGAGGCTTAACCTTACAACGCCGAAGATGTTTTGGCGGATGAGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGATGaaa                                                                                                                                                   >  2:797802/1‑149 (MQ=34)
      ggTGTGTAAGCGCAGCGATGCGTTGAGCTAACCGGTACTAATGAACCGTGAGGCTTAACCTTACAACGCCGAAGATGTTTTGGCGGATGAGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGATGAAACa                                                                                                                                                 <  2:507472/149‑1 (MQ=35)
            tAAGCGCAGCGATGCGTTGAGCTAACCGGTACTAATGAACCGTGAGGCTTAACCTTACAACGCCGAAGATGTTTTGGCGGATGAGAGAAGATTTTCAGCCTGGTACAGATTAAATCAGAACGCAGAAGCGGTCTGATGAAACAGAAttt                                                                                                                                           <  1:406304/149‑1 (MQ=21)
                 gcAGCGATGCGTTGAGCTAACCGGTACTAATGAACCGTGAGGCTTAACCTTACAACGCCGAAGATGTTTTGGCGGATGAGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGATGAAACAGAATTTGCCTg                                                                                                                                      >  1:451434/1‑149 (MQ=35)
                   aGCGATGCGGTGAGCTAACCGGTACTAATGAACCGTGAGGCTTAAGCTGACAACGCCGAAGATGTTTTGGCGGATGAGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAACGCACACGCGGCCTGCTGAAAAAGAATTTACATggc                                                                                                                                    >  2:138259/1‑149 (MQ=11)
                         gCGTTGAGCTAACCGGTACTAATGAACCGTGAGGCTTAACCTTACAACGCCGAAGATGTTTTGGCGGATGAGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGATGAAACAGAATTTGCCTGGCGGCCGt                                                                                                                              <  1:447151/149‑1 (MQ=35)
                               aGCTAACCGGTACTAATGAACCGTGAGGCTTAACCTTACAACGCCGAAGATGTTTTGGCGGATGAGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGATGAAACAGAATTTGCCTGGCGGCCGTAGcgcg                                                                                                                        >  1:600177/1‑149 (MQ=35)
                                      cGGTACTAATGAACCGTGAGGCTTAACCTTACAACGCCGAAGATGTTTTGGCGGATGAGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGATGAAACAGAATTTGCCTGGCGGCCGTAGCGCGGTGGTcc                                                                                                                 >  1:66206/1‑149 (MQ=35)
                                      cGGTACTAATGAACCGTGAGGCTTAACCTTACAACGCCGAAGATGTTTTGGCGGATGAGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGATGAAACAGAATTTGCCTGGCGGCCGTAGCGCGGTGGTcc                                                                                                                 <  1:334185/149‑1 (MQ=35)
                                       ggTACTAATGAACCGTGAGGCTTAACCTTACAACGCCGAAGGTGTTTTGGCGGATGAGAGAAGATTTTCAGCCTGCTACAGATTAAATCAGAACGCAGAACCGGCCCGAGGAAACAGAATTTGCCTGGCGGCCGAAGCGCGGTGGTccc                                                                                                                >  1:147357/1‑149 (MQ=11)
                                        gTACTAATGAACCGTGAGGCTTAACCTTACAACGCCGAAGATGTTTTGGCGGATGAGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGATGAAACAGAATTTGCCTGGCGGCCGTAGCGCGGTGGTCCCa                                                                                                               >  2:260674/1‑149 (MQ=35)
                                         tACTAATAAACCGTGAGGCTTAACCTTACAACGCCGAAGATGTTTTGGCGGATGAGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGATGAAACAGAATTTGCCTGGCGGCCGTAGAGCGGCGGTcacag                                                                                                              <  2:559403/149‑2 (MQ=255)
                                          aCTAATGAACCGTGAGGCTTAACCTTACAACGCCGAAGATGTTTTGGCGGATGAGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGATGAAACAGAATTTGCCTGGCGGCCGTAGCGCGGTGGTCCCAcc                                                                                                             <  2:66206/149‑1 (MQ=35)
                                           cTAATGAACCGTGAGGCTTAACCTTACAACGCCGAAGATGTTTTGGCGGATGAGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGATGAAACAGAATTTGCCTGGCGGCCGTAGCGCGGTGGTCCCACCt                                                                                                            <  1:558823/149‑1 (MQ=35)
                                             aaTGAACCGTGAGGCTTAACCTTACAACGCCGAAGATGTTCTGGCGGATGAGAGAAGATTTTCAGGCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGATGAAACAGAAGTTGCCTGGCGGCCGTAGCGCGGTGGTCCCACCTGa                                                                                                          <  2:410396/149‑1 (MQ=255)
                                              aTGAACCGTGAGGCTTAACCTTACAACGCCGAAGATGTTTTGGCGGATGAGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGATGAAAAAGAATTTGCCTGGCGGCCGTAGCGCGGTGGTCCCACCTGAc                                                                                                         >  2:567890/1‑149 (MQ=21)
                                                   ccGTGAGGCTTAACCTTACAACGCCGAAGATGTTTTGGCGGATGAGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGATGAAACAGAATTTGCCTGGCGGCCGTAGCGCGGTGGTCCCACCTGACCACAt                                                                                                    >  1:164480/1‑149 (MQ=25)
                                                      tGAGGCTTAACCTTACAACGCCGAAGATGTTTTGGCGGATGAGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGATGAAACAGAATTTGCCTGGCGGCCGTAGCGCGGTGGTCCCACCTGACCCCATGcc                                                                                                 >  1:608497/1‑149 (MQ=35)
                                                         ggCTTAACCTTACAACGCCGAAGATGTTTTGGCGGATGAGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTGTGATGAAACAGAATTTGCCTGGCGGCCGTAGCGCGGTGGTCCCACCTGACCCCATGCCGaa                                                                                              <  1:626580/149‑1 (MQ=21)
                                                         ggCTTAACCTTACAACGCCGAAGATGTTTTGGCGGATGAGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGATGAAACAGAATTTGCCTGGCGGCCGTAGCGCGGTGGTCCCACCTGACCCCATGCCGaa                                                                                              <  1:578167/149‑1 (MQ=35)
                                                           cTTACCCTCACAACCTCTAATATTTTTTGGCGGATGAGAGAAGCTTTTCAGCCTGATACAGATGAAATCAGAACGCAGAAGCGCTCTGATGAAACAGAATTCGCCTGGCGGCCGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACt                                                                                            <  1:617618/149‑1 (MQ=11)
                                                           cTTAACCTTACAACGCCGAAGATGTTTTGGCGGATGAGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGATGAAACAGACTTTGCCTGGCGGCCGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACt                                                                                            >  2:498101/1‑149 (MQ=25)
                                                            ttAACCTTACAACGCCGAAGATGTTTTGGCGGATGAGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGATGAAACAGAATTTGCCTGGCGGCCGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTc                                                                                           >  2:61285/1‑149 (MQ=35)
                                                                      aaCGCCGAAGATGTTTTGGCGGATGAGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGATGAAACAGAATTTGCCTGGCGGCCGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAG‑‑TGaaa                                                                               >  1:204731/1‑149 (MQ=35)
                                                                      aaCGCCGAAGATGTTTTGGCGGATGAGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGATGAAACAGAATTTGCCTGGCGGCCGTAGCGCGGTGGGCCCACCTGAACCCATGCCGAACTCAGGAG‑‑TGaaa                                                                               >  2:445156/1‑149 (MQ=255)
                                                                        cGCCGAAGATGTTTTGGCGGATGAGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGATGAAACAGAATTTGCCTGGCGGCCGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAG‑‑TGAAACg                                                                             >  1:689810/1‑149 (MQ=35)
                                                                          ccGAACTTGTTTTGGCGGATGAGAGAATATCTGCAGCCTGATACAGATTACATCAAAACGCAGAAGCGGTCTGATGAAACATCATCTGCTTTGTGCCTGTATAAAGCGGCTCTCACCTGACTCCTTGCCGAACTCAGAAG‑‑TGAAACGcc                                                                           <  2:369070/149‑1 (MQ=2)
                                                                            gAAGATGTTTTGGCGGATGAGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGATGAAACAGAATTTGCCTGGCGGCCGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAG‑‑TGAAACGCCGt                                                                         >  1:486224/1‑149 (MQ=35)
                                                                                  gTTTTGGCGGATAGGAGAAGGTTTTCAGCCTGATACAGTTTAAATTAGAACGCAGTAGCGGTCTGATGAAACAGAATTTGCCTGGCGGCCGTAGCGCGGTGGTCCCACCTGTCCCCATGCCGATCTCTGAAG‑‑TGAACCGCCGTAGCGcc                                                                   <  2:419047/149‑1 (MQ=11)
                                                                                   ttttGGCGGATGAGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGATGAAACAGAATTTGCCTGGCGGCCGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAG‑‑TGAAACGCCGTAGCGCCg                                                                  >  2:304143/1‑149 (MQ=32)
                                                                                     ttGGCGGATGAGAGAAGATTTTCAGCCTGATACAGATTAAATCAGACCGCAGAAGCGGTCTGATGAAACATAATTTGCCTGGCGGCCGTAGCGCGGTGGTCCCTCCTGACCCCATGCCGAACTCAGAAG‑‑TGAAACGCCGTAGCGCCGAt                                                                <  1:633970/149‑1 (MQ=12)
                                                                                         cGGATGAGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGATGAAACAGAATTTGCCTGGCGGCCGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAG‑‑TGAAACGCCGTAGCGCCGATGGTa                                                            <  2:204731/149‑1 (MQ=32)
                                                                                                 agaagaTTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGATGAAACAGAATTTGCCTGGCGGCCGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAG‑‑TGAAACGCCGTAGCGCCGATGGTAGTGTgggg                                                    >  1:780487/1‑149 (MQ=32)
                                                                                                  gaagaTTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGATGAAACAGAATTTGCCTGGCGGCCGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAG‑‑TGAAACGCCGTAGCGCCGATGGTAGTGTggagc                                                   <  2:614483/149‑4 (MQ=17)
                                                                                                             gCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGATGAAACAGAATTTGCCTGGCGGCCGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAG‑‑TGAACCGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCCCATGCg                                        >  1:665086/1‑149 (MQ=17)
                                                                                                              ccTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGATGAAACAGAATTTGCCTGGCGGCCGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAG‑‑TGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCCCATGCga                                       >  1:574954/1‑149 (MQ=32)
                                                                                                              ccTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGATGAAACAGAATTTGCCTGGCGGCCGTAGCACGGTGGTCCCACATGACCCCATGCCGAACTCAGAAG‑‑TGAAACGCCGTAGCGCCGACGGTAGGGTGGGGTCTCCCAATGCga                                       >  2:760324/1‑149 (MQ=11)
                                                                                                                tGATACAGATTAAATCAGACCGCAGAAGCGGTCTGATGAAACATAATTTGCCTGGCGGCCGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAG‑‑TGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCCCATGCgaga                                     <  2:723425/149‑1 (MQ=17)
                                                                                                                  aTACAGATTAAATCAGAACGCAGAAGCGGTCAGATGAAAAAGAATTTGCCTGGAGGCCGTAACGCGGTGGTCCCACCTGACCCCATCCTATTCGGAAGATTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCCCATGCGAGAGt                                   >  1:183940/1‑149 (MQ=11)
                                                                                                                    aCAGATTAAATCAGAACGCAGAAGCGGTCTGATGAAACAGAATTTGCCTGGCGGCCGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAG‑‑TGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCCCATGCGAGAGTAg                                 <  1:797802/149‑1 (MQ=32)
                                                                                                                        aTTAAATCAGAACGCAGAAGCGGACTAATGAAACAGAATTAGCCTGGCGACCGTAGCGCGGTGGTCCCACATGACCCCATGCCGAACTCAGAAG‑‑TGAAATGCCGTAGCGCCGACGGTAGTGGTGGGTGTCCCCATGGCAGAGTAGGGaa                             >  2:173266/1‑149 (MQ=11)
                                                                                                                            aaTCAGAACGCAGAAGCGGTCGGAGGAAACAGAATTTGCCTGGCGGCCGCAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAG‑‑TGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGc                         <  1:570350/149‑1 (MQ=12)
                                                                                                                              tCAGAACGCAGAAGCGGTCTGATGAAACAGAATTTGCCTGGAGGCCGTAGCGCGGAGGTCCAAGCTGACCCCATGCCGAACTCAGAAG‑‑TGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCACCCCATGCGAGAGTAGGGAACTGCCa                       >  2:255637/1‑149 (MQ=11)
                                                                                                                                        gAAGCGGTCTGATGAAACAGAATTTGCCTGGCGGCCGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAG‑‑TGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAAt             >  2:450620/1‑149 (MQ=32)
                                                                                                                                           gCGGTCTGATGAAACCGAATTTGCCTGGCGGCCGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAG‑‑TGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATTaa          <  1:590113/149‑1 (MQ=17)
                                                                                                                                               tCTGATGAAACAGAATTTGCCTGGCGGCCGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAG‑‑TGAAACGCCGTAGCGCCGATGGTAGTGTGGGGGCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATTAAGCAg      >  1:758258/1‑149 (MQ=17)
                                                                                                                                                  gatgaAACAGAATTTGCCTGGCGGCCGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAG‑‑TGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATTAAGCAGTaa   <  2:14059/149‑1 (MQ=32)
                                                                                                                                                   atgaAACAGAATTTGCCTGGCGGCCGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAG‑‑TGAAACGCCGTAGCGCCGACGGTAGTGTCGGGTCTCCCCATGCGAGAGTAAGGAACTGCCAGGCATCAAATTAAGCAGTAAg  >  2:18497/1‑149 (MQ=12)
                                                                                                                                                    |                                                                                                                                                     
AGGCCGGGTGTGTAAGCGCAGCGATGCGTTGAGCTAACCGGTACTAATGAACCGTGAGGCTTAACCTTACAACGCCGAAGATGTTTTGGCGGATGAGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGATGAAACAGAATTTGCCTGGCGGCCGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAG‑‑TGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATTAAGCAGTAAG  >  NC_000913_3_bme_pgi/3946540‑3946835

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.