New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913_3_bme_pgi	= 1529619	15 (1.040)	4 (0.290)	4/250	1.4	21.4%	pseudogene (1718/2037 nt)	rhsE	pseudogene, Rhs family
	?	NC_000913_3_bme_pgi	= 1529633	15 (1.080)					pseudogene (1732/2037 nt)	rhsE	pseudogene, Rhs family

AGGTGGACTAATTTATGCCGTATCTGGTGTACCTGGATTGATTGCTGCAAACAGCATTACTAACAGTGCTTACCAGTTCGGTTATGATATGGATGCTATTGTTGGCGGAGCTCATAATGGGGCCGCCGATGCAATGAGACATTGTTACTTGATGTGTCGAATGACTAAGACATTTGGATCAACAATAGCTGACGTGATAGGTAAAAATCATGAGGCGGCAGGGGATAGACAAGGTCAGCCAGCTAAAGAAAGAATCA  >  NC_000913_3_bme_pgi/1529512‑1529768                                                                                                                          |                                                                                                                                        aGGTGGACTAATTTATGCCGTATCTGGTGTACCTGGATTGATTGCTGCAAACAGCATTACTAACAGTGCTTACCAGTTCGGTTATGATATGGATGCTATTGTTGGCGGAGCTCATAATGGGGCCGCCGATGCAATg                                                                                                                           <  1:162102/136‑1 (MQ=255) aGGTGGACTAATTTATGCCGTATCTGGTGTACCTGGATTGATTGCTGCAAACAGCATTACTAACAGTGCTTACCAGTTCGGTTATGATATGGATGCTATTGTTGGCGGAGCTCATAATGGGGCCGCCGATGCAATg                                                                                                                           >  2:162102/1‑136 (MQ=255)                                  tGGATTGATTGCTGCAAACAGCATTACTAACAGTGCTTACCAGTTCGGTTATGATATGGATGCTATTGTTGGCGGAGCTCATAATGGGGCCGCCGATGCAATGAGACATTGTTACTTGATGTGTCGAATGACTAAGACATTTGGATcaac                                                                            >  3:59645/1‑150 (MQ=255)                                      ttgattgCTGCAAACAGCATTACTAACAGTGCTTACCAGTTCGGTTATGATATGGATGCTATTGTTGGCGGAGCTCATAATGGGGCCGCCGATGCAATGAGACATTGTTACTTGATGTGTCGAATGACTAAGACATTTGGATCAACAATa                                                                        >  1:133825/1‑150 (MQ=255)                                           tgctgcAAACAGCATTACTAACAGTGCTTACCAGTTCGGTTATGATATGGATGCTATTGTTGGCGGAGCTCATAATGGGGCCGCCGATGCAATGAGACATTGTTACTTGATGTGTCGAATGACTAAGACATTTGGATCAACAATAGCTGa                                                                   <  4:61308/150‑1 (MQ=255)                                              tgcAAACAGCATTACTAACAGTGCTTACCAGTTCGGTTATGATATGGATGCTATTGTTGGCGGAGCTCATAATGGGGCCGCCGATGCAATGAGACATTGTTACTTGATGTGTCGAATGACTAAGACATTTGGATCAACAATAGCTGACGt                                                                >  4:98724/1‑150 (MQ=255)                                                           aCTAACAGTGCTTACCAGTTCGGTTATGATATGGATGCTATTGTTGGCGGAGCTCATAATGGGGCCGCCGATGCAATGAGACATTGTTACTTGATGTGTCGAATGACTAAGACATTTGGATCAACAATAGCTGACGTGATAGGTAAAAAt                                                   >  4:100439/1‑150 (MQ=255)                                                             tAACAGTGCTTACCAGTTCGGTTATGATATGGATGCTATTGTTGGCGGAGCTCATAATGGGGccgccaacaatagcatccatatcataaccgaactggtaagcactgtta                                                                                         <  1:51906‑M2/110‑49 (MQ=255)                                                             tAACAGTGCTTACCAGTTCGGTTATGATATGGATGCTATTGTTGGCGGAGCTCATAATGGGGccgccaacaatagcatccatatcataaccgaactggtaagcactgtta                                                                                         >  2:51906‑M2/1‑62 (MQ=255)                                                              aaCAGTGCTTACCAGTTCGGTTATGATATGGATGCTATTGTTGGCGGAGCTCATAATGGGGCCGCCGATGCAATGAGACATTGTTACTTGATGTGTCGAATGACTAAGACATTTGGATCAACAATAGCTGACGTGATAGGTAAAAATCAt                                                >  2:47604/1‑150 (MQ=255)                                                                                 gTTATGATATGGATGCTATTGTTGGCGGAGCTCATAATGGGGccgccaacaatagcatccatatcataaccgaactggtaagcactgttagtaatgctgtttgcagcaatcaatccaggtacaccagatacggcataaattagtc                                  <  2:147849‑M2/145‑104 (MQ=255)                                                                                 gTTATGATATGGATGCTATTGTTGGCGGAGCTCATAATGGGGccgccaacaatagcatccatatcataaccgaactggtaagcactgttagtaatgctgtttgcagcaatcaatccaggtacaccagatacggcataaattagtc                                  >  1:147849‑M2/1‑42 (MQ=255)                                                                                       atatGGATGCTATTGTTGGCGGAGCTCATAATGGGGCCGCCGATGCAATGAGACATTGTTACTTGATGTGTCGAATGACTAAGACATTTGGATcaacaa                                                                          <  3:60025/99‑1 (MQ=255)                                                                                       atatGGATGCTATTGTTGGCGGAGCTCATAATGGGGCCGCCGATGCAATGAGACATTGTTACTTGATGTGTCGAATGACTAAGACATTTGGATcaacaa                                                                          >  4:60025/1‑99 (MQ=255)                                                                                            gATGCTATTGTTGGCGGAGCTCATAATGGGGCCGCCGATGCAATGAGACATTGTTACTTGATGTGTCGAATGa                                                                                               <  2:140354/73‑1 (MQ=255)                                                                                            gATGCTATTGTTGGCGGAGCTCATAATGGGGCCGCCGATGCAATGAGACATTGTTACTTGATGTGTCGAATGa                                                                                               >  1:140354/1‑73 (MQ=255)                                                                                            gATGCTATTGTTGGCGGAGCTCATAATGGGGCCGCCGATGCAATGAGACATTGTTACTTGATGTGTCGAATGACTAAGACATTTGGATCAACAATAGCTGACGTGATAGGTAAAAATCATGAGGCGGCAGGGGATAGACAAGGTcagcca                  <  1:47604/150‑1 (MQ=255)                                                                                                  aTTGTTGGCGGAGCTCATAATGGGGCCGCCGATGCAATGAGACATTGTTACTTGATGTGTCGAATGACTAAGACATTTGGATCAACAATAGCTGACGTGATAGGTAAAAATCATGAGGCGGCAGGGGATAGACAAGGTCAGCCAGCTaaa            <  4:59645/150‑1 (MQ=255)                                                                                                            gAGCTCATAATGGGGCCGCCGATGCAATGAGACATTGTTACTTGATGTGTCGAATGACTAAGACATTTGGATCAACAATAGCTGACGTGATAGGTAAAAATCATGAGGCGGCAGGGGATAGACAAGGTCAGCCAGCTAAAGAAAGAATCa  >  2:68562/1‑150 (MQ=255)                                                                                                                          |                                                                                                                                        AGGTGGACTAATTTATGCCGTATCTGGTGTACCTGGATTGATTGCTGCAAACAGCATTACTAACAGTGCTTACCAGTTCGGTTATGATATGGATGCTATTGTTGGCGGAGCTCATAATGGGGCCGCCGATGCAATGAGACATTGTTACTTGATGTGTCGAATGACTAAGACATTTGGATCAACAATAGCTGACGTGATAGGTAAAAATCATGAGGCGGCAGGGGATAGACAAGGTCAGCCAGCTAAAGAAAGAATCA  >  NC_000913_3_bme_pgi/1529512‑1529768

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