New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913 = 87106216 (0.820)4 (0.230) 3/118 1.6 20.9% coding (96/921 nt) gsiC glutathione ABC transporter permease
?NC_000913 = 871080 16 (0.920)coding (114/921 nt) gsiC glutathione ABC transporter permease

GGTTGATTCCGACGCTGTTTATCGTCTCGGTGCTGGTGTTTTTATTTGTCCATATGCTGCCCGGCGATCCGGCGCGATTGATTGCCGGGCCCGAAGCTGATGCGCAGGTTATAGAACTGGTGCGTCAG  >  NC_000913/870998‑871125
                                                                |                                                               
ggTTGATTCCGACGCTGTTTATCGTCTCGGTGCTGGTGTTTTTATTTGTCCATATGCTGCCCGGCGAt                                                              >  2:483810/1‑68 (MQ=255)
   tGATTCCGACGCTGTTTATCGTCTCGGTGCTGGTGTTTTTATTTGTCCATATGCTGCCCGGCGATCCg                                                           >  2:583212/1‑68 (MQ=255)
          gACGCTGTTTATCGTCTCGGTGCTGGTGTTTTTATTTGTCCATATGCTGCCCGGCGATCCGGCGCgat                                                    <  2:491567/68‑1 (MQ=255)
                    aTCGTCTCGGTGCTGGTGTTTTTATTTGTCCATATGCTGCCCGGCGATCCGGCGCGATTGATTgccgg                                          >  2:131685‑M2/1‑63 (MQ=255)
                      cGTCTCGGTGCTGGTGTTTTTATTTGTCCATATGCTGCCCGGCGATCCGGCGCGATTGATTgccgggc                                        >  1:270526‑M2/1‑61 (MQ=255)
                       gTCTCGGTGCTGGTGTTTTTATTTGTCCATATGCTGCCCGGCGATCCGGCGCGATTGATTGCCGGGcc                                       >  1:245495/1‑68 (MQ=255)
                         ctcGGTGCTGGTGTTTTTATTTGTCCATATGCTGCCCGGCGATCCGGCGCGATTGATTgtcgggcagc                                     >  2:292424‑M2/1‑58 (MQ=255)
                         ctcGGTGCTGGTGTTTTTATTTGTCCATATGCTGCCCGGCGATCCGGCGCGATTGATTgccgggcagc                                     >  1:37768‑M2/1‑58 (MQ=255)
                               gCTGGTGTTTTTATTTGTCCATATGCTGCCCGGCGATCCGGCGCGATTGATTGCCGGGCCCGAAGCTg                               >  1:358126/1‑68 (MQ=255)
                                    tgtTTTTATTTGTCCATATGCTGCCCTGCGATCCGGCGCGATTGATTgccgggcagcatatggacaaa                          >  1:516662‑M2/1‑47 (MQ=255)
                                    tgtTTTTATTTGTCCATATGCTGCCCGGCaatcaatcgcgccggatcgccgggcagcatatggacaaa                          >  1:429025‑M1/1‑29 (MQ=255)
                                        tttATTTGTCCATATGCTGCCCGGCGATCCGGCGCGATTGATTGCCGGGCCCGAAGCTGATGCGCAgg                      <  1:633164/68‑1 (MQ=255)
                                         ttATTTGTCCATATGCTGCCCGGCGATCCGGCGCGATTGATTGCCGGGCCCGAAGCTGATGCGCAGGt                     <  1:549204/68‑1 (MQ=255)
                                          tATTTGTCCATATGCTGCCCGGCGATCCGGCGCGATTGATTGCCGGGCCCGATGCTGATGCGCAGGtt                    <  1:392223/68‑1 (MQ=255)
                                          tATTTGTCCATATGCTGCCCGGCGATCCGGCGCGATTGATTGCCGGGCCCGAAGCTGATGCGCAGGtt                    <  1:149467/68‑1 (MQ=255)
                                          tATTTGTCCATATGCTGCCCGGCGATCCGGCGCGATTGATTGCCGGGCCCGAAGCTGATGCGCAGGtt                    <  1:164539/68‑1 (MQ=255)
                                          tATTTGTCCATATGCTGCCCGGCGATCCGGCGCGATTGATTGCCGGGCCCGAAGCTGATGCGCAGGtt                    <  1:269939/68‑1 (MQ=255)
                                            tttGTCCATATGCTGCCCGGCaatcaatcgcgccggatcgcagggcagcatatggacaaataaaaaca                  >  1:516662‑M1/1‑21 (MQ=255)
                                            tttGTCCATATGCTGCCCGGCGATCCGGCGCGATTGATTgccgggcagcatatggacaaataaaaaca                  >  1:429025‑M2/1‑39 (MQ=255)
                                             ttGTCCATATGCTGTCCGGCGATCCGGCGCGATTGATTGCCGGGCCCGAAGCTGATGCGCAGGGtata                 <  1:18082/68‑1 (MQ=255)
                                               gTCCATATGCTGCCCGGCGATCCGGCGCGATTGATTGCCGGGCCCGAAGCTGATGCGCAGGTTATAGa               <  1:217998/68‑1 (MQ=255)
                                                   atatGCTGCCCGGCGAGCCGGCGCGATTGATTGCCGGGCCCGAAGCTGATGCGCAGGGTATAGAACTg           <  1:141806/68‑1 (MQ=255)
                                                      tgctgcCCGGCGATCCGGCGCGATTGATTGCCGGGCCCGAAGCTGATGCGCAGGTTATAGAACTGGTg        >  1:22196/1‑68 (MQ=255)
                                                       gctgcCCGGCaatcaatcgcgccggatcgccgggcagcatatggacaaataaaaacaccagcaccgag       >  1:37768‑M1/1‑10 (MQ=255)
                                                       gctgcCCGACaatcaatcgcgccggatcgccgggcagcatatggacaaataaaaacaccagcaccgag       >  2:292424‑M1/1‑10 (MQ=255)
                                                       gctgcACGGCGATCCGGCGAGATTGATTGCCGGGCCCGAAGCTGATGCGCTGGTTATAGAACTGGTGc       >  2:254517/1‑68 (MQ=255)
                                                          gcCCGGCaatcaatcgcgccggatcgccgggcagcatatggacaaataaaaacaccagcaccgagacg    >  1:270526‑M1/1‑7 (MQ=255)
                                                            ccGGCaatcaatcgcgccggatcgccgggcagcatatggacaaataaaaacaccagcaccgagacgat  >  2:131685‑M1/1‑5 (MQ=255)
                                                                |                                                               
GGTTGATTCCGACGCTGTTTATCGTCTCGGTGCTGGTGTTTTTATTTGTCCATATGCTGCCCGGCGATCCGGCGCGATTGATTGCCGGGCCCGAAGCTGATGCGCAGGTTATAGAACTGGTGCGTCAG  >  NC_000913/870998‑871125

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