New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913	= 4481898	15 (0.570)	3 (0.120)	3/282	2.9	17.5%	coding (1940/2856 nt)	valS	valyl‑tRNA synthetase
	?	NC_000913	= 4481970	14 (0.560)					coding (1868/2856 nt)	valS	valyl‑tRNA synthetase
Rejected: Frequency below cutoff threshold.

AGCGGAAGCTGTCCAGCGCTTCGCGGTACGCTTTGATGGTCTGGTTGAACTCCGCCAGAATCCAGCGGTCCGCCAGCGACAGCGTCATTTCGCCGCCGTTGAAGCCGCAATCCTGACCTTCTGTGTTCATCAGCACAAAGCGGCTGGCGTTCCACAGCTTGTTACAGAAGTTACGGTAACCTTCCAGACGCTTCATATCCCAGTTGATGTCACGACCGGTAGACGCCAGCGCCGCCAGGGTGAAGCGCAGCGCGTCAGTACCGTGCGGCTCAATACCGTTCGGGAACTGCTTCTCGG  >  NC_000913/4481749‑4482045                                                                                                                                                      |                                                                                                                                                    aGCGGAAGCTGTCCAGCGCTTCGCGGTACGCTTTGATGGTCTGGTTGAACTCCGCCAGAATCCAGCGGTCCGCCAGCGACAGCGTCATTTCGCCGCCGTTGAAGCCGCAATCCTGACCTTCTGTGTTCATCAGCACAAAGCGGCTGGCGtt                                                                                                                                                    >  2:156141/1‑151 (MQ=255)      aaGCTGTCCAGCGCTCCGCGGTCCGCTTTGATGGTCTGGTTGAACTCCGCCAGAATCCCGCGGTCCGCCAGCGACAGCGTCATTTCGCCGCCGTTGAAGCCGCAATCCTGACCTTCTGTGTTCATCAGCACAAAGCGGCTGGCGTTCcac                                                                                                                                                <  2:366133/150‑1 (MQ=255)                                           ggTTGAACTCCGCCAGAATCCAGCGGTCCGCCAGCGACAGCGTCATTTCGCCGCCGTTGAAGCCGCAATCCTGACCTTCTGTGTTCATCAGCACAAAGCGGCTGGCGTTCCACAGCTTGTTACAGAAGTTACGGTAACCTTCCAGACGCtt                                                                                                          >  2:38105/1‑151 (MQ=255)                                                 aCTCCGCCAGAATCCAGCGTCCCGCCAGCGACAGCGTCATTTCGCCGCCGTTGAAGCCGCAATCCTGGCCTTCTGTGTTCATCAGCACAAAGCGGCTGGCGTTCCACAGCTTGTTACAGAAGGTACGGTAACCTTCCAGACGCTTCATATc                                                                                                    <  2:296166/151‑1 (MQ=255)                                                          gAATCCAGCGGTCCGCCAGCGACAGCGTCATTTCGCCGCCGTTGAAGCCGCAATCCTGACCTTCTGTGTTCATCAGCACAAAGCGGCTGGCGTTCCACAGCTTGTTACAGAAGTTACGGTAACCTTCCAGACGCTTCATATCCCAGTTGAt                                                                                           >  1:51961/1‑151 (MQ=255)                                                          gAATCCAGCGGTCCGCCAGCGACAGCGTCATTTCGCCGCCGTTGAAGCCGCAATCCTGACCTTCTGTGTTCATCAGCACAAAGCGGCTGGCGGctaccggtcctgacatcaactgggatatgaagcgtctggaaggttaccgtaactt                                                                                              >  2:232655‑M1/1‑93 (MQ=255)                                                            aTCCAGCGGTCCGCCAGCGACAGCGTCATTTCGCCGCCGTTGAAGCCGCAATCCTGACCTTCTGTGTTCATCAGCACAAAGCGGCTGGCGTTCCACAGCTTGTTACAGAAGTTACGGTAACCTTCCAGACGCTTCATATCCCAGTTGATGt                                                                                         >  2:244967/1‑151 (MQ=255)                                                                  cGGTCCGCCAGCGACAGCGTCATTCCGCCGCCGTTGAAGCCGCAATCCTGACCTTCTGTGTTCATCAGCACAAAGCGGCTGGCGTTCCACAGCTTGTTACAGAAGTTACGGTAACCTTCCAGACGCTTCATATCCCAGTTGATGTCACGAc                                                                                   <  1:123463/151‑1 (MQ=255)                                                                             cGACAGCGTCATTTCGCCGCCGTTGAAGCCGCAATCCTCACCTTCTGTGTTCATCAGCACAAAGCGGCTGGCGTctaccggtcgtgacctcaactgcgatatgatgcgtcaggaaggttaccgtaactttcttaacaacctgttgcacgc                                                                         >  1:220242‑M1/1‑74 (MQ=255)                                                                                       aTTTCGCCGCCGTTGAAGCCGCAATCCTGACCTTCTGTGTTCATCAGCACAAAGCGGCGGGCGTTCCACAGCTTGTTACAGAAGTTACGGGAACCTTCCAGACGCTTCATATCCCAGTTGATGTCACGACCGGTAGacgccagccgctttg                                                              >  2:43783‑M2/1‑136 (MQ=255)                                                                                                              aTCCTGACCTTCTGTGTTCATCAGCACAAAGCGGCTGGCGTTCCACAGCTTGTTACAGAAGTTACGGTAACCTTCCAGACGCTTCATATCCCAGTTGATGTCACGACCGGTAGACGCCAGCGCCGCCAGGGTGAAGCGCAGCGCGTCAGTa                                       >  1:327816/1‑151 (MQ=255)                                                                                                              aTCCTGACCTTCTGTGTTCATCAGCACAAAGCGGCTGGCGTTCCACAGCTTGTTACAGAAGTTACGGTAACCTTCCAGACGCTTCATATCCCAGTTGATGTCACGACCGGTAGACGCCAGCGCCGCCAGGGTGAAGCGCAGCGCGTCAGTa                                       >  1:288605/1‑151 (MQ=255)                                                                                                                ccTGACCTTCTGTGTTCATCAGCACAAAGCGGCTGGCGTTCCACAGCTTGTTACAGAAGTTACGGTAACCTTCCAGACGCTTCATATCCCAGTTGATGTCACGACCGGTAGACTCCAGCGCCGCCAGGGTGAAGCGCAGCGCGTCAGTAc                                      >  2:88810/1‑150 (MQ=255)                                                                                                                        tCTGTGTTCATCAGCACAAAGCGGCTGGCGTTCCACAGCTTGTTACAGAAGTTACGGTAAACTTCCAGACCCTTCATATCCCAGTTGCTGTCACGACCGGTAGACGCCAGCGCCGCCAGGGTGAAGCGCAGCGCGTCAGTACCGTGCgg                               >  2:271072/1‑149 (MQ=255)                                                                                                                         cTGTGTTCATCAGCACAAAGCGGCTGGCGTTCCACAGCTTGTTACAGAAGTTACGGTAACCTTCCAGACGCTTCATATCCCAGTTGATGTCACGACCGGTAGACGCCAGCGCCGCCAGGGTGAAGCGCAGCGCGTCAGTACCGTGCGGCt                             >  2:18682/1‑150 (MQ=255)                                                                                                                                     gaacaAAGCGGCTGGCGTTCCACAGCTTGTTACAGAAGTTACGGTAACCTTCCAGACGCTTCATATCCCAGTTGATGTCACGACCGGTAGACGCCAGCGCCGCCAGGGTGAAGCGCAGCGCGTCAGTACCGTGCGGCTCAATACCGTTCgg                >  2:55614/3‑151 (MQ=255)                                                                                                                                        caAAGCGGCTGGCGTctaccggtcgtgacatcaactgggatatgaagcgtctggaaggttcccgtaacttctgtaacaagctgtggaacgcccgccgctttgtgctgatgaacacagaaggtcaggattgcggcttcaacggcggcgaaat             >  2:43783‑M1/1‑15 (MQ=255)                                                                                                                                                   gCGTGCAACAGGTTGTTAAGAAAGTTACGGTAACCTTCCTGACGCATCATATCGCAGTTGAGGTCACGACCGGTAGacgccagccgctttgtgctgatgaacacagaaggtgaggattgcggcttcaacggcggcgaaatgacgctgtcg   >  1:220242‑M2/1‑76 (MQ=255)                                                                                                                                                    cGTTCCACAGCTTGTTACAGAAGTTACGGTAACCTTCCAGACGCTTCATATCCCAGTTGATGTCACGACCGGTAGACGCCAGCGCCGCCAGGGTGAAGCGCAGCGCGTCAGTACCGTGCGGCTCAATACCGTTCGGGAACTGCTTCTCgg  <  1:55608/150‑1 (MQ=255)                                                                                                                                                    cGTTCCAAAGCTTGTTACAGAAGTTACGGTAACCTTCCAGACGCTTCATATCCCAGTTGATGTCACGACCGGTAGACGCCAGCGCCGCCAGGGTGAAGCGCAGCGCGTCAGTACCGTGCGGCTCAATACCGTTCGGGAACTGCTTCTCgg  <  2:10114/150‑1 (MQ=255)                                                                                                                                                      |                                                                                                                                                    AGCGGAAGCTGTCCAGCGCTTCGCGGTACGCTTTGATGGTCTGGTTGAACTCCGCCAGAATCCAGCGGTCCGCCAGCGACAGCGTCATTTCGCCGCCGTTGAAGCCGCAATCCTGACCTTCTGTGTTCATCAGCACAAAGCGGCTGGCGTTCCACAGCTTGTTACAGAAGTTACGGTAACCTTCCAGACGCTTCATATCCCAGTTGATGTCACGACCGGTAGACGCCAGCGCCGCCAGGGTGAAGCGCAGCGCGTCAGTACCGTGCGGCTCAATACCGTTCGGGAACTGCTTCTCGG  >  NC_000913/4481749‑4482045

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