New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913 3423665 =NA (NA)6 (0.140) 5/280 5.6 10.8% intergenic (‑10/+3) thrV/rrfD tRNA‑Thr/5S ribosomal RNA of rrnD operon
?NC_000913 = 4040683 50 (1.120)intergenic (+47/+223) rrfA/mobB 5S ribosomal RNA of rrnA operon/molybdopterin‑guanine dinucleotide biosynthesis protein B
Rejected: Coverage evenness skew score above cutoff.
Rejected: Frequency below/above cutoff threshold.

TCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCTCATGCGAGAGTAGGGAACTGCCAGGCATCAAATAAAACGAAAGGCTCAGTCGGAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGACAAATCCGCCGGGAGCGGATTTGAACGTTGCGAAGCAACGGCCCGGAGGGTGGCGGGCAGGACGCCCGCCATAAACTGCCAG  >  NC_000913/4040541‑4040800
                                                                                                                                              |                                                                                                                     
gcccACCTGACCCCATGCCGAACGCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCTCATGCGAGAGTAGGGAACTGCCAGGCATCAAATAAAACGAAAGGCTCAGTCGGAAGACTGGGCCTTTCGTTTTAt                                                                                                                       <  1:100673/142‑1 (MQ=25)
   cACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCTCATGCGAGAGTAGGGAACTGCCAGGCATCAAATAAAACGAAAGGCTCAGTCGGCAGCCTGGGCCTTTCGTTTTATCtg                                                                                                                    >  2:268726/1‑143 (MQ=21)
   cACCTGACCCCATGCCGAACGCAGAAGTGAAAAGCCGTAGTGCAGATGGTAGTGTGGGGGCTCCTCATGCGAGAGTAGGGAACTACCAGGCATCAAATAAAACGAAAGGCTCAGTCGGAAGACGGGGCCGTTGGGATTATCtg                                                                                                                    >  1:102315/1‑143 (MQ=11)
    accgggccccCTGCCGAACTCAGAAGGCACACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCTCATGCGAGAGTATGGAACTGCCAGGCATCAAATAAAACGAAAGGCTCAGTCGGAAGACTGGGCCTTTCGTTTTATCtgt                                                                                                                   <  2:151874/137‑1 (MQ=11)
    aCCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCTCATGCGAGAGTAGGGAACTGCCAGGCATCAAATAAAACGAAAGGCTCCGTCGGCAGACTGAGCCTTTCGGTTAATCtgt                                                                                                                   >  2:794965/1‑143 (MQ=11)
           cccATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCTCATGCGAGAGTAGGGAACTGCCAGGCATCAGATAAAACGAAAGGCTCAGTCGGAAGACTGGGCCTTTCGTTTTATCTGttgtttgt                                                                                                            <  1:156611/143‑1 (MQ=25)
           cccATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCTCATGCGAGAGTAGGGAACTGCCAGGCATCAAATAAAACGAAAGGCTCAGTCGGAAGACTGGGCCTTTCGTTTTATCTGttgtttgt                                                                                                            <  1:191093/143‑1 (MQ=35)
           cccATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCTCATGCGAGAGTAGGGAACTGCCAGGCATCAAATAAAACGAAAGGCTCAGTCGGAAGACTGGGCCTTTCGTTTTATCTGttgtttgt                                                                                                            >  1:495403/1‑143 (MQ=35)
            ccATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCGCATGCGAGAGTAGGGAACTGCCAGGCATCAAATAAAACGAAAGGCTCAGTCGGAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTc                                                                                                           <  1:387470/143‑1 (MQ=17)
             cATGCCGAGCTCGGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCTCATGCGAGAGTAGGGAACTGCCAGGCATCAAATAAAACGAAAGGCTCAGTCGGAAGACTGGGCCTTGCGTTTTATCTGTTGTTTGTCg                                                                                                          <  2:238191/143‑1 (MQ=255)
                gCCGAACTCGGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCTCATGCGAGAGTAGGGAACTGCCAGGCATCAAATAAAACGAAAGGCTCAGTCGGAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTg                                                                                                       <  2:109914/143‑1 (MQ=25)
                   gAACTCAGAAGTCAAACGCCGTAGCGACGATGGTAGGGTGGGGTCTCGTCAGGCGAGAGTAGGGAACTGCCAGGCATCTAATAAAACGAAAGGCTCAGTCGGAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAAc                                                                                                    <  2:197312/143‑1 (MQ=11)
                       tCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCTCATGCGAGAGTAGGGAACTGCCAGGCATCAAATAAAACGAAAGGCTCAGTCGGAAGACTGGGCCGTTCGTTTTATCTGTTGTTTGTCGGTGAACGctc                                                                                                >  1:523008/1‑143 (MQ=25)
                           aaGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCTCATGCGAGAGTAGGGAACTGCGAGGCATCAAATAAACCGAAAGGCTCAGTCGGAAGACTGGGCCTTTCGTTTTATCTGTTGTCTGTCGGTGAACGCTCTCCt                                                                                            >  2:528541/1‑143 (MQ=255)
                            aGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCTCATGCGAGAGTAGGGAACTGCCAGGCATCAAATAAAACGAAAGGCTCAGTCGGAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTg                                                                                           >  2:683685/1‑143 (MQ=35)
                            aGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCTCATGCGAGAGTAGGGAACTGCCAGGCATCAAATAAAACGAAAGGCTCAGTCGGAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTg                                                                                           >  2:720807/1‑143 (MQ=35)
                              tGAAACGCCGTAGCACCGATTGTAGTGTGGGGCCTCTTCATGCGAGAGTAGGGAACTGCCAGGCATCTAATAAAACGAAAGGCTCAGTCGGAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGATCGCTCTCCTGAg                                                                                         <  1:616318/143‑1 (MQ=11)
                                 aaCGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCTCATGCGAGCGTAGGGAACTGCCAGGCATCAAATACAACGAAAGGCTCAGTCGGAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAg                                                                                      >  2:145734/1‑143 (MQ=21)
                                           cgcCGATGGTAGTGTGGGGTCTCCTCATGCGAGAGTAGGGAACTGCCAGGCATCAAATAAAACGAAAGGCTCAGTCGGAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGACAAATccg                                                                            <  1:73935/143‑1 (MQ=34)
                                               gATGGTAGTGTGGGGTCTCCTCATGCGAGAGTAGGGAACTGCCAGGCATCAAATAAAACGAAAGGCTCAGTCGGCAGACTGTGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGACAACTCCGCCgg                                                                        >  2:8587/1‑143 (MQ=255)
                                                     aGTGTGGGGTCTCCTCATGCGAGAGTAGGGAACTGCCAGGCATCAAATAAAACGAAAGGCTCAGTCGGAAGACTGGGCCTTTCGGTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGACAAATCCGCCGGGAGCgg                                                                  >  1:643809/1‑143 (MQ=25)
                                                        gtgGGGTCTCCTCATGCGAGAGTAGGGAACTGCCAGGCATCAAATAAAACGAAAGGCTCAGTCGGAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGACAAATCCGCCGGGGGCGGagt                                                               >  1:375737/1‑141 (MQ=25)
                                                        gtgGGGTCTCCTCATGCGAGAGTAGGGAACTGCCAGGCATCAAATAAAACGAAAGGCTCAGTCGGAAGACAGGGCCTTTCGTTTTATCTGTCGCTACTCAGTGTACGCTCTCCTGTGTAACACAAATGCGCCGGGAGCGGCtt                                                               >  2:693687/1‑143 (MQ=11)
                                                         tgGGGTCTGCTCATGCGAGAGTAGGGAACTGCCAGGCATCAAATAAAACGAAAGGCTCAGTCGGAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGACAAATCCGCCGGGAGCGGAttt                                                              <  1:239521/143‑1 (MQ=25)
                                                                  ctcATGCGAGAGTAGGGAACTGCCAGGCATCAAATAAAACGAAAGGCTCAGTCGGAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGACAAATCCGCCGGGAGCGGATTTGAACGTTGc                                                     >  1:785792/1‑143 (MQ=34)
                                                                   tcATGCGAGAGTAGGGAACTGCCAGGCATCAAATAAAACGAAAGGCTCAGTCGGAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGACAAATCCGCCGGGAGCGGATTTGAACGTTGCg                                                    <  1:720807/143‑1 (MQ=33)
                                                                   tcATGCGAGAGTAGGGAACTGCCAGGCATCAAATAAAACGAAAGGCTCAGTCGGAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGACAAATCCGCCGGGAGCGGATTTGAACGTTGCg                                                    <  1:753361/143‑1 (MQ=33)
                                                                   tcATGCGAGAGTAGGGAACTGCCAGGCATCAAATAAAACGAAAGGCTCAGTCGGAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGACAAATCCGCCGGGAGCGGATTTGAACGTTGCg                                                    >  2:255718/1‑143 (MQ=33)
                                                                    cATGCGAGAGTAGGGAACTGCCAGGCATCAAATAAAACGAAAGGCTCAGTCGGAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGACAAATCCGCCGGGAGCGGATTTGAACGTTGCGa                                                   >  1:780653/1‑143 (MQ=32)
                                                                            atgaGGGAACAGCAAGACGACAAATAAAACCAAACGCCCAGTCGGAAGACTGGGCCTTTCGTTATATCCGTTGTCTGCCGGTGAACGCTCTCCTGAGTAGGACAAATCCGCCGGGAGCGGATTTGAACGTTGCGAAGCAACgg                                           <  1:291328/140‑1 (MQ=11)
                                                                               aGGGAACTGCCAGGCATCAAATAAAACGAAAGGCTCAGTCGGAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGACAAATCCGCCGGGAGCGGATTTGAACGTTGCGAAGCAACGGccc                                        <  1:529190/143‑1 (MQ=32)
                                                                                 ggAACTGCCAGGCATCAAATAAACCGAAAGGCTCAGTCGGGAGACTGTGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGACAAATCCGCCGGGAGCGGATTTGAACGTTGCGAAGCAACGGCCCgg                                      <  2:276320/143‑1 (MQ=12)
                                                                                       gCCAGGCATCAAATAAAACGAAAGGCTCAGTCGGAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGCCAAATCCGCCGGGAGCGGATTTGAACGCTGCGAAGCAACGGCCCGGATGGTg                                >  1:237039/1‑143 (MQ=12)
                                                                                        ccAGGCATCAAATAAAACGAAAGGCTCAGTCGGAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGACAAATCCGCCGGGAGCGGATTTGAACGTTGCGAAGCAACGGCCCGGAGGGTgg                               >  1:478138/1‑143 (MQ=32)
                                                                                        ccAGGCATCAAATAAAACGAAAGGCTCAGTCGGAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGACAAATCCGCCGGGAGCGGATTTGAACGTTGCGAAGCAACGGCCCGGAGGGTgg                               <  1:745964/143‑1 (MQ=32)
                                                                                          aGGCATTAAATAAATCGAAAGGCTCAGTCGGAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGACAAATCCGCCGGGAGCGGATTTGAACGTTGCGAAGCAACGGCCCGGAGGGTGGCg                             <  1:569775/143‑1 (MQ=17)
                                                                                          aGGCATCAAATAAAACGAAAGGCTCAGTCGGAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGACAAATCCGCCGGGAGCGGATTTGAACGTTGCGAAGCAACGGCCCGGAGGGTGGCg                             <  1:683685/143‑1 (MQ=32)
                                                                                                cAACTAAACCGAAAGGCTCAGGCGGAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGCGTAGGACAAATCCGCCGGGAGCGGATTTGAACGTTGCGAAGCAACGGCCCGGAGGGTGGCGGGCAgg                       <  1:801528/143‑1 (MQ=255)
                                                                                                cAAATAAAACGAAAGGCTCAGTCGGAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGACAAATCCGCCGGGAGCGGATTTGAACGTTGCGAAGCAACGGCCCGCAGGGTGGCGGGCAgg                       >  1:426030/1‑143 (MQ=17)
                                                                                                 aaaTAAAACGAAAGGCTCAGGCGGAAGACTGGGCCGTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGACAAATCCGCCGGGAGCGGATTTGAACGTTGCGAAGCAACGGCCCGGAGGGTGGCGGGCAGGa                      <  1:268726/143‑1 (MQ=17)
                                                                                                   aTAAAACGAAAGGCTCAGTCGGAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGACAAATCCGCCGGGAGCGGATTTGAACGTTGCGAAGCAACGGCCCGGAGGGTGGCGGGCAGGAcg                    >  2:277715/1‑143 (MQ=32)
                                                                                                        ccGAAAGGCGCCGGCGGAAGACTGGGCCGTTCGTGTGATCTGTTCTTCCTCGGTGAACGCTCTCCTGAGTAGGACAAATCCGCCGGGAGCGGATTTGAACGTTGCGAAGCAACGGCCCGGAGGGTGGCGGGCAGGAcgcccgc               <  2:693737/142‑1 (MQ=11)
                                                                                                        aCGAACGGCTCGGTCGGAAGCCTGGGCCTTTCGTGTTATCTGTTGTTTGTCGGGGGAGGCGCTCGTGAGTAGGACAAATCCGGCGGGGGCGGATGTGAACGTTGCGAAGCAACGGCCCGGAGGGTGGCGGGCAGGAcgcccgc               <  2:681434/143‑1 (MQ=11)
                                                                                                           agaGGCTCAGTCTGAAGACTGGGTTTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGACAAATCCGCCGGGAGCGGATTTGAACGTTGCGAAGCAACGGCCCGGAGGGTGGCGGGCAGGCCGCCCGCCAt            <  1:561344/141‑1 (MQ=255)
                                                                                                             aGGCTCAGTCGGAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGACAAATCCGCCGGGAGCGGATTTGAACGTTGCGAAGCAACGGCCCGGAGGGTGGCGGGCAGGACGCCCGCCATaa          >  1:91763/1‑143 (MQ=32)
                                                                                                               gCTCAGTCGGAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGACCAATCCGCCGGGAGCGGATTTGAACGTTGCGAAGCAACGGCCCGGAGGGCGGCGGGCAGGACGCCCGCCATAAAc        >  2:136562/1‑143 (MQ=17)
                                                                                                               gCTCAGTCGGAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGACAAATCCGCCGGGAGCGGATTTGAACGTTGCGAAGCAACGGCCCGGAGGGTGGCGGGCAGGACGCCCGCCATAAAc        >  2:95002/1‑143 (MQ=32)
                                                                                                               gCTCAGTCGGAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGACAAATCCGCCGGGAGCGGATTTGAACGTTGCGAAGCAACGGCCCGGAGGGTGGCGGGCAGGACGCCCGCCATAAAc        <  1:699365/143‑1 (MQ=32)
                                                                                                               gCTCAGTCGGAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGACAAATCCGCCGGGAGCGGATTTGAACGTTGCGAAGCAACGGCCCGGAGGGTGGCGGGCAGGACGCCCGCCATAAAc        >  1:706589/1‑143 (MQ=32)
                                                                                                                     tCGGAGGACTGGGCGTGTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGACAAATCCGCCGGGAGCGGATTTGAACGTTGCGAAGCAACGGCCCGGAGGGTGGCGGGCAGGACGCCCGCCATAAACTGCCAg  <  1:620975/143‑1 (MQ=12)
                                                                                                                     tCGGAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGACAAATCCGCCGGGAGCGGATTTGAACGTTGCGAAGCAACGGCACGGAGGGTGGCGGGCAGGACGCCCGCCATAAACTGCCAg  >  1:229/1‑143 (MQ=17)
                                                                                                                                              |                                                                                                                     
TCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCTCATGCGAGAGTAGGGAACTGCCAGGCATCAAATAAAACGAAAGGCTCAGTCGGAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGACAAATCCGCCGGGAGCGGATTTGAACGTTGCGAAGCAACGGCCCGGAGGGTGGCGGGCAGGACGCCCGCCATAAACTGCCAG  >  NC_000913/4040541‑4040800

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.