New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913_3_pae_pgi = 4040816NA (NA)9 (0.180) 7/262 6.0 14.3% intergenic (+180/+90) rrfA/mobB 5S ribosomal RNA of rrnA operon/molybdopterin‑guanine dinucleotide biosynthesis protein B
?NC_000913_3_pae_pgi = 4040840 58 (1.050)intergenic (+204/+66) 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.

CGCTCTCCTGAGTAGGACAAATCCGCCGGGAGCGGATTTGAACGTTGCGAAGCAACGGCCCGGAGGGTGGCGGGCAGGACGCCCGCCATAAACTGCCAGGCATCAAATTAAGCAGAAGGCCATCCTGACGGATGGCCTTTTTGCATTGGCGCAGAAAAAAATGCCGGATGCGACGCTGGCGCGTCTTATCCAGCCTACTCTTGCTTATCCGTTCTGCTTTTGCATCCACTCCACCACAAAATCAGCCAGCCCCTCAACATCA  >  NC_000913_3_pae_pgi/4040702‑4040963
                                                                                                                  |                                                                                                                                                   
cGCTCTCCTGAGTAGGACAAATCCGCCGGGAGCGGATTTGAACGTTGCGAAGCAACGGCCCGGAGGGTGGCGGGCAGGACGCCCGCCATAAACTGCCAGGCATCAAATTAAGCAGAAGGCCAGCCTGACGGATGGCCTTTTTGAATTgt                                                                                                                   >  2:865471/1‑148 (MQ=14)
 gCTCTCCTGAGTAGGACAAATCCGCCGGGAGCGGATTTGAACGTTGCGAAGCAACGGCCCGGAGGGTGGCGGGCAGGACGCCCGCCATAAACTGCCAGGCATCAAATTAAGCAGAAGGCCATCCTGACGGATGGCCTTTTTGCATTGgc                                                                                                                  >  1:496856/1‑149 (MQ=34)
   tctcCTGAGTAGGACAAATCCGCCGGGAGCGGATTTGAACGTTGCGAAGCAACGGCCCGGAGGGTGGCGGGCAGGACGCCCGCCATAAACTGCCAGGCATCAAATTAAGCAGAAGGCCATCCTGACGGATGGCCTTTTTGCATTGgcgc                                                                                                                >  2:69524/1‑149 (MQ=35)
    ctcCTGAGTAGGACAAATCCGCCGGGAGCGGATTTGAACGTTGCGAAGCAACGGCCCGGAGGGTGGCGGGCAGGACGCCCGCCATAAACTGCCAGGCACCAAATTAAGCAGAAGGCCATCCTGACGGATGGCCTTTTTGCATTGGCGCa                                                                                                               <  1:551423/149‑1 (MQ=37)
         gAGTAGGACAAATCCGCCGGGAGCGGATTTGAACGTTGCGAAGCAACGGCCCGGAGGGTGGCGGGCAGGACGCCCGCCATAAACTGCCAGGCATCAAATTAAGCAGAAGGCCATCCTGACGGATGGCCTTTTTGCATTGGCGCAGaaaa                                                                                                          <  1:285389/149‑1 (MQ=40)
           gTAGGACAAATCCGCCGGGAGCGGATTTGAACGTTGCGAAGCAACGGCCCGGAGGGTGGCGGGCCGGACGCCCGCCATAAACTGCCAGGCATCAAATTCACTCAAATGTCATCCTGACGGATGGCCTTTTTGCATTGGCGCAGcaaaaa                                                                                                        >  1:564190/1‑149 (MQ=14)
            tAGGACAAATCCGCCGGGAGCGCATTTGAACGTTGCGAACCAACGGCCCGGAGGGTGGCGGGCAGGACGCCCGCCATACACTGCCAGGCATCAAATTAAGCAGAAGGCCATCCTGACGGATGGCCTTTTTGCATTGGCGCAGaaaaaaa                                                                                                       <  1:231819/149‑1 (MQ=255)
                              aGCGGATTTGAACGTTGCGAAGCAACGGCCCGGAGGGTGGCGGGCAGGACGCCCGCCATAAACTGCCAGGCATCAAATTAAGCAGAAGGCCATCCTGACGGATGGCCTTTTTGCATTGGCGCAGAAAAAAATGCCGGATGCGACGCTgg                                                                                     <  1:207730/149‑1 (MQ=255)
                                   aTTTGAACGTTGCGAAGCAACGGCCCGGAGGGTGGCGGGCAGGACGCCCGCCATAAACTGCCAGGCATCAAATTAAGCAGAAGGCCATCCTGACGGATGGCCTTTTTGCATTGGCGCAGAAAAAAATGCCGGATGCGACGCTGACGCGt                                                                                >  2:888952/1‑149 (MQ=255)
                                    tttGAACGTTTCGAAGCAACGGCCCGGAGGGTGGCGGGCAGGACGCCCGCCATAAACTGCCAGGCATCAAATTAAGCAGAAGGCCATCCTGACGGATGGCCTTTTTGCATTGGCGCAGAAAAAAATGCCGGATGCGACGCTGGCGCGTc                                                                               <  2:30831/149‑1 (MQ=255)
                                        aaCGTTGCGAAGCAACGGCCCGGAGGGTGGCGGGCAGGACGCCCGCCATAAACTGCCAGGCATCAAATTAAGCAGAAGGCCATCCTGACGGATGGCCTTTTTGCATTGGCGCAGAAAAAAATGCCGGATGCGACGCTGGCGCGTCTTAt                                                                           >  2:454727/1‑149 (MQ=255)
                                           gTTGCGAAGCAACGGCCCGGAGGGTGGCGGGCAGGACGCCCGCCATAAACTGCCAGGCATCAAATTAAGCAGAAGGCCATCCTGACGGATGGCCTTTTTGCATTGGCGCAGAAAAAAATGCCGGATGCGACGCTGGCGCGTCTTATCCa                                                                        <  2:399348/149‑1 (MQ=255)
                                           gTTGCGAACCAACGGCCCGGAGGGTGGCGGGCAGGACGCCCGCCATAAACTGCCAGGCATCAAATTAAGCAGAAGGCCATCCTGACGGATGGCCTTTTTGCATTGGCGCAGAAAAAAATGCCGGATGCGACGCTGGCGCGTCTTATCCa                                                                        <  2:635306/149‑1 (MQ=255)
                                             tGCGAAGCAACGGCCCGGAGGGTGGCGGGCAGGACGCCCGCCATAAACTGCCAGGCATCAAATTAAGCAGAAGGCCATCCTGACGGATGGCCTTTTTGCATTGGCGCAGAAAAAAATGCCGGATGCGACGCTGGCGCGTCTTATCCAGc                                                                      <  2:501488/149‑1 (MQ=255)
                                                    cAACGGCCCGGAGGGTGGCGGGCAGGACGCCCGCCATAAACTGCCAGGCATCAAATTAAGCAGAAGGCCATCCTGACGGATGGCCTTTTTGCATTGGCGCAGAAAAAAATGCCGGATGCGACGCTGGCGCGTCTTATCCAGCCTActct                                                               <  1:69524/149‑1 (MQ=255)
                                                     aaCGGCCCGGAGGGTGGCGGGCAGGACGCCCGCCATAAACTGCCAGGCATCAAATTAAGCAGAAGGCCATCCTGACGGATGGCCTTTTTGCATTGGCGCAGAAAAAAATGCCGGATGCGACGCTGGCGCGTCTTATCCAGCCTACTCtt                                                              <  1:437906/149‑1 (MQ=255)
                                                      aCGGCCCGGAGGGTGGCGGGCAGGACGCCCGCCATAAACTGCCAGGCATCAAATTAAGCAGAAGGCCATCCTGACGGATGGCCTTTTTGCATTGGCGCAGAAAAAAATGCCGGATGCGACGCTGGCGCGTCTTATCCAGCCTACTCTTg                                                             >  2:300203/1‑149 (MQ=255)
                                                      aCGGCCCGGAGGGTGGCGGGCAGGACGCCCGCCATAAACTGCCAGGCATCAAATTAAGCAGAAGGCCATCCTGACGGATGGCCTTTTTGCATTGGCGCAGAAAAAAATGCCGGATGCGACGCTGGCGCGTCTTATCCAGCCTACTCTTg                                                             >  1:126499/1‑149 (MQ=255)
                                                       cGGCCCGGAGGGTGGCGGGCAGGACGCCCGCCATAAACTGCCAGGCATCAAATTAAGCAGAAGGCCATCCTGACGGATGGCCTTTTTGCATTGGCGCAGAAAAAAATGCCGGATGCGACGCTGGCGCGTCTTATCCAGCCTACTCTTGc                                                            <  2:126499/149‑1 (MQ=255)
                                                          cccGGAGGGTGGCGGGCAGGACGCCCGCCATAAACTGCCAGGCATCAAATTAAGCAGAAGGCCATCCTGACGGATGGCCTTTTTGCATTGGCGCAGAAAAAAATGCCGGATGCGACGCTGGCGCGTCTTATCCAGCCTACTCTTGCTTa                                                         <  1:584981/149‑1 (MQ=255)
                                                                 ggTGGCGGGCAGGACGCCCGCCATAAACTGCCAGGCATCAAATTAAGCAGAAGGCCATCCTGACGGATGGCCTTTTTGCATTGGCGCAGAAAAAAATGCCGGATGCGACGCTGGCGCGTCTTATCCAGCCTACTCTTGCTTATCCGTTc                                                  <  1:395243/149‑1 (MQ=255)
                                                                       gggCAGGACGCCCGCCATAAACTGCCAGGCATCAAATTAAGCAGAAGGCCATCCTGACGGATGGCCTTTTTGCATTGGCGCAGAAAAAAATGCCGGATGCGACGCTGGCGCGTCTTATCCAGCCTACTCTTGCTTATCCGTTCTGCttt                                            <  1:250147/149‑1 (MQ=255)
                                                                       gggCAGGACGCCCGCCATAAACTGCCAGGCATCAAATTAAGCAGAAGGCCATCCTGACGGATGGCCTTTTTGCATTGGCGAAGAAAAAAATGCCGGATGCGACGCTGGCGCGTCTTATCCAGCCTACTCTTGCTTATCCGTTCTGCttt                                            <  2:333522/149‑1 (MQ=255)
                                                                        ggCAGGCCTCCCGCCATAAACTGCCAGGCATCAAATTAAGCAGAAGGCCATCCTGACGGATGGCCTTTTTGCATTGGTGTAGAAAAAAATGCCGGATGCGACGCTGGCGCGTCTTATCCAGCCTACTCTTGCTTATCCGTTCTGCtttt                                           <  1:631566/149‑1 (MQ=255)
                                                                             taCGCCCGCCATAAACTGCCAGGCATCAAATTAAGCAGAAGGCCATCCTGACGGATGGCCTTTTTGCATTGGCGCAGAAAAAAATGCCGGATGCGACGCTGGCGCGTCTTATCCAGCCTACTCTTGCTTATCCGTTCTGCTTTTGCAtc                                      >  1:199423/2‑149 (MQ=255)
                                                                             gACGCCCGCCATAAACTGCCAGGCATCAAATTAAGCAGAAGGCCATCCTGACGGATGGCCTTTTTGCATTGGCGCAGAAAAAAATGCCGGATGCGACGCTGGCGCGTCTTATCCAGCCTACTCTTGCTTATCCGTTCTGCTTTTGCAtc                                      <  1:300203/149‑1 (MQ=255)
                                                                               tgcccgccAAAAACGGCCAGGGATCAAATTAAGCAGAAGGGCATCCTGACGGATGGCTTTTTTGCATTGGCGCAGAAAAAAATGCCGGATGCGACGCTGGCGCGTCTTATCCAGCCTACTCTTGCTTATCCGTTCTGCTTTTGCAtcca                                    <  2:763215/148‑1 (MQ=25)
                                                                                 cccgccATAAACTGCCAGGCATCAAATTAAGCAGAAGGCCATCCTGACGGATGGCCTTTTTGCATTGGCGCAGAAAAAAATGCCGGATGCGACGCTGGCGCGTCTTATCCAGCCTCCTCTTGCTTATCCCTTCTGCTTTTGCAtccact                                  >  2:142267/1‑149 (MQ=255)
                                                                                  gcgcgATAACGTGCCAGGCATCAAATGAAGAAGAAGGCCATCCTGACGGATGGCTTTTTTGCATTGGCGCAGAAAAAAACGCCGGATGCGACGCTGGCGCGTCTTATCCAGCCTACTCTTGCTTATCCGTTCTGCTTTTGCAtccactc                                 <  2:302950/148‑1 (MQ=25)
                                                                                    gccATAAACTGCCAGGCATCAAATTAAGCAGAAGGCCATCCTGACGGATGGCCTTTTTGCATTGGCGCAGAAAAAAATGCCGGATGCGACGCTGGCGCGTCTTATCCAGCCTACTCTTGCTTATCCGTTCTGCTTTTGCAtccactcca                               >  1:609771/1‑149 (MQ=255)
                                                                                     ccATAAACTGCCAGGCATCAAATTAGGCAGAAGGCCATCCTGACGGATGGCCTTTTTGCATTGGCGCAGAAAAAAATGCCGGATGCAACGCTGGCGCGTCTTATCCAGCCTACTCTTGCTTATCCGTTCTGCTTTTGCAtccactccac                              <  1:953426/149‑1 (MQ=255)
                                                                                      cATAAACTGCCAGGCATCAAATTAAGCAGAAGGCCATCCTGACGGATGGCCTTTTTGCATTGGCGCAGAAAAAAATGCCGGATGCGACGc                                                                                        <  1:845256/90‑1 (MQ=255)
                                                                                      cATAAACTGCCAGGCATCAAATTAAGCAGAAGGCCATCCTGACGGATGGCCTTTTTGCATTGGCGCAGAAAAAAATGCCGGATGCGACGc                                                                                        >  2:845256/1‑90 (MQ=255)
                                                                                           aCTGCCAGGCATCAAATTAAGCAGAAGGCCATCCTGACGGATGGCCTTTTTGCATTGGCGCAGAAAAAAATGCCGGATGCGACGCTGGCGCGTCTTATCCAGCCTACTCTTGCTTATCCGTTCTGCTTTTGCATCCACCCCACCAaaaa                        >  2:130665/1‑149 (MQ=255)
                                                                                             tGCCAGGCATCAAATTAAGCAGAAGGCCATCCTGACGGATGGCCTTTTTGCATTGGCGCAGAAAAAAATGCCGGATGCGACGCTGGCGCGTCTTATACAGCCTACTCTTGCTTATCCGTTCTGCTTTTGCATCCACTCCACCACAAAAt                      >  2:875096/1‑149 (MQ=255)
                                                                                              gCCAGGCATCAAATTAAGCAGAAGGCCATCCTGACGGATGGCCTTTTTGCATTGGCGCAGAAAAAAATGCCGGATGCGACGCTGGCGCGTCTTATCCAGCCTACTCTTGCTTATCCGTTCTGCTTTTGCATCCACTCCACCACAAAATc                     <  1:46906/149‑1 (MQ=255)
                                                                                               ccAGGCATCAAATTAAGCAGAAGGCCATCCTGACGGATGGCCTTTTTGCATTGGCGCAGAAAAAAATGCCGGATGCGACGCTGGCGCGTCTTATCCAGCCTACTCTTGCTTATCCGTTCTGCTTTTGCATCCACTCCACCACAAAATca                    <  1:454727/149‑1 (MQ=255)
                                                                                               ccAGGCATCAAATTAAGCAGAAGGCCATCCTGACGGATGGCCTTTTTGCATTGGCGCAGAAAAAAATGCCGGATGCGACGCTGGCGCGTCTTATCCAGCCTACTCTTGCTTATCCGTTCTGCTTTTGCATCCACTCCACCACAAAATca                    >  2:840882/1‑149 (MQ=255)
                                                                                                   gCATCAAATTAAGCAGAAGGCCATCCTGACGGATGGCCTTTTTGCATTGGCGCAGAAAAAAATGCCGGATGCGACGCTGGCGCGTCTTATCCAGCCTACTCTTGCTTATCCGTTCTGCTTTTGCATCCACTCCACCACCAAATcaacca                >  1:708091/1‑149 (MQ=255)
                                                                                                         aaTTAAGCAGAAGGACCTCTTGAGGTATGGCCTTTTTGCATTGGCGCAGAAAAAAATGCCGGATGCGACGCTGGCGCGTCTTATCCAGCCTACTCTTGCTTATCCGTTCTGCTTTTGCATCCACTCCACCACAAAATCAGCCAGCCCCt          >  2:460883/1‑149 (MQ=38)
                                                                                                            tAAGCAGAAGGCCATCCTGACGGATGGCCTTTTTGCATTGGCGCAGAAAAAAATGCCGGATGCGACGCTGGCGCGTCTTATCCAGCCTACTCTTGCTTATCCGTTCTGCTTTTGCATCCACTCCACCACAAAATCAGCCAGCCCCTCaa       >  1:270976/1‑149 (MQ=255)
                                                                                                                 aGAAGGCCATCCTGACGGATGGCCTTTTTGCATTGGCGCAGAAAAAAATGCCGGATGCGACGCTGGCGCGTCTTATCCAGCCTACTCTTGCTTATCCGTTCTGCTTTTGCATCCAATCCACCCAACAATAAGCCACCCCCGAAAcataa  >  2:374581/1‑147 (MQ=25)
                                                                                                                  |                                                                                                                                                   
CGCTCTCCTGAGTAGGACAAATCCGCCGGGAGCGGATTTGAACGTTGCGAAGCAACGGCCCGGAGGGTGGCGGGCAGGACGCCCGCCATAAACTGCCAGGCATCAAATTAAGCAGAAGGCCATCCTGACGGATGGCCTTTTTGCATTGGCGCAGAAAAAAATGCCGGATGCGACGCTGGCGCGTCTTATCCAGCCTACTCTTGCTTATCCGTTCTGCTTTTGCATCCACTCCACCACAAAATCAGCCAGCCCCTCAACATCA  >  NC_000913_3_pae_pgi/4040702‑4040963

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