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

AGGGAACTGCCAGGCATCAAATAAAACGAAAGGCTCAGTCGGAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGACAAATCCGCCGGGAGCGGATTTGAACGTTGCGAAGCAACGGCCCGGAGGGTGGCGGGCAGGACGCCCGCCATAAACTGCCAGGCATCAAATTAAGCAGAAGGCCATCCTGACGGATGGCCTTTTTGCATTGGCGCAGAAAAAAATGCCGGATGCGACGCTGGCGCGTCTTATCCAGCCTACTCTTGCTTATCCGTTC  >  NC_000913_3_bme_pgi/4040620‑4040915                                                                                                                                                      |                                                                                                                                                   aGGGAACTGCCAGGCATCAAATAAAACGAAAGGCTCAGTCGGAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGACAAATCCGCCGGGAGCGGATTTGAACGTTGCGAAGCAACGGCCCGGAGGGt                                                                                                                                                    >  2:185137/1‑150 (MQ=32)   ggtggagtggatgcaaaagcagaacggataagcaagagtaggctggataagacgcgccagcgtcgcatccggcatttttttctgcgccaatgcaaaaaggccatccgtcaggatggccttctgcttaatttgatgcctggcagtttaTgg                                                                                                                                                  >  3:160045‑M1/148‑150 (MQ=255)       gagtggatgcaaaagcagaacggataagcaagagtaggctggataagacgcgccagcgtcgcatccggcatttttttctgcgccaatgcaaaaaggccatccgtcaggatggccttctgcttaatttgatgcctggcagtttaTGGCggg                                                                                                                                              <  4:77506‑M1/7‑1 (MQ=255)       gagtggatgcaaaagcagaacggataagcaagagtaggctggataagacgcgccagcgtcgcatccggcatttttttctgcgccaatgcaaaaaggccatccgtcaggatggccttctgcttaatttgatgcctggcagtttaTGGCggg                                                                                                                                              <  1:268679‑M1/7‑1 (MQ=255)       cTGCCAGGCATCAAATAAAACGAAAGGCTCAGTCGGAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGACCAATCCGCCGGGAGCGGATTTGAACGTTGCGAAGCAACGGCCCGGAGGGTGGCggg                                                                                                                                              <  2:120645/150‑1 (MQ=17)             atgcaaaagcagaacggataagcaagagtaggctggataagacgcgccagcgtcgcatccggcatttttttatgcgtcaatgcaaaaaggccatccgtcaggatggccttgtgcttaatttgatgcctggcagttcaTGgggggccttct                                                                                                                                        <  2:251538‑M1/13‑6 (MQ=255)                              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                                                                                                             <  2:172465‑M1/40‑1 (MQ=255)                                           ggctggataagacgcgccagcgtcgcatccggcatttttttctgcgccaatgcaaaaaggccatccgtcaggatggccttctgcttaatttgatgcctggcagtttaTGGCGGGCAGGACGCCCGCCATAAACTGCCAGGCATCaaa                                                                                                             >  1:172465‑M1/108‑147 (MQ=255)                                                          gccagcgtcgcatccggcatttttttctgcgccaatgcaaaaaggccatccgtcaggatggccttctgcttaatttgatgcctggcagtttaTGGCGGGCAGGAcgcccgcc                                                                                                                                 >  2:293137‑M1/93‑112 (MQ=255)                                                          gccagcgtcgcatccggcatttttttctgcgccaatgcaaaaaggccatccgtcaggatggccttctgctgaatttgatgcctggcagtttaTGGCGGGCAGGAcgcccgcc                                                                                                                                 <  1:293137‑M1/20‑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                                                                                                                        >  2:407141‑M1/32‑60 (MQ=255)                                                                                                                       ccttctgcttaatttgatgcctggcagtttaTGGCGGGCAGGACGCCCGCCATAAACTGc                                                                                                                        <  1:407141‑M1/29‑1 (MQ=255)                                                                                                                           aaCGTTGCGAAGCAACGGCCCGGAGGGTGGCGGGCAGGACGCCCGCCATAAACTGCCAGGCATCAAATTAAGCAGAAGGCCATCCTGACGGATGGCCTTTTTGCATTGGCGCAGAAAAAAATGCCGGATGCGACGCTGGCGCGTCTTATc                          >  4:135954/1‑150 (MQ=255)                                                                                                                                   gAAGCAACGGCCCGGAGGGTGGCGGGCAGGACGCCCGCCATAAACTGCCAGGCATCAAATTAAGCAGAAGGCCATCCTGACGGATGGCCTTTTTGCATTGGCGCAGAAAAAAATGCCGGATGCGACGCTGGCGCGTCTTATCCAGCCTAc                  >  1:360087/1‑150 (MQ=255)                                                                                                                                   gAAGCAACGGCCCGGAGGGTGGCGGGCAGGACGCCCGCCATAAACTGCCAGGCATCAAATTAAGCAGAAGGCCATCCTGACGGATGGCCTTTTTGCATTGGCGCAGAAAAAAATGCCGGATGCGACGCTGGCGAGTCTTATCCAGCCTAc                  >  1:251538/1‑150 (MQ=255)                                                                                                                                              ccGGAGGGTGGCGGGCAGGACGCCCGCCATAAACTGCCAGGCATCAAATTAAGCAGAAGGCCATCCTGACGGATGGCCTTTTTGCATTGGCGCAGAAAAAAATGCCGGATGCGACGCTGGCGCGTCTTATCCAGCCTACTCTTGCTTATc       >  1:264847/1‑150 (MQ=255)                                                                                                                                              ccGGAGGGTGGCGGGCAGGACGCCCGCCATAAACTGCCAGGCATCAAATTAAGCAGAAGGCCATCCTGACGGATGGCCTTTTTGCATTGGCGCAGAAAAAAATGCCGGATGCGACGCTGGCGCGTCTTATCCAGCCTACTCTTGCTTATc       <  4:132485/150‑1 (MQ=255)                                                                                                                                                  aGGGTGGCGGGCAGGACGCCCGCCATAAACTGCCAGGCATCAAATTAAGCAGAAGGCCATCCTGACGGATGGCCTTTTTGCATTGGCGCAGAAAAAAATGCCGGATGCGa                                           <  3:273522/110‑1 (MQ=255)                                                                                                                                                  aGGGTGGCGGGCAGGACGCCCGCCATAAACTGCCAGGCATCAAATTAAGCAGAAGGCCATCCTGACGGATGGCCTTTTTGCATTGGCGCAGAAAAAAATGCCGGATGCGa                                           >  4:273522/1‑110 (MQ=255)                                                                                                                                                   gggTGGCGGGCAGGACGCCCGCCATAAACTGCCAGGCATCAAATTAAGCAGAAGGCCATCCTGACGGATGGCCTTTTTGCATTGGCGCAGAAAAAAATGCCGGATGCGACGCTGGCGCGTCTTATCCAGCCTACTCTTGCTTATCCGTTc  >  4:172668/1‑150 (MQ=255)                                                                                                                                                      |                                                                                                                                                   AGGGAACTGCCAGGCATCAAATAAAACGAAAGGCTCAGTCGGAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGACAAATCCGCCGGGAGCGGATTTGAACGTTGCGAAGCAACGGCCCGGAGGGTGGCGGGCAGGACGCCCGCCATAAACTGCCAGGCATCAAATTAAGCAGAAGGCCATCCTGACGGATGGCCTTTTTGCATTGGCGCAGAAAAAAATGCCGGATGCGACGCTGGCGCGTCTTATCCAGCCTACTCTTGCTTATCCGTTC  >  NC_000913_3_bme_pgi/4040620‑4040915

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.