New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913 = 403495434 (0.940)3 (0.080) 3/272 3.8 7.5% coding (347/546 nt) hemG protoporphyrinogen oxidase
?NC_000913 = 4035048 41 (1.160)coding (441/546 nt) hemG protoporphyrinogen oxidase
Rejected: Coverage evenness skew score above cutoff.
Rejected: Frequency below/above cutoff threshold.

AACAGCTACGCGCGCAAGTTTCTGATGAACTCGCAATGGCGTCCCGATCGCTGCGCGGTCATTGCCGGGGCGCTGCGTTACCCACGTTATCGCTGGTACGACCGTTTTATGATCAAGCTGATTATGAAGATGTCAGGCGGTGAAACGGATACGCGCAAAGAAGTTGTCTATACCGATTGGGAGCAGGTGGCGAATTTCGCCCGAGAAATCGCCCATTTAACCGACAAACCGACGCTGAAATAAGCATAAAGAATAAAAAATGCGCGGTCAGAAAAT  >  NC_000913/4034911‑4035186
                                                                                                                                         |                                                                                                                                          
aaCAGCTACGCGCGCAAGTTTCTGATGAACTCGCAATGGCGTCCCGATCGCTGCGCGGTCATTGCCGGGGCGCTGCGTTACCCACGTTATCGCTGGTACGACCGTTTTATGATCAAGCTGATTATGAAGATGTCAGGCgg                                                                                                                                          <  2:570631/140‑1 (MQ=255)
  cAGCTACGCGCGCAAGTTTCTGATGAACTCGCAATGGCGTCCCGATCGCTGCGCGGTCATTGCCGGGGCGCTGCGTTACCCACGTTATCGCTGGTACGACCGTTTTATGATCAAGCTGATTATGAAGATGTCAGGCggag                                                                                                                                        >  2:357174/1‑138 (MQ=255)
  cAGCTACGCGCGCAAGTTTCTGATGAACTCGCAATGGCGTCCCGATCGCTGCGCGGTCATTGCCGGGGCGCTGCGTTACCCACGTTATCGCTGGTACGACCGTTTTATGATCAAGCTGATTATGAAGATGTCAGGCGGTg                                                                                                                                        >  2:93527/1‑140 (MQ=255)
  cAGCTACGCGCGCAAGTTTCTGATGAACTCGCAATGGCGTCCCGATCGCTGCGCGGTCATTGCCGGGGCGCTGCGTTACCCACGTTATCGCTGGTACGACCGTTTTATGATCAAGCTGATTATGAAGATGTCAGGCGGTg                                                                                                                                        >  2:104269/1‑140 (MQ=255)
     cTACGCGCGCAAGTTTCTGATGAACTCGCAATGGCGTCCCGATCGCTGCGCGGTCATTGCCGGGGCGCTGCGTTACCCACGTTATCGCTGGTACGACCGTTTTATGATCAAGCTGATTATGAAGATGTCAGGCGGTGaaa                                                                                                                                     <  2:156958/140‑1 (MQ=255)
         gcgcgcAAGTTTCTGATGAACTCGCAATGGCGTCCCGATCGCTGCGCGGTCATTGCCGGGGCGCTGCGTTACCCACGTTATCGCTGGTACGACCGTTTTATGATCAAGCTGATTATGAAGATGTCAGGCGGTGAAACGGa                                                                                                                                 >  1:107982/1‑140 (MQ=255)
                      tgatgaGCTCGCAATGGCGTCCCGATCGCTGCGCTGTCATTGCCGGGGCGCTGCGTTTCCCACGTTATCGCTGGTACGACCGTTTTATGATCAAGATGATTATGAAGATGTCAGGCGGTGAAACGGATACGCGCAaagaa                                                                                                                    <  2:584958/140‑1 (MQ=255)
                      tgatgaACTCGCAATGGCGTCCCGATCGCTGCGCGGTCATTGCCGGGGCGCTGCGTTACCCACGTTATCGCTGGTACGACCGTTTTATGATCAAGCTGATTATGAAGATGTCAGGCGGTGAAACGGATACGCGCAaagaa                                                                                                                    <  1:581588/140‑1 (MQ=255)
                       gatgaACTCGCAATGGCGTCCCGATCGCTGCGCGGTCATTGCCGGGGCGCTGCGTTACCCACGTTATCGCTGGTACGACCGTTTTATGATCAAGCTGATTATGAAGATGTCAGGCGGTGAAACGGATACGCGCAaagaag                                                                                                                   <  2:107982/140‑1 (MQ=255)
                        atgaACTCGCAATGGCGTCCCGATCGCTGCGCGGTCATTGCCGGGGCGCTGCGTTACCCACGTTATCGCTGGTACGACCGTTTTATGATCAAGCTGATTATGAAGATGTCAGGCGGTGAAACGGATACGCGCAAAGAAGt                                                                                                                  <  2:345707/140‑1 (MQ=255)
                           aaCTCGCAATGGCGTCCCGATCGCTGCGCGGTCATTGCCGGGGCGCTGCGTTACCCACGTTATCGCTGGTACGACCGTTTTATGATCAAGCTGATTATGAAGATGTCAGGCGGTGAAACGGATAcg                                                                                                                             <  2:537919/126‑1 (MQ=255)
                           aaCTCGCAATGGCGTCCCGATCGCTGCGCGGTCATTGCCGGGGCGCTGCGTTACCCACGTTATCGCTGGTACGACCGTTTTATGATCAAGCTGATTATGAAGATGTCAGGCGGTGAAACGGATAcg                                                                                                                             >  1:537919/1‑126 (MQ=255)
                             cTCGCAATGGCGTCCCGATCGCTGCGCGGTCATTGCCGGGGCGCTGCGTTACCCACGTTATCGCTGGTACGACCGTTTTATGATCAAGCTGATTATGAAGATGTCAGGCGGTGAAACGGATACGCGCAAAGAAGt                                                                                                                  <  1:353051/135‑1 (MQ=255)
                             cTCGCAATGGCGTCCCGATCGCTGCGCGGTCATTGCCGGGGCGCTGCGTTACCCACGTTATCGCTGGTACGACCGTTTTATGATCAAGCTGATTATGAAGATGTCAGGCGGTGAAACGGATACGCGCAAAGAAGt                                                                                                                  >  2:353051/1‑135 (MQ=255)
                             cTCGCAATGGCGTCCCGATCGCTGCGCGGTCATTGCCGGGGCGCTGCGTTACCCACGTTATCGCTGGTACGACCGTTTTATGATCAAGCTGATTATGAAGATGTCAGGCGGTGAAACGGATACGCGCAAAGAAGTTGTCt                                                                                                             >  1:183772/1‑140 (MQ=255)
                                gCAATGGCGTCCCGATCGCTGCGCGGTCATTGCCGGGGCGCTGCGTTACCCACGTTATCGCTGGTACGACCGTTTTATGATCAAGCTGATTATGAAGATGTCAGGCGGTGAAACGGATACGCGCAAAGAAGTTGTCtata                                                                                                          >  2:29765/1‑140 (MQ=255)
                                                           cATTGCCGGGGCGCTGCGTTACCCACGTTATCGCTGGTACGACCGTTTTATGATCAAGCTGATTATGAAGATGTCAGGCGGTGAAACGGATACGCGCAAAGAAGTTGTCTATACCGATTGGGAGCAGGTGGCGAATTTCg                                                                               <  1:489971/140‑1 (MQ=255)
                                                                   gggCGCTGCGTTACCCACGTTATCGCTGGTACGACCGTTTTATGATCAAGCTGATTATGAAGATGTCAGGCGGTGAAACGGATACGCGCAAAGAAGTTGTCTATACCGATTGGGAGCAGGTGGCGAATTTCGCCCGAGaa                                                                       >  1:416193/1‑140 (MQ=255)
                                                                       gcTGCGTTACCCACGTTATCGCTGGTACGACCGTTTTATGATCAAGCTGATTATGAAGATGTCAGGCGGTGAAACGGATACGCGCAAAGAAGTTGTCTATACCGATTGGGAGCAGGTGGTGAATTTCGCCCGAGAAATCg                                                                   <  2:32983/140‑1 (MQ=255)
                                                                                     gTTATCGCTGGTACGACCGTTTTATGATCAAGCTGATTATGAAGATGTCAGGCGGTGAAACGGATACGCGCAAAGAAGTTGTCTATACCGATTGGGAGCAGGTGGCGAATTTCGCCCGAGAAATCGCCCATTTAACCGAc                                                     >  1:144417/1‑140 (MQ=255)
                                                                                        aTCGCTGGTACGACCGTTTTATGATCAAGCTGATTATGAAGATGTCAGGCGGTGAAACGGATACGCGCAAAGAAGTTGTCTATACCGATTGGGAGCAGGTGGCGAATTTCGCCCGAGAAATCGCCCATTTAACCGACaaa                                                  >  1:106677/1‑140 (MQ=255)
                                                                                         tCGCTGGTACGACCGTTTTATGATCAAGCTGATTATGAAGATGTCAGGCGGTGAAACGGATACGCGCAAAGAAGTTGTCTATACCGATTGGGAGCAGGTGGCGAATTTCGCCCGAGAAATCGCCCATTTAACCGACAAAc                                                 <  2:183772/140‑1 (MQ=255)
                                                                                          cGCTGGTACGACCGTTTTATGATCAAGCTGATTATGAAGATGTCAGGCGGTGAAACGGATACGCGCAAAGAAGTTGTCTATACCGATTGGGAGCAGGTGGCGAATTTCGCCCGAGAAATCGCCCATTTAACCGACAAAcc                                                >  1:326580/1‑140 (MQ=255)
                                                                                          cGCTGGTACGACCGTTTTATGATCAAGCTGATTATGAAGATGTCAGGCGGAGAAACGGATACGCGCAAAGAAGTTGTCTATACCGATTGGGAGCAGGTGGCGAATTTCGCCCGAGAAATCGCCCATTTAACCGACAAAcc                                                >  2:436681/1‑140 (MQ=255)
                                                                                                tACGACCGTTTTATGATCAAGCTGATTATGAAGATGTCAGGCGGTGAAACGGATACGCGCAAAGAAGTTGTCTATACCGATTGGGAGCAGGTGGCGAATTTCGCCCGAGAAATCGCCCATTTAACCGACAAACCGACGCt                                          <  1:563794/140‑1 (MQ=255)
                                                                                                    aCCGTTTTATGATCAAGCTGATTATGAAGATGTCAGGCGGTGAAACGGATACGCGCAAAGAAGTTGTCTATACCGATTGGGAGCAGGTGGCGAATTTCGCCCGAGAAATCGCCCATTTAACCGACAAACCGACGCTGaaa                                      <  1:93527/140‑1 (MQ=255)
                                                                                                     ccGTTTTATGATCAAGCTGATTATGAAGATGTCAGGCGGTGAAACGGATACGCGCAAAGAAGTTGTCTATACCGATTGGGAGCAGGTGGCGAATTTCGCCCGAGAAATCGCCCATTTAACCGACAAACCGACGCTGAAAt                                     <  2:326580/140‑1 (MQ=255)
                                                                                                       gTTTTATGATCAAGCTGATTATGAAGATGTCAGGCGGTGAAACGGATACGCGCAAAGAAGTTGTCTATACCGATTGGGAGCAGGTGGCGAATTTCGCCCGAGAAATCGCCCATTTAACCGACAAACCGACGCTGAAATaa                                   <  1:388184/140‑1 (MQ=255)
                                                                                                         tttATGATCAAGCTGATTATGAAGATGTCAGGCGGTGAAACGGATACGCGCAAAGAAGTTGTCTATACCGATTGGGAGCAGGTGGCGAATTTCGCCCGAGAAATCGCCCATTTAACCGACAAACCGACGCTGAAATAAGc                                 >  1:192697/1‑140 (MQ=255)
                                                                                                             tGATCAAGCTGATTATGAAGATGTCAGGCGGAGAAACGGATACGCGCAAAGAAGTTGTCTATACCGATTGGGAGCAGGTGGCGAATTTCGCCCGAGAAATCGCCCATTTAACCGACAAACCGACGCTGAAATAAGCATaa                             <  1:357174/140‑1 (MQ=255)
                                                                                                                  aaGCTGATTATGAAGATGTCAGGCGGTGAAACGGATACGCGCAAAGAAGTTGTCTATACCGATTGGGAGCAGGTGGCGAATTTCGCCCGAGAAATCGCCCATTTAACCGACAAACCGACGCTGAAATAAGCATAAAGAAt                        >  1:448029/1‑140 (MQ=255)
                                                                                                                   aGCTGATTATGAAGATGTCAGGCGGTGAAACGGATACGCGCAAAGAAGTTGTCTATACCGATTGGGAGCAGGTGGCGAATTTCGCCCGAGAAATCGCCCATTTAACCGACAAACCGACGCTGaaa                                      <  1:550230/125‑1 (MQ=255)
                                                                                                                   aGCTGATTATGAAGATGTCAGGCGGTGAAACGGATACGCGCAAAGAAGTTGTCTATACCGATTGGGAGCAGGTGGCGAATTTCGCCCGAGAAATCGCCCATTTAACCGACAAACCGACGCTGaaa                                      >  2:550230/1‑125 (MQ=255)
                                                                                                                   aGCTGATTATGAAGATGTCAGGCGGTGAAACGGATACGCGCAAAGAAGTTGTCTATACCGATTGGGAGCAGGTGGCGAATTTCGCCCGAGAAATCGCCCATTTAACCGACAAACCGACGCTGAAATAAGCATAAAGAATa                       >  1:406157/1‑140 (MQ=255)
                                                                                                                            tGAAGATGTCAGGCGGTGAAACGGATACGCGCAAAGAAGTTGTCTATACCGATTGGGAGCAGGTGGCGAATTTCGCCCGAGAAATCGCCCATTTAACCGACAAACCGACGCTGAAATAAGCATAAAGAATAAAAAATgcg              <  2:376572/140‑1 (MQ=255)
                                                                                                                              aaGATGTCAGGCGGTGAAACGGATACGCGCAAAGAAGTTGTCTATACCGATTGGGAGCAGGTGGCGAATTTCGCCCGAGAAATCGCCCATTTAACCGACAAACCGACGCTGAAATAAGCATAAAGAATAAAAAATGcgcg            >  2:298630/1‑140 (MQ=255)
                                                                                                                                  tGTCAGGCGGTGAAACGGATACGCGCAAAGAAGTTGTCTATACCGATTGGGAGCAGGTGGCGAATTTCGCCCGAGAAATCGCCCATTTAACCGACAAACCGACGCTGAAATAAGCATAAAGAATAAAAAATgcg              <  1:86506/134‑1 (MQ=255)
                                                                                                                                  tGTCAGGCGGTGAAACGGATACGCGCAAAGAAGTTGTCTATACCGATTGGGAGCAGGTGGCGAATTTCGCCCGAGAAATCGCCCATTTAACCGACAAACCGACGCTGAAATAAGCATAAAGAATAAAAAATgcg              >  2:86506/1‑134 (MQ=255)
                                                                                                                                  tGTCAGGCGGTGAAACGGATACGCGCAAAGAAGTTGTCTATACCGATTGGGAGCAGGTGGCGAATTTCGCCCGAGAAATCGCCCATTTAACCGACAAACCGACGCTGAAATAAGCATAAAGAATAAAAAATGCGCGGTCa        >  1:128813/1‑140 (MQ=255)
                                                                                                                                   gTCAGGCGGTGAAACGGATACGCGCAAAGAAGTTGTCTATACCGATTGGGAGCAGGTGGCGAATTTCGCCCGAGAAATCGCCCATTTAACCGACAAACCGACGCTGAAATAAGCATAAAGAATAAAAAATGCGCGGTCAg       >  1:484957/1‑140 (MQ=255)
                                                                                                                                   gTCAGGCGGTGAAACGGATACGCGCAAAGAAGTTGTCTATACCGATTGGGAGCAGGTGGCGAATTTCGCCCGAGAAATCGCCCATTTAACCGACAAACCGACGCTGAAATAAGCATAAAGAATAAAAAATGCGCGGTCAg       >  1:524311/1‑140 (MQ=255)
                                                                                                                                      aGGCGGTGAAACGGATACGCGCAAAGAAGTTGTCTATACCGATTGGGAGCAGGTGGCGAATTTCGCCCGAGAAATCGCCCATTTAACCGACAAACCGACGCTGAAATAAGCATAAAGAATAAAAAATGCGCGGTCAGaaa    >  2:534310/1‑140 (MQ=255)
                                                                                                                                        gCGGTGAAACGGATACGCGCAAAGAAGTTGTCTATACCGATTGGGAGCAGGTGGCGAATTTCGCCCGAGAAATCGCCCATTTAACCGACAAACCGACGCTGAAATAAGCATAAAGAATAAAAAATGCGCGGTCAGAAAat  >  1:265560/1‑140 (MQ=255)
                                                                                                                                         |                                                                                                                                          
AACAGCTACGCGCGCAAGTTTCTGATGAACTCGCAATGGCGTCCCGATCGCTGCGCGGTCATTGCCGGGGCGCTGCGTTACCCACGTTATCGCTGGTACGACCGTTTTATGATCAAGCTGATTATGAAGATGTCAGGCGGTGAAACGGATACGCGCAAAGAAGTTGTCTATACCGATTGGGAGCAGGTGGCGAATTTCGCCCGAGAAATCGCCCATTTAACCGACAAACCGACGCTGAAATAAGCATAAAGAATAAAAAATGCGCGGTCAGAAAAT  >  NC_000913/4034911‑4035186

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