New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913 381258 =NA (NA)4 (0.110) 3/206 3.9 NA pseudogene (3/417 nt) yaiX pseudogene, interrupted by IS2A, hexapeptide transferase superfamily, putative transferase
?NC_000913 = 382589 NA (NA)noncoding (1330/1331 nt) IS2 repeat region
Rejected: Coverage evenness skew score above cutoff.

TTCGAGCATTATAACGAATGGCATCCGCATAGTGCGCTGGGTTATCGCTCGCCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAGGGGCAAATCCAATTACCTATCAGGCAGTTTGCGCCAATCACTACCGGACCCTGAATAACGGCACCATGACAAATACGGGTATTCGCGCCGATAACAACTTCACCTGCAGTTTCATCAATAATTACCTGGTCTGCAATCTGAA  >  NC_000913/382463‑382721
                                                                                                                              |                                                                                                                                    
ttCGAGCATTATAACGAATGGCATCCGCATAGTGCGCTGGGTTATCGCTCGCCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAGGGGCAAATCCAATTACCTATCa                                                                                                                          >  1:871206/1‑139 (MQ=37)
           taACGAATGGCATCCGCATAGTGCGCTGGGTTATCGCTCGCCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAGGGGCAAATCCAATTACCTATCAGGCAGTTTgcg                                                                                                               <  2:36086/139‑1 (MQ=255)
           taACGAATGGCATCCGCATAGTGCGCTGGGTTATCGCTCGCCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAGGGGCAAATCCAATTACCTATCAGGCAGTTTgcg                                                                                                               <  1:26246/139‑1 (MQ=255)
                                       ggTTATCGCTCGCCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAGGGGCAAATCCAATTACCTATCAGGCAGTTTGCGCCAATCACTACCGGACCCTGAATAACgg                                                                                   >  2:866326/1‑139 (MQ=255)
                                          tATCGCTCGCCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAGGGGCAAATCCAATTACCTATCAGGCAGTTTGCGCCAATCACTACCGGACCCTGAATAACGGCAc                                                                                <  2:871206/139‑1 (MQ=255)
                                          tATCGCTCGCCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAGGGGCAAATCCAATTACCTATCAGGCAGTTTGCGCCAATCACTACCGGACCCTGAATAACGGCAc                                                                                >  2:671078/1‑139 (MQ=255)
                                                                 gCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAGGGGCAAATCCAATTACCTATCAGGCAGTTTGCGCCAATCACTACCGGACCCTGAATAACGGCACCATGACAAATACGGGTATTcgcg                                                         <  1:866326/139‑1 (MQ=255)
                                                                        cGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAGGGGCAAATCCAATTACCTATCAGGCAGTTTGCGCCAAt                                                                                                          >  2:813756/1‑83 (MQ=255)
                                                                        cGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAGGGGCAAATCCAATTACCTATCAGGCAGTTTGCGCCAAt                                                                                                          <  1:813756/83‑1 (MQ=255)
                                                                        cGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAGGGGCAAATCCAATTACCTATCAGGCAGTTTGCGCCAATCACTACCGGACCCTGAATAACGGCACCATGACAAAt                                                                      >  2:677283/1‑119 (MQ=255)
                                                                        cGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAGGGGCAAATCCAATTACCTATCAGGCAGTTTGCGCCAATCACTACCGGACCCTGAATAACGGCACCATGACAAAt                                                                      <  1:677283/119‑1 (MQ=255)
                                                                                   tGGGTTAAGTGATAACAGATGTCTGGAAATATAGGGGCAAATCCAATTACCTATCAGGCAGTTTGCGCCAATCACTACCGGACCCTGAATAACGGCACCATGACAAATACGGGTATTCGCGCCGATAACAACTTCACCt                                       <  1:671078/139‑1 (MQ=255)
                                                                                     ggTTAAGTGATAACAGATGTCTGGAAATATAGGGGCAAATCCAATTACCTATCAGGCAGTTTGCGCCAATCACTACCGGAccc                                                                                             >  2:169552/1‑83 (MQ=255)
                                                                                     ggTTAAGTGATAACAGATGTCTGGAAATATAGGGGCAAATCCAATTACCTATCAGGCAGTTTGCGCCAATCACTACCGGAccc                                                                                             <  1:169552/83‑1 (MQ=255)
                                                                                       ttAAGTGATAACAGATGTCTGGAAATATAGGGGCAAATCCAATTACCTATCAGGCAGTTTGCGCCAATCACTACCg                                                                                                  <  2:59966/76‑1 (MQ=255)
                                                                                       ttAAGTGATAACAGATGTCTGGAAATATAGGGGCAAATCCAATTACCTATCAGGCAGTTTGCGCCAATCACTACCg                                                                                                  >  1:59966/1‑76 (MQ=255)
                                                                                          aGTGATAACAGATGTCTGGAAATATAGGGGCAAATCCAATTACCTATCAGGCAGTTTGCGCCAATCACTACCGGa                                                                                                <  1:453407/75‑1 (MQ=255)
                                                                                          aGTGATAACAGATGTCTGGAAATATAGGGGCAAATCCAATTACCTATCAGGCAGTTTGCGCCAATCACTACCGGa                                                                                                >  2:453407/1‑75 (MQ=255)
                                                                                          aGTGATAACAGATGTCTGGAAATATAGGGGCAAATCCAATTACCTATCAGGCAGTTTGCGCCAATCACTACCGGAcc                                                                                              >  2:119474/1‑77 (MQ=255)
                                                                                          aGTGATAACAGATGTCTGGAAATATAGGGGCAAATCCAATTACCTATCAGGCAGTTTGCGCCAATCACTACCGGAcc                                                                                              <  1:119474/77‑1 (MQ=255)
                                                                                          aGTGATAACAGATGTCTGGAAATATAGGGGCAAATCCAATTACCTATCAGGCAGTTTGCGCCAATCACTACCGGACCCTGAATaa                                                                                      <  1:46717/85‑1 (MQ=255)
                                                                                          aGTGATAACAGATGTCTGGAAATATAGGGGCAAATCCAATTACCTATCAGGCAGTTTGCGCCAATCACTACCGGACCCTGAATaa                                                                                      >  2:46717/1‑85 (MQ=255)
                                                                                          aGTGATAACAGATGTCTGGAAATATAGGGGCAAATCCAATTACCTATCAGGCAGTTTGCGCCAATCACTACCGGACCCTGAATAACGGCAc                                                                                <  1:597117/91‑1 (MQ=255)
                                                                                          aGTGATAACAGATGTCTGGAAATATAGGGGCAAATCCAATTACCTATCAGGCAGTTTGCGCCAATCACTACCGGACCCTGAATAACGGCAc                                                                                >  2:597117/1‑91 (MQ=255)
                                                                                          aGTGATAACAGATGTCTGGAAATATAGGGGCAAATCCAATTACCTATCAGGCAGTTTGCGCCAATCACTACCGGACCCTGAATAACGGCACCATGa                                                                           >  1:27327/1‑96 (MQ=255)
                                                                                          aGTGATAACAGATGTCTGGAAATATAGGGGCAAATCCAATTACCTATCAGGCAGTTTGCGCCAATCACTACCGGACCCTGAATAACGGCACCATGa                                                                           <  2:27327/96‑1 (MQ=255)
                                                                                                aaCAGATGTCTGGAAATATAGGGGCAAATCCAATTACCTATCAGGCAGTTTGCGCCAATCACTACCGGACCCTGAATAACGGCACCATGACAAATACGGGTATTCGCGCCGATAACAACTTCACCTGCAGTTTCATCaa                          >  1:379602/1‑139 (MQ=255)
                                                                                                 aCAGATGTCTGGAAATATAGGGGCAAATCCAATTACCTATCAGGCAGTTTGCGCCAATCACTACCGGACCCTGAATAACGGCACCATGACAAATACGGGTATTCGCGCCGATAACAACTTCACCTGCa                                    <  2:43859/128‑1 (MQ=255)
                                                                                                 aCAGATGTCTGGAAATATAGGGGCAAATCCAATTACCTATCAGGCAGTTTGCGCCAATCACTACCGGACCCTGAATAACGGCACCATGACAAATACGGGTATTCGCGCCGATAACAACTTCACCTGCa                                    >  1:43859/1‑128 (MQ=255)
                                                                                                 aCAGATGTCTGGAAATATAGGGGCAAATCCAATTACCTATCAGGCAGTTTGCGCCAATCACTACCGGACCCTGAATAACGGCACCATGACAAATACGGGTATTCGCGCCGATAACAACTTCACCTGCAGTTTCATCaat                         >  1:856780/1‑139 (MQ=255)
                                                                                                     aTGTCTGGAAATATAGGGGCAAATCCAATTACCTATCAGGCAGTTTGCGCCAATCACTACCGGACCCTGAATAACGGCACCATg                                                                            <  1:754353/84‑1 (MQ=255)
                                                                                                     aTGTCTGGAAATATAGGGGCAAATCCAATTACCTATCAGGCAGTTTGCGCCAATCACTACCGGACCCTGAATAACGGCACCATg                                                                            >  2:754353/1‑84 (MQ=255)
                                                                                                       gTCTGGAAATATAGGGGCAAATCCAATTACCTATCAGGCAGTTTGCGCCAATCACTACCGGACCCTGAATAACGGCACCATGACAAATACGGGTATTCGCGCCGATAACAACTTCACCTGCAGTTTCATCAATAATTAc                   >  1:197328/1‑139 (MQ=255)
                                                                                                                  taGGGGCAAATCCAATTACCTATCAGGCAGTTTGCGCCAATCACTACCGGAcc                                                                                              >  1:731841/1‑53 (MQ=255)
                                                                                                                  taGGGGCAAATCCAATTACCTATCAGGCAGTTTGCGCCAATCACTACCGGAcc                                                                                              <  2:731841/53‑1 (MQ=255)
                                                                                                                     aggCAAATCCAATGACCTATCAGGCAGCTGGCGTCAATCACTACCGTACCCTGAATAACGGCACCATGACAAATGCGGGTATTCGCGCCGATAACAACTTCACCTGCAGTTTCATCAATAATTACCTGGTCTGCAATCt     <  1:687294/138‑1 (MQ=255)
                                                                                                                        cAAATCCAATTACCTATCAGGCAGTTTGCGCCAATCACTACCGGACCCTGAATAACGGCACCATGACAAATACGGGTATTCGCGCCGATAACAACTTCACCTGCAGTTTCATCCATAATTACCTGGTCTGCAATCTGaa  >  1:217828/1‑139 (MQ=255)
                                                                                                                         aaaTCCAATTACCTATCAGGCAGTTTGCGCCAATCACTACCGGACCCTGAATAACGGCACCATGACAAATACGGGTATTCGCGCCGATAACAACTTCACCTGCAGTTTCATCaa                          >  2:11453/1‑114 (MQ=255)
                                                                                                                         aaaTCCAATTACCTATCAGGCAGTTTGCGCCAATCACTACCGGACCCTGAATAACGGCACCATGACAAATACGGGTATTCGCGCCGATAACAACTTCACCTGCAGTTTCATCaa                          <  1:11453/114‑1 (MQ=255)
                                                                                                                              |                                                                                                                                    
TTCGAGCATTATAACGAATGGCATCCGCATAGTGCGCTGGGTTATCGCTCGCCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAGGGGCAAATCCAATTACCTATCAGGCAGTTTGCGCCAATCACTACCGGACCCTGAATAACGGCACCATGACAAATACGGGTATTCGCGCCGATAACAACTTCACCTGCAGTTTCATCAATAATTACCTGGTCTGCAATCTGAA  >  NC_000913/382463‑382721

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