New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913 3955668 =35 (0.610)9 (0.170) 6/306 4.5 26.0% coding (338/1545 nt) ilvA l‑threonine dehydratase, biosynthetic; also known as threonine deaminase
?NC_000913 3955897 = 18 (0.330)coding (567/1545 nt) ilvA l‑threonine dehydratase, biosynthetic; also known as threonine deaminase
Rejected: Coverage evenness skew score above cutoff.

AGTCGACGCGGTGCGCGGCTTCGGCGGCGAAGTGCTGCTCCACGGCGCGAACTTTGATGAAGCGAAAGCCAAAGCGATCGAACTGTCACAGCAGCAGGGGTTCACCTGGGTGCCGCCGTTCGACCATCCGATGGTGATTGCCGGGCAAGGCACGCTGGCGCTGGAACTGCTCCAGCAGGACGCCCATCTCGACCGCGTATTTGTGCCAGTCGGCGGCGGCGGTCTGGCTGCTGGCGTGGCGGTGCTGATCAAACAACTGATGCCGCAAATCAAAGTGATCGCCGTAGAAGCGGAAGACTCCGCCTGCCTGAAAGCAGCGCTGGATGCGGGTCATCCGGTTGATCTGCCGCGCGTAGGGCTATTTGCTGAAGGCGTAGCGGTAAAACGCATCGGTGACGAAACCTTCCGTTTATGCCAGGAGTATCTCGACGACATCATCACCGTCGATA  >  NC_000913/3955681‑3956129
                                                                                                                                                                                                                        |                                                                                                                                                                                                                                        
aggcGACGCGGTGGGGGGCTTCGGCGGCGAAGTGCTGCTCCACGGCGCGAACTTTGATGAAGCGAAAGCCAAAGCGCTCGAACTGTCACAGCAGCAGGGGTTCACCTGGGTGCCGCCGTTCGACCATCCGATGGTGATTGCCGGGCAAGGCACGCTGGCGCTGGAACTGCTCCAGCAGGACGCCCATCTCGACCGCGTATTTGTGCCAGTCGGCGGCGGCGGTCTGGCTGCtggcgtgg                                                                                                                                                                                                                    <  2:402551/236‑1 (MQ=255)
                     cggcggCGAAGTGCTGCTCCACGGCGCGAACTTTGATGAAGCGAAAGCCAAAGCGATCGAACTGTCACAGCAGCAGGGGTTCACCTGGGTGCCGCCGTTCGACCATCCGATGGTGATTGCCGGGCAAGGCACGCTGGCGCTGGAACTGCTCCAGCAGGACGCCCATCTCGACCGCGTATTTGTGCCAGTCGGCGGCGGCGGTCTGGCTGCTGGCGTGGCGGTGCTGATCAAACAACTGa                                                                                                                                                                                               <  1:637201/239‑1 (MQ=255)
                                        cACGGCGCGAACTTTGATGAAGCGAAAGCCAAAGCGATCGAACTGTCACAGCAGCAGGGGTTCACCTGGGTGCCGCCGTTCGACCATCCGATGGTGATTGCCGGGCAAGGCACGCTGGCGCTGGAACTGCTCCAGCAGGACGCCCATCTCGACCGCGTATTTGTGCCAGTCGGCGGCGGCGGTCTGGCTGCTGGCGTGGCGGTGCTGATCAAACAACTGATGCCGCAAATCAAAGTGAt                                                                                                                                                                            <  2:640091/239‑1 (MQ=255)
                                          cGGCGCGAACTTTGATGAAGCGAAAGCCAAAGCGATCGAACTGTCACAGCAGCAGGGGTTCACCTGGGTGCCGCCGTTCGACCATCCGATGGTGATTGCCGGGCAAGGCACGCTGGCGCTGGAACTGCTCCAGCAGGACGCCCATCTCGACCGCGTATTTGTGCCAGTCGGCGGCGGCGGTCTGGCTGCTGGCGTGGCGGTGCTGATCAAACAACTGATGCCGCAAATCAAAGTGAt                                                                                                                                                                            >  1:518208/1‑237 (MQ=255)
                                          cGGCGCGAACTTTGATGAAGCGAAAGCCAAAGCGATCGAACTGTCACAGCAGCAGGGGTTCACCTGGGTGCCGCCGTTCGACCATCCGATGGTGATTGCCGGGCAAGGCACGCTGGCGCTGGAACTGCTCCAGCAGGACGCCCATCTCGACCGCGTATTTGTGCCAGTCGGCGGCGGCGGTCTGGCTGCTGGCGTGGCGGTGCTGATCAAACAACTGATGCCGCAAATCAAAGTGAt                                                                                                                                                                            <  2:518208/237‑1 (MQ=255)
                                                            aGCGAAAGCCAAAGCGATCGAACTGTCACAGCAGCAGGGGTTCACCTGGGTGCCGCCGTTCGACCATCCGATGGTGATTGCCGGGCAAGGCACGCTGGCGCTGGAACTGCTCCAGCAGGACGCCCATCTCGACCGCGTATTTGTGCCAGTCGGCGGCGGCGGTCTGGCTGCTGGCGTGGCGGTGCTGATCAAACAACTGATGCCGCAAATCAAAGTGATCGCCGTAGAAGCGGAAGAc                                                                                                                                                         >  1:31644/1‑238 (MQ=255)
                                                                                                            ggTGCCGCCGTTCGACCATCCGATGGTGATTGCCGGGCAAGGCACGCTGGCGCTGGAACTGCTCCAGCAGGACTCCCATCTCGACCGCGTATTTGAGCCAGTCGGCGGCGGCGGTCTGGCTGCTGGCGTGGCGGTGCTGATCAAACAACTGATGCCGCAAATCCAAGTGATCGCCGTAGAAGCGGAAGACTCCCCCTGCCTGAAAGCAGCGCTGGATGCGGGGCCTCCGGTTGa                                                                                                             >  1:717894/1‑234 (MQ=255)
                                                                                                                              tCCGATGGTGATTGCCGGGCAAGGCACGCTGGCGCTGGAACTGCTCCAGCAGGACGCCCATCTCGACCGCGTATTTGTGCCAGTCGGCGGCGGCGGTCTGGCTGCTGGCGTGGCGGTGCTGATCAAACAACTGATGCCGCAAATCAAAGTGATCGCCGTAGAAGCGGAAGACTCCgcc                                                                                                                                                   <  1:692060/178‑1 (MQ=255)
                                                                                                                                        aTTGCCGGGCAAGGCACGCTGGCGCTGGAACTGCTCCAGCAGGACGCCCATCTCGACCGCGTATTTGTGCCAGTCGGCGGCGGCGGTCTGGCTGCTGGCGTGGCGGGGCTGCTCAAACAACTGAAGCGCCAAATCAAAATGATTGCCGTAGAAACGGAAAACTCCgcctgcctttaag                                                                                                                                         >  2:692060/1‑173 (MQ=255)
                                                                                                                                         ttGCCGGGCAAGGCACGCTGGCGCTGGAACTGCTCCAGCAGGACGCCCATCTCGACCGCGTATTTGTGCCAGTCGGCGGCGGCGGTCTGGCTGCTGGCGTGGCGGTGCTGATCAAACAACTGATGCCGCAAATCAAAGTGATCGCCGTAGAAGCGGAAGACTCCGCCTGCCTGAAAGCAGCGCTGGATGCGGGTCATCCGGTTGATCTGCCGCGCGTAGGGCTATTTGCTGAAGGCGTa                                                                           <  2:31644/239‑1 (MQ=255)
                                                                                                                                                gCAAGGCACGCTGGCGCTGGAACTGCTCCAGCAGGACGCCCATCTCGACCGCGTATTTGTGCCAGTCGGCGGCGGCGGTCTGGCTGCTGGCGTGGCGGTGCTGATCAAACAACTGATGCCGCAAATCAAAGTGATCGCCGTAGAAGCGGAAGACCCCGCCCGCCCGAAAGCAGCGCTGGATGCGGGGCATCCGGGTGTTCTTCCGCGCGTAGGGCGATTTGCTGAAGGCGGTACGGTaa                                                                    >  2:595440/1‑239 (MQ=255)
                                                                                                                                                    ggCACGCTGGCGCTGGAACTGCTCCAGCAGGACGCCCATCTCGACCGCGTATTTGTGCCAGTCGGCGGCGGCGGTCTGGCTGCTGGCGTGGCGGTGCTGATCAAACAACTGATGCCGCAAATCAAAGTGATCGCCGTAGAAGCGGAAGACTCCGCCTGCCTGAAAGCAGCGCTGGATGCGGGTCATCCGGTTg                                                                                                              >  1:101999/1‑193 (MQ=255)
                                                                                                                                                    ggCACGCTGGCGCTGGAACTGCTCCAGCAGGACGCCCATCTCGACCGCGTATTTGTGCCAGTCGGCGGCGGCGGTCTGGCTGCTGGCGTGGCGGTGCTGATCAAACAACTGATGCCGCAAATCAAAGTGATCGCCGTAGAAGCGGAAGACTCCGCCTGCCTGAAAGCAGCGCTGGATGCGGGTCATCCGGTTg                                                                                                              <  2:101999/193‑1 (MQ=255)
                                                                                                                                                               gctggAACTGCTCCAGCAGGACGCCCATCTCGACCGCGTATTTGTGCCAGTCGGCGGCGGCGGTCTGGCTGCTGGCGTGGCGGTGCTGATCAAACAACTGATGCCGCAAATCAAAGTGATCGCCGTAGAAGCGGAAGACTCCGCCTGCCTGAAAGCAGCGCTGGATGCGGGTCATCCGGTTGAt                                                                                                            <  1:72313/184‑1 (MQ=255)
                                                                                                                                                                         cTCCAGCAGGACGCCCATCTCGACCGCGTATTTGTGCCAGTCGGCGGCGGCGGTCTGGCTGCTGGCGTGGCGGTGCTGATCAAACAACCGATGCCGCCAATCAAAGAGATTGCCCTAGAAGCGGAAAACTCCCCCTGCCCGAAAGCCGACCTGGCTGCGGGGCATCCGGGTGGTCTTCcccgcg                                                                                                  >  2:72313/1‑184 (MQ=255)
                                                                                                                                                                                    cGCCCATCTCGACCGCGTATTTGTGCCAGTCGGCGGCGGCGGTCTGGCTGCTGGCGTGGCGGTGCTGATCAAACAACTGATGCCGCAAATCAAAGTGATCGCCGTAGAAGCGGAAGACTCCgc                                                                                                                                                    >  1:19362/1‑123 (MQ=255)
                                                                                                                                                                                    cGCCCATCTCGACCGCGTATTTGTGCCAGTCGGCGGCGGCGGTCTGGCTGCTGGCGTGGCGGTGCTGATCAAACAACTGATGCCGCAAATCAAAGTGATCGCCGTAGAAGCGGAAGACTCCgc                                                                                                                                                    <  2:19362/123‑1 (MQ=255)
                                                                                                                                                                                     cccaagcgcggccGCGTTTTTGTGCCAGTCGGCGGCGGAGGTCGGGCTGCTGGCGTGGCGGTGCTGATCAAACAACTGATGCCGCAAAACAAAGTGATCGCCGTAGAGGCGGAAGACTCCGCCTGCCTGAAAGCAGCGCTGGATGCGGGTCATCCGGTTGATCTGCCGCGCGTAGGGCTATTTGCTGAAGGCGTAGCGGTAAAACGCATCGGTGACGAAACCTTCCGTTTATGCCAgg                                <  2:717894/227‑1 (MQ=255)
                                                                                                                                                                                                                gTCGGCGGCGGGGGTCTGGCTGATGGCGGCGCGGTGGTGATCAAACAACGAGTGGCGCAAATCAAAGTGATCGCCGTAGAAGCGGAAGACTCCGCCTGCCTGAAAGCAGCGGTGGATGCGGGTCATCCGGTTGATCTGCCGCGCGTAGGGCTATTTGCTGAAGGCGTAGCGGTAAAACGCATCGGTGACGAAACCTTCCGTTTATGCCAGGAGTATCTCGACGACATCATCACCGTCGa    <  2:867243/239‑1 (MQ=255)
                                                                                                                                                                                                                  cggcggcggcggTCTGGCTGCTGGCGTGGCGGTGCTGATCAAACAACTGATGCCGCAAATCAAAGTGATCGCCGTAGAAGCGGAAGGCTCCGCCTGCCTGAAAGCAGCGCTGGATGCGGGTCATCCGGTTGATCTGCCGCGCGTAGGGCTATTTGCTGAAGGCGTAGCGGTAAAACGCATCGGTGACGAAACCTTCCGTTTATGCCAGGAGTATCTCGACGACATCATCACCGTCGATa  <  1:595440/239‑1 (MQ=255)
                                                                                                                                                                                                                  cggcggcggcggTCTGGCTGCTGGCGTGGCGGTGCTGATCAAACAACTGATGCCGCAAATCAAAGTGATCGCCGTAGAAGCGGAAGACTCCGCCTGCCTGAAAGCAGCGCTGGATGCGGGTCATCCGGTTGATCTGCCGCGCGTAGGGCTATTTGCTGAAGGCGTAGCGGTAAAACGCATCGGTGACGAAACCTTCCGTTTATGCCAGGAGTATCTCGACGACATCATCACCGTCGATa  <  1:234182/239‑1 (MQ=255)
                                                                                                                                                                                                                        |                                                                                                                                                                                                                                        
AGTCGACGCGGTGCGCGGCTTCGGCGGCGAAGTGCTGCTCCACGGCGCGAACTTTGATGAAGCGAAAGCCAAAGCGATCGAACTGTCACAGCAGCAGGGGTTCACCTGGGTGCCGCCGTTCGACCATCCGATGGTGATTGCCGGGCAAGGCACGCTGGCGCTGGAACTGCTCCAGCAGGACGCCCATCTCGACCGCGTATTTGTGCCAGTCGGCGGCGGCGGTCTGGCTGCTGGCGTGGCGGTGCTGATCAAACAACTGATGCCGCAAATCAAAGTGATCGCCGTAGAAGCGGAAGACTCCGCCTGCCTGAAAGCAGCGCTGGATGCGGGTCATCCGGTTGATCTGCCGCGCGTAGGGCTATTTGCTGAAGGCGTAGCGGTAAAACGCATCGGTGACGAAACCTTCCGTTTATGCCAGGAGTATCTCGACGACATCATCACCGTCGATA  >  NC_000913/3955681‑3956129

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