New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913_3_bme_pgi	2737519 =	34 (0.740)	4 (0.090)	4/268	6.2	10.7%	intergenic (+24/‑80)	raiA/pheL	cold shock protein associated with 30S ribosomal subunit/pheA gene leader peptide
	?	NC_000913_3_bme_pgi	2737539 =	35 (0.810)					intergenic (+44/‑60)	raiA/pheL	cold shock protein associated with 30S ribosomal subunit/pheA gene leader peptide
Rejected: Coverage evenness skew score above cutoff.
Rejected: Frequency below/above cutoff threshold.

AAGCTGGAACGGCAGCTCAATAAACTGCAGCACAAAGGCGAAGCACGTCGTGCCGCAACATCGGTGAAAGACGCCAACTTCGTCGAAGAAGTTGAAGAAGAGTAGTCCTTTATATTGAGTGTATCGCCAACGCGCCTTCGGGCGCGTTTTTTGTTGACAGCGTGAAAACAGTACGGGTACTGTACTAAAGTCACTTAAGGAAACAAACATGAAACACATACCGTTTTTCTTCGCATTCTTTTTTACCTTCCCCTGAATGGGAGGCGTTTCGTCGTGTGAAACAGAATGCGAAGACG  >  NC_000913_3_bme_pgi/2737391‑2737686                                                                                                                                                     |                                                                                                                                                    aaGCTGGAACGGCAGCTCAATAAACTGCAGCACAAAGGCGAAGCACGTCGTGCCGCAACATCGGTGAAAGACGCCAACTTCGTCGAAGAAGTTGAAGAAGAGTAGTCCTTTATATTGAGTGTATCGCCAACGCGCCTTCGGGCGCGttt                                                                                                                                                     <  1:138809/149‑1 (MQ=255) aaGCTGGAACGGCAGCTCAATAAACTGCAGCACAAAGGCGAAGCACGTCGTGCCGCAACATCGGTGAAAGACGCCAACTTCGTCGAAGAAGTTGAAGAAGAGTAGTCCTTTATATTGAGTGTATCGCCAACGCGCCTTCGGGCGCGttt                                                                                                                                                     <  2:97307/149‑1 (MQ=255)  aGCCGGAACGGCAGCCCAATAAACTGCAGCACAAAGGCGAAGCACGTCGTGCCGCAACATCGGTGAAAGACGCCAACTTCGTCGAAGAAGTTGAAGAAGAGTAGTCCTTTATATTGAGTGTATCGCCAACGCGCCTTCGGGCGCGtttt                                                                                                                                                    <  1:335190/149‑1 (MQ=255)           ggCAGCTCAATAAACTGCAGCACAAAGGCGAAGCACGTCGTGCCGCAACATCGGTGAAAGACGCCAACTTCGTCGAAGAAGTTGAAGAAGAGTAGTCCTTTATATTGAGTGTATCGCCAACGCGCCTTCGGGCGCGTTTTTTGTTGACa                                                                                                                                           <  1:233876/149‑1 (MQ=255)            gCAGCTCAATAAACTGCAGCACAAAGGCGAAGCACGTCGTGCCGCAACATCGGTGAAAGACGCCAACTTCGTCGAAGAAGTTGAAGAAGAGTAGTCCTTTATATTGAGTGTATCGCCAACGCGCCTTCGGGCGCGTTTTTTGTTGACAg                                                                                                                                          >  2:38408/1‑149 (MQ=255)                    aTAAACTGCAGCAGAAAGGCGAAGCAGGTCGTGCCGCAACATCGGTGAAAAACGCCAACTTCGTCGAAGAAGTTGAAGAAGAGTAGTCCTTTATATTGGGTGTATCGCCAACGCGCCTTCGGGCGGGTTTTTTGTTGATAGCGTGaaac                                                                                                                                  >  1:664044/1‑148 (MQ=255)                    aTAAACTGCAGCACAAAGGCGGAGCACGTCGTGCCGCAACATCGGTGAAAGACGCCAACTTCGTCGAAGAAGTTGAAGAAGAGTAGTCCTTTATATTGAGTGTATCGCCAACGCGCCTTCGGGCGCGTTTTTTGTTGACAGCGTGaaaa                                                                                                                                  <  2:558223/149‑1 (MQ=255)                     tAAACTGCAGCACAAAGGCGAAGCACGTCGTGCCGCAACATCGGTGAAAGACGCCAACTTCGTCGAAGAAGTTGAAGAAGAGTAGTCCTTTATATTGAGTGTATCGCCAACGCGCCTTCGGGCGCGTTTTTTGTTGACAGCGTGAAAAc                                                                                                                                 >  1:249179/1‑149 (MQ=255)                            cagcaCAAAGGCGAGGCACGTCGTGCCGCAACATCGGTGAAAGACGCCAACTTCGTCGAAGAAGTTGAAGAAGAGTAGTCCTTTATATTGAGTGTATCGCCAACGCGCCTTCGGGCGCGTTTTTTGTTGACAGCGTGAAAACAGTACgg                                                                                                                          <  1:390016/149‑1 (MQ=255)                             agcaCAAAGGCGAAGCACGTCGTGCCGCAACATCGGTGAAAGACGCCAACTTCGTCGAAGAAGTTGAAGAAGAGTAGTCCTTTATATTGAGTGTATCGCCAACGCGCCTTCGGGCGCGTTTTTTGTTGACAGCGTGAAAACAGTACggg                                                                                                                         <  2:235402/149‑1 (MQ=255)                                    aGGCGAAGCACGTCGTGCCGCAACATCGGTGAAAGACGCCAACTTCGTCGAAGAAGTTGAAGAAGAGTAGTCCTTTATATTGAGTGTATCGCCAACGCGCCTGCGGGCGCGTTTTTTGTTGACAGCGTGAAAACAGTACGGgtactgta                                                                                                                  >  2:38024/1‑149 (MQ=255)                                        gAAGCACGTCGTGCCGCAACATCGGTGAAAGACGCCAACTTCGTCGAAGAAGTTGAAGAAGAGTAGTCCTTTATATTGAGTGTATCGCCAACGCGCCTTCGGGCGCGTTTTTTGTTGACAGCGTGAAAACAGTACGGGTACTGTACTaa                                                                                                              <  1:327372/149‑1 (MQ=255)                                         aaGCACGTCGCGCCACAACATCGGTGAAAGACGCCAACTTCGTCGAAGAAGTTGAAGAAGAGTAGTCCTTTATATTGAGTGTATCGCCAACGCGCCTTCGGGCGCGTTTTTTGTTGACAGCGTGAAAACAGTACGGGTACTGTACTaaa                                                                                                             <  2:329883/149‑1 (MQ=255)                                           gCACGTCGTGCCGACACATCGGTGAAAGACGCCAACTTCGTCGAAGAAGTTGAAGAAGAGTAGTCCTTTATATTGAGTGTATCGCCAACGCGCCTTCGGGCGCGTTTTTTGTTGACAGCGTGAAAACAGTACGGGTACTGTACTAAAGt                                                                                                           <  1:627703/149‑1 (MQ=255)                                                            aTCGGTGAAAGACGCCAACTTCGTCGAAGAAGTTGAAGAAGAGTAGTCCTTCATATTGAGTGTATCGCCAACGCGCCTTCGGGCGCGTTTTTTGTTGACAGCGTGAAAACAGTACGGGTACTGTACTAAAGTCACTTAAGGaaacaaac                                                                                          <  1:603880/149‑1 (MQ=255)                                                             gctgTGAAAGACGCCAACTTCGTCGAAGAAGTTGAAGAAGAGTAGTCCTTTATATTGAGTGTATCGCCAACGCGCCTTCGGGCGCGTTTTTTGTTGACAGCGTGAAAACAGTACGGGTACTGTACTAAAGTCACTTAAGGAaacaaaca                                                                                         <  1:181425/146‑1 (MQ=255)                                                                gTGAAAGACGCCAACTTCGTCGAAGAAGTTGAAGAAGAGTAGTCCTTTATATTGAGTGTATCGCCAACGCGCCTTCGGGCGCGTTTTTTGTTGACAGCGTGAAAACAGTACGGGTACTGTACTAAAGTCACTTAGGGAAACAAACATGa                                                                                      >  1:192181/1‑149 (MQ=255)                                                                     aGACGCCAACTTCGTCGAAGAAGTTGAAGAAGAGTAGTCCTTTATATTGAGTGTATCGCCAACGCGCCTTCGGGCGCGTTTTTTGTTGACAGCGTGAAAACAGTACGGGTACTGTACTAAAGTCACTTAAGGAAACAAACATGAAacac                                                                                 >  1:383895/1‑149 (MQ=255)                                                                         gCCAACTTCGTCGAAGAAGTTGAAGAAGAGTAGTCCTTTATATTGAGTGTATCGCCAACGAGCCTTCGGGCGCGTTTTTTGTTGACAGCGTGAAAACAGTACGGGTACTGTACTCAAGTCCCTTAAGGAAACAAACATGAACCACATAc                                                                             >  2:306671/1‑149 (MQ=255)                                                                                tcgtcgAAGAAGTTGAAGAAGAGTAGTCCTTTATATTGAGTGTATCGCCAACGCGCCTTCGGGCGCGTTTTTTGTTGACAGCGTGAAAACAGTACGGGTACTGTACTAAAGTCACTTAAGGAAACAAACATGAAACACATACCGttttt                                                                      <  1:676696/149‑1 (MQ=255)                                                                                    cgAAGAAGTTGAAGAAGAGTAGTCCTTTATATTGAGTGTATCGCCAACGCGCCTTCGGGCGCGTTTTTTGTTGACAGCGTGAAAACAGTACGGGTACTGTACTAAAGTCACTTAAGGAAACAAACATGAAACACATACCGTTTttcttc                                                                  <  1:232434/149‑1 (MQ=255)                                                                                                  aagaGTAGCCCTTGATATTGATTGTATCGGCAACGCGCCTTCGGGCGCGTTGTTTGTTGACAGCGTTAAAACAGTACGGGTACTGTACTACAGTCACTTAAGGTAACAAACATGAATCACATATCGCTTCTCTTCGCATTCTTTTTTAc                                                    <  2:354410/149‑1 (MQ=255)                                                                                                         gtCCTTTATATTGAGTGTATCGCCAACGCGCCTTCGGGCGCGTTTTTTGTTGAGAGCGTGAAAACAGTACGGGTACTGTACTAAAGTCACTTAAGGAAACAAACATGAAACACATACCGTTTTTCTTCGCATTCTTTTTTACCTTcccc                                             >  1:529674/1‑149 (MQ=255)                                                                                                          tCCTTTATATTGAGTGTATCGCCAACGCGCCTTCGGGCGCGTTTTTTGTTGACAGCGTGAAAACAGTACGGGTACTGTACTAAAGTCACTTAAGGAAACAAACATGAAACACATACCGTTTTTCTTCGCATTCTTTTTTACCTTCCCCt                                            <  1:173557/149‑1 (MQ=255)                                                                                                            cTTTATATTGAGTGTATCGCCAACGCGCCTTCGGGCGCG‑TTTTTGTTGACAGCGTGAAAACAGTACGGGTACTGTACTAAAGTCACTTAAGGAAACAAACATGAAACACATACCGTTTTTCTTCGCATTCTTTTTTACCTTCCCCTGaa                                         <  2:422830/149‑1 (MQ=255)                                                                                                                 tatTGAGTGTACCGCCAACGCGCCTTCGGGCGCGTTTTTTGTTGACAGCGTGAAAACAGTACGGGTACTGTACTAAAGTCACGTAAGGAAACAAACATGAAACACATACCGTTTTTCTTCGCATTCTTTTTTACCTTCCCCTGAATggg                                     <  1:141172/149‑1 (MQ=255)                                                                                                                     taGTGTATCGCCAACGCGCCTTCGGGCGCGTTTTTTGTTGACAGCGTGAAAACAGTACGGGTACTGTACTAAAGTCACTTAAGGAAACAAACATGAAACACATACCGTTTTTCTTCGCATTCTTTTTTACCTTCCCCTGAATGGGAGGc                                 <  1:541913/148‑1 (MQ=255)                                                                                                                          tATCGCCAACGCGCCTTCGGGTGTGTTTTTTGTTGACAGCGTGAAACCAGGACGGGTACTGTACTTAAGTCTCTTAACGAAACAAACATGAATCACATCCCGTTTTGCGTCGCATTTTTTTTTACCTTTTCCTGAATGGGAGGCGTTtc                            <  2:432221/149‑1 (MQ=255)                                                                                                                            tCGCCAACGCGCCTTCGGGCGCGTTTTTTGTTGACAGCGTGAAAACAGTACGGGTACTGTACTAAAGTCACTTAAGGAAACAAACATGAAACACATACCGTTTTTCTTCGCATTCTTTTTTACCTTCCCCTGAATGGGAGGCGTTtcgt                          >  1:740672/1‑149 (MQ=255)                                                                                                                               ccAACGCGCCTTCGGGCGCGTTTTTTGTTGACAGCGTGAAAACAGTACGGGTACTGTACTAAAGTCACTTAAGGAAACAAACATGAAACACATACCGTTTTTCTTCGCATTCTTTTTTACCTTCCCCTGAATGGGAGGCGTTTcgtcgt                       <  1:26280/149‑1 (MQ=255)                                                                                                                                  aCGCGCCTTCGGGCGCGTTTTTTGTTGACAGCGTGAAAACAGTACGGGTACTGTACTAAAGTCACTTAAGGAAACAAACATGAAACACATACCGTTTTTCTTCGCATTCTTTTTTACCTTCCCCTGAATGGGCGGCGTTTCGTCGtgtg                    >  1:129780/1‑149 (MQ=255)                                                                                                                                   cgcgCCTTCGGGCGCGTTTTTTGTTGACAGCGTGAAAACAGTACGGGTACTGTACTAAAATCACGTAAGGACACAAACATGAAACACATACCGTTTTTCTTCGCATTCTTTTTTACCTTCCCCTGAATGTGAGGCTTTTCGTCgtgtta                   >  2:348099/1‑147 (MQ=255)                                                                                                                                       ccTTCGGGCGCGTTTTTTGTTGACAGCGTGAAAACAGTACGGGTACTGTACTAAAGTCACTTAAGGAAACAAACATGAAACACATACCGTTTTTCTTCGCATTCTTTTTTACCTTCCCCTGAATGGGAGGCGTTTCGTCGTGTGAAACa               >  1:329152/1‑149 (MQ=255)                                                                                                                                        cTTCGGGCGCGTTTTTTGTTGACAGCGTGAAAACAGTACGGGTACTGTACTAAAGTCACTTAAGTAAACAAACATGAAACACATACCGTTTTTCTTCGCATTCTTTTTTACCTTCCCCTGAATGGGAGGCGTTTCGTCGTGTGAAACAg              >  1:703086/1‑149 (MQ=255)                                                                                                                                           cGGGCGCGTTTTTTGTTGACAGCGTGAAAACAGTACGGGTACTGTACTAAAGTCACTTAAGGAAACAAACATGAAACACATACCGTTTTTCTTCGCATTCTTTTTTACCTTCCCCTGAATGGGAGGCGTTTCGTCGTGTGAAACAGAAt           <  2:122255/149‑1 (MQ=255)                                                                                                                                                    tttttGTTGACAGCGTGAAAACAGTACGGGTACTGTACTAAAGTCACTTAAGGAAACAAACATGAAACACATACCGTTTTTCTTCGCATTCTTTTTTACCTTCCCCTGAATGGGAGGCGTTTCGTCGTGTGAAACAGAATGCGAAGACg  <  1:505701/149‑1 (MQ=255)                                                                                                                                                     |                                                                                                                                                    AAGCTGGAACGGCAGCTCAATAAACTGCAGCACAAAGGCGAAGCACGTCGTGCCGCAACATCGGTGAAAGACGCCAACTTCGTCGAAGAAGTTGAAGAAGAGTAGTCCTTTATATTGAGTGTATCGCCAACGCGCCTTCGGGCGCGTTTTTTGTTGACAGCGTGAAAACAGTACGGGTACTGTACTAAAGTCACTTAAGGAAACAAACATGAAACACATACCGTTTTTCTTCGCATTCTTTTTTACCTTCCCCTGAATGGGAGGCGTTTCGTCGTGTGAAACAGAATGCGAAGACG  >  NC_000913_3_bme_pgi/2737391‑2737686

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