New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913	3783137 =	24 (0.800)	4 (0.140)	4/208	2.9	18.0%	coding (525/1020 nt)	gpsA	glycerol‑3‑phosphate dehydrogenase (NAD+)
	?	NC_000913	3783176 =	14 (0.500)					coding (486/1020 nt)	gpsA	glycerol‑3‑phosphate dehydrogenase (NAD+)
Rejected: Frequency below/above cutoff threshold.

AAACCGATACCGTCGGACATCCCCGCACCAATGGCAATAACGTTTTTCACCGCGCCGCCAAGCTGCACGCCAATGAAATCCGGATTGCTGTAAACGCGGAAACTTTTGCCGCAGTGCAGCAGCTGCTGGAGATCATCGGCAAAGGTCTGATCGGTCGAGGCCAGCGAAATAGCTGTCGGTAAACCTGCCGCCAGTTCTTTCGCAAACGTTGGGCCAGAGATAACCGCCAGCGGAATTTGATCGCCTAAG  >  NC_000913/3783033‑3783281                                                                                                                                                |                                                                                                          aaaCCGATACCGTCGGACATCCCCGCACCAATGGCAATAACGTTTTTCACCGCGCCGCCAAGCTGCACGCCAATGAAATCCGGATTGCTGTAAACGCGGAAACTTTTGCCGCAGTGCAGCAGCTGCTGGAGATCATCGGCAAAGGTCt                                                                                                       >  5:73943/1‑148 (MQ=255)                             cAATGGCAATAACGTTTTTCACCGCGCCGCCAAGCTGCACGCCAATGAAATCCGGATTGCTGTAAACGCGGAAACTTTTGCCGCAGTGCAGCAGCTGCTGGAGATCATCGGCAAAGGTCTGATCGGTCGAGGCCAGCGAAATAGCTGTc                                                                          >  4:103677/1‑149 (MQ=255)                                              ttCACCGCGCCGCCAAGCTGCACGCCAATGAAATCCGGATTGCTGTAAACGCGGAAACTTTTGCCGCAGTGCAGCAGCTGCTGGAGATCATCGGCAAAGGTCTGATCGGTCGAGGCCAGCGAAATAGCTGTCGGTAAACCTGCCGCCa                                                          <  3:103677/148‑1 (MQ=255)                                                                      cgacagctatttcgctggcctcgaccgatcagaccTTTGCCGCAGTGCAGCAGCTGCTGGAGATCATCGGCAAAGGTc                                                                                                        >  5:4689‑M1/36‑78 (MQ=255)                                                                      cgacagctatttcgctggcctcgaccgatcagaccTTTGCCGCAGTGCAGCAGCTGCTGGAGATCATCGGCAAAGGTc                                                                                                        <  6:4689‑M1/43‑1 (MQ=255)                                                                                        cTGTAAACGCGGAAACTTTTGCCGCAGTGCAGCAGCTGCTGGAGATCATCGGCAAAGGTCTGATCGGTCGAGGCCAGCGAAATAGCTGGCGGTAAACCTGCCGCCAGTTCTTTCGCAAACGTTGGGCCAGAGATAACCGCCAGCGGaa                <  6:73943/148‑1 (MQ=255)                                                                                         tGTAAACGCGGAAACTTTTGCCGCAGTGCAGCAGCTGCTGGAGATCATCGGCAAAGGTCTGATCGGTCGAGGCCAGCGAAATAGCTGTCGGTAAACCTGCCGCCAGTTCTTTCGCAAACGTTGGGCCAGAGATAACCGCCAGCGGAAtt              >  1:45562/1‑149 (MQ=255)                                                                                                 cgGAAACTTTTGCCGCAGTGCAGCAGCTGCTGGAGATCATCGGCAAAGGTCTGATCGGTCGAGGCCAGCGAAATAGCTGTCGGTAAACCTGCCGCCAGTTCTTTCGCAAACGTTGGGCCAGAGATAACCGCCAg                     >  2:135496/1‑134 (MQ=255)                                                                                                 cgGAAACTTTTGCCGCAGTGCAGCAGCTGCTGGAGATCATCGGCAAAGGTCTGATCGGTCGAGGCCAGCGAAATAGCTGTCGGTAAACCTGCCGCCAGTTCTTTCGCAAACGTTGGGCCAGAGATAACCGCCAg                     <  1:135496/134‑1 (MQ=255)                                                                                                  ggAAACTTATGCCGCAGTGCAGCAGCTGCTGGAGATCATCGGCAAAGGTCTGATCGGTCGAGGCCAGCGAAATAGCTGTCGGTAAACCTGCCGCCATTTCTTTCGCAAACGTTGGGCCAGAGATAACCGCCAGCGGAATTTGATCGcc      >  5:144030/1‑148 (MQ=255)                                                                                                    aaaCTTTTGCCGCAGTGCAGCAGCTGCTGGAGATCATCGGCAAAGGTCTGATCGGTCGAGGCCAGCGAAATAGCTGTCGGTAAACCTGCCGCCAGTTCTTTCGCAAACGTTGGGCCAGAGATAACCGCCAGCGGAATTTGATCGCCTaa   <  2:45562/149‑1 (MQ=255)                                                                                                     gacctttgccgatgatctccagcagctgctgcactgcggcaaaGGTCTGATCGGTCGAGGCCAGCGAAATAGCTGTCg                                                                         <  6:4689‑M2/35‑1 (MQ=255)                                                                                                     gacctttgccgatgatctccagcagctgctgcactgcggcaaaGGTCTGATCGGTCGAGGCCAGCGAAATAGCTGTCg                                                                         >  5:4689‑M2/44‑78 (MQ=255)                                                                                                     aaCGTTTGCCGCAGTGCAGCAGCTGCGGGAGATCATCGGCAAAGGTCTGATCGGGCGAGGCCAGCGAAATAGCTGTCGGTAAACCTGCCGCCAGTTCTTTCGCAAACGTTGGGCCAGAGATAACCGCCAGCGGAATTTGATCGCCTAAg  <  6:144030/149‑1 (MQ=255)                                                                                                        ttttGCCGCAGTGCAGCAGCTGCTGGAGATCATCGGCAAAGGTCTGATCGGTCGAg                                                                                            >  8:39127/1‑56 (MQ=255)                                                                                                        ttttGCCGCAGTGCAGCAGCTGCTGGAGATCATCGGCAAAGGTCTGATCGGTCGAg                                                                                            <  7:39127/56‑1 (MQ=255)                                                                                                        ttttGCCGCAGTGCAGCAGCTGCTGGAGATCATCGGCAAAGGTCTGATCGGTCGAGGCCAGCGAAATAGCTGTCGGTAAACCTGCCGCCAGTTCttt                                                   >  6:147714/1‑97 (MQ=255)                                                                                                        ttttGCCGCAGTGCAGCAGCTGCTGGAGATCATCGGCAAAGGTCTGATCGGTCGAGGCCAGCGAAATAGCTGTCGGTAAACCTGCCGCCAGTTCttt                                                   <  5:147714/97‑1 (MQ=255)                                                                                                         tttgccgatgatctccagcagctgctgcactgcggcaaaGGTCTGATCGGTCGAGGCCAGCGAAATAGCTGTCGGTAAACCTGCCGCCAGTTCTTTCGCAAACGTTGGGCCAGAGATAACCGCCAGCGGAATTTGa           >  4:81910‑M2/40‑136 (MQ=255)                                                                                                         tttgccgatgatctccagcagctgctgcactgcggcaaaGGTCTGATCGGTCGAGGCCAGCGAAATAGCTGTCGGTAAACCTGCCGCCAGTTCTTTCGCAAACGTTGGGCCAGAGATAACCGCCAGCGGAATTTGa           <  3:81910‑M2/97‑1 (MQ=255)                                                                                                                                            caaaGGTCTGATCGGTCGAGGCCAGCGAAATAGCTGTCGGTAAACCTGCCGCCAGTTCTTTCGCAAACGTTggg                                      >  4:64056‑M2/5‑74 (MQ=255)                                                                                                                                                |                                                                                                          AAACCGATACCGTCGGACATCCCCGCACCAATGGCAATAACGTTTTTCACCGCGCCGCCAAGCTGCACGCCAATGAAATCCGGATTGCTGTAAACGCGGAAACTTTTGCCGCAGTGCAGCAGCTGCTGGAGATCATCGGCAAAGGTCTGATCGGTCGAGGCCAGCGAAATAGCTGTCGGTAAACCTGCCGCCAGTTCTTTCGCAAACGTTGGGCCAGAGATAACCGCCAGCGGAATTTGATCGCCTAAG  >  NC_000913/3783033‑3783281

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