New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913_3_pae_pgi	= 4192177	6 (0.370)	4 (0.260)	4/242	1.6	35.1%	coding (459/771 nt)	thiG	thiamine biosynthesis ThiGH complex subunit
	?	NC_000913_3_pae_pgi	= 4192203	9 (0.580)					coding (433/771 nt)	thiG	thiamine biosynthesis ThiGH complex subunit

AGCGCCTGCGCGGCATGGCTGGGAACGCCGATGCCAGCATCGACAACCACCGGCACTGTGGCCTGCTGGATAATAATCTCCAGCATGGCGCGGGTTTCCAGTCCCTGATTCGAGCCAATCGGCGCGCCGAGCGGCATCACCGCTGCACAGCCGACTTCTTCCAGACGTTTACACAATACCGGATCGGCCCCGCAGTAAGGCAGCACGACAAATCCCTGTTGTACCAGCGTTTCGGCGGCTTTCAGGGTTTCGATGGGATCGGGCAACAGCCA  >  NC_000913_3_pae_pgi/4192055‑4192326                                                                                                                           |                                                                                                                                                      aGCGCCTGCGCGGCATGGCTGGGAACGCCGATGCCAGCATCGACAACCACCGGCACTGTGGCCTGCTGGATAATAATCTCCAGCATGGCGCGGGTTTCCAGTCCCTGATTCGAGCCAATCGGCGCGCCGAGCGGCATCACCGCTGCACAg                                                                                                                            <  1:83849‑M2/150‑2 (MQ=255)                           gCCGATGCCAGCATCGACAACCACCGGCACTGTGGCCTGCTGGATAATAATCTCCAGCATGGCGCGGGTTTCCAGTCCCTGATTCGAGCCAATCGGCGCGCCGAGCGGCATCACCGCTGCACAGCCGACTTCTTCCAGACGTTTACACaa                                                                                                  >  2:6009/1‑150 (MQ=255)                                 gCCAGCATCGACAACCACCGGCACTGTGGCCTGCTGGATAATAATCTCCAGCATGGCGCGGGTTTCCAGTCCCTGATTCGAGCCAATCGGCtgtgcagcggtgatgccgctcggcgcgccgattggctc                                                                                                                 <  1:169889‑M1/129‑39 (MQ=255)                                 gCCAGCATCGACAACCACCGGCACTGTGGCCTGCTGGATAATAATCTCCAGCATGGCGCGGGTTTCCAGTCCCTGATTCGAGCCAATCGGCtgtgcagcggtgatgccgctcggcgcgccgattggctc                                                                                                                 >  2:169889‑M1/1‑91 (MQ=255)                                          gACAACCACCGGCACTGTGGCCTGCTGGATAATAATCTCCAGCATGGCGCGGGTTTCCAGTCCCTGATTCGAGCCAATCGGCtgtgcagcggtgatgccgctcggcgcgccg                                                                                                                         <  3:60315‑M1/112‑31 (MQ=255)                                          gACAACCACCGGCACTGTGGCCTGCTGGATAATAATCTCCAGCATGGCGCGGGTTTCCAGTCCCTGATTCGAGCCAATCGGCtgtgcagcggtgatgccgctcggcgcgccg                                                                                                                         >  4:60315‑M1/1‑82 (MQ=255)                                                         tgtgGCCTGCTGGATAATAATCTCCAGCATGGCGCGGGTTTCCAGTCCCTGATTCGAGCCAATCGGCGCGCCGAGCGGCATCACCGCTGCACAGCCGACTTCTTCCAGACGTTTACACAATACCGGATCGGCCCCGCAGTAAGGCAGCAc                                                                    >  2:147316/1‑150 (MQ=255)                                                                     gATAATAAACGCCAGCATGGCGCGCGTTTCCAGACTCTGATTCGAGCCAAACGGCGCGCCGAGCGGCATAACCGCTGCCCAGCCGACCTCATCCAGACGATTCCCCAATACCGGAACGACCCCACAGCAAGACAGCACGACAAATccccg                                                        >  2:138399/1‑148 (MQ=255)                                                                                     aTGGCGCGGGTTTCCAGTCCCTGATTCGAGCCAATCGGCGCGCCGAGCGGCATCACCGCTGCACAGCCGACTTCTTCCAGACGTTTACACAATACCGGATCGGCCCCGCAGTAAGGCAGCACGACAAATCCCTGTTGTACCAGCGTTTcg                                        >  4:108002/1‑150 (MQ=255)                                                                                                         cTGATTCGACCCAATCGGCGCGCCGAGCGGCATCACCGCTGCACAGCCGACTTCTTCCAGACGTTTACACAATACCGGATCGGCCCCGCAGTAAGGCAGCACGACAAATCCCTGTTGTACCAGCGTTTCGGCGGCTTTCAGGGTTTCGAt                    >  1:137038/1‑150 (MQ=255)                                                                                                                gAGCCAATCGGCGCGCCGAGCGGCATCACCGCTGCACAgccgattggctcgaatcagggactggaaacccgcgccatgctggagattattatccagcaggccacagtgccggtggttgtcgatgctggc                                  <  1:169889‑M2/129‑92 (MQ=255)                                                                                                                gAGCCAATCGGCGCGCCGAGCGGCATCACCGCTGCACAgccgattggctcgaatcagggactggaaacccgcgccatgctggagattattatccagcaggccacagtgccggtggttgtcgatgctggc                                  >  2:169889‑M2/1‑38 (MQ=255)                                                                                                                        cGGCGCGCCGAGCGGCATCACCGCTGCACAgccgattggctcgaatcagggactggaaacccgcgccatgctggagattattatccagcaggccacagtgccggtggttgtc                                           <  3:60315‑M2/112‑83 (MQ=255)                                                                                                                        cGGCGCGCCGAGCGGCATCACCGCTGCACAgccgattggctcgaatcagggactggaaacccgcgccatgctggagattattatccagcaggccacagtgccggtggttgtc                                           >  4:60315‑M2/1‑30 (MQ=255)                                                                                                                          gcgcgcCGAGCGGCATCACCGCTGCACAGCCGACTTCTTCCAGACGTTTACACAATACCGGATCGGCCCCGCAGTAAGGCAGCACGACAAATCCCTGTTGTACCAGCGTTTCGGCGGCTTTCAGGGTTTCGATGGGATCGGGCAAcagcc   <  1:147316/150‑1 (MQ=255)                                                                                                                           ctgtgcagcggtgatgccgctcggcgcgccgattggctcgaatcagggactggaaacccgcgccatgctggagattattatccagcaggccacagtgccggtggttgtcgatgctggcatcggcgttcccagccatgccgcgcaggcgct  <  1:83849‑M1/150‑150 (MQ=255)                                                                                                                           |                                                                                                                                                      AGCGCCTGCGCGGCATGGCTGGGAACGCCGATGCCAGCATCGACAACCACCGGCACTGTGGCCTGCTGGATAATAATCTCCAGCATGGCGCGGGTTTCCAGTCCCTGATTCGAGCCAATCGGCGCGCCGAGCGGCATCACCGCTGCACAGCCGACTTCTTCCAGACGTTTACACAATACCGGATCGGCCCCGCAGTAAGGCAGCACGACAAATCCCTGTTGTACCAGCGTTTCGGCGGCTTTCAGGGTTTCGATGGGATCGGGCAACAGCCA  >  NC_000913_3_pae_pgi/4192055‑4192326

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