New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913_3_bme_pgi 4040769 =20 (0.540)8 (0.230) 8/228 2.8 29.8% intergenic (+133/+137) rrfA/mobB 5S ribosomal RNA of rrnA operon/molybdopterin‑guanine dinucleotide biosynthesis protein B
?NC_000913_3_bme_pgi 4040790 = 19 (0.550)intergenic (+154/+116) rrfA/mobB 5S ribosomal RNA of rrnA operon/molybdopterin‑guanine dinucleotide biosynthesis protein B

AAGGCTCAGTCGGAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGACAAATCCGCCGGGAGCGGATTTGAACGTTGCGAAGCAACGGCCCGGAGGGTGGCGGGCAGGACGCCCGCCATAAACTGCCAGGCATCAAATTAAGCAGAAGGCCATCCTGACGGATGGCCTTTTTGCATTGGCGCAGAAAAAAATGCCGGATGCGACGCTGGCGCGTCTTATCCAGCCTACTCTTGCTTATCCGTTCTGCTTTTGCATCCACTCCACC  >  NC_000913_3_bme_pgi/4040649‑4040936
                                                                                                                                             |                                                                                                                                                  
aaGGCTCAGTCGGAAGACTGGGCCTTTCGTTTTATCTGTTGTTTCTCGGTGAACGCTCTCCTGAGTAGGACAAATCCGCCGGGAGCGGATTTGAACGTTGCGAAGCAACGGCCCGGAGGGTGGCGGGCAGGACGCCCGCCATAAACTGcc                                                                                                                                            >  2:327369/1‑150 (MQ=17)
  ggCTCAGTCGGAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGACAAATCCGCCGGGAGCGGATTTGAACGTTGCGAAGCAACGGCCCGGAGGGTGGCGGGCAGGACGCCCGCCATAAACTGCCAg                                                                                                                                          <  1:327369/150‑1 (MQ=32)
   gCTCAGTCGGAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGACAAATCCGCCGGGAGCGGATTTGAACGTTGCGAAGCAACGGCCCGGAGGGTGGCGGGCAGGACGCCCGCCATAAACTGCCAgg                                                                                                                                         >  1:76390/1‑150 (MQ=32)
     tCAGTCGGAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGACAAATCCGCCGGGAGCGGATTTGAACGTTGCGAAGCAACGGCCCGGAGGGTGGCGGGCAGGACGCCCGCCATAAACTGCCAGGCa                                                                                                                                       >  3:77506/1‑150 (MQ=32)
        gTCGGAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGACAAATCCGCCGGGAGCGGATTTGAACGTTGCGAAGCAACGGCCCGGAGGGTGGCGGGCAGGACGCCCGCCATAAACTGCCAGGCATCa                                                                                                                                    >  3:132485/1‑150 (MQ=32)
             ggctggataagacgcgccagcgtcgcatccggcatttttttctgcgccaatgcaaaaaggccatccgtcaggatggccttctgcttaatttgatgcctggcagtttaTGGCGGGCAGGACGCCCGCCATAAACTGCCAGGCATCaaa                                                                                                                                  >  1:172465‑M1/108‑147 (MQ=255)
             ggctggataagacgcgccagcgtcgcatccggcatttttttctgcgccaatgcaaaaaggccatccgtcaggatggccttctgcttaatttgatgcctggcagtttaTGGCGGGCAGGACGCCCGCCATAAACTGCCAGGCATCaaa                                                                                                                                  <  2:172465‑M1/40‑1 (MQ=255)
                                                            cTGAGTAGGACAAATCCGCCGGGAGCGGATTTGAACGTTGCGAAGCAACGGCCCGGAGGGTGGCGGGCAGGACGCCCGCCATAAACTGCCAGGCATCAAATTAAGCAGAAGGCCATCCTGACGGATGGCCTTTTTGCATTGGCGCAGaaa                                                                                <  3:128627/150‑1 (MQ=40)
                                                                  aaaaaggccatccgtcaggatggccttctgcttaatttgatgcctggcagtttaTGGCGGGCAGGACGCCCGCCATaa                                                                                                                                                  <  1:284210‑M1/24‑1 (MQ=255)
                                                                  aaaaaggccatccgtcaggatggccttctgcttaatttgatgcctggcagtttaTGGCGGGCAGGACGCCCGCCATaa                                                                                                                                                  >  2:284210‑M1/55‑78 (MQ=255)
                                                                    gACAAATGCGCCGGGAGCGGATTTGAACGTTGCGAAGCAACGGCCCGGAGGGTGGCGGGCAGGACGCCCGCCATAAACTGCCAGGCATCAAATTAAGCAGAAGGCCATCCTGACGGATGGCCTTTTTGCATTGGCGCAGAAAAAAATGcc                                                                        >  2:472122/1‑150 (MQ=255)
                                                                    gACAAATCCGCCGGGAGCGGATTTGAACGTTGCGAAGCAACGGCCCGGAGGGTGGCGGGCAGGACGCCCGCCATAAACTGCCAGGCATCAAATTAAGCAGAAGGCCATCCTGACGGATGGCCTTTTTGCATTGGCGCAGAAAAAAATGcc                                                                        >  2:268679/1‑150 (MQ=255)
                                                                                     cGGATTTGAACGTTGCGAAGCAACGGCCCGGAGGGTGGCGGGCAGGACGCCCGCCATAAACTGCCAGGCATCAAATTAAGCAGAAGGCCATCCTGACGGATGGCCTTTTTGCATTGGCGCAGAAAAAAAt                                                                           <  2:25364/130‑1 (MQ=255)
                                                                                     cGGATTTGAACGTTGCGAAGCAACGGCCCGGAGGGTGGCGGGCAGGAAGCCCGCCATAAACTGCCAGGCATCAAATTAAGCAGAAGGCCATCCTGACGGATGGCCTTTTTGCATTGGCGCAGAAAAAAAt                                                                           >  1:25364/1‑130 (MQ=255)
                                                                                         ccttctgcttaatttgatgcctggcagtttaTGGCGGGCAGGACGCCCGCCATAAACTGc                                                                                                                                             <  1:407141‑M1/29‑1 (MQ=255)
                                                                                         ccttctgcttaatttgatgcctggcagtttaTGGCGGGCAGGACGCCCGCCATAAACTGc                                                                                                                                             >  2:407141‑M1/32‑60 (MQ=255)
                                                                                             aaCGTTGCGAAGCAACGGCCCGGAGGGTGGCGGGCAGGACGCCCGCCATAAACTGCCAGGCATCAAATTAAGCAGAAGGCCATCCTGACGGATGGCCTTTTTGCATTGGCGCAGAAAAAAATGCCGGATGCGACGCTGGCGCGTCTTATc                                               >  4:135954/1‑150 (MQ=255)
                                                                                                     tttgatgcctggcagtttatggcgggcgtcctgcccgccaTAAACTGCCAGGCATCAAATTAAGCAGAAGGCCATCCTGACGGATGGCCTTTTTGCATTGGCGCAGAAAAAAATGCCGGATGCGACGCTGGCGCGTCTTATCCAGcc                                          >  1:172465‑M2/41‑147 (MQ=255)
                                                                                                     tttgatgcctggcagtttatggcgggcgtcctgcccgccaTAAACTGCCAGGCATCAAATTAAGCAGAAGGCCATCCTGACGGATGGCCTTTTTGCATTGGCGCAGAAAAAAATGCCGGATGCGACGCTGGCGCGTCTTATCCAGcc                                          <  2:172465‑M2/107‑1 (MQ=255)
                                                                                                     gAAGCAACGGCCCGGAGGGTGGCGGGCAGGACGCCCGCCATAAACTGCCAGGCATCAAATTAAGCAGAAGGCCATCCTGACGGATGGCCTTTTTGCATTGGCGCAGAAAAAAATGCCGGATGCGACGCTGGCGCGTCTTATCCAGCCTAc                                       >  1:360087/1‑150 (MQ=255)
                                                                                                     gAAGCAACGGCCCGGAGGGTGGCGGGCAGGACGCCCGCCATAAACTGCCAGGCATCAAATTAAGCAGAAGGCCATCCTGACGGATGGCCTTTTTGCATTGGCGCAGAAAAAAATGCCGGATGCGACGCTGGCGAGTCTTATCCAGCCTAc                                       >  1:251538/1‑150 (MQ=255)
                                                                                                                gcagtttatggcgggcgtcctgcccgccaTAAACTGCCAGGCATCAAATTAAGCAGAAgg                                                                                                                      <  1:407141‑M2/31‑1 (MQ=255)
                                                                                                                gcagtttatggcgggcgtcctgcccgccaTAAACTGCCAGGCATCAAATTAAGCAGAAgg                                                                                                                      >  2:407141‑M2/30‑60 (MQ=255)
                                                                                                                ccGGAGGGTGGCGGGCAGGACGCCCGCCATAAACTGCCAGGCATCAAATTAAGCAGAAGGCCATCCTGACGGATGGCCTTTTTGCATTGGCGCAGAAAAAAATGCCGGATGCGACGCTGGCGCGTCTTATCCAGCCTACTCTTGCTTATc                            >  1:264847/1‑150 (MQ=255)
                                                                                                                ccGGAGGGTGGCGGGCAGGACGCCCGCCATAAACTGCCAGGCATCAAATTAAGCAGAAGGCCATCCTGACGGATGGCCTTTTTGCATTGGCGCAGAAAAAAATGCCGGATGCGACGCTGGCGCGTCTTATCCAGCCTACTCTTGCTTATc                            <  4:132485/150‑1 (MQ=255)
                                                                                                                    aGGGTGGCGGGCAGGACGCCCGCCATAAACTGCCAGGCATCAAATTAAGCAGAAGGCCATCCTGACGGATGGCCTTTTTGCATTGGCGCAGAAAAAAATGCCGGATGCGa                                                                >  4:273522/1‑110 (MQ=255)
                                                                                                                    aGGGTGGCGGGCAGGACGCCCGCCATAAACTGCCAGGCATCAAATTAAGCAGAAGGCCATCCTGACGGATGGCCTTTTTGCATTGGCGCAGAAAAAAATGCCGGATGCGa                                                                <  3:273522/110‑1 (MQ=255)
                                                                                                                     ttatggcgggcgtcctgcccgccaTAAACTGCCAGGCATCAAATTAAGCAGAAGGCCATCCTGACGGATGGCCttttt                                                                                               <  1:284210‑M2/54‑1 (MQ=255)
                                                                                                                     ttatggcgggcgtcctgcccgccaTAAACTGCCAGGCATCAAATTAAGCAGAAGGCCATCCTGACGGATGGCCttttt                                                                                               >  2:284210‑M2/25‑78 (MQ=255)
                                                                                                                     gggTGGCGGGCAGGACGCCCGCCATAAACTGCCAGGCATCAAATTAAGCAGAAGGCCATCCTGACGGATGGCCTTTTTGCATTGGCGCAGAAAAAAATGCCGGATGCGACGCTGGCGCGTCTTATCCAGCCTACTCTTGCTTATCCGTTc                       >  4:172668/1‑150 (MQ=255)
                                                                                                                         ggcgggcgtcctgcccgccaTAAACTGCCAGGCATCAAATTCAGCAGAAGGCCATCCTGACGGATGGCCTTTTTGCATTGGCGCAGAAAAAAATGCCGGATGCGACGCTGgc                                                         <  1:293137‑M2/92‑1 (MQ=255)
                                                                                                                         ggcgggcgtcctgcccgccaTAAACTGCCAGGCATCAAATTAAGCAGAAGGCCATCCTGACGGATGGCCTTTTTGCATTGGCGCAGAAAAAAATGCCGGATGCGACGCTGgc                                                         >  2:293137‑M2/21‑112 (MQ=255)
                                                                                                                         GGCGGGCAGGACGCCCGCCATAAACTGCCAGGCATCAAATTAAGCAGAAGGCCATCCTGACGGATGGCCTTTTTGCATTGGCGCAGAAAAAAATGCCGGATGCGACGCTGGCGCGTCTTATCCAGCCTACTCTTGCTTATCCGTTc                       >  2:420626‑M1/1‑146 (MQ=255)
                                                                                                                         GGCGGGCAGGACGCCCGCCATAAACTGCCAGGCATCAAATTAAGCAGAAGGCCATCCTGACGGATGGCCTTTTTGCATTGGCGCAGAAAAAAATGCCGGATGCGACGCTGGCGCGTCTTATCCAGCCTACTCTTGCTTATCCGTTc                       <  1:420626‑M1/146‑1 (MQ=255)
                                                                                                                           CGGGCAGGACGCCCGCCATAAACTGCCAGGCATCAAATTAAGCAGAAGGCCATCCTGACGGATGGCCTTTTTGCATTGGCGCAGAAAAAAATGCCGGATGCGACGCTGGCGCGTCTTATCCAGCCTACTCTTGCTTATCCGTTCTGCttt                 >  3:200220‑M1/1‑150 (MQ=255)
                                                                                                                            GGGCAGGACGCCCGCCATAAACTGCCAGGCATCAAATTAAGCAGAAGGCCATCCTGACGGATGGCCTTTTTGCATTGGCGCAGAAAAAAATGCCGGATGCGACGCTGGCGCGTCTTATCCAGCCTACTCTTGCTTATCCGTTCTGCtttt                >  4:35724‑M1/1‑150 (MQ=255)
                                                                                                                                agaaggccccccaTGAACTGCCAGGCATCAAATTAAGCACAAGGCCATCCTGACGGATGGCCTTTTTGCATTGACGCATAAAAAAATGCCGGATGCGACGCTGGCGCGTCTTATCCAGCCTACTCTTGCTTATCCGTTCTGCTTTTGCAt            <  2:251538‑M2/137‑1 (MQ=255)
                                                                                                                                    cgcccgccATAAACTGCCAGGCATCAAATTAAGCAGAAGGCCATCCTGACGGATGGCCTTTTTGCATTGGCGCAGAAAAAAATGCCGGATGCGACGCTGGCGCGTCTTATCCAGCCTACTCTTGCTTATCCGTTCTGCTTTTGCAtccac        <  3:172668/150‑1 (MQ=255)
                                                                                                                                      cccgccaTAAACTGCCAGGCATCAAATTAAGCAGAAGGCCATCCTGACGGATGGCCTTTTTGCATTGGCGCAGAAAAAAATGCCGGATGCGACGCTGGCGCGTCTTATCCAGCCTACTCTTGCTTATCCGTTCTGCTTTTGCAtccactc      <  1:268679‑M2/143‑1 (MQ=255)
                                                                                                                                      cccgccaTAAACTGCCAGGCATCAAATTAAGCAGAAGGCCATCCTGACGGATGGCCTTTTTGCATTGGCGCAGAAAAAAATGCCGGATGCGACGCTGGCGCGTCTTATCCAGCCTACTCTTGCTTATCCGTTCTGCTTTTGCAtccactc      <  4:77506‑M2/143‑1 (MQ=255)
                                                                                                                                          ccaTAAACTGCCAGGCATCAAATTAAGCAGAAGGCCATCCTGACGGATGGCCTTTTTGCATTGGCGCAGAAAAAAATGCCGGATGCGACGCTGGCGCGTCTTATCCAGCCTACTCTTGCTTATCCGTTCTGCTTTTGCATCCACTccacc  >  3:160045‑M2/4‑150 (MQ=255)
                                                                                                                                             tAAACTGCCAGGCATCAAATTAAGCAGAAGGCCATCCTGACGGATGGCCTTTTTGCATTGGCGCAGAAAAAAATGCCGGATGCGACGCTGGCGCGTCTTATCCAGcc                                          <  2:13238/107‑1 (MQ=255)
                                                                                                                                             tAAACTGCCAGGCATCAAATTAAGCAGAAGGCCATCCTGACGGATGGCCTTTTTGCATTGGCGCAGAAAAAAATGCCGGATGCGACGCTGGCGCGTCTTATCCAGcc                                          >  1:13238/1‑107 (MQ=255)
                                                                                                                                             |                                                                                                                                                  
AAGGCTCAGTCGGAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGACAAATCCGCCGGGAGCGGATTTGAACGTTGCGAAGCAACGGCCCGGAGGGTGGCGGGCAGGACGCCCGCCATAAACTGCCAGGCATCAAATTAAGCAGAAGGCCATCCTGACGGATGGCCTTTTTGCATTGGCGCAGAAAAAAATGCCGGATGCGACGCTGGCGCGTCTTATCCAGCCTACTCTTGCTTATCCGTTCTGCTTTTGCATCCACTCCACC  >  NC_000913_3_bme_pgi/4040649‑4040936

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