New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913_3_pae_pgi = 404068350 (1.120)13 (0.300) 7/274 4.7 21.6% intergenic (+47/+223) rrfA/mobB 5S ribosomal RNA of rrnA operon/molybdopterin‑guanine dinucleotide biosynthesis protein B
?NC_000913_3_pae_pgi = 4213162 45 (1.030)intergenic (+3/‑72) rrfE/yjaA 5S ribosomal RNA of rrnE operon/stress‑induced protein
Rejected: Coverage evenness skew score above cutoff.

TCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCTCATGCGAGAGTAGGGAACTGCCAGGCATCAAATAAAACGAAAGGCTCAGTCGGAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGACAAATCCGCCGGGAGCGGATTTGAACGTTGCGAAGCAACGGCCCGGAGGGTGGCGGGCAGGACGCCCGCCATAAACTGCCAGGCATCAAA  >  NC_000913_3_pae_pgi/4040541‑4040808
                                                                                                                                              |                                                                                                                             
gcccACCTGACCCAATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCTCAGGCGAGAGTAGGGAACTGCCAGGCATCTAATAATACGACAGGCACAGTCGGAAGACTGGGCCTTTCGGTGTATcgtttg                                                                                                                         <  2:684747/148‑6 (MQ=11)
   cACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCTCATGCGAGAGTAGGGAACTGCCAGGCATCAAATAAAACGAAAGGCTCAGTCGGAAGACTGGGCCTTTCGTTTTATCTGttgttt                                                                                                                      <  1:47130/149‑1 (MQ=35)
       tGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCTCATGCGAGAGTAGGGAACTGCCAGGCATCAAATAAAACGAAAGGCTCAGTCGGAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCg                                                                                                                  >  1:535716/1‑149 (MQ=35)
       tGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCTCATGCGAGAGTAGGGAACTGCCAGGCATCAAATAAAACGAAAGGCTCAGTCGGAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCg                                                                                                                  <  1:443512/149‑1 (MQ=35)
         aCCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCTCATGCGAGAGTAGGGAACTGCCAGGCATCAAATAAAACGAAAGGCTCAGTCGGAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGt                                                                                                                <  2:535716/149‑1 (MQ=35)
         aCCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCTCATGCGAGAGTAGGGAACTGCCAGGCATCAAATAAAACGAAAGGCTCAGTCGGAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGt                                                                                                                <  1:251837/149‑1 (MQ=35)
           cccATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCTCATGCGAGAGTAGGGAACTGCCAGGCATCAAATAAAACGAAAGGCTCAGTCGGAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGa                                                                                                              >  2:233142/1‑149 (MQ=35)
            ccATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCTCATGCGAGAGTAGGGAACTGCCAGGCATCAAATAAAACGAAAGGCTCAGTCGGAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGaa                                                                                                             <  2:358155/149‑1 (MQ=35)
                gCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCTCATGCGAGAGTAGGGAACTGCCAGGCATCAAATAAAACGAAAGGCTCAGTCGGAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGct                                                                                                         >  2:37801/1‑149 (MQ=35)
                   gaccgcAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCTCATGCGAGAGTAGGGAACTGCCAGGCATCAAATAAAACGAAAGGCTCAGTCGGAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCtctc                                                                                                      <  1:247124/144‑1 (MQ=25)
                   gAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCTCATGCGAGAGTAGGGAACTGCCAGGCATCAAATAAAACTAAAGGCTCAGTCGGAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCtctc                                                                                                      <  1:265135/149‑1 (MQ=21)
                        cAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCTCATGCGAGAGTAGGGAACTGCCAGGCATCAAATAAAACGAAAGGCTCAGTCGGAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAg                                                                                                 >  1:42306/1‑149 (MQ=35)
                                 aaCGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCTCATGCGAGAGTAGGGAACTGCCAGGCATCAAATAAAACGAAAGGCTCAGTCGGAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGACaaa                                                                                        >  2:22961/1‑149 (MQ=35)
                                              cGATGGTAGTGTGGGGTCTCCTCATGCGAGAGTAGGGAACTGCCAGGCATCAAATAAAACGAAAGGCTCAGTCGGAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGACAAATCGGCCGGGAGCg                                                                           >  1:430669/1‑149 (MQ=25)
                                              cGATGGTAGTGTGGGGTCTCCTCATGCGAGAGTAGGGAACTGCCAGGCATCAAATAAAACGAAAGGCTCAGTCGGAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGACAAATCCGCCGGGAGCg                                                                           <  1:313894/149‑1 (MQ=35)
                                                aTGGTAGTGTGGGGTCTCCTCATGCGAGAGTAGGGAACTGCCAGGCATCAAATAAAACGAAAGGCTCAGTCGGAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGACAAATCCGCCGGGAGCGGa                                                                         <  2:573671/149‑1 (MQ=35)
                                                 tGGTAGTGTGGGGTCTCCTCATGCGAGAGTAGGGAACTGCCAGGCATCAAATAAAACGAAAGGCTCAGTCGGAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGACAAATCCGCCGGGAGCGGAt                                                                        <  2:42306/149‑1 (MQ=35)
                                                  ggtggtGAGAGGTCTCCTCATGCGAGAGTAGGGAACGGCCAGGGAGCAGATAAAACGAAAGGCTCAGTCGGAAGACTGGGCCTTTCGTTGTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGACAAATCCGCCGGGAGCGGAtt                                                                       <  2:9364/149‑1 (MQ=11)
                                                  ggtggggtgGGGTCTGCTCATGCGAGAGTAGGGAACTGCCAGGCATCAAATAAAGCGAAAGGCTCAGTCGGAAGACTGGGGCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGACAAATCCGCCGGGAGCGGAtt                                                                       <  2:13079/143‑1 (MQ=11)
                                                       tgtgGGGTCTCCTCATGCGAGAGTAGGGAACTGCCAGGCATCAAATAAAACGAAAGTCGCAGTCGGAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGACAAATCCGCCGGGAGCGGGTTTGAAc                                                                  >  1:548266/1‑149 (MQ=255)
                                                       tgaggggTCTCCTCATGCGCGAGTAGGGAACTGCCAGGCATCAAATAAAACGCAAGGCTCAGTCGGAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTATACTGCGTAGGCAGCCTGATCAGGGAGCGGATTTGAAc                                                                  >  2:749645/4‑149 (MQ=11)
                                                           gaggctcctcATGCGAGAGTAGGGAACTGCCAGGCATCAAATAAAACGAAAGGCTCAGTCGGAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGACAAATCCGCCGGGAGCGGATTTGAACGTTg                                                              <  2:657052/145‑1 (MQ=25)
                                                              tctcCTCATGCGAGAGTAGGGAACTGCCAGGCATCAAATAAAACGAAAGGCTCAGTCGGAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGACAAATCCGCCGGGAGCGGATTTGAACGTTGCGa                                                           >  2:492257/1‑149 (MQ=35)
                                                                tcctcATGCGAGAGTAGGGAACTGCCAGGCATCAAATAAAACGAAAGGCTCAGTCGGAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTAAGGAGGACAAACCCGCCGGGCGCGGATATGAACGTTGCGagg                                                         >  2:399955/1‑147 (MQ=11)
                                                                 cctcATGCGAGAGTAGGGAACTGCCAGGCATCAAATAAAACGAAAGGCTCAGTCGGAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGACAAATCCGCCGGGAGCGGATTTGAACGTTGCGAAGc                                                        <  1:37801/149‑1 (MQ=35)
                                                                  ctcATGCGAGAGTAGGGAACTGCCAGGCATCCAATTAAACGAAAGGCTCAGTCGGAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGACAAATCCGCCGGGAGCGGATTTGAACGTTGCGAAGCa                                                       >  2:170609/1‑149 (MQ=21)
                                                                   tcATGCGGGAGTAGGGAACGGCCAGGCATCAAATAAAACGAAAGGCTCAGTCGGAAGACTGGGCCTTTCGTTTTATCTGTGGTTTGTCGGTGAACGCTCTCCTGAGTAGGACAAATCCGACGGGAGCGGATTTGAACGTTGCGAAGCaa                                                      <  2:404996/149‑1 (MQ=255)
                                                                     aTGCGAGAGTAGGGAACTGCCAGGCATCAAATAAAACGAAAGGCTCAGTCGGAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGACAAATCCGCCGGGAGCGGATTTGAACGTTGCGAAGCAACg                                                    >  1:45180/1‑149 (MQ=32)
                                                                        cGAGAGTAGGGAACTGCCAGGCATCAAATAAAACGAAAGGCTCAGTCGGAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGACAAATCCGCCGGGAGCGGATTTGAACGTTGCGAAGCAACGGcc                                                 >  1:317761/1‑149 (MQ=32)
                                                                         gagaGTAGGGAACTGCCAGGCATCAAATAAAACGAAAGGCTCAGTCGGAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGACAAATCCGCCGGGAGCGGATTTGAACGTTGCGAAGCAACGGccc                                                >  2:402465/1‑149 (MQ=32)
                                                                         gacagTAGGTAATTGCCATGCATCAAAGAAGACTAAAGGCTCAGTCGGAAGTGTGTGCCTTTCGGTTTAGCTGTTGTTTGTCGGTGAACTCTCTCCTGAGTAGGACAAATCCGCCGGGAGCGGATTTGAACGTTGCGAAGCAACGGccc                                                <  1:749645/146‑1 (MQ=11)
                                                                            agTAGGGAACTGCCAGGCATCAAATAAAACGAAAGGCTCAGTCGGAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGACAAATCCGCCGGGAGCGGATTTGAACGTTGCGAAGCAACGGCCCGGa                                             <  2:520092/149‑1 (MQ=32)
                                                                                   aaCTGCCAGGCATCAAATAAAACGAAAGGCTCAGTCGGAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGACAAATACGCCGGGAGCGGATTTGAACGTTGCGAAGCAACGGCCCGGCGGGTGGc                                      >  2:757531/1‑149 (MQ=17)
                                                                                     cTGCCAGGCAGCAAATAAAACGAAAGGCTCAGTCGGAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGACAAATCCGCCGGGAGCGGATTTGAACGTTGCGAAGCAACGGCCCGGAGGGTGGCgg                                    <  1:317831/149‑1 (MQ=17)
                                                                                        ccAGGCATCAAATAAAACGAAAGGCTCAGTCGGAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGACAAATCCGCCGGGAGCGGATTTGAACGTTGCGAAGCAACGGCCCGGAGGGTGGCGGGCa                                 >  1:763110/1‑149 (MQ=32)
                                                                                           ggCATCAAATAAACCGAAAGGATCCGGCGGAAGACTGGGCCTTTCGTTGTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGACAAATCCGCCGGGAGCGGATTTGAACGTTGCGAAGCAACGGCCCGGAGGGTGGCGGGCAGGa                              <  2:253055/149‑1 (MQ=11)
                                                                                                 aaaTAAAACGAAAGGCTCAGTCGGAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGACCGCTCTCCTGAGTAGGACAAATCCGCCGGGAGCGGATTTGAACGTTGCGAAGCACCGGCCCGGAGGGTGGCGGGCAGGAcgcccg                        >  1:212205/1‑149 (MQ=17)
                                                                                                 aaaTAAAACGAAAGGCTCAGTCGGAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGACAAATCCGCCGGGAGCGGATTTGAACGTTGCGAAGCAACGGCCCGGAGGGTGGCGGGCAGGAcgcccg                        >  1:674111/1‑149 (MQ=32)
                                                                                                  aaTAAAACGAAAGGCTCAGTCGGAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGACAAATCCGCCGGGAGCGGATTTGAACGTTGCGAAGCAACGGCCCGGAGGGTGGCGGGGAGGAcggccgc                       >  1:274944/1‑149 (MQ=17)
                                                                                                  aaTAAAACGAAAGGCTCAGTCGGAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGACAAATCCGCCGGGAGCGGATTTGAACGTTGCGAAGCAACGGCCCGGAGGGTGGCGGGCAGGAcgcccgc                       >  2:77506/1‑149 (MQ=32)
                                                                                                      aaaCGAAAGGCTCAGTCGGAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGACAAATCCGCCGGGAGCGGATTTGAACGTTGCGAAGCAACGGCCCGGAGGGTGGCGGGCAGGACGCCCGCCATa                   >  2:756329/1‑149 (MQ=32)
                                                                                                        aCGAAAGGCTCAGTCGGAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGACAAATCCGCCGGGAGCGGATTTGAACGTTGCGAAGCAACGGCCCGGAGGGTGGCGGGCAGGACGCCCGCCATaaa                 <  1:233142/149‑1 (MQ=32)
                                                                                                         cGAAAGGCTCAGTCGGAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGACAAATCCGCCGGGAGCGGATTTGAACGATACGAAGCAACGTCCCGGAGGGTGGCGGGCAGGACGACCGCCATAAAc                >  1:457623/1‑149 (MQ=255)
                                                                                                              ggCTCAGTCGGAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGACAAATCCGCCGGGAGCGGATTTGAACGTTGCGAAGCAACGGCCCGTAGGGTGGCGGGCAGGACGCCCGCCATAAACTGCCa           <  2:317761/149‑1 (MQ=17)
                                                                                                                  cAGTCGGAAACCGGGGCCTAACGTATTATCTGTTGTTTGTCGGTGAACGCCGCCCTGAGGAGGACAAAACAGTCCGGAGCGGATTTGAACGTTGCGAAGCAACAGCCCGGAGGGTGGCAGGAAGGCGGACAGCCATAAAATgccagcca       <  2:722658/149‑1 (MQ=2)
                                                                                                                       ggAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGACAAATCCGCCGGGAGCGGATTTGAACGTTGCGAAGCAACGGCCCGGAGGGTGGCGGGCAGGAc                             <  2:15841/122‑1 (MQ=31)
                                                                                                                       ggAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGACAAATCCGCCGGGAGCGGATTTGAACGTTGCGAAGCAACGGCCCGGAGGGTGGCGGGCAGGAc                             >  1:15841/1‑122 (MQ=31)
                                                                                                                       ggAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGACAAATCCGCCGGGAGCGGATTTGAACGTTGCGAAGCAACGGCCCGGAGGGTGGCGGGCAGGACGCCCGCCATAAACTGCCAGGCATCaaa  >  1:489289/1‑149 (MQ=31)
                                                                                                                       ggAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGACAAATCCGCCGGGAGCGGATTTGAACGTTGCGAAGCAACGGCCCGGAGGGTGGCGGGCAGGACGCCCGCCATAAACTGCCAGGCATCaaa  >  1:182274/1‑149 (MQ=31)
                                                                                                                       ggAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCCCCGGAGTAGGACAAATCCGCCGGGAGCGGATTTGAACGCTGCGAAGCAACGGCCCGGAGGGTGGCGGGCAGGACGCCCGCCATAAACTGCCAGGCATCaaa  <  1:19582/149‑1 (MQ=9)
                                                                                                                                              |                                                                                                                             
TCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCTCATGCGAGAGTAGGGAACTGCCAGGCATCAAATAAAACGAAAGGCTCAGTCGGAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGACAAATCCGCCGGGAGCGGATTTGAACGTTGCGAAGCAACGGCCCGGAGGGTGGCGGGCAGGACGCCCGCCATAAACTGCCAGGCATCAAA  >  NC_000913_3_pae_pgi/4040541‑4040808

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