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.

CAACTACCGAGTCGGCAATAAAACATTGCGGACCAATCGTTGCTTCCGCTTCAATAACCGCATTTTTAATCTCGGTGGCAAAACCAATTTTTACGCCATTGCTGATTATTGTGCCAGGACGAATAAACGCATAATTACTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGTTTTTAGACTAGTCATTGGAGAACAGATGATTGATGTCTTAGGGCCGGAGAAACGCAGA  >  NC_000913/381122‑381383
                                                                                                                                        |                                                                                                                             
cAACTACCGAGTCGGCAATAAAACATTGCGGACCAATCGTTGCTTCCGCTTCAATAACCGCATTTTTAATCTCGGTGGCAAAACCAATTTTTACGCCATTGCTGATTATTGTGCCAGGACGAATAAACGCATAATTACt                                                                                                                             >  1:409598/1‑139 (MQ=255)
              gCAATAAAACATTGCGGACCAATCGTTGCTTCCGCTTCAATAACCGCATTTTTAATCTCGGTGGCAAAACCAATTTTTACGCCATTGCTGATTATTGTGCCAGGACGAATAAACGCATAATTACTAGACTGGCCCCCTg                                                                                                               >  2:277894/1‑139 (MQ=255)
                       cATTGCGGACCAATCGTTGCTTCCGCTTCAATAACCGCATTTTTAATCTCGGTGGCAAAACCAATTTTTACGCCATTGCTGATTATTGTGCCAGGACGAATAAACGCATAATTACTAGACTGGCCCCCTGAATCTCCAg                                                                                                      <  2:564417/139‑1 (MQ=255)
                         ttGCGGACCAATCGTTGCTTCCGCTTCAATAACCGCATTTTTAATCTCGGTGGCAAAACCAATTTTTACGCCATTGCTGATTATTGTGCCAGGACGAATAAACGCATAATTACTAGACTGGCCCCCTGAATCTCCAGAc                                                                                                    <  2:43115/139‑1 (MQ=255)
                                              cGCTTCAATAACCGCATTTTTAATCTCGGTGGCAAAACCAATTTTTACGCCATTGCTGATTATTGTGCCAGGACGAATAAACGCATAATTACTAGACTGGccc                                                                                                                   <  2:802145/103‑1 (MQ=255)
                                              cGCTTCAATAACCGCATTTTTAATCTCGGTGGCAAAACCAATTTTTACGCCATTGCTGATTATTGTGCCAGGACGAATAAACGCATAATTACTAGACTGGccc                                                                                                                   >  1:802145/1‑103 (MQ=255)
                                                           gCATTTTTAATCTCGGTGGCAAAACCAATTTTTACGCCATTGCTGATTATTGTGCCAGGACGAATAAACGCATAATTACTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAAt                                                                  >  1:560999/1‑139 (MQ=255)
                                                              tttttAATCTCGGTGGCAAAACCAATTTTTACGCCATTGCTGATTATTGTGCCAGGACGAATAAACGCATAATTACTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACt                                                               <  2:409598/139‑1 (MQ=255)
                                                                         ggTGGCAAAACCAATTTTTACGCCATTGCTGATTATTGTGCCAGGACGAATAAACGCATAATTACTAGACTGGCCCCCTGAATCTcc                                                                                                        >  1:464307/1‑87 (MQ=255)
                                                                         ggTGGCAAAACCAATTTTTACGCCATTGCTGATTATTGTGCCAGGACGAATAAACGCATAATTACTAGACTGGCCCCCTGAATCTcc                                                                                                        <  2:464307/87‑1 (MQ=255)
                                                                             gCAAAACCAATTTTTACGCCATTGCTGATTATTGTGCCAGGACGAATAAACGCATAATTACTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGTTTTTAGACTAGt                                                <  1:277894/139‑1 (MQ=255)
                                                                                                 aTTGCTGATTATTGTGCCAGGACGAATAAACGCATAATTACTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGAt                                                                           <  2:76026/92‑1 (MQ=255)
                                                                                                 aTTGCTGATTATTGTGCCAGGACGAATAAACGCATAATTACTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGAt                                                                           >  1:76026/1‑92 (MQ=255)
                                                                                                          tattGTGCCAGGACGAATAAACGCATAATTACTAGACTGGCCCCCTGAATCTCCAGACAACCAat                                                                                             <  1:617075/65‑1 (MQ=255)
                                                                                                          tattGTGCCAGGACGAATAAACGCATAATTACTAGACTGGCCCCCTGAATCTCCAGACAACCAat                                                                                             >  2:617075/1‑65 (MQ=255)
                                                                                                             tgtgCCAGGACGAATAAACGCATAATTACTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGttt                                                          >  1:757917/1‑97 (MQ=255)
                                                                                                             tgtgCCAGGACGAATAAACGCATAATTACTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGttt                                                          <  2:757917/97‑1 (MQ=255)
                                                                                                               tgCCAGGACGAATAAACGCATAATTACTAGACTGGCCCCCTGAATCTCCAGa                                                                                                     >  2:454271/1‑52 (MQ=255)
                                                                                                               tgCCAGGACGAATAAACGCATAATTACTAGACTGGCCCCCTGAATCTCCAGa                                                                                                     >  1:902237/1‑52 (MQ=255)
                                                                                                               tgCCAGGACGAATAAACGCATAATTACTAGACTGGCCCCCTGAATCTCCAGa                                                                                                     <  2:902237/52‑1 (MQ=255)
                                                                                                               tgCCAGGACGAATAAACGCATAATTACTAGACTGGCCCCCTGAATCTCCAGa                                                                                                     <  1:454271/52‑1 (MQ=255)
                                                                                                               tgCCAGGACGAATAAACGCATAATTACTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAAt                                                                                  >  1:200489/1‑71 (MQ=255)
                                                                                                               tgCCAGGACGAATAAACGCATAATGACTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAAt                                                                                  <  2:200489/71‑1 (MQ=255)
                                                                                                                    ggACGAATAAACGCATAATTACTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGTTTTTAGACTAGTCATTGGAGAACAGATGATTGATGTCTTAGGGCCGGAGaa         >  1:125118/1‑139 (MQ=255)
                                                                                                                     gACGAATAAACGCATAATTACTAGACTGGCCCCCTGAATCTCCAGACAACCAat                                                                                             <  1:68886/54‑1 (MQ=255)
                                                                                                                     gACGAATAAACGCATAATTACTAGACTGGCCCCCTGAATCTCCAGACAACCAat                                                                                             >  2:68886/1‑54 (MQ=255)
                                                                                                                      aCGAATAAACGCATAATTACTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGTTTTTAGACTAGTCATTGGAGAACAGATGATTGATGTCTTAGGGCCGGAGAAAc       >  1:195583/1‑139 (MQ=255)
                                                                                                                         aaTAAACGCATAATTACTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGTTTTTAGACTAGTCATTGGAGAACAGATGATTGATGTCTTAGGGCCGGAGAAACGCa    >  2:910294/1‑139 (MQ=255)
                                                                                                                           tAAACGCATAATTACTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGTTTTTAGACTAGTCATTGGAGAACAGATGATTGATGTCTTAGGGCCGGAGAAACGCAGa  <  2:195583/139‑1 (MQ=40)
                                                                                                                           tAAACGCATAATTACTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGTTTTTAGACTAGTCATTGGAGAACAGATGATTGATGTCTTAGGGCCGGAGAAACGCAGa  <  2:125118/139‑1 (MQ=40)
                                                                                                                           tAAACGCATAATTACTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGTTTTTAGACTAGTCATTGGAGAACAGATGATTGATGTCTTAGGGCCGGAGAAACGCAGa  <  1:337775/139‑1 (MQ=40)
                                                                                                                                        |                                                                                                                             
CAACTACCGAGTCGGCAATAAAACATTGCGGACCAATCGTTGCTTCCGCTTCAATAACCGCATTTTTAATCTCGGTGGCAAAACCAATTTTTACGCCATTGCTGATTATTGTGCCAGGACGAATAAACGCATAATTACTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGTTTTTAGACTAGTCATTGGAGAACAGATGATTGATGTCTTAGGGCCGGAGAAACGCAGA  >  NC_000913/381122‑381383

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