New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913_3_bme_pgi	3950444 =	23 (0.620)	9 (0.260)	9/230	2.6	31.2%	intergenic (+24/‑63)	ilvL/ilvX	ilvG operon leader peptide/uncharacterized protein
	?	NC_000913_3_bme_pgi	3950470 =	18 (0.520)					intergenic (+50/‑37)	ilvL/ilvX	ilvG operon leader peptide/uncharacterized protein

TTCGAACAAGATGCAAGAAAAGACAAAATGACAGCCCTTCTACGAGTGATTAGCCTGGTCGTGATTAGCGTGGTGGTGATTATTATCCCACCGTGCGGGGCTGCACTTGGACGAGGAAAGGCTTAGAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGGGTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAATGAATAACAGCACAAAATTCTGTTTCTCAAGATTCAGGACGGGGAACTAACTATGAATGGCGCACAGTGGGTGGTAC  >  NC_000913_3_bme_pgi/3950295‑3950584                                                                                                                                                      |                                                                                                                                             ttCGAACAAGATGCAAGAAAAGACAAAATGACAGCCCTTCTACGAGTGATTAGCCTGGTCGTGATTAGCGTGGTGGTGATTATTATCCCACCGTGCGGGGCTGCACTTGGACGAGGAAAGGCTTAGAGATCAAGCCTTAACGAACTAAGa                                                                                                                                              >  2:386698/1‑150 (MQ=255) aaaacggtgttcacaccctgtgcccgcaacgcatgtaccacccactgtgcgccattcatagttagttccccgtcctgaatcttgagaaacagaattttgtgctgttattcattgtctgctcctcggttatgtttttaaggtcaaaaaaaa                                                                                                                                              <  3:236946‑M1/1‑1 (MQ=255)  tCGAACAAGATGCAAGAAAAGACAAAATGACAGCCCTTCTACGAGTGATTAGCCTGGTCGTGATTAGCGTGGTGGTGATTATTATCCCACCGTGCGGGGCTGCACTTGGACGAGGAAAGGCTTAGAGATCAAGCCTTAACGAACTAAGAc                                                                                                                                             <  2:468301/150‑1 (MQ=255)   cGAACAAGATGCAAGAAAAGACAAAATGACAGCCCTTCTACGAGTGATTAGCCTGGTCGTGATTGGCGTGGTGGTGATTATTATCCCACCGTGCGGGGCTGCACTTGGACGAGGAAAGGCTTAGAGATCAAGCCTTAACGAACTAAGAcc                                                                                                                                            <  2:331569/150‑1 (MQ=255)    gAACAAGATGCAAGAAAAGACAAAATGACAGCCCTTCTACGAGTGATTAGCCTGGTCGTGATTAGCGTGGTGGTGATTATTATCCCACCGTGCGGGGCTGCACTTGGACGAGGAAAGGCTTAGAGATCAAGCCTTAACGAACTAAGAccc                                                                                                                                           <  2:439888/150‑1 (MQ=255)     cggtgttcacaccctgtgcccgcaacgcatgtaccacccactgtgcgccattcatagttagttccccgtcctgaatcttgagaaacagaattttgtgctgttattcattgtctgctcctcggttatgtttttaaggtcaaaaaaaAcccc                                                                                                                                          >  1:346103‑M1/146‑150 (MQ=255)                ccctgtgcccgcaacgcatgtaccacccactgtgcgccattcatagttagttccccgtcctgaatcttgagaaacagaattttgtgctgttattcattgtctgctcctcggttatgtttttaaggtcaaaaaagACCCCCGCACCGAAAg                                                                                                                               <  2:212318‑M1/16‑1 (MQ=255)                aGAAAAGACAAAATGACAGCCCTTCTACGAGTGATTAGCCTGGTCGTGATTAGCGTGGTGGTGATTATTATCCCACCGTGCGGGGCTGCACTTGGACGAGGAAAGGCTTAGAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAg                                                                                                                               >  4:119463/1‑150 (MQ=255)                    gtgcccgcaacgcatgtaccacccactgtgcgccattcatagttagttccccgtcctgaatcttgagaaacagaattttgtgctgttattcattgtctgctcctcggttatgtttttaaggtcaaaaaaaACCCCCGCACCGAAAGGTcc                                                                                                                           >  2:64130‑M1/131‑150 (MQ=255)                      gACAAAATGACAGCCCTTCTACGAGTGATTAGCCTGGTCGTGATTAGCGTGGTGGTGATTATTATCCCACCGTGCGGGGCTGCACTTGGACGAGGAAAGGCTTAGAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCgg                                                                                                                         >  1:421684/1‑150 (MQ=255)                       aCAAAATGACAGCCCTTCTACGAGTGATTAGCCTGGTCGTGATTAGCGTGGTGGTGATTATTATCCCACCGTGCGGGGCTGCACTTGGACGAGGAAAGGCTTAGAGATCAAGCCTTAACGAACTAAGACCCACGCACCGAAAGGTCCggg                                                                                                                        <  1:386698/150‑1 (MQ=255)                          gcaacgcatgtaccacccactgtgcgccattcatagttagttccccgtcctgaatcttgagaaacagaattttgtgctgttattcattgtctgctcctcggctatgtttttaaggccaaaaaaaACACCCGCACCGAAAGGTCcgggggg                                                                                                                     <  1:128206‑M1/26‑2 (MQ=255)                                 cgCCCTTCTACGAGTGATTAGCCTGGTCGTGATTAGCGTGGTGGTGATTATTATCCCACCGTGCGGGGCTGCACTTGGACGAGGAAAGGCTTAGAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGGGttattttt                                                                                                              >  4:236946/2‑150 (MQ=255)                                   gtaccacccactgtgcgccattcatagttagttccccgtcctgaatcttgagaaacagaattttgtgctgttattcattgtctgctcctcggttatgtttttaaggtcaaaaaagACCCCCGCACCGAAAGGTCCGGGGGTTTTTTTTGa                                                                                                            >  1:212318‑M1/116‑150 (MQ=255)                                      ttCTACGAGTGATTAGCCTGGTCGTGATTAGCGTGGTGGTGATTATTATCCCACCGTGCGGGGCTGCACTTGGACGAGGAAAGGCTTAGAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGGGTTTTTTTTGACCt                                                                                                         <  3:119463/150‑1 (MQ=255)                                                                      ccgtcctgaatcttgagaaacagaattttgtgctgttattcattgtctgctcctcggttatgtttttaaggtcaaaaaaaACCCCCGCACCGAAAGGTCCGGGGGtttttttt                                                                                                              >  4:24179‑M1/81‑113 (MQ=255)                                                                      ccgtcctgaatcttgagaaacagaattttgtgctgttattcattgtctgctcctcggttatgtttttaaggtcaaaaaaaACCCCCGCACCGAAAGGTCCGGGGGtttttttt                                                                                                              <  3:24179‑M1/33‑1 (MQ=255)                                                                       cgtcctgaatgttgaggaacggaattttgtgctgttattcattgtctgctcctcggttatgtttttaaggtcaaaaaaaACCCCCGCACCGAAAGGTCCGGGGGTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAATGAATaa                                                                        >  2:128206‑M1/80‑150 (MQ=255)                                                                          tggtgATTATTATCCCACCGTGCGGGGCTGCACTTGGACGAGGAAAGGCTTAGAGATCAAGCCTTAACGAACTAAGAc                                                                                                                                             >  2:24275/1‑78 (MQ=255)                                                                          tggtgATTATTATCCCACCGTGCGGGGCTGCACTTGGACGAGGAAAGGCTTAGAGATCAAGCCTTAACGAACTAAGAc                                                                                                                                             <  1:24275/78‑1 (MQ=255)                                                                                                cGGGGCTGCACTTGGACGAGGAAAGGCTTAGAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGGGTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAATGAATAACAGCACAAAATTCTGTTTCTCAAGa                                               <  2:421684/150‑1 (MQ=255)                                                                                                     cTGCACTTGGACGAGGAAAGGCTTAGAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGGGTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAATGAATAACAGCACAAAATTCTGttt                                                      >  2:26803/1‑138 (MQ=255)                                                                                                     cTGCACTTGGACGAGGAAAGGCTTAGAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGGGTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAATGAATAACAGCACAAAATTCTGttt                                                      <  1:26803/138‑1 (MQ=255)                                                                                                        cACTTGGACGAGGAAAGGCTTAGAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGGGTTTTTTTTGACCTTAAAAACATAACCga                                                                                          >  2:77800/1‑99 (MQ=255)                                                                                                        cACTTGGACGAGGAAAGGCTTAGAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGGGTTTTTTTTGACCTTAAAAACATAACCga                                                                                          <  1:77800/99‑1 (MQ=255)                                                                                                                ttgtctgctcctcggttatgtttttaaggtcaaaaaaaACCCCCGCACCGAAAGGTCCgggg                                                                                                                       <  2:214634‑M1/24‑1 (MQ=255)                                                                                                                ttgtctgctcctcggttatgtttttaaggtcaaaaaaaACCCCCGCACCGAAAGGTCCgggg                                                                                                                       >  1:214634‑M1/39‑62 (MQ=255)                                                                                                                cGAGGAAAGGCTTAGAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGGGTGTTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAATGAATAACAGCACAAAATTCTGTTTCTCAAGATTCAGGACGGGGaact                               >  1:475958/1‑150 (MQ=255)                                                                                                                              gagaTCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGGGTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAATGAATAACAGCACAAAATTCTGTTTCTCAAGATTCAGGACGGGGAACTAACTAt                         <  1:458603/142‑1 (MQ=255)                                                                                                                              gagaTCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGGGTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAATGAATAACAGCACAAAATTCTGTTTCTCAAGATTCAGGACGGGGAACTAACTAt                         >  2:458603/1‑142 (MQ=255)                                                                                                                                 aTCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGGGTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAATGAATAACAGCACAAAATTCTGTTTCTCAAGATTCAGGACGGGGAACTAACTATGAATGGCGcac              >  4:91734/1‑150 (MQ=255)                                                                                                                                           aaCGAACTAAGACCCCCGCACCGAAAGGTCCGGGGGTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAATGAATAACAGCACAAAATTCTGTTTCTCAAGATTCAGGACGGGGAACTAACTATGAATGGCGCACAGTGggtggt    >  2:88851/1‑150 (MQ=255)                                                                                                                                           aaCGAACTAAGACCCCCGCACCGAAAGGTCCGGGGGTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAATGAATAACAGCACAAAATTCTGTTTCTCAAGATTCAGGACGGGGAACTAACTATGAATGGCGCACAGTGggtggt    >  2:194086/1‑150 (MQ=255)                                                                                                                                             cGAACTAAGACCCCCGCACCGAAAGGTCCGGGGGTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAATGAATAACAGCACAAAATTCTGTTTCTCAAGATTCAGGACGGGGAACTAACTATGAATGGCGCACAGTGGGTGGTAc  >  4:14598/1‑150 (MQ=255)                                                                                                                                                      |                                                                                                                                             TTCGAACAAGATGCAAGAAAAGACAAAATGACAGCCCTTCTACGAGTGATTAGCCTGGTCGTGATTAGCGTGGTGGTGATTATTATCCCACCGTGCGGGGCTGCACTTGGACGAGGAAAGGCTTAGAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGGGTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAATGAATAACAGCACAAAATTCTGTTTCTCAAGATTCAGGACGGGGAACTAACTATGAATGGCGCACAGTGGGTGGTAC  >  NC_000913_3_bme_pgi/3950295‑3950584

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.