New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913_3_bme_pgi	2095931 =	23 (0.490)	3 (0.070)	3/268	NT	9.4%	coding (88/777 nt)	hisF	imidazole glycerol phosphate synthase, catalytic subunit with HisH
	?	NC_000913_3_bme_pgi	2096003 =	35 (0.760)					coding (160/777 nt)	hisF	imidazole glycerol phosphate synthase, catalytic subunit with HisH

AATCCCATGTCTCGACGTTCGTGATGGTCAGGTGGTGAAAGGCGTACAGTTTCGCAACCATGAAATCATTGGCGATATCGTGCCGCTGGCAAAACGCTACGCTGAAGAAGGCGCTGACGAACTGGTGTTCTACGA‑TATCACCGCTTCCAGCGATGGCCGTGTGGTAGATAAAAGCTGGGTATCTCGCGTGGCGGAAGTGATCGACATTCCGTTTTGTGTGGCGGGTGGGATTAAGTCTCTGGAAGATGCCGCGAAAATTCTTTCCTTTGGCGCGGATAAAA  >  NC_000913_3_bme_pgi/2095861‑2096141                                                                                                                                                |                                                                                                                                           acggaatgtcgatcacttccgccacgcgagatacccagcttttatctaccacacggccatcgctggaagcGGCGATATCGTGCCGCTGGCAAAACGCTACGCTGAAGAAGGCGCTGACGAACTGGTGTTCTACGA‑TATCACCGCTTCCa                                                                                                                                      >  2:776884‑M1/71‑149 (MQ=255)      gacgTCGCGACGTTCGTGATGGTCAGGTTGTGAAAGGCGTACTGTTTCGCAACCATGAAATCATTGGCGATATCGTGCCGCTGGCAAAACGCTACGCTGAAGAAGGCGCTGACGAACTGGTGTTCTACGA‑TATCACCGCTTCCAGCGAt                                                                                                                                 <  2:392851/146‑1 (MQ=255)            tcGACGTTCGTGATGGTCAGGTGGTGGAAGGCGTACAGTTTCGCAACCATGAAAGCATTGGCGATATCGTGCCGCTGGCCAAACGCTACGCTGGAGAAGGCGCTGACGAACTGGTGTTACACGA‑TATCCCGGATTCCAGCGATGGCCgt                                                                                                                           >  1:295862/1‑149 (MQ=255)                      tGATGGTCAGGTGGTGAAAGGCGTACAGTTTCGCAACCATGAAATCATTGGCGATATCGTGCCGCTGGCAAGACGCTACGCTGAAGAAGGCGCTGACGAACTGGTGTTCTACGA‑TATCACCGCCTCCAGCACTGGCCGTGTgggctgta                                                                                                                 >  1:87494/1‑143 (MQ=255)                      tGATGGTCAGGTGGTGAAAGGCGTACAGTTTCGCAACCATGAAATCATTGGCGATATCGTGCCGCTGGCAAAACGCTACGCTGAAGAAGGCGCTGACGAACTGGTGTTCTACGA‑TGTCACCGCTTCCAGCGCTGGCCGTTTGGTCGATa                                                                                                                 >  1:441930/1‑149 (MQ=255)                        aTGGTCAGGTGGTGAAAGGCGTACAGTTTCGCAACCATGAAATCATTGGCGATATCGTGCCGCTGGCAAAACGCTACGCTGAAGAAGGCGCTGACGAACTGGTGTTCTACGA‑TATCACAGCTTCCAGTGATGCTCGTATGGTAGATaaa                                                                                                               >  1:470860/1‑149 (MQ=255)                              aGGTGGTGAAAGGCGTAAAGTTTCGCAACCATGAAATCATTGGCGATATCGTGCCGCTGGCAAAACGCTACGCTGAAGAAGGCGCTGACGAACTGGTGTTCTACGA‑TATCACCGCTTCCAGCGATGGCCGTGTGGTAGATAAAAGCTgg                                                                                                         <  2:754470/149‑1 (MQ=255)                                     gAAAGGCGTACAGTTTCGCAACCATGAAATCATTGGCGATATCGTGCCGCTGGCAAAACGCTACGCTGAAGAAGGCGCTGACGAACTGGTGTTCTACGA‑TATCACCGCTTCCAGCGATGGCCGGGTGGTAGATAAAAGCCGGGTAtctc                                                                                                  >  1:820968/1‑149 (MQ=255)                                            gTACAGTTTCGCAACCATGAAATCATTGGCGATATCGTGCCGCTGGCAAAACGCTACGCTGAAGAAGGCGCTGACGAACTGGTGTTCTACGACAACCACCG‑TGCTATCGATGGCCGTGTGGTAGATAAAAGCTGGGTATCTCGCGTGGc                                                                                           <  2:597409/149‑1 (MQ=255)                                             tACAGTTTCGCAACCATGAAATCATTGGCGATATCGTGCCGCTGGCAAAACGCTACGCTGAAGAAGGCGCTGACGAACTGGTGTTCTACGA‑TATCACCGCTTCCCGCGATGGCCGTGTGGTAGATTAAAGCTGGGATTCTCGCATTGCg                                                                                          >  1:413847/1‑149 (MQ=255)                                             gaccgTTTCGCAGCCATGAAATCATTGGCGATATCGTGCCGCTGGCAAAACGCTACGCTGAAGAAGGCGCTGACGAACTGGTGTTCTACGA‑TATCACCGCTGCGAGGGATGGCCGTGTGGTAGATAAAAGCTGGGTATCTCGCGTGGCg                                                                                          <  1:484925/145‑1 (MQ=255)                                                 gTTGCGGAAGCATGAAATCATGGGCGATATCGAGCCGCGGGCAAAACGCTACGCTGAAGAAGGCGCTGACGAACTGGTGTTCTACGA‑TATCACCGCTTCCAGCGATGGCCGTGTGGTAGATAAAAGCTGGGTATCTCGCGTGGCGGAAg                                                                                      <  2:335945/149‑1 (MQ=255)                                                    tCGCAACCATGAAATCATTGGCGATATCGTGCCGCTGGCAAAACGCTACGCTGAAGAAGGCGCTGACGAACTGGTGTTCTACGA‑TATCACCGCTTCCAGCGATGGCCGTGTGGTAGATAAAAGCTGGGTATCTCGCGTGGCGGAAGTGa                                                                                   >  2:479185/1‑149 (MQ=255)                                                      ggagcgaTGAAATGGTTGGCGATATCGTGCCGCTGGCAAAACGCTACGCTGAAGAAGGCGCTGACGAACTGGTGTTCTACGA‑TATCACCGCTTCCAGCGATGGCCGTGTGGTAGATAAAAGCTGGGTATCTCGCGTGGCGGAAGTGATc                                                                                 <  1:649109/143‑1 (MQ=255)                                                           cATGAAATCATTGGCGATATCGTGCCGCTGGCAAAACGCTACGCTGAAGAAGGCGCTGACGAACTGGTGTTCTACGA‑TATCACCGCTTCCAGCGATGGCCGTGTGGTAGATAAAAGCTGGGTATCCCTCGTGGCGGCTATGATCTACAt                                                                            >  1:381431/1‑149 (MQ=255)                                                                 tggaagcggtgatatcgtagaacaccagttcgtcagcgccttcttcagcgtagcgttttgccagcggcacgatatcgccGCTTCCAGCGATGGCCGTGTGGTAGATAAAAGCTGGGTATCTCGCGTGGCGGAAGTGATCGACATTCCGt                                                                       >  2:776884‑M2/80‑149 (MQ=255)                                                                     ttGGCGATATCGTGCCACTGGCAAAACGCCACGCTGAAGAAGGCGGTGACGAACTTGTGTTCTACGA‑TATCACCGCTGCCAGCGATGGCCGTGTGGTAGATAAAAGCTGGGTATCTCGCGTGGCGGAAGTGATCGACATTCCGTTTtgt                                                                  <  2:87494/149‑1 (MQ=255)                                                                         cGATATCGTGCCGCTGGCAAAACGCTACGCTGAAGAAGGCGCTGACGAACTGGTGTTCTACGA‑TATCACCGCTTCCAGCGCTGGCCGTGTGGTAGATAGAAGCTGGGTATCTCGCGTGGCGGACGTGATCCACACTCCTTTCTGTGTgg                                                              >  1:215190/1‑149 (MQ=255)                                                                                  gCCGGTGGCAAAACGTCACGCTGAAGAAGGCGCTAACGAACTGGTGTTGTACGA‑TTTCACCGCTTCCAGCGATGGCCGTGTGGTAGATAAAAGCTGGGTATCTCGCGTGGCGGAAGTGATCGACATTCCGTTTTGTGTGGCGGGTGGGa                                                     <  1:225432/149‑1 (MQ=255)                                                                                   ccGCTGGCAAAACGCTACGCGGAAGAAGGCGCTGACGAAATGGTTTTCTCCGA‑TATTAGCGCTTCCAGCGATGGCCGTGTGGGAGATAAAAGCTGGGTATCTCGCGTGGCGGAAGTGATCGACATTCCGTTTTGTGTGGCGGGTGGGAt                                                    <  2:453402/149‑1 (MQ=255)                                                                                     gCTGGCAAAACGCTACGCTGAAAAAGAGGCATACTAACTGGTGTTCTACGA‑TATCACCGCTTCCAGCGATGGCCGTGTGGTAGATAAAAGCTGGGTATCTCGCGTGGCGGAAGTGATCGACATTCCGTTTTGTGTGGCGGGTGGGATTa                                                  >  2:870888/1‑149 (MQ=255)                                                                                      gtGGCAAAACGCGACGCGGAAGAAGGCGCTGACGATCTGGTGGTCTTCGA‑TATTACCGCTTCCAGCGATGGCCGTGTGGTAGATAAAAGCTGGGTATCTCGCGTGGGGGAAGTGATCGACATTCCGTTTTGTGTGGCGGGTGGGATTaa                                                 <  1:292083/148‑1 (MQ=255)                                                                                      cTGGCAAAACGCTACGCTGAAGAAGGCGCTGACGAACTGGTGTTCTACGA‑TATCACCGCTTCCAGCGATGGCCGTGTGGTAGATAAAAGCTGGGTATCTCGCGTGGCGGAAGTGATCGACAGTCAGTGTTGTGTGGCGGGTGGGATTaa                                                 >  1:446806/1‑149 (MQ=255)                                                                                          cAAAACGCTACGCTGAAGAAGGCGCTGACGAACTGGTGTTCTACGA‑TATCACCGCTTCCAGCGATGGCCGTGTGGTAGATAAAAGCTGGGTATCTCGCGTGGCGGAAGTGATCGACATTCCATTTTGTGTGGCGGGTGGGATTAAGact                                             >  1:806696/1‑149 (MQ=255)                                                                                           aaaaCGCTACGCTGAAGAAGGCGCTGACGAACTGGTGTTCTACGA‑TATCACCGCTTCCAGCGATGGCCGCGTGGTAGATAAAAGCTGGGGATCTCGCGTGGCGGAAGTGATCGCCATTCCGTTTTTGTTTGCGGGTGGGATAAAAtctc                                            >  1:285012/1‑149 (MQ=255)                                                                                            aaaCGCGACGCTGAACAAGCGTCTGACGAACTGGTGTTCTACGA‑TATGCCCACTGCCAGCGATGGCCGTGTGGTAGATAAAAGCTGGGTATCTCGCGTGGCGGAAGTGATCGACATTCCGTTTTGTGTGGCGGGTGGGATTAAGTctct                                           <  1:618945/149‑1 (MQ=255)                                                                                                    gcgccttcttcagcgtagcgttttgccagcggcacgatatcgccGCTTCCAGCGATGGCCGTGTGGTAGATAAAAGCTGGGTATCTCGCGTGGCGGAAGTGATCGACATTCCGTTTTGTGTGGCAGGTGGGATTAAGTCTCTGGAAGAt                                    <  1:776884‑M2/105‑1 (MQ=255)                                                                                                     ggcgaagaaGGCGGTGAGGAAATGGTGCTCTACGA‑TATCACCGCTTGCAGGGATTGCCGTGTGGTAGAATAAAGCCGGGTATCTAGCGTATCGGAAGTGATCGACATTCCGTTTTGTGTGGCGGGTGGGATTAAGTCTCTGGAAGATgg                                  <  2:680291/146‑2 (MQ=255)                                                                                                       tGAAGAAGGCGCTGACGAACTGGTGTTCTACGA‑TATCACCGCTTCCAGCGATGGCCGTGTGGTAGATACAAGCTGGGTATCTCGCGTGGCGGGAGCGATCGACATTCCGATTTGTATAACGGGTGGGATTTAGTCGCTGAAAGATACcg                                >  1:686372/1‑149 (MQ=255)                                                                                                       tGAAGAAGGCGCTGACGAACTGGTGTTCTACGA‑TATCACCGCTTCCAGCGATGGCCGTGTGGTAGATAAAATCTGGGTATCTCGCGTGTCGGAAGTGATCGACATTCGGTTTTGTGTGGCGAATAGGATTAAGTGTCTGCAAGATGcag                                >  1:19914/1‑147 (MQ=255)                                                                                                          agaagGCGCTGACGAACTGGTGTTCTACGA‑TATAACCGCTTCCAGCGATGGCCGTGTGGTAGATAAAAGCTGGGTATCTCGCGTGGCGGAAGTGATCGACATTCCGTTTTGTGTGGCGGGTGGGATTAAGTCTCTGGAAGATGcaccga                             >  2:869965/1‑144 (MQ=255)                                                                                                               gcgcCGACGAACTGGTTTTCTACGA‑TATCACCGCCTCCAGCGATCGCCGTATTCTATACAAAAGCTGGGTATCTAGCGTGGCGGAAGTGATCGACAGTCCGTTTTGTGTGGCGGGTGGGATTAAGTCTCTGGAAGATGCCGCGAAAAtt                        <  1:119713/149‑1 (MQ=255)                                                                                                                          cTGGTGTTCTACGA‑TATCACCGCTTCCAGCGATGGCCGTGTGGTAGATAAAAGCTGGGTATCTCGCGTGGCGGAAGTGATCGACATTCCGTTTTGTGTGGCGGGTGGGATTAAGTATCTGGAAGAGGACGCGAAAATTCTTTCATTTgg             >  1:723159/1‑149 (MQ=255)                                                                                                                           tGGTGTTCTACGA‑TATCACCGCTTCCAGCGATGGCCGTGTGGTAGATAAAAGCTGGGTATCTCGCGTGGCGGAAGTGATCGACATTCCGtttt                                                                    <  1:779686/93‑1 (MQ=255)                                                                                                                           tGGTGTTCTACGA‑TATCACCGCTTCCAGCGATGGCCGTGTGGTAGATAAAAGCTGGGTATCTCGCGTGGCGGAAGTGATCGACATTCCGtttt                                                                    >  2:779686/1‑93 (MQ=255)                                                                                                                             gggtTCTTCGA‑TATCACCGCTTCCAGCGATGGCCGTGTGGTAGATAAAAGCTGGGTATCTCGCGTGGCGGAAGTGATCGACATTCCGTTTTGTGTGGCGGGTGGGATTAAGTCTCTGGAAGATGCCGCGAAAATTCTTTCCTTTGgcgc          <  1:456892/147‑1 (MQ=255)                                                                                                                                 gtctacgatatcgccGCTTCCAGCGATGGCCGTGTGGTAGATAAAAGCTGGGTATCTCGCGTGGCGGAAGTGATCGACATTCCGTTTTGTGTGGTGGGTGGGGTTAAGTGTCTGGAAGATGCCGCGAAAATTCTTTCCTTTGGCGCGGa       <  1:306753‑M2/134‑1 (MQ=255)                                                                                                                                   tACGA‑TATCACCGCTTCCAGCGATGGCCGTGTGGTAGATAAAAGCTGGGTATCTCGCGTGGCGGAAGTGATCGACATTCCGTTTTGTGTGGCGGGTGGGATTAAGTCTCTGGAAGATGCCGCGAAAATTCTTTCCTTTGGCGCGGATaa    >  2:889833/1‑149 (MQ=255)                                                                                                                                     cGA‑TACCATCGCTGCCAGCGATGGCTCTGTGGGAGATAAAAGCTTGGTATCTCGCGTGGCGGAAGTGATCGACATTCCGTTTTGTGGGGCGGGTGGGATTAAGTCTCTGGAAGATGCCGCAACAATTCTTTCCTTTAGCggggcgtaaa  <  2:872905/149‑7 (MQ=255)                                                                                                                                                |                                                                                                                                           AATCCCATGTCTCGACGTTCGTGATGGTCAGGTGGTGAAAGGCGTACAGTTTCGCAACCATGAAATCATTGGCGATATCGTGCCGCTGGCAAAACGCTACGCTGAAGAAGGCGCTGACGAACTGGTGTTCTACGA‑TATCACCGCTTCCAGCGATGGCCGTGTGGTAGATAAAAGCTGGGTATCTCGCGTGGCGGAAGTGATCGACATTCCGTTTTGTGTGGCGGGTGGGATTAAGTCTCTGGAAGATGCCGCGAAAATTCTTTCCTTTGGCGCGGATAAAA  >  NC_000913_3_bme_pgi/2095861‑2096141

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