New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913_3_bme_pgi 3950444 =35 (0.790)3 (0.070) 3/272 6.5 8.1% intergenic (+24/‑63) ilvL/ilvX ilvG operon leader peptide/uncharacterized protein
?NC_000913_3_bme_pgi 3950470 = 35 (0.830)intergenic (+50/‑37) ilvL/ilvX ilvG operon leader peptide/uncharacterized protein
Rejected: Coverage evenness skew score above cutoff.
Rejected: Frequency below/above cutoff threshold.

CAAGATGCAAGAAAAGACAAAATGACAGCCCTTCTACGAGTGATTAGCCTGGTCGTGATTAGCGTGGTGGTGATTATTATCCCACCGTGCGGGGCTGCACTTGGACGAGGAAAGGCTTAGAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGGGTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAATGAATAACAGCACAAAATTCTGTTTCTCAAGATTCAGGACGGGGAACTAACTATGAATGGCGCACAGTGGGTGGTACATGCGTT  >  NC_000913_3_bme_pgi/3950301‑3950591
                                                                                                                                               |                                                                                                                                                   
cAAGATGCAAGAAAAGACAAAATGACAGCCCTTCTACGAGTGATTAGCCTGGTCGTGATTAGCGTGGTGGTGATTATTATCCCACCGTGCGGGGCTGCACTTGGACGAGGAAAGGCTTAGAGATCAAGCCTTAACGAACTAAGAccccc                                                                                                                                                >  2:246515/1‑149 (MQ=255)
      gCAAGAAAAGACAAAATGACAGCCCTTCTACGAGTGATTAGGATGGCCGTTATTAGCGTAGAGGGGATTATTATCCCACCGTGCGGGGCTGCACTTGGACGAGGAAAGGCTTAGAGATCAAGCCTTAACGAACTAAGACCCCCGCACCg                                                                                                                                          >  1:571036/1‑149 (MQ=255)
                aCAAAATGACAGCCCTTCTACGAGTGATTAGCCTGGTCGTGATTAGCGTGGTGGTGATTATTATCCCACCGTGCGGGGCTGCACTTGGACGAGGAAAGGCTTAGAGATCAAGCCTTAACGAACTAAGACCCCCGCCCCGAAAGGTCCgg                                                                                                                                >  2:302095/1‑149 (MQ=255)
                   aaaGGACAGACCTTTTACGAGTGATTAGCGTTGTCGGGATTAGCGTGGTGGTGATTATTATCCCACCGTGCGGGGCTGCACTTGGACGAGGAAAGGCTTAGAGATCAAGCGTTAACGAACTAAGACCCCCGCACCTAAAGGTCCggggg                                                                                                                             <  2:678261/149‑1 (MQ=255)
                      tGACAGCCCTTCTACGAGTGATTAGCCTGGTCGTGATTAGCGTGGTGGTGATTATTATCCCACCGTGCGGGGCTGCACTTGGACGAGGAAAGGCTTAGAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGGGttt                                                                                                                          <  1:704433/149‑1 (MQ=255)
                              cTTCTACGAGTGATTAGCCTGGTCGTGATTAGCGTGGTGGTGATTATTATCCCACCGTGCGGGGCTGCACTTGGACGAGGAAAGGCTTAGAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGGGTTTTTTTTGAc                                                                                                                  <  1:524040/149‑1 (MQ=255)
                              cTTCTACGAGTGATTAGCCTGGTCGTGATTAGCGTGGGGGTGATTATTATCCCACCGTGCGGGGCTGCACTTGGACGAGGAAAGGCTTAGAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGGGATTTTTTTGAc                                                                                                                  >  1:596920/1‑149 (MQ=255)
                                     gAGTGATTAGCCTGGTCGTGATTAGCGTGGTGGTGATTATTATCCCACCGTGCGGGGCTGCACTTGGACGAGGAAAGGCTTAGAGATCAAGCCTTAACGAACTAGGGTTCCCGCACCGTAAGGTGCGGGGGTTTTTTTTTTGCTTaaaa                                                                                                           >  1:34930/1‑149 (MQ=255)
                                     gAGTGATTAGCCTGGTCGTGATTAGCGTGGTGGTGATTATTATCCCACCGTGCGGGGCTGCACTTGGACGAGGAAAGGCTTAGAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAACGGTCCGGGGGATTTTTTTGGCCTTaaaa                                                                                                           >  2:689594/1‑149 (MQ=255)
                                     gAGTGATTAGCCTGGTCGTGATTAGCGTGGTGGTGATTATTATCCCACCGTGCGGGGCTGCACTTGGACGAGGAAAGGCTTAGAGATCAAGCCTTAACGAACTAAGACCCCCGAACCGAAAGGTCCGGGGGTTTTTTTTGACCTTaaaa                                                                                                           >  2:608582/1‑149 (MQ=255)
                                        tGATTAGCCTGGTCGTGATTAGCGTGGTGGTGATTATTATCCCACCGTGCGGGGCTGCACTTGGACGAGGAAAGGCTTAGAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGGGTTTTTTTTGACCTTAAAAACa                                                                                                        >  2:613277/1‑149 (MQ=255)
                                           ttAGCCTGGTCGTGATTAGCGTGGTGGTGATTATTATCCCACCGTGCGGGGCTGCACTTGGACGAGGAAAGGCTTAGAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGGGTTTTTTTTGACCTTAAAAACATaa                                                                                                     <  1:361576/149‑1 (MQ=255)
                                              gCCTGGTCGTGATTAGCGTGGTGGTGATTATTATCCCACCGTGCGGGGCTGCACTTGGACGAGGAAAGGCTTAGAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGGGTTTTTTTTGACCTTAAAAACATAACCg                                                                                                  <  1:246515/149‑1 (MQ=255)
                                                 tGGTCGTGATTAGCGTGGTGGTGATTATTATCCCACCGTGCGGGGCTGCACTTGGACGAGGAAAGGCTTAGAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGGGTTTTTTTTGACCTTAAAAACATACCGgagg                                                                                               >  2:150637/1‑149 (MQ=255)
                                                 tGGTCGTGATTAGCGTGGTGGTGATTATTATCCCACCGTGCGGGGCTGCACTTGGACGAGGAAAGGCTTAGAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGGGTTTTTTTTGACCTTAAAAACATAACCgagg                                                                                               >  1:400484/1‑149 (MQ=255)
                                                    tCGTGATTAGCGTGGTGGTGATTATTATCCCACCGTGCGGGGCTGCACTTGGACGAGGAAAGGCTTAGAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGGGTTTTTTTTGACCTTAAAAACATAACCGAGGAGc                                                                                            >  2:604560/1‑149 (MQ=255)
                                                       tGATTAGCGTGGTGGTGATTATTATCCCACCGTGCGGGGCTGCACTTGGACGAGGAAAGGCTTAGAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGGGTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGa                                                                                         <  1:604560/149‑1 (MQ=255)
                                                                                          ggggCTGCACTTGGACGAGGAAAGGCTTAGAGATCAAGCCTTAACGAACTAAGAc                                                                                                                                                    <  1:370736/55‑1 (MQ=255)
                                                                                          ggggCTGCACTTGGACGAGGAAAGGCTTAGAGATCAAGCCTTAACGAACTAAGAc                                                                                                                                                    >  2:370736/1‑55 (MQ=255)
                                                                                                   cTTGGACGAGGAAAGGCTTAGAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGGGTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAATGAATAACAGCACAAAATTc                                                                  >  1:350904/1‑128 (MQ=255)
                                                                                                   cTTGGACGAGGAAAGGCTTAGAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGGGTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAATGAATAACAGCACAAAATTc                                                                  <  2:350904/128‑1 (MQ=255)
                                                                                                        aCGAGGAAAGGCTTAGAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGGGTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAATGAATAACAGCACAAAATTCTGTTTCTCAAGATTCAGGACGGGGaa                                        <  1:644921/149‑1 (MQ=255)
                                                                                                        aCGAGGAAAGGCTTAGAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGGGTTTTTTTTGACCTTAAAAACATAACCGAGCAGCAAACAATGAATAACAGCAGAAAATTCTGTTTCCCAAGATTCCGGAAgggcac                                        >  2:695472/1‑146 (MQ=255)
                                                                                                            ggAAAGGCTTAAAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGGTCGGAGGTAATTTTAGACCTTCAAAACATAACCGAGGAGCAGAAGATGAATAAAAGCCCAAAATTCTTTTTCTCAAGATTCAGCACGGGAaactaa                                    >  2:560321/1‑149 (MQ=255)
                                                                                                             gAAAGGCTTAGAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGGGTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAATGAATAACa                                                                             <  2:95112/107‑1 (MQ=255)
                                                                                                             gAAAGGCTTAGAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGGGTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAATGAATAACa                                                                             >  1:95112/1‑107 (MQ=255)
                                                                                                                    ttAGAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGGGTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAATGAATAACAGCACAAAATTCTGTTTCTCAAGATTCAGGACGGGGAACTAACTATGAAt                            >  1:423778/1‑149 (MQ=255)
                                                                                                                     tAGAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGGTTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAATGAATAACAGCACAAAATTCTGTTTCTCAAGATTGAGGACGGGGAACTAACTATGCATg                           >  1:742269/1‑149 (MQ=255)
                                                                                                                     tAGAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGGGTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAATGAATAACAGCACAAAATTCTGTTTCTCAAGATTCAGGACGGGGAACTAACTATGAATgg                          <  1:385790/149‑1 (MQ=255)
                                                                                                                      agagATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGGGTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAATGAATAACAGCACAAAATTCTGTTTCTCAAGATTCCGGACGGGGAACTAACTATGAATgg                          >  2:15407/1‑149 (MQ=255)
                                                                                                                          aTCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTGCGGGGGTTTTTTTTGACCGTAAAAACATAACCGCGGAGCAGACAATGAATAACGGCACAAAATGCTGTTTCTCAGGATTCAGGAGGGGGAACTGACTATGAAGGgcgcg                      >  1:50806/1‑148 (MQ=255)
                                                                                                                            cAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGGGTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAATGAATAACAGCACAAAATTCTGTTTCTCAAGATTCAGGACGGGGAACTAACTATGAATGGCGcaca                    <  1:110453/149‑1 (MQ=255)
                                                                                                                                ccTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGGGTTTTTTTTGACCTTAAAAACATTACCGAGGGGCAGACAATGAATAACAGCACAAAGTTCTGTTTCTCAAGATTGTGGGCGGGGTACTAACTATGCGTGgcgcggggtgg                >  2:415400/1‑142 (MQ=255)
                                                                                                                                            aaGACCCCCGCACCGAAAGGTCCGGGGGTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAATGAATAACAGCACAAAATTCTGTTTCTCAAGATTCAGGACGGGGAACTAACTATGAATGGCGCACAGTGGGTGGTACATGCg    >  1:81119/1‑149 (MQ=255)
                                                                                                                                              gACCCCCGCACCGAAAGGTACGGGGGTTTTTTTTGACATTAAAAACCTAACCGGGGAGCAGACAATGAATAATGGCACGAAATTCTGTGTCTCAAGTTTCAGGGCGGGGATCTAACTATGAATGGCGCACAGTGGGTGGTACATGCGtt  <  2:571036/149‑1 (MQ=255)
                                                                                                                                               |                                                                                                                                                   
CAAGATGCAAGAAAAGACAAAATGACAGCCCTTCTACGAGTGATTAGCCTGGTCGTGATTAGCGTGGTGGTGATTATTATCCCACCGTGCGGGGCTGCACTTGGACGAGGAAAGGCTTAGAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGGGTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAATGAATAACAGCACAAAATTCTGTTTCTCAAGATTCAGGACGGGGAACTAACTATGAATGGCGCACAGTGGGTGGTACATGCGTT  >  NC_000913_3_bme_pgi/3950301‑3950591

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