New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913 = 403524426 (0.220)7 (0.060) 4/168 NT 22.2% intergenic (+91/‑287) hemG/rrsA protoporphyrin oxidase, flavoprotein/16S ribosomal RNA of rrnA operon
?NC_000913 = 4035271 NA (NA)intergenic (+118/‑260) hemG/rrsA protoporphyrin oxidase, flavoprotein/16S ribosomal RNA of rrnA operon

CTGAAATAAGCATAAAGAATAAAAAATGCGCGGTCAGAAAATTATTTTAAATTTCCTCTTGTCAGGCCGGAATAACTCCCTATAATGCGCCACCACTGACACGGAACAACGGCAAACACGCCGCCGGGTCAGCGGGGTTCTCCTGAGAACTCCGGCAGAGAAAGCAAAAATAAATGCTTGACTCTGTAGCGGGAAGGCG  >  NC_000913/4035145‑4035343
                                                                                                   |                                                                                                   
cTGAAATAAGCATAAAGAATAAAAAATGCCCGGTCAGAAAATTATATTAAATTTCCTCTTGTCACGCCGGAATAACTCCCTATAATGCGCCACCACTGAC                                                                                                     >  2:2268659‑M1/1‑100 (MQ=255)
  gAAATAAGCATAAAGAATAAAAAATGCGCGGTCAGAAAATTATTTTAAATTTCCTCTTGTCAGGCCGGAATAACTCCCTATAATGCGCCACCACTGacac                                                                                                   >  2:3228383/1‑100 (MQ=255)
  gAAATAAGCATAAAGAATAAAAAATGCGCGGTCAGAAAATTATTTTAAATTTCCTCTTGTCAGGCCGGAATAACTCCCTATAATGCGCCACCACTGacac                                                                                                   >  2:2490896/1‑100 (MQ=255)
    aaTAAGCATAAAGAATAAAAAATGCGCGGTCAGAAAATTATTTTAAATTTCCTCTTGTCAGGCCGGAATAACTCCCTATAATGCGCCACCACTGACACgg                                                                                                 <  1:2579254/100‑1 (MQ=255)
     aTAAGCATAAAGAATAAAAAATGCGCTGTCAGAAAATTATTTTAAATTTCCTCTTGTCAGGCCGGAATAACTCCCTATAATGCGCCACCACTGACACGGa                                                                                                <  1:3074465/100‑1 (MQ=255)
     aTAAGCATAAAGAATAAAAAATGCGCGGTCAGAAAATTATTTTAAATTTCCTCTTGTCAGGCCGGAATAACTCCCTATAATGCGCCACCACTGACACGGa                                                                                                <  1:1565450/100‑1 (MQ=255)
     aTAAGCATAAAGAATAAAAAATGCGCGGTCAGAAAATTATTTTAAATTTCCTCTTGTCAGGCCGGAATAACTCCCTATAATGCGCCACCACTGACACGGa                                                                                                >  2:493091/1‑100 (MQ=255)
     aTAAGCATAAAGAATAAAAAATGCGCGGTCAGAAAATTATTTTAAATTTCCTCTTGTCAGGCCGGAATAACTCCCTATAATGCGCCACCACTGACACGGa                                                                                                <  1:3003894/100‑1 (MQ=255)
      tAAGCATAAAGAATAAAAAATGCGCGGTTAGAAAATTATTTTAAATTTCCTCTTGTCAGGCCGGAATAACTCCCTATAATGCGCCACCACTGACACGGaa                                                                                               <  1:691284/100‑1 (MQ=255)
       aaGCATAAAGAATAAAAAATGCGCGGTCAGAAAATTATTTTAAATTTCCTCTTGTCAGGCCGGAATAACTCCCTATAATGCGCCACCACTGACACGGaac                                                                                              <  1:1448043/100‑1 (MQ=255)
       aaGCATAAAGAATAAAAAATGCGCGGTCAGAAAATTATTTTAAATTTCCTCTTGTCAGGCCGGAATAACTCCCTATAATGCGCCACCACTGACACGGaac                                                                                              <  1:2024207/100‑1 (MQ=255)
         gCATAAAGAATAAAAAATGCGCGGTCAGAAAATTATTTTAAATTTCCTCTTGTCAGGCCGGAATAACTCCCTATAATGCGCCACCACTGACACGGaacaa                                                                                            >  2:2430248/1‑100 (MQ=255)
          cATAAAGAATAAAAAATGCGCGGTCAGAAAATTATTTTAAATTTCCTCTTGTCAGGCCGGAATAACTCCCTATAATGCGCCACCACTGACACGGaacaac                                                                                           >  2:1407465/1‑100 (MQ=40)
          cATAAAGAATAAAAAATGCGCGGTCAGAAAATTATTTTAAATTTCCTCTTGTCAGGCCGGAATAACTCCCTATAATGCGCCACCACTGACACGGaacaac                                                                                           >  1:345714/1‑100 (MQ=40)
          cATAAAGAATAAAAAATGCGCGGTCAGAAAATTATTTTAAATTTCCTCTTGTCAGGCCGGAATAACTCCCTATAATGCGCCACCACTGACACGGaacaac                                                                                           >  2:169458/1‑100 (MQ=40)
          cATAAAGAATAAAAAATGCGAGGTCAGAAAATTATTTTAAATTTCCTCTTGTCAGGCCGGAATAACCCCCTATAATGCGCCACCACTGACACGGaacaac                                                                                           >  2:1987356/1‑100 (MQ=255)
           aTAAAGAATAAAAAATGCGCGGTCAGAAAATTATTTTAAATTTCCTCTTGTCAGGCCGGAATAACTCCCTATAATGCGCCACCACTGACACGGAACAACg                                                                                          >  1:2947466/1‑100 (MQ=40)
           aTAAAGAATAAAAAATGCGCGGTCAGAAAATTATTTTAAATTTCCTCTTGTCAGGCCGGAATAACTCCCTATAATGCGCCACCACTGACACGGAACAACg                                                                                          <  1:1211413/100‑1 (MQ=40)
           aTAAAGAATAAAAAATGCGCGGTCAGAAAATTATTTTAAATTTCCTCTTGTCAGGCCGGAATAACTCCCTATAATGCGCCACCACTGACACGGAACAACg                                                                                          <  1:2577955/100‑1 (MQ=40)
            tAAAGAATAAAAAATGCGCGGTCAGAAAATTATTTTAAATTTCCTCTTGTCAGGCCGGAATAACTCCCTATAATGCGCCACCACTGACACGGAACAACgg                                                                                         >  1:482643/1‑100 (MQ=38)
            tAAAGAATAAAAAATGCGCGGTCAGAAAATTATTTTAAATTTCCTCTTGTCAGGCAGGAATAACTCCCTATAATGCGCCACCACTGACACGGAACAACgg                                                                                         >  2:244778/1‑100 (MQ=255)
             aaaGAATAAAAAATGCGCGGTCAGAAAATTATTTTAAATTTCCTCTTGTCAGGCCGGAATAACTCCCTATAATGCGCCACCACTGACACGGAACAACGGc                                                                                        >  1:2526242/1‑100 (MQ=38)
              aaGAATAAAAAATGCGCGGTCAGAAAATTATTTTAAATTTCCTCTTGTCAGGCCGGAATAACTCCCTATAATGCGCCACCACTGACACGGAACAACGGCa                                                                                       <  1:1461663/100‑1 (MQ=35)
              aaGAATAAAAAATGCGCGGTCAGAAAATTATTTTAAATTTCCTCTTGTCAGGCCGGAATAACTCCCTATAATGCGCCACCACTGACACGGAACAACGGCa                                                                                       <  2:45982/100‑1 (MQ=35)
              aaGAATAAAAAATGCGCGGTCAGAAAATTATTTTAAATTTCCTCTTGTCAGGCCGGAATAACTCCCTATAATGCGCCACCACTGACACGGAACAACGGCa                                                                                       <  2:108583/100‑1 (MQ=35)
              aaGAATAAAAAATGCGCGGTCAGAAAATTATTTTAAATTTCCTCTTGTCAGGCCGGAATAACTCCCTATAATGCGCCACCACTGACACGGAACAACGGCa                                                                                       <  1:2234002/100‑1 (MQ=35)
              aaGAATAAAAAATGCGCGGTCAGAAAATTATTTTAAATTTCCTCTTGTCAGGCCGGAATAACTCCCTATAATGCGCCACCACTGACACGGAACAACGGCa                                                                                       <  1:2234013/100‑1 (MQ=35)
                gAATAAAAAATGCGCGGTCAGAAAATTATTTTAAATTTCCTCTTGTCAGGCCGGAATAACTCCCTATAATGCGCCACCACTGACACGGAACAACGGCaaa                                                                                     >  2:471209/1‑100 (MQ=33)
                 aaTAAAAAATGCGCGGTCAGAAAATTATTTTAAATTTCCTCTTGTCAGGCCTGAATATCTCCCTATAATGCGCCACCACTTACACGGAACAACGGCAAac                                                                                    <  2:1933322/100‑1 (MQ=255)
                 aaTAAAAAATGCGCGGTCAGAAAATTATTTTAAATTTCCTCTTGTCAGGCCGGAATAACTCCCTATAATGCGCCACCACTGACACGGAACAACGGCAAac                                                                                    <  1:1667979/100‑1 (MQ=32)
                 aaTAAAAAATGCGCGGTCAGAAAATTATTTTAAATTTCCTCTTGTCAGGCCGGAATAACTCCCTATAATGCGCCACCACTGACACGGAACAACGGCAAac                                                                                    <  1:2236639/100‑1 (MQ=32)
                  aTAAAAACTGCGCGGTCAGAAAATTATTTTAAATTTCCTCTTGTCAGGCCGGAATAACTCCCTATGATGCGCCACCACTGACACGGAACAACGGCAAaca                                                                                   >  2:1705840/1‑100 (MQ=9)
                   tAAAAAATGCGCGGTCAGAAAATTATTTTAAATTTCCTCTTGTCAGGCCGGAATAACTCCCTATAATGCGCCACCACTGACACGGAACAACGGCAAac                                                                                    <  2:144587/98‑1 (MQ=11)
                   tAAAAAATGCGCGGTCAGAAAATTATTTTAAATTTCCTCTTGTCAGGCCGGAATAACTCCCTATAATGCGCCACCACTGACACGGAAAAAAAGCaaaaa                                                                                   >  1:32373/1‑97 (MQ=11)
                               ggTCAGAAAATTATTTTAAATTTCCTCTTGTCAGGCCGGAATAACTCCCTATAATGCGCCACCACTGACccggcggcgtgtttgccgttgttccgtgtca                                                                      <  1:2222261‑M1/100‑32 (MQ=255)
                                 tCAGAAAATTATTTTAAATTTCCTCTTGTCAGGCCGGAATAACTCCCTATAATGCGCCACCACTGACccggcggcgtgtttgccgttgttccgtgtcagt                                                                    <  1:394384‑M1/100‑34 (MQ=255)
                                 tCAGAAAATTATTTTAAATTTCCTCTTGTCAGGCCGGAATAACTCCCTATAATGCGCCACCACTGACccggcggcgtgtttgccgttgttccgtgtcagt                                                                    <  2:1150597‑M1/100‑34 (MQ=255)
                                 tCAGAAAATTATTTTAAATTTCCTCTTGTCAGGCCGGAATAACTCCCTATAATGCGCCACCACTGACccggcggcgtgtttgccgttgttccgtgtcagt                                                                    <  1:393631‑M1/100‑34 (MQ=255)
                                  cAGAAAATTATTTTAAATTTCCTCTTGTCAGGCCGGAATAACTCCCTATAATGCGCCACCACTGACccggcggcgtgtttgccgttgttccgtgtcagtg                                                                   <  1:1660364‑M1/100‑35 (MQ=255)
                                  cAGAAAATTATTTTAAATTTCCTCTTGTCAGGCCCGAATAACTCCCTATAATGCGCCACCACTGACccggcggcgtgtttgccgttgttccgtgtcagtg                                                                   <  2:1710355‑M1/100‑35 (MQ=255)
                                   aGAAAATTATTTTATATTTCCTCTTGTCAGGCCGGAATAACTCCCTATAATGCGCCACCACTGACccggcggcgtgtttgccgttgttccgtgtcagtgg                                                                  <  2:2209228‑M1/100‑36 (MQ=255)
                                                                                           gatcttgacccggcggcgtgtttgccgttgttccgtgtcagtggtggcgcattatagggagttattccggcctgacaagaggaaatttaaaataattttc          <  2:3213769‑M1/95‑92 (MQ=255)
                                                                                               cTGACccggcggcgtgtttgccgttgttccgtgtcagtggtggcgcattatagggagttattccggcctgacaagaggaaatttaaaataattttctgac      >  1:1566877‑M1/1‑5 (MQ=255)
                                                                                                tGACccggcggcgtgtttgccgttgttccgtgtcagtggtggcgcattatagggagttattccggcctgacaagaggaaatttaaaataattttcagatc     >  1:3213797‑M1/1‑4 (MQ=255)
                                                                                                 gACccggcggcgtgtttgccgttgttccgtgtcagtggtggcgcattatagggagttattccggcctgacaagaggaaatttaaaataattttttgaccg    >  1:2992246‑M1/1‑3 (MQ=255)
                                                                                                 gACccggcggcgtgtttgccgttgttccgtgtcagtggtggcgcattatagggagttattccggcctgacaagaggaaatttaaaataattttctgaccg    >  2:1066555‑M1/1‑3 (MQ=255)
                                                                                                 gACccggcggcgtgtttgccgttgttccgtgtcagtggtggcgcattatagggagttattccggcctgacaagaggaaatttaaaataattttctgaccg    >  1:2886867‑M1/1‑3 (MQ=255)
                                                                                                  aCccggcggcgtgtttgccgttgttccgtgtcagtggtggcgcattatagggagttattccggcctgacaagaggaaatttaaaataattttctgaccgc   >  2:848663‑M1/1‑2 (MQ=255)
                                                                                                   cccgggggcgtgtttgccgttgttccgtgtcagtggtggcgcattatagggagttattccggcctgacaagaggaaatttaaaataattttctgaccgcg  <  2:2684730‑M1/100‑100 (MQ=255)
                                                                                                   cccggcggcgtgtttgccgttgttccgtgtcagtggtggcgcattatagggagttattccggcctgacaagaggaaatttagaataattttctgaccgcg  >  1:3192196‑M1/1‑1 (MQ=255)
                                                                                                   cccggcggcgtgtttgccgttgttccgtgtcagtggtggcgcattatagggagttattccggcctgacaagaggaaatttaaaataattttctgaccgcg  >  1:444681‑M1/1‑1 (MQ=255)
                                                                                                   cccggcggcgtgtttgccgttgttccgtgtcagtggtggcgcattatagggagttattccggcctgacaagaggaaatttaaaataattttctgaccgcg  >  2:232570‑M1/1‑1 (MQ=255)
                                                                                                   cccggcggcgtgtttgccgttgttccgtgtcagtggtggcgcattatagggagttattccggcctgacaagaggaaatttaaaataattttctgaccgcg  <  2:195902‑M1/100‑100 (MQ=255)
                                                                                                   cccggcggcgtgtttgccgttgttccgtgtcagtggtggcgcattatagggagttattccggcctgacaagaggaaatttaaaataattttctgaccgcg  <  2:1902462‑M1/100‑100 (MQ=255)
                                                                                                   cccggcggcgtgtttgccgttgttccgtgtcagtggtggcgcattatagggagttattccggcctgacaagaggaaatttaaaataattttctgaccgcg  >  1:495357‑M1/1‑1 (MQ=255)
                                                                                                   cccggcggcgtgtttgccgttgttccgtgtcagtggtggcgcattatagggagttattccggcctgacaagaggaaatttaaaataattttctgaccgcg  >  1:431082‑M1/1‑1 (MQ=255)
                                                                                                   cccggcggcgtgtttgccgttgttccgtgtcagtggtggcgcattatagggagttattccggcctgacaagaggaaatttaaaataattttctgaccgcg  >  1:2923080‑M1/1‑1 (MQ=255)
                                                                                                   cccggcggcgtgtttgccgttgttccgtgtcagtggtggcgcattatagggagttattccggcctgacaagaggaaatttaaaataattttctgaccgcg  >  1:1809526‑M1/1‑1 (MQ=255)
                                                                                                   cccggcggcgtgtttgccgttgttccgtgtcagtggtggcgcattatagggagttattccggcctgacaagaggaaatttaaaataattttctgaccgcg  <  1:1217625‑M1/100‑100 (MQ=255)
                                                                                                   cccggcggcgtgtttgccgttgttccgtgtcagtggtggcccattatagggagttattccggcctgacaagaggaaatttaaaataattttctgaccgcg  >  2:2824073‑M1/1‑1 (MQ=255)
                                                                                                   |                                                                                                   
CTGAAATAAGCATAAAGAATAAAAAATGCGCGGTCAGAAAATTATTTTAAATTTCCTCTTGTCAGGCCGGAATAACTCCCTATAATGCGCCACCACTGACACGGAACAACGGCAAACACGCCGCCGGGTCAGCGGGGTTCTCCTGAGAACTCCGGCAGAGAAAGCAAAAATAAATGCTTGACTCTGTAGCGGGAAGGCG  >  NC_000913/4035145‑4035343

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