New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913 3423821 =NA (NA)6 (0.140) 4/272 6.0 11.3% intergenic (‑36/+57) rrfD/rrlD 5S ribosomal RNA of rrnD operon/23S ribosomal RNA of rrnD operon
?NC_000913 = 3946687 49 (1.100)intergenic (+81/‑12) rrlC/rrfC 23S ribosomal RNA of rrnC operon/5S ribosomal RNA of rrnC operon
Rejected: Coverage evenness skew score above cutoff.
Rejected: Frequency below/above cutoff threshold.

TAAGCGCAGCGATGCGTTGAGCTAACCGGTACTAATGAACCGTGAGGCTTAACCTTACAACGCCGAAGATGTTTTGGCGGATGAGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGATGAAACAGAATTTGCCTGGCGGCCGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATTAAGC  >  NC_000913/3946551‑3946828
                                                                                                                                        |                                                                                                                                             
tAAGCGCAGCGATGCGTTGAGCTAACCGGTACTAATGAACCGTGAGGCTTAACCTTACAACGCCGAAGATGTTTTGGCGGATGAGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGATGAAACa                                                                                                                                         >  1:641628/1‑143 (MQ=34)
 aaGCGCAGCGATGCGTTGAGCTAACCGGTACTAATGAACCGTGAGGCTTAACCTTACAACGCCGAAAGTGTTTTGGCGGATGAGAGAAGATTTTCAGCCTGATACAGATTAAATCCGAACGCAGAAGCGGTCTGATGACACAg                                                                                                                                        >  2:310993/1‑143 (MQ=255)
  aGCGGAGCGATGCGTGGAGCCAACCGGGACTAAGGAACCGTGAGGCTTAACCTTACAACGCCGAAGATGTTTTGGCGGATGAGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGATGAAACAGa                                                                                                                                       <  2:777970/143‑1 (MQ=11)
        ggggtGCGTTGAGCTAACGGGTTCTATTGGTCCGCGAGGATTAACGTTACAACGCCGAAGATGTTTTGGCGGATGAGAGAAGATTTTCAGCCTGATACAGATTAAAGCAGAACGCAGAAGCGGTCTGATGAAACAGAATTTGc                                                                                                                                 <  2:597066/139‑1 (MQ=11)
         cGATGCGTTGAGCTAACCGGTACTAATGAACCGTGAGGCTTAACCTTACAACGCCGAAGATGTTTTGGCGGATGAGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAACGGTCTGATGAAACAGAATTTGcc                                                                                                                                >  2:504225/1‑143 (MQ=21)
           aTGCGTTGAGCTAACCGGTACTAATGAACCGTGAGGCTTGACCTTACAACGCCGAAGATGTTTTGGCGGATGAGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGGCTGATGAAACAGAATTTGCCTg                                                                                                                              >  1:743399/1‑143 (MQ=21)
             gCGTTGGACTAACCGGTACTAATTAGCCGTGAGACTTAACCTTACAACGCCGAAGATGTTTTGGCGGATGAGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGATGAAACAGAATTTGCCTggc                                                                                                                            <  2:353114/143‑1 (MQ=11)
                ttGAGCTAACCGGTACTAATGAACCGTGAGGCTTAACCTTACAACGCCGAAGATGTTTTGGCGGATGAGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGATGAAACAGAATTTGCCTggcggc                                                                                                                         >  1:111512/1‑143 (MQ=34)
                  gAGCTAACCGGTACTAATGTACCGTGAGGCTTATCCTTACAACGCCGAAGATGTTTTGGCGGATGAGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGATGAAACAGAATTTGCCTGGCGGCCg                                                                                                                       <  1:509926/143‑1 (MQ=21)
                                tAATGAACCGTGAGGCTTAACCTTACAACGCCGAAGATGTTTTGGCGGATGAGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGATGACACAGAATTTGCCTGGCGGCCGTAGCGCGGTGGTcc                                                                                                         >  1:583657/1‑143 (MQ=25)
                                 aaTGAACCGTGAGGCTTAACCTTACAACGCCGAAGATGTTTTGGCGGATGAGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGATGAAACAGAATTTGCCTGGCGGCCGTAGCGCGGCGGGccc                                                                                                        >  1:602432/1‑143 (MQ=25)
                                  aTGAACCGTGAGGCTTAACCTTACAACGCCGAAGATGTTTTGGCGGATGAGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGATGAAACAGAATTTGCCTGGCGGCCGTAGCGCGGTGGTCCCa                                                                                                       >  1:166016/1‑143 (MQ=35)
                                        cGTGAGGCTTAACCTTACAACGCCGGAGATGGTGTGGCGGATGAGAGAAGATTTTCAGCCTGATACAGAGTAAATCAGAGCGCAGAAGCGGTCTGATGAAACAGAATTTGCCTGGCGGCCGTAGCGCGGTGGTCCCACCTGAc                                                                                                 <  2:200940/143‑1 (MQ=11)
                                            aGGCTTAACCTTACAACGCCGAAGATGTTTTGGCGGATGAGAGAAGATTTTCAGCCTGATACAGGTTAAATCAGAACGCAGAAGCGGTCTGATGAAACAGAAATTCTCGGTCTTCGGTAGCGCGGTGGTCCCACCTGACCCCa                                                                                             >  2:549382/1‑143 (MQ=11)
                                            aGGCTTAACCTTACAACGCCGAAGATGTTTTGGCGGATGAGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGATGAAACAGAATTTGCCTGGCGGCCGTAGCGCGGTGGTCCCACCTGACCCCa                                                                                             >  2:595831/1‑143 (MQ=35)
                                              gCTCAACCGTACAACGCCGAAGATGTTTTGGCGGATGAGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGATGAAACAGAATTTGCCTGGCGGCCGTAGCGCGGTGGTCCCACCTGACCCCATg                                                                                           <  1:271858/143‑1 (MQ=25)
                                               cgagaccgtagagcGCCGAAGAGGTTGTGGCGGATGAGAGAAGATTTTCAGCCTGATACAGATTAAAGCAGAGCGCAGAAGCGGTGTGATGAAACAGAATTTGCCTGGCGGCCGTAGCGCGGTGGTCCCACCTGACCCCATGc                                                                                          <  1:345757/130‑1 (MQ=11)
                                                ttAACCTTACAACGCCGAAGATGTTTTGGCGGATGAGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGGCTGATGAAACAGAACTTGCCTGGCGGCCGTAGCACGGTGGTCCCACCTGACCCCATGca                                                                                         >  2:501595/1‑142 (MQ=255)
                                                   aCCTTACAACGCCGAAGATGTTTTGGCGGATGAGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGATGAAACAGAATTTGCCTGGCGGCCGTAGCGCGGTGGTCCCACCTGACCCCATGCCgcc                                                                                      >  2:103631/1‑141 (MQ=25)
                                                       tACAACGCCGAAGATGTTTTGGCGGATGAGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGATGAAACAGAATTTGCCTGGCGGCCGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCa                                                                                  <  1:164155/143‑1 (MQ=35)
                                                         cAACGCCGAAGATGTTTTGGCGGATGAGAGAAGATTTCCAGCCTGATACGGATTAAATCAGAACGCAGAAGCGGTCTGATGAAACAGAATTTGCCTGGCGGCCGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGa                                                                                <  2:673790/143‑1 (MQ=21)
                                                          aaCGGCGTAGATGGTTTGGCGGATGAGAGAAGATTTTCAGCCTGATACATATTAAATCAGAACGCAGAAGCGGTCTGATGAAACAGAATTTGCCTGGCGGCCGTAGCGCGGTTGTCCCACCTGACCCCATGCCGAACTCAGaa                                                                               <  2:595749/143‑1 (MQ=11)
                                                          aaCGCTGAAGATGTTTTGGCGGATGAGAGAAGATTTGCGGCCGGATACAGATTGAATCAGAACGCAGAGGCGGTCTGATGAAACAGAATTTGCCTGGCGGCCGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGaa                                                                               <  2:397890/143‑1 (MQ=11)
                                                              cccaaGATGTTTTGACGGAGGAGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGATGAAACAGAATTTGCCTGGCGGCCGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGa                                                                           <  1:183622/140‑1 (MQ=255)
                                                              ccGAAGATGTTTTGGCGGATGAGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGATGAAACAGAATTTGCCTGGCGGCCGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGa                                                                           <  1:33461/143‑1 (MQ=35)
                                                                gAAGATGTTTTGGCGGATGAGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGATGAAACAGAATTTGCCTGGCGGCCGTAGCGCGGTGGACCCGCCTGACCCCATGCCGAACCAAGAAGTGcaa                                                                         >  2:628929/1‑143 (MQ=11)
                                                                 aaGATGTTTTGGCGGATGAGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGATGAAACAGAATTTGCCTGGCGGCCGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAAc                                                                        >  1:551264/1‑143 (MQ=35)
                                                                 aaGATGTTTTGGCGGATGAGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGATGAAACAGAATTTGCCTGGCGGCCGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGAGAAAc                                                                        >  1:389887/1‑143 (MQ=21)
                                                                           ggCGGATGCGAGAAAGTTTTATGCCGGAGAGAGATTAAATCAGAACGCAGAAGCGGTCTGATGAAACAGAATTTGCCTGGCGGCCGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAgcgc                                                              <  2:546993/143‑1 (MQ=11)
                                                                                 tGAGAGAAGTTGTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCCGATGAAACAGAATTTGCCTGGCTGCCGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATgg                                                        <  2:391437/143‑1 (MQ=255)
                                                                                  gagagaAGATTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGATGAAACAGAATTTGCCTGGCGGCCGTAGCGCGGTGGTCCCACCTGACCCCATGCCGCACTCACACGTGAAACGCCGTAGCGCCGATGGt                                                       >  2:447729/1‑143 (MQ=12)
                                                                                     agaagaTTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGATGAAACAGAATTTGCCTGGCGGCCGTAGCGCGGTGGTACCACCTGACCCCATGCCGAACTCAGAACGGAACAGCGGTAGCGCCGATGGtggc                                                    >  2:80122/1‑140 (MQ=11)
                                                                                                  ccTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGATGAAACAGAATTTGCCTGGCGGCCGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGTATCGTTGTTTGTGTggggaaggccc                                       >  1:565076/1‑136 (MQ=11)
                                                                                                      aTACAGATTAAATCAGAACGCAGAAGCGGTCTGATGAAACAGAATTTGCCTGGCGGCCGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCCCATGc                                   >  2:212905/1‑143 (MQ=32)
                                                                                                        aCAGATTAGATCAGAAGGCAGAAGCGGTCTGATGAAACAGAATTTGCCTGGGGGCCGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCAGTAGCGCCGATGGGAGTGCGGGGTCTCACCAGGCgc                                 >  2:537246/1‑142 (MQ=11)
                                                                                                        aCAGATTAAATCAGAACGCACAAGCGGTCTGAAGAAACAGAATTTGCCTGGCGGCCGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGCGGGGTCTCCCCATGCga                                 >  2:407936/1‑143 (MQ=12)
                                                                                                          cgATTAAATCAGAACGCAGAAGCGGTCTGATGAAACAGAATTTGCCTGGCGGCCGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCCCATGCgaga                               <  1:335218/142‑1 (MQ=17)
                                                                                                           gagtAAATCAGAACGCAGAAGCGGTCTGATGAAACAGAATTGGCCTGGCGGCCGTAGCGCGGTGGTGGCTCCGTACCCCATGCCGAACGCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCCCATGCGagag                              <  1:572865/140‑1 (MQ=11)
                                                                                                              tAAATCAGAACGCAGAAGCGGTCTGATGAAACAGAATTTGCCTGGCGGCCGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCCCATGCGAGAGTAg                           <  1:595831/143‑1 (MQ=32)
                                                                                                                  tCAGAACGCAGAGGCTGTCTGCTGAAACGCAACTTGCCGCGCGGCCGGAGCGCGGTGGTCGCACCGGGCCCCATGCCGAACTCAGAAGTGAAACTCCTTAGCGCCGCTGGTAGTGTGGGGTCTCCCCATGCGAGAGTAGGGaa                       <  2:427074/143‑1 (MQ=11)
                                                                                                                  tCAGAACGCAGAAGCGGTCTGATGAAACAGAATTTGCCTGGCGGCCGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCCCATGCGAGAGTAGGGaa                       >  2:739879/1‑143 (MQ=32)
                                                                                                                      gacggggAAGCGGGCTGATGAAACAGGGTGTGCCGGGCGGGCGGAGCGCGGTGGTCCCGCCGGGCGCCATGCAGGACTCCGAAGTGAAGCGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGc                   <  2:11253/137‑1 (MQ=11)
                                                                                                                        cGCAGAAGCGGTCTGAGGAAACAGAATTTGCGTGGCGGCCGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCa                 <  1:125715/143‑1 (MQ=17)
                                                                                                                         gCAAAAGCGGTCCGATGAAACAGAATGTGGCTGGCGGCCGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAg                <  1:274443/143‑1 (MQ=255)
                                                                                                                            gAATCGGTCTGATGAAACAGAATTTGCCTGGCGGCCGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCCAATGCGAGGGGAGGGAAGTGccggccg             >  2:641782/1‑138 (MQ=11)
                                                                                                                             aggCGGTCTGATGAAACAGAATTTGCCGGGCGGCCGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCAt            <  2:295698/141‑1 (MQ=17)
                                                                                                                                cGGTCGGAGGAACCAGAACTTCCCTGCCGGCCGTAGCGCTGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCaa         <  2:323267/143‑1 (MQ=11)
                                                                                                                                     tgatgaAACAGAATTTGCCTGGCGGCCGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATTaa    <  1:631221/143‑1 (MQ=32)
                                                                                                                                       atgaAACAGAATTTGCCTGGCGGCCGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATTAAGc  >  1:303029/1‑143 (MQ=32)
                                                                                                                                        |                                                                                                                                             
TAAGCGCAGCGATGCGTTGAGCTAACCGGTACTAATGAACCGTGAGGCTTAACCTTACAACGCCGAAGATGTTTTGGCGGATGAGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGATGAAACAGAATTTGCCTGGCGGCCGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATTAAGC  >  NC_000913/3946551‑3946828

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