New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913_3_bme_pgi	= 3390511	49 (1.240)	6 (0.160)	5/276	4.9	11.4%	intergenic (+62/+72)	aaeR/tldD	transcriptional regulator for aaeXAB operon/putative peptidase
	?	NC_000913_3_bme_pgi	= 3390524	46 (1.200)					intergenic (+75/+59)	aaeR/tldD	transcriptional regulator for aaeXAB operon/putative peptidase
Rejected: Coverage evenness skew score above cutoff.
Rejected: Frequency below/above cutoff threshold.

AGGTCGTGATCAACTCGCTGACGGATTATTTTGTTGAGGTCGGTAAATTGTTTCAGGAGATGCACGGGCGCGGGAAAGAGAAGTAATTCATATTGTACTGTTACGTTGTACAAACCTGTGCCAACGGGTTCCCCTCACCCTAACCCTCACCCCAAAGGGGCGAGGGGACCGTTCGTGCACGTAGAAAGATTAATTATCCTTCTGAAAATAGTGAAATTATTACGCAGTACCGCCAACAGTCAGGTTATCGACTTTCAACGTTGGCTGGCCCACGCCAACCGGCAAACTTTGCC  >  NC_000913_3_bme_pgi/3390364‑3390656                                                                                                                                                    |                                                                                                                                                  aGGTCGTGATCAACTCGCTGACGGATTATTTTGTTGAGGTCGGTAAATTGTTTCAGGAGATGCACGGGCGCGGGAAAGAGAAGTAATTCATATTGTACTGTTACGTTGTACAAACCTGTGCCAACGGGTTCCCCTCACCCTAACCCTcg                                                                                                                                                  >  2:346959/1‑148 (MQ=255)   gTCGTGATCAACTCGCTGACGGATTATTTTGTTGAGGTCGGTAAATTGTTTCAGGAGATGCACGGGCGCGGGAAAGAGAAGTAATTCATATTGTACTGTTACGTTGTACAAACCTGTGCCAACGGGTTCCCCTCACCCTAACCCTCAcc                                                                                                                                                >  2:300856/1‑149 (MQ=255)         attaaCGCGCTGACGTATTATTTTGTTGAGGGCGGTAAATTGTTGTAGGAGATTCACGGTCGCGGGAAAGAGAAGTAATACATATTGTACTGTTACGTTGTACAAACCTGTGCCAACGTGGTCCCCTCCCCCTAAGCCTCACCCCAAAg                                                                                                                                          <  2:196285/146‑1 (MQ=255)          tCAACTCGCTGACGGATTATTTTGTTGAGGTCGGTAAATTGTTTCAGGAGATGCACGGGCGCGGGAAAGAGAAGTAATTCATATTGTACTGTTACGTTGTACAAACCTGTGCCAACGGGTTCCCCTCACCCTAACCCTCACCCCAAAgg                                                                                                                                         <  2:58140/149‑1 (MQ=255)                cGCTGCGGGATTATTTTGTTGAGGTCGGTAAATTGTTTCAGGAGATGCACGGGCGCGGGAAAGAGAAGTAATTCATATTGTACTGTTACGTTGTACAAACCTGTGCCAACGGGTTCCCCTCACCCTAACCCTCACCCCAAAGGGGCGAg                                                                                                                                   <  2:308519/149‑1 (MQ=255)                  cTGACGGATTATTTTGTTGAGGTCGGTAAATTGTTTCAGGAGATGCACGGGCGCGGGAAAGAGAAGTAATTCATATTGTACTGTTACGTTGTACAAACCTGTGCCAACGGGTTCCCCTCACCCTAACCCTCACCCCAAAGGGGCGAggg                                                                                                                                 <  2:364809/149‑1 (MQ=255)                     aCGGATTATTTTGTTGAGGTCGGTAAATTGTTTCAGGAGATGCACGGGCGCGGGAAAGAGAAGTAATTCATATTGTACTGTTACGTTGTACAAACCTGTGCCAACGGGTTCCCCTAACCCTAACCCTCACCCCAAAGGGCCGAGGGGAc                                                                                                                              >  2:160021/1‑149 (MQ=255)                      cGGATTATTTTGTTGAGGTCGGTAAATTGTTTCAGGAGATGCACGGGCGCGGGAAAGAGAAGTAATTCATATTGTACTGTTACGTTGTACAAACCTGTGCCAACGGGTTCCCCTCACCCTAACCCTCGCCCCTTTGGGGTGAgggttag                                                                                                                             <  1:346959/149‑5 (MQ=255)                        gATTATTTTGTTGAGGTCGGTAAATTGTTTCAGGAGATGCACGGGCGCGGGAAAGAGAAGTAATTCATATTGTACTGTTACGTTGTACAAACCTGTGCCAACGGGTTCCCCTCACCCTAACCCTCACCCCAAAGGGGCGAGGGGACCGt                                                                                                                           >  1:150497/1‑149 (MQ=255)                             ttttGTTGAGGTCGGTAAATTGTTTCAGGAGATGCACGGGCGCGGGAAAGAGAAGTAATTCATATTGTACTGTTACGTTGTACAAACCTGTGCCAACGGGTTCCCCTCACCCTAACCCTCACCCCAAAGGGGCGAGGGGACCGTTCGTg                                                                                                                      >  1:409012/1‑149 (MQ=255)                             ttttGTTGAGGTCGGTAAATTGTTTCAGGAGATGCACGGGCGCGGGAAAGAGAAGTAATTCATATTGTACTGTTACGTTGTACAAACCTGTGCCAACGGGTTCCCCTCACCCTAACCCTCACCCCAAAGGGGCGAGGGGACCGTTCGTg                                                                                                                      >  1:621887/1‑149 (MQ=255)                                tgttgAGGTCGGTAAATTGTTTCAGGAGATGCACGGGCGCGGGAAAGAGAAGTAATTCATATTGTACTGTTACGTTGTACAAACCTGTGCCAACGGGTTCCCCTCACCCTAACCCTCACCCCAAAggg                                                                                                                                        <  1:31136/128‑1 (MQ=255)                                tgttgAGGTCGGTAAATTGTTTCAGGAGATGCACGGGCGCGGGAAAGAGAAGTAATTCATATTGTACTGTTACGTTGTACAAACCTGTGCCAACGGGTTCCCCTCACCCTAACCCTCACCCCAAAggg                                                                                                                                        >  2:31136/1‑128 (MQ=255)                                tgttgAGGTCGGTAAATTGTTTCAGGAGATGCACGGGCGCGGGAAAGAGAAGTAATTCATATTGTACTGTTACGTTGTACAAACCTGTGCCAACGGGTTCCCCTCACCCTAACCCTCACCCCAAAGGGGCGAGGGGACCGTTCGTGCAc                                                                                                                   <  2:164151/149‑1 (MQ=255)                                   tgATGTCGGCAAGTTGTTTCAGGAGATGCACGGGCGTGGGAAAAAGAAGGAATTCATATTGTACTGTTACGTTGTACAAACCTGTGCCAACGGGTCCCCCTCACCCTAACCCTCACCCCAAAGGGGCGAGGGGACCGTTCGTGCACGTa                                                                                                                <  1:317793/149‑1 (MQ=38)                                       gTCGGTAAATTGTTTCAGGAGATGCACGGGCGCGGGAAAGAGAAGTAATTCATATTGTACTGTTACGTTGTACAAACCTGTGCCAACGGGTTCCCCTCACCCTAACCCTCACCCCAAAGGGGCGAGGGGACCGTTCGTGCACGTAGaaa                                                                                                            >  2:494730/1‑149 (MQ=255)                                         cGGTAAAGTGTTTCAGGAGATGCACGGGCGCGGGAAAGAGAAGTAATTCATATTGTACTGTTACGTTGTACAAACCTGTGCCAACGGGTTCCCCTCACCCTAACCCTCACCCCAAAGGGGCGAGGGGACCGTTCGTGCACGTAGAAAGa                                                                                                          <  2:587785/149‑1 (MQ=255)                                            tAAATTGTTTCAGGAGATGCACGGGCGCGGGAAAGAGAAGTAATTCATATTGTACTGTTACGTTGTACAAACCTGTGCCAACGGGTTCCCCTCACCCTAACCCTCACCCCAAAGGGGCGAGGGGACCGTTCGTGCACGTAGAAAGatta                                                                                                       >  1:227474/1‑149 (MQ=255)                                             aaaTTGTTTCAGGAGATGCACGGGCGCGGGAAAGAGAAGTAATTCATATTGTACTGTTACGTTGTACAAACCTGTGCCAACGGGTTCCCCTCACCCTAACCCTCACCCCAAAGGGGCGAGGGGACCGTTCGTGCACGTAGAAAGattaa                                                                                                      <  1:615943/149‑1 (MQ=255)                                                ttGTTTCAGGAGATGCACGGGCGCGGGAAAGAGAAGTAATTCATATTGTACTGTTACGTTGTACAAACCTGTGCCAACGGGTTCCCCTCACCCTAACCCTCACCCCAAAGGGGCGAGGGGCCCGTTCGTGCACGTAGAAAGattaatta                                                                                                   >  2:460034/1‑149 (MQ=255)                                                ttGTTTCAGGAGATGCACGGGCGCGGGAAAGAGAAGTAATTCATATTGTACTGTTACGTTGTACAAACCTGTGCCAACAGGTTCCCCGCAGCCTAACCCTCACCCCAAAGGGGCAAGGGGACCGTTCGTGCACGTAGAAAGattaatta                                                                                                   >  2:124398/1‑149 (MQ=255)                                                  gTTTCAGGAGATCCACGGGCGCGGGAAAGAGAAGTAATTCATATTGTACTGTTACGTTTTACAAACCTGTGCCAACGGGTTCCCCTCACCCTCACCCTAACCACATAAGCCCGAGGGGACCGTTCGTGCACGTAGAAAGATTAATTATc                                                                                                 <  2:213838/149‑1 (MQ=255)                                                     tCAGGAGATGCACGGGCGCGGGAAAGAGAAGTAATTCATATTGTACTGTTACGTTGTACAAACCTGTGCCAACGGGTTCCCCTCACCCTAACCCTCACCCCAAAGGGGCGAGGGGACCGTTCGTGCACGTAGAAAGATTAATTATCCtt                                                                                              <  1:300856/149‑1 (MQ=255)                                                      cAGGAGATGCACGGGCGCGGGAAAGAGAAGTAATTCATATTGTACTGTTACGTTGTACAAACCTGTGCCAACGGGTTCCCCTCACCCTAACCCTCACCCCAAAGGGGCGAGGGGACCGTTGGTGCACGTAGAAAGATTAATTATCCTTc                                                                                             >  1:54103/1‑149 (MQ=255)                                                       aGGAGATGCACGGGCGCGGGAAAGAGAAGTAATTCATATTGTACTGTTACGTTGTACAAACCTGTGCCAACGGGTTCCCCTCACCCTAACCCTCACCCCAAAGGGGCGAGGGGACCGTTCGTGCACGTAGAAAGATTAATTATCCTTCt                                                                                            <  2:285276/149‑1 (MQ=255)                                                                   ggCGCGGGAAAGAGAAGTAATTCATATTGTACTGTTACGTTGTACAAACCTGTGCCAACGGGTTCCCCTCACCCTAACCCTCACCCCAAAGGGGCGAGGGGACCGTTCGTGCACGTAGAAAGATTAATTATCCTTCTGAAAATAGTGaa                                                                                <  2:227474/149‑1 (MQ=255)                                                                        gggAAAGAGAAGTAATTCATATTGTACTGTTACGTTGTACAAACCTGTGCCAACGGGTTCCCCTCACCCTAACCCTCACCCCAAAGGGGCGAGGGGACCGTTCGTGCACGTAGAAAGATTAATTATCCTTCTGAAAATAGTGAAattat                                                                           >  2:331546/1‑149 (MQ=255)                                                                               agaAGTAATTCATATTGTACTGTTACGTTGTACAAACCTGTGCCAACGGGTTCCCCTCACCCTAACCCTCACCCCAAAGGGGCGAGGGGACCGTTCGTGCACGTAGaaa                                                                                                            >  2:621678/1‑109 (MQ=255)                                                                               agaAGTAATTCATATTGTACTGTTACGTTGTACAAACCTGTGCCAACGGGTTCCCCTCACCCTAACCCTCACCCCAAAGGGGCGAGGGGACCGTTCGTGCACGTAGaaa                                                                                                            <  1:621678/109‑1 (MQ=255)                                                                                gaAGTAATACATATTGTACTGTTACGTTGTACAAACCTGTGCCAACGGGTTCCCCTCACCCTAACCCTCACCCCAAAGGGGCGAGGGGACCGTTCGTGCACGTAGAAAGATTAATTATGCTTTTGAGAATAGTGAAATTATTTCGAAGt                                                                   >  1:181768/1‑149 (MQ=38)                                                                                      atgcATATTGTACTGTTACGGTGTACAAACCTGCGCCAACGGGTCCCCCGCACCCTAACCCTCACCCCAAAGGGGCGAGGGAACCGTTCGTGCACGTAGAAAGATTAATTATCCTTCTGAAAATAGTGAAATTATTACGCAGTACCGcc                                                             <  1:12268/146‑1 (MQ=255)                                                                                        tCATATTGTACTGTTACGTTGTACAAACCTGTGCCAACGGGGTCCCCTCGCCCTAACCCGCACCCGAGAAGGGCGAGGGGACCGTTCGTGCAAGTAAAAAGATTTATTATCCTTATCAACATAGTGAAATTATTACACAGTaatcccca                                                           >  2:282830/1‑142 (MQ=255)                                                                                        tCATATTGTACTGTTACGATGTACAAACCTGTGACAACGGGTTCCCCTCACCCTAACCCTCACCCCAAAGGGGCAAGGGGACCGTTCGTGCACGTAGAAAGATTAATTATCCTTCTGAGAATAGTTAAATTATTACGCAGTACCGCCaa                                                           >  2:340609/1‑149 (MQ=39)                                                                                        tCACATTGTACGGTTACGTTGTCCAAACCTGTGCCAACGGGTTCCCCCCACCCTAACCCTCACCCCAAAGGGGCGAGGGGTGCATTCGTGCACGGAGAAAGATTAATTATCCTTCTGAAAATAGTTAAATTATTACGCAGTACCGCCaa                                                           <  1:307469/149‑1 (MQ=255)                                                                                          atatTGTACTGTTACGTTGTACAAACCTGTGCCAACGGGTTCCCCTCACCCTAACCCTCACCCCAAAGGGGCGAGGGGACCGTTCGTGCACGTAGAAAGATTAATTATCCTTCTGAAAATAGTGAAATTATTACGCAGTACCGCCAACa                                                         <  2:621887/149‑1 (MQ=255)                                                                                               gTACTGTTACGTTGTACAAACCTGTGCCAACGGGTTCCCCTCACCCTAACCCTCACCCCAAAGGGGCGAGGGGACCGTTCGTGCACGTAGAAAGATTAATTATCCTTCTGAAAATAGTGAAATTATTACGCAGTACCGCCAACAGTCAg                                                    <  1:376308/149‑1 (MQ=255)                                                                                                     ttACGTTGTACAAACCTGTGCCAACGGGTTCCCCTCACCCTAACCCTCACCCCAAAGGGGCGAGGGGACCGTTCGTGCACGTAGAAAGATTAATTATCCTTCTGAAAATAGTGAAATTATTACGCAGTACCGCCAACAGTCAGGTTATc                                              <  1:331546/149‑1 (MQ=255)                                                                                                       aCGTTGTACAAACCTGTGCCAACGGGTTCCCCTCACCCTAACCCTCACCCCAAAGGGGCGAGGGGACCGTTCGTGCACGTAGAAAGATTAATTATCCTTCTGAAAATAGTGAAATTATTACGCAGTACCGCCAACAGTCAGGTTATCGa                                            <  2:409012/149‑1 (MQ=255)                                                                                                             tACAAACCTGTGCCAACGGGTTCCCCTCACCCTAACCCTCACCCCAAAGGGGCGAGGGGACCGTTCGTGCACGTAGAAAGATTAATTATCCTTCTGAAAATAGTGAAATTATTACGCCGTACCGCCAACAGTCAGGTTATCGACTTTCa                                      >  2:57940/1‑149 (MQ=255)                                                                                                             tACAAACCTGTGCCAACGGGTTCCCCTCACCCTAACCCTCACCCCAAAGGGGCGAGGGGACCGTTCGTGCACGTAGAAAGATTAATTATCCTTCTGAAAATAGTGAAATTATTACGCAGTACCGCCAACAGTCAGGTTATCGACt                                          <  1:45737/145‑1 (MQ=255)                                                                                                             tACAAACCTGTGCCAACGGGTTCCCCTCACCCTAACCCTCACCCCAAAGGGGCGAGGGGACCGTTCGTGCACGTAGAAAGATTAATTATCCTTCTGAAAATAGTGAAATTATTACGCAGTACCGCCAACAGTCAGGTTATCGACt                                          >  2:45737/1‑145 (MQ=255)                                                                                                             tACAAACCTGTGCCAACGGGTTCCCCTCACCCTAACCCTCACCCCAAAGGGGCGAGGGGACCGTTCGTGCACGTAGAAAGATTAATTATCCTTCTGAAAATAGTGAAATTATTACGCAGTACCGCCAACAGTCAGGTTATCGACTTTCa                                      >  2:68402/1‑149 (MQ=255)                                                                                                                       tgCCAACGGGTTCCCCTCACCCTAACCCTCACCCCAAAGAGGCGAGGGGACCGTTCGTGCCCGTAGAAAGATTAATTATCCTGCTGAAAATAGCGAAATGAGGACGCAGTACCGCCAACAGTCAGGTTATCGACTTTCAACGTcggctg                            >  2:145877/1‑149 (MQ=255)                                                                                                                       tgCCAACGGGTCCCCCTCACCCTAACCCTCACCCCAAAGGGGCGAGGGGACCGTTCGTGCACGTAGAAAGATTAATTACCCTTCTGAAAATAGTGAAATTATTACGCAGTACCGCCAACAGTCAGGTTATCGACTTTCAACGTtggctg                            <  1:221172/149‑1 (MQ=255)                                                                                                                            aCGGGTTCCCCTCACCCTAACCCTCACCCCAAAGGGGCGAGGGGACCGTTCGTGTACGTAGAAAGATTAATTATCCTTCTGAAAATAGTGAAATTATTACGCAGTACCGCCAACAGTCAGGTTATCGACTTTCAACGTTGGCTGGCCCa                       <  1:595455/149‑1 (MQ=255)                                                                                                                                         aCCCTAACCCTCACCCCAAAGGGGCGAGGGGACCGTTCGTGCACGTAGAAAGATTAATTATCCTTCTGAAAGTAGTGAAATTATTGCGCAGTACCGCCAACAGTCAGGTTATCGACTTTCAACGTTGGCTGGCCCACGCCAACCGGCaa          >  1:613880/1‑149 (MQ=255)                                                                                                                                           ccTAACCCTCACCCCAAAGGGGCGAGGGGACCGTTCGTGCACGTAGAAAGATTAATTATCCTTCTGAAAATAGTGAAATTATTACGCAGTACCGCCAACAGTCAGGTTATCGACTTTCAACGTTGGCTGGCCCACGCCAACCGGCAAAc        >  1:402536/1‑149 (MQ=255)                                                                                                                                            cTAACCCTCACCCCAAAGGGGCGAGGGGACCGTTCGTGCACGTAGAAAGATTAATTATCCTTCTGAAAATAGTGAAATTATTACGCAGTACCACCAACAGTCAGGTTATCGACTTTCAACGTTGGCTGGCCCACGCCAACCGGCAAACt       >  1:253617/1‑149 (MQ=255)                                                                                                                                             tAACCCTCACCCCAAAGGGGCGAGGGGACCGTTCGTGCACGTAGAAAGATTAAGTATCCTTCTGAAAATAGTGAAATTATTACGCAGTACCGCCAAAAGTCAGGTTATCGACTTTCAACGCTGGATGCCCCACGCCAACCGCCAAACtt      >  2:326532/1‑149 (MQ=255)                                                                                                                                                 ccTCACCCCAAAGGGGCGAGGGGACCGTTCGTGCACGTAGAAAGATTAATTATCCTTCTGAAAATAGTGAAATTATTACGCAGTACCGCCAACAGTCAGGTTATCGACTTTCAACGTTGGCTGGCCCACGCCAACCGGCAAACTTTGcc  >  1:347974/1‑149 (MQ=255)                                                                                                                                                 ccTCACCCCAAAGGGGCGAGGGGACCGTTCGTGCACGTAGAAAGATTAATTATCCTTCTGAAAATAGTGAAATTATTACGCAGTACCGCCAACAGTCAGGTTATCGACTTTCAACGTTGGCTGCCCCACGCCAACCGGCAAACTTTGcc  >  1:600929/1‑149 (MQ=255)                                                                                                                                                    |                                                                                                                                                  AGGTCGTGATCAACTCGCTGACGGATTATTTTGTTGAGGTCGGTAAATTGTTTCAGGAGATGCACGGGCGCGGGAAAGAGAAGTAATTCATATTGTACTGTTACGTTGTACAAACCTGTGCCAACGGGTTCCCCTCACCCTAACCCTCACCCCAAAGGGGCGAGGGGACCGTTCGTGCACGTAGAAAGATTAATTATCCTTCTGAAAATAGTGAAATTATTACGCAGTACCGCCAACAGTCAGGTTATCGACTTTCAACGTTGGCTGGCCCACGCCAACCGGCAAACTTTGCC  >  NC_000913_3_bme_pgi/3390364‑3390656

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.