New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913_3_pae_pgi 3423823 =NA (NA)16 (0.300) 10/270 4.8 29.5% intergenic (‑36/+57) rrfD/rrlD 5S ribosomal RNA of rrnD operon/23S ribosomal RNA of rrnD operon
?NC_000913_3_pae_pgi = 3946688 40 (0.710)intergenic (+81/‑12) rrlC/rrfC 23S ribosomal RNA of rrnC operon/5S ribosomal RNA of rrnC operon
Rejected: Coverage evenness skew score above cutoff.

GCCGGGTGTGTAAGCGCAGCGATGCGTTGAGCTAACCGGTACTAATGAACCGTGAGGCTTAACCTTACAACGCCGAAGATGTTTTGGCGGATGAGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGATGAAACAGAATTTGCCTGGCGGCCGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATTAAGCAGTAA  >  NC_000913_3_pae_pgi/3946542‑3946834
                                                                                                                                                  |                                                                                                                                                  
gCCGGGTGTGTAAGCGCAGCGATGCGTTGAGCTAACCGGTACTAATGAACCGTGAGGCTTAACCTTACAACGCCGAAGATGTTTTGGCGGATGAGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAACGCAGACGCGGTCTGATGaa                                                                                                                                                 >  2:40050/1‑150 (MQ=25)
 ccGGGTGTGTAAGCGCAGCGATGCGTTGAGCTAACCGGTACTAATGAACCGTGAGGCTTAACCTTACAACGCCGAAGATGTTTTGGCGGATGAGAGAAGATTTTCAGCCTGATACAGATTAAATCAGTACGCAGAAGCGGTCTGATGaaa                                                                                                                                                <  2:216330/150‑1 (MQ=21)
 ccGGGTGTGTAAGCGCAGCGATGCGTTGAGCTAACCGGTACTAATGAACCGTGAGGCTTAACCTTACAACGCCGAAGATGTTTTGGCGGATGAGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGATGaaa                                                                                                                                                >  2:893901/1‑150 (MQ=34)
          tAAGCGCAGCGATGCGTTGAGCTAACCGGTACTAATGAACCGTGAGGCTTAACCTTACAACGCCGAAGATGTTTTGGCGGATGAGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGATGAAACAGAATTTg                                                                                                                                       >  1:874708/1‑150 (MQ=34)
             gcgcAGCGATGCGTTGAGCTAACCGGTACTAATGAACCGTGAGGCTTAACCTTACAACGCCGAAGATGTTTTGGCGGATGAGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGATGAAACAGAATTTGCCt                                                                                                                                    <  2:617521/150‑1 (MQ=35)
                       gCGTTGAGCTAACCGGTACTAATGAACCGTGAGGCTTAACCTTACAACGCCGAAGATGTTTTGGCGGATGAGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGATGAAACAGAATTTGCCTGGCGGCCGTa                                                                                                                          <  1:679692/150‑1 (MQ=35)
                          ttGAGCTAACCGGTACTAATGAACCGTGAGGCTTAACCTTACAACGCCGAAGATGTTTTGGCGGATGAGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCtgat                                                                                                                                                    <  2:709371/121‑1 (MQ=32)
                          ttGAGCTAACCGGTACTAATGAACCGTGAGGCTTAACCTTACAACGCCGAAGATGTTTTGGCGGATGAGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCtgat                                                                                                                                                    >  1:709371/1‑121 (MQ=33)
                             aGCTAACCGGTACTAATGAACCGTGAGGCTTAACCTTACAACGCCGAAGATGTTTTGGCGGATGAGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGATGAAACAGAATTTGCCTGGCGGCCGTAGCGCgg                                                                                                                    <  1:292388/150‑1 (MQ=35)
                              gCTAACCGGTACTAATGAACCGTGAGGCTTAACCTTACAACGCCGAAGATGTTTTGGCGGATGAGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGAGGAAACAGAATTTGCCTGGCGGCCGTAGCGCggg                                                                                                                   >  1:294096/1‑149 (MQ=21)
                                    cGGTACTAATGAACCGTGAGGCTTAACCTTACAACGCCGAAGATGTTTTGGCGGATGAGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGATGAAACAGAATTTGCCTGGCGGCCTGAGCGCGGTGGTccc                                                                                                             >  1:120704/1‑150 (MQ=25)
                                    cGGTACTAATGAACCGTGAGGCTTAACCTTACAACGCCGAAGATGTTTTGGCGGATGAGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGATGAAACAGAATTTGCCTGGCGGCCGTAGCGCGGTGGTccc                                                                                                             >  2:818705/1‑150 (MQ=35)
                                     ggTACTAATGAACCGTGAGGCTTAACCTTACAACGCCGAAGATGTTTTGGCGGATGAGAGAAGATTTTCAGCCTGATACAGAGTAATTCAGAACGCAGAAGCGGTCTGATGAAACAGAATTTGCCTGGCGGCCGTAGCGCGGTGGTCCCa                                                                                                            <  1:526206/150‑1 (MQ=21)
                                       aaCTAATGAACCGTGAGGCTTAACCTTACAACGCCGAAGATGTTTTGGCGGATGAGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGATGAAACAGAATTTGCCTGGCGGCCGTAGCGCGGTGGTCCCAcc                                                                                                          <  1:593487/149‑1 (MQ=25)
                                               aaCCGTGAGGCTTAACCTTACAACGCCGAAGATGTTTTGGCGGATGAGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGATGAAACAGAATTTGCCTGGCGGCCGTAGCGCGGTGGTCCCACCTGACCCCa                                                                                                  >  2:952838/1‑150 (MQ=35)
                                                actgTGAGGCTTAACCTTACAACGCCGAAGATGTTTTGGCGGATGAGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGATAAAACAGAATTTGCCTGGCGGCCGTAGCGCGGTGATCCCACCTGACCCCAt                                                                                                 <  2:57375/147‑1 (MQ=0)
                                                           tAACCTTACAACGCCGAAGATGTTTTGGCGGATGAGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGATGAAACAGAATTTGCCTGGCGGCCGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAg                                                                                      <  1:893901/150‑1 (MQ=35)
                                                             aCCTTACAACGCCGAAGATGTTTTGGCGGATGAGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGATGAAACAGAATTTGCCTGGCGGCCGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGaa                                                                                    <  1:616272/150‑1 (MQ=35)
                                                               cTTACAACGCCGAAGATGTTTTGGCGGATGAGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGATGAAACAGAATTTGCCTGGCGGCCGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGt                                                                                  >  1:705686/1‑150 (MQ=35)
                                                                          gAAGATGTTTTGGCGGATGAGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGATGAAACAGAATTTGCCTGGCGGCCGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGaa                                                                                    >  2:17402/1‑137 (MQ=37)
                                                                          gAAGATGTTTTGGCGGATGAGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGATGAAACAGAATTTGCCTGGCGGCCGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGaa                                                                                    <  1:17402/137‑1 (MQ=35)
                                                                               tGTTTTGGCGGATGAGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGATGAAACAGAATTTGCCTGGCGGCCGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGcc                                                                  <  1:945187/150‑1 (MQ=32)
                                                                                   ttGGCGGATGAGAGCAGATTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGATGAAACAGAATTTGCCTGGCGGCCGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATg                                                              <  1:333226/150‑1 (MQ=17)
                                                                                      ggttaTGAGCGAAGATTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGATGAAACAGAATTTGCCTGGCGGCCGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTa                                                           <  1:471285/146‑1 (MQ=17)
                                                                                             agagaAGATTTTCATCCTGATACAGATTAAACCAGAACGCAGAAGCGGTCTGATGAAACAGAATTTGCCTGGCGGCCGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTggg                                                    <  2:7619/150‑1 (MQ=17)
                                                                                              gagaAGATTTTCAACCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGATGAAACAGAATTTGCCTGGCGGCCGTAg                                                                                                                         <  1:132874/80‑1 (MQ=255)
                                                                                              gagaAGATTTTCAACCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGATGAAACAGAATTTGCCTGGCGGCCGTAg                                                                                                                         >  2:132874/1‑80 (MQ=255)
                                                                                                gaagaTTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGATGAAACAGAATTTGCCTGGCGGCCGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGtc                                                 >  2:583252/1‑150 (MQ=32)
                                                                                                 aagaTTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGATGAAACAGAATTTGCCTGGCGGCCGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGtgtg                                                      <  2:859522/144‑1 (MQ=32)
                                                                                                 aagaTTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGATGAAACAGAATTTGCCTGGCGGCCGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGtct                                                >  1:564822/1‑150 (MQ=32)
                                                                                                 aagaTTTTCAGCCTGAAACAGATTAAATCAGAACGCAGAAGCGGTCTGATGAAACAGAATTTGCCTGGCGGCCGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGtgtg                                                      >  1:859522/1‑144 (MQ=17)
                                                                                                   gaTGTTCAGGCTGAGACAGATTAAATCAGAACGCAGAAGCGGTCTGATGAAACAGAATTTGCCTGGCGGCCGTAGCGCGGTGGGCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTcc                                              <  2:403183/150‑1 (MQ=255)
                                                                                                       ttCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGATGAAACAGAATTTGCCTGGCGGCCGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCCCAt                                          <  1:513199/150‑1 (MQ=32)
                                                                                                           gCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGATGAAACAGAATTTGCCTGGCGGCCGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCCCATGCga                                      <  1:934778/150‑1 (MQ=32)
                                                                                                              tGATACAGATTAAATCAGAACGCAGAAGCGGTCTGATGAAACAGAATTTGCCTGGCGGCCGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCCCATGCGagag                                   <  1:642134/150‑1 (MQ=32)
                                                                                                                 aaCAGATTAAATCAGAACGCAGAAGCGGTCTGATGAAACAGAATTTGCCTGGCGGCCGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCCCATGCGAGAGTAg                                <  2:835014/149‑1 (MQ=17)
                                                                                                                  aCAGATTAAATCAGAACGCAGAAGCGGTCTGATGAAACAGAATTTGCCTGGCGGCCGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCCCATGCGAGAGTAgg                               <  2:874708/150‑1 (MQ=32)
                                                                                                                    aGATTAAATCAGAACGCAGAAGCGGTCTGATGAAACAGAATTTGCCTGGCGGCCGTAGCGCGGTGGTCCCACCTGACCCCATTCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCCCATGCGAGAGTAGGGa                             <  1:952838/150‑1 (MQ=17)
                                                                                                                      aTTAAATCAGAACGCAGAAGCGGTCTGATGAAACAGAATTTGCCTGGCGGCCGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGCTAGTGTGGGGTCTCCCCGTGCGAGCGTCGGTAAc                           >  1:412451/1‑150 (MQ=11)
                                                                                                                          aaTCAGAACGCAGAAGCGGTCTGATGAAACAGAATTTGCCTGGCGGCCGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGcc                       >  2:842394/1‑150 (MQ=32)
                                                                                                                           aTCAGAACGCAGAAGCGGTCTGATGAAACAGAATTTGCCTGGCGGCCGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCa                      <  1:89460/150‑1 (MQ=32)
                                                                                                                                               tgatgaAACAGAATTTGCCTGGCGGCCGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCGGGCATGCAAGTAAGCGGtta  >  2:439759/1‑148 (MQ=11)
                                                                                                                                                  |                                                                                                                                                  
GCCGGGTGTGTAAGCGCAGCGATGCGTTGAGCTAACCGGTACTAATGAACCGTGAGGCTTAACCTTACAACGCCGAAGATGTTTTGGCGGATGAGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGATGAAACAGAATTTGCCTGGCGGCCGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATTAAGCAGTAA  >  NC_000913_3_pae_pgi/3946542‑3946834

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