New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913_3_bme_pgi	2729475 =	38 (0.850)	3 (0.070)	3/266	6.4	10.8%	intergenic (‑31/+141)	gltW/rrsG	tRNA‑Glu/16S ribosomal RNA of rrnG operon
	?	NC_000913_3_bme_pgi	= 4168363	14 (0.340)					intergenic (+163/‑9)	rrsB/gltT	16S ribosomal RNA of rrnB operon/tRNA‑Glu
Rejected: Coverage evenness skew score above cutoff.
Rejected: Frequency below/above cutoff threshold.

TGAGATATTGAGGTCGGCGACTTTCACTCACAAACCAGCAAGTGGCGTCCCCTAGGGGATTCGAACCCCTGTTACCGCCGTGAAAGGGCGGTGTCCTGGGCCTCTAGACGAAGGGGACACGAAAATTGCTTATCACGCG‑TTGCGTGATATTTTCGTGTAGGGTGAGCTTTCATTAATAGAAAGCGAACGGCCTTATTCTCTTCAGCCTCACTCCCAACGCGTAAACGCCTTGCTATTCACTTTTCATCAGACAATCTGTGTGAGCACTGCAAAGTACGCTTCTTTAAGGTAAGGAGG  >  NC_000913_3_bme_pgi/2729327‑2729623                                                                                                                                                      |                                                                                                                                                     tGAGATATTGAGGTCGGCGACTTTCACTCACAAACCAGCAAGTGGCGTCCCCTAGGGGATTCGAACCCCTGTTACCGCCGTGAAAGGGCGGTGTCCTGGGCCTCTAGACGAAGGTGACACGAAAATTGCTTATCACGCG‑TTGCGTGata                                                                                                                                                      <  2:745106/149‑1 (MQ=17)  gagaTATTGAGGTCGGCGACTTTCACTCACAAACCAGCAAGTGGCGTCCCCTAGGGGATTCGAACCCCTGTTACCGCCGTGAAAGGGCGGTGTCCTGGGCCTCTAGACGAAGGGGACACGAAAATTGCTTATCACGCG‑TTGCGTGatat                                                                                                                                                     >  1:422518/1‑149 (MQ=32)  gagaTATTGAGGTCGGCGACTTTCACTCACAAACCAGCAAGTGGCGTCCCCTAGGGGATTCGAACCCCTGTTACCGCCGTGAAAGGGCGGTGTCCTGGGCCTCTAGACGAAGGGGACACGAAAATTGCTTATCACGCG‑TTGCGTGatat                                                                                                                                                     >  1:666771/1‑149 (MQ=32)  gagaTATTGAGGTCGGCGACTTTCACTCACAAACCAGCAAGTGGCGTCCCCTAGGGGATTCGAACCCCTGTTACCGCCGTGAAAGGGCGGTGTCCTGGGCCTCTAGACGAAGGGGACACGAAAATTGCTTATCACGCG‑TTGCGGGatat                                                                                                                                                     >  2:640268/1‑149 (MQ=17)   agaTATTGAGGGCGGCGACTTTCGGTCCTGAACCTGCGAGTGGCGTCCGCTAGGGGGTTCGATCCACTGTTACCGCCGTGAAAGGGCGGTCTCCTGGGCCTCGAGACGAAGGGGACACGAAAATTGCATATCACGCA‑ACGCGTGataag                                                                                                                                                    <  2:83101/149‑3 (MQ=11)    gaTATTGAGGTCGGCGACTTTCACCCACAAACCAGCAAGTGGCGTCCCCTAGGGGATTCGAACCCCTGTTACCGCCGTGAAAGGGCGGTGTCCTGGGCCTCTAGACGAAGGGGACACGAAAATTGCTTATCACGCG‑TTGCGTGATAttt                                                                                                                                                   >  2:366514/1‑149 (MQ=17)      tatTGAGGTCGGCGACTTTCACGCACAAACCAGCAAGTGGCGTCCCCTAGGGGATTCGAACCCCTGTTACCGCCGTGAAAGGGCGGTGTCCTGGGCCTCTAGACGAAGGGGACACGAAAATTGCTTATCACGCG‑TTGCGTGATATTTTc                                                                                                                                                 <  1:319197/149‑1 (MQ=17)       atTGAGGTCGGCGACTTTCACTCACAAACCAGCAAGTGGCGTCCCCTAGGGGATTCGAACCCCTGTTACCGCCGTGAAAGGGCGGTGTCCTGGGCCTCTAGACGAAGGGGACACGAAAATTGCTTATCACGCG‑TTGCGTGATATTTTCg                                                                                                                                                >  2:468382/1‑149 (MQ=32)        ttaagcacGGCGACTCTCGCGCACAAACCCGCAAGGGGCGGCCCCTAGGGGATTCGCACCCCTGTTCCCGCCGTGAAAGCGCCGTGTCCTGGGCCTCTAGACGCTGGGGTCACGAAAATTGCTTATCCTGCG‑TTGCGTGATATTTTCgt                                                                                                                                               <  2:625166/142‑1 (MQ=11)              tctgCGACTTTCACTCACAAACCAGCAACGGGGGTCCGCGAGGGGATTAGAACCCCTGTTACCGCGGTGAAAGGGCGGTGTCCTGGGCCTCTAGACGAAGGTGACACGAAAATTGCTTATCACGCG‑TTGCGTGATATTTTCGTGTAggg                                                                                                                                         <  2:191071/146‑1 (MQ=1)                                                                                       ggCGGTGTCCTGGGCCTCTAGACGAAGGGGACACGAAAATTGCTTATCACGCG‑TTGCGTGATATTTTCGTGTAGGGTGAGCTTTCATTAATAGAAAGCGAACGGCCTTATTCTCTTCAGCCTCACTCCCAACGCGTAAACGCCTTGCTa                                                                >  1:616025/1‑149 (MQ=11)                                                                                           gtgtCCTGGGCCTCTAGACGAAGGGGACACGAAAATTGCTTATCACGCG‑TTGCGTGATATTTTCGTGTAGGGTGAGCTTTCATTAATAGAAAGCGAACGGCCTTAGTCTCTTCAGCCTCACTCCCAACGCGTAAACGCCTTGCTATTCa                                                            <  1:640268/149‑1 (MQ=17)                                                                                              tCCTGGGCCTCTAGACGAAGGGGACACGAAAATTGCTTATCACGCG‑TTGCGTGATATTTTCGTGTAGGGTGAGCTTTCATTAATAGAAAGCGAACGGCCTTATTCTCTTCAGCCTCACTCCCAACGCGTAAACGCCTTGCTATTCACtt                                                         <  2:565027/149‑1 (MQ=32)                                                                                                cTGGGCCTCTAGACGAAGGGGACACGAAAATTGCTTATCACGCG‑TTGCGTGATATTTTCGTGTAGGGTGAACTTTCATTAATAAAAAACGAACGGCATTATTCTCTTCAGCCTCACTCCCAACGCGTAAACGACTTGCTATTCACtttt                                                       >  1:56915/1‑149 (MQ=11)                                                                                                 tGGGCCTCTAGACGAAGGGGACACGAAAATTGCTTATCACGCG‑TTGCGTGATATTTTCGTGTAGGGTGAGCTTTCATTAATAGAAAGCGAACGGCCTTATTCTCTTCAGCCTCACTCCCAACGCGTAAACGCCTTGCTATTCACTTTtc                                                      <  1:468382/149‑1 (MQ=34)                                                                                                     ccTCTAGACGAAGGGGACACGAAAATTGCTTATCACGCG‑TTGCGTGATATTTTCGTGTAGGGTGAGCTTTCATTAATAGAAAGCGAACGGCCTTATTCTCTTCAGCCTCACTCCCAACGCGTAAACGCCTTGCTATTCACTTTtcatca                                                  <  1:436123/149‑1 (MQ=38)                                                                                                      ctctAGACGAAGGGGACACGAAAATTGCTTATCACGCG‑TTGCGTGATATTTTCGTGTAGGGTGAGCTTTCATTAATAGAAAGCGAACGGCCTTATTCTCTTCAGCCTCACTCCCAACGCGTAAACGCCTTGCTATTCACTTTTCATCAg                                                 <  2:530700/149‑1 (MQ=38)                                                                                                       tctAGACGAAGGGGACACGAAAATTGCTTATCACGCG‑TTGCGTGATATTTTCGTGTAGGGTGAGCTTTCATTAATAGAAAGCGAACGGCCTTATTCTCTTCAGCCTCACTCCCAACGCGTAAACGCCTTGCTATTCACTTTTCATCAGa                                                >  2:153301/1‑149 (MQ=38)                                                                                                           gACGGAGGGGACACGAAAATTGCGTATCACGCG‑TTGCGGGATATTCTCGTGTAGGGTGAGCTGTCATTAATAGAAAGCTAACGGCCTCATTTTCTTCATCCTCAGTCCCACCGCGTAAAGGCCTTGCTATTCACTTTTCATCAGACAAt                                            <  1:652731/149‑1 (MQ=255)                                                                                                               caGGGGAGACAAAGACTGCGTATCACGCGCTTGC‑TGATATTTTCGTGTAGGGTGAGCTTTCATGAATAGGAAGCGATCGCCCTTATCCTCTTCAGCCTCACTCGCAACGCGTAAACGCCTTGCTATTCACTTTTCATCAGACAATCtgt                                        <  1:333103/148‑1 (MQ=11)                                                                                                                aGGGGAGACGAGAATTGCTTATCACTCG‑TTGCGTGATATTTTCGTGTAGGTTGAGCTTTGATTAATAGAAAGCTAACGGCCTTATTCTCTTCAGCCTCATTCCCAACGCGTAAACGCCTTGCGATTCACTTTTCATCAGACAATCtgtg                                       <  1:125101/149‑1 (MQ=11)                                                                                                                 tgggATACGAAAATTGCTTATCACGCG‑TTGCGTGATATTTTCGTGTAGGGTGAGCTTTCATTAATAGAAAGCGAACGGCCTTATTCTCTTCAGCCTCACTCCCAACGCGTAAACGCCTTGCTATTCACTTTTCATCAGACAATCtgtgt                                      <  2:737780/148‑1 (MQ=37)                                                                                                                       acGAAAATTGCTTATCACGCG‑TTGCGTGATATTTTCGTGTAGGGTGAGCTTTCATTAATAGAAAGCGAACGGCCTTATTCTCTTCAGCCTCACTCCCAACGCGTAAACGCCTTGCTATTCACTTTTCATCAGACAATCTGTGTGAGctc                                >  2:242556/1‑147 (MQ=39)                                                                                                                       acGAAAATTGCTTATCACGCG‑TTGCGTGATATTTTCGTGTAGGGTGAGCTTTCATTAATAAAAAGCGCACGGCATTATTCTCTTCAGCCTCACTCCCAACGCGTAAACGCCTTGCTATTCCCTTTTAATCAGACAATTTTTGGGAGCAc                                >  2:179285/1‑149 (MQ=11)                                                                                                                       acGAAAAGTGCTTATCACGCG‑TTGCGTGATATTTTCGTGTAGGGTGAGCTTTCATTAATAGAAAGCGAACGGCCTTATTCTCTTCAGCCTCACTCCCAACGCGTAAACGCCTTGCTATTCACTTTTCATCAGACAATCTGTGTGAGCAc                                <  1:292334/149‑1 (MQ=255)                                                                                                                        cGAAAATTGCTTATCACGCG‑TTGCGTGATATTTTCGTGTAGGGTGAGCTTTCATTAATAGAAAGCGAACGGCCTTATTCTCTTCAGCCTCACTCCCAACGCGTAAACGCCTTGCTATTCACTTTTCATCAGACAATCTGTGTGAGCACt                               >  2:200671/1‑149 (MQ=38)                                                                                                                         gAAAATTGCTTATCACGCG‑TTGCGTGATATTTTCGTGTAGGGTGAGCTTTCATTAATAGAAAGCGAACGGCCTTATTCTCTTCAGCCTCACTCCCAACGCGTAAACGCCTTGCTATTCACTTTTCATCAGACAATCTGTGTGAGCACTg                              >  2:572124/1‑149 (MQ=38)                                                                                                                         gAAAATTGCTTATCACGCG‑TTGCGTGATATTTTCGTGTAGGGTGAGCTTTCATTAATAGAAAGCGAAAGGCCTTATTCTCTTCAGCCTCACTCCCAACGCGTAAACGCCTTGATATTCACTTTTCATCAGACAATCTGTGTCAGAACTg                              >  2:233942/1‑149 (MQ=255)                                                                                                                          aaaaTTGCTTATCACGCG‑TTGCGTGATATTTTCGTGTAGGGTGAGCTTTCATTAATAGAAAGCGAACGGCCTTATTCTCTTCAGCCTCACTCCCAACGCGTAAACGCCTTGCTATTCACTTTTCATCAGACAATCTGTGTGAGCACTGc                             >  1:701412/1‑149 (MQ=40)                                                                                                                           aaaTTGCTTATCACGCG‑TTGCGTGATATTTTCGTGTAGGGTGAGCTTTCATTAATAGAAAGCGAACGGCCTTATTCTCTTCAGCCTCACTCCCAACGCGTAAACGCCTTGCTATTCACTTTTCATCAGACAATCTGTGTGAGCACTGCa                            <  1:153301/149‑1 (MQ=255)                                                                                                                                 tttATCAGGCA‑TTGTGTGATACTGTCGTGTCGGTGGATCTTTTATGAATAGAAAGCGAACGGCCTTATTTTCTTCAGCTTTACTCTTAACGCGTAAACGCCTCTCTATACACTTTGTACCAGTCAATCTCTGTGCGCCCTGCAACGTAc                      <  2:64956/148‑1 (MQ=2)                                                                                                                                 cTTATCACGCG‑TTGCGTGATATTTTCGTGTAGGGTGAGCTTTCATTAATAGAAAGCGAACGGCCTTATTCTCTTCAGCCTCACTCCCAACGCGTAAACGCCTTGCTATTCACCTTTCATCAGACAATCTGTGTGAGCACTGCAAAGTAc                      <  1:79888/149‑1 (MQ=255)                                                                                                                                  ttATCACGCG‑TTGCGTGATATTTTCGTGTAGGGTGAGCTTTCATTAATAGAAAGCGAACGGCCTTATTCTCTTCAGCCTCACTCCCAACGCGTAAACGCCTTGCTATTCACTTTTCATCAGACAATCTGTGTGAGCACTGCAAAGTACg                     >  1:44197/1‑149 (MQ=40)                                                                                                                                     tCACGCG‑TTGCGTGATATTTTCGTGTAGGGTGAGCTTTCATTAATAGAAAGCGAACGGCCTTATTCTCTTCAGCCTCACTCCCAACGCGTAAACGCCTTGCTATTCACTTTTCATCAGACAATCTGTGTGAGCACTGCAAAGTACGctt                  >  2:410355/1‑149 (MQ=40)                                                                                                                                           g‑TTGCGTGATATTTTCGTGTAGGGTGAGCTTTCATTAATAGAAAGCGAACGGCCTTATTCTCTTCAGCCTCACTCCCAACGCGTAAACGCCTTGCTATTCACTTTTCATCAGACAATCTGTGTGAGCACTGCAAAGTACGCTTCTTtaa            <  2:702731/149‑1 (MQ=255)                                                                                                                                           g‑TTGCGTGATATTTTCGTGTAGGGTGAGCTTTCATTAATAGAAAGCGAACGGCCTTATTCTCTTCAGCCTCACTCCCAACGCGTAAACGCCTTGCTATTCACTTTTCATCAGACAATCTGTGTGAGCACTGCAAAGTACGCTTCTTtaa            <  1:200671/149‑1 (MQ=255)                                                                                                                                             ttGCGTGATATTTTCGTGTAGGGTGAGCTTTCATTAATAGAAAGCGAACGGCCTTATTCTCTTCAGCCTCACTCCCAACGCGTAAACGCCTTGCTATTCACTTTTCATCAGACAATCTGTGTGAGCACTGCAAAGTACGCTTCTTtaag           <  2:427245/149‑1 (MQ=255)                                                                                                                                                     tatTTTCGTGTAGGGTGAGCTTTCATTAATAGAAAGCGAACGGCCTTATTCTCTTCAGCCTCACTCCCAACGCGTAAACGCCTTGCTATTCACTTTTCATCAGACAATCTGTGTGAGCACTGCAAAGTACGCTTCTTTAAGGTAaggag   <  1:622783/149‑1 (MQ=255)                                                                                                                                                      atTTTCGTGTAGGGTGAGCTTTCATTAATAGAAAGCGAACGGCCTTATTCTCTTCAGCCTCACTCCCAACGCGTAAACGCCTTGCTATTCACTTTTCATCAGACAATCTGTGTGAGCACCGCAAAGTACGCTTCTTTAAGGTAaggagg  <  1:473574/149‑1 (MQ=255)                                                                                                                                                      |                                                                                                                                                     TGAGATATTGAGGTCGGCGACTTTCACTCACAAACCAGCAAGTGGCGTCCCCTAGGGGATTCGAACCCCTGTTACCGCCGTGAAAGGGCGGTGTCCTGGGCCTCTAGACGAAGGGGACACGAAAATTGCTTATCACGCG‑TTGCGTGATATTTTCGTGTAGGGTGAGCTTTCATTAATAGAAAGCGAACGGCCTTATTCTCTTCAGCCTCACTCCCAACGCGTAAACGCCTTGCTATTCACTTTTCATCAGACAATCTGTGTGAGCACTGCAAAGTACGCTTCTTTAAGGTAAGGAGG  >  NC_000913_3_bme_pgi/2729327‑2729623

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.