New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913 = 2539674NA (NA)7 (0.190) 6/248 3.8 NA intergenic (‑91/+43) cysM/cysA cysteine synthase B (O‑acetylserine sulfhydrolase B)/sulfate/thiosulfate transporter subunit
?NC_000913 = 2539688 NA (NA)intergenic (‑105/+29) cysM/cysA cysteine synthase B (O‑acetylserine sulfhydrolase B)/sulfate/thiosulfate transporter subunit
Rejected: Coverage evenness skew score above cutoff.

GCTGCAACTTCACCAGAGGCGTATTGCCTATTGTTTGTTCTAATGTACTCACGATCTCTATCCATACGTGGTGTTGCCTGATGCGACGCTTGCGCGTCTTATCAGGTCTACAGGTTACAAACCTTGCCATTAAAAAGCCCGGATCGCGGGAAGCGCCTCCGGGCGTTTAACATTCACTCAACCTATCAGGCGCTTTGTGCGAGAGCAAGTTCCTCATCGCGGGTTTCGATACGCTCGTCGCCGTTATACAGCCGCGCATGTTGCAGACCAACGAATAAACGCTC  >  NC_000913/2539532‑2539815
                                                                                                                                              |                                                                                                                                             
gCTGCAACTTCACCAGAGGCGTATTGCCTATTGTTTGTTCTAATGTACTCACGATCTCTATCCATACGTGGTGTTGCCTGATGCGACGCTTGCGCGTCTTATCAGGTCTACAGGTTACAAACCTTGCCATTAAAAAGCCCGGATCGCgg                                                                                                                                         >  7:47187/1‑149 (MQ=255)
gCTGCAACTTCACCAGAGGCGTATTGCCTATTGTTTGTTCTAATGTACTCACGATCTCTATCCATACGTGGTGTTGCCTGATGCGACGCTTGCGCGTCTTATCAGGTCTACAGGTTACAAACCTTGCCATTAAAAAGCCCGGATCGCgg                                                                                                                                         <  7:66801/149‑1 (MQ=255)
  tGCAACTTCACCAGAGGCGTATTGCCTATTGTTTGTTCTAATGTACTCACGATCTCTATCCATACGTGGTGTTGCCTGATGCGACGCTTGCGCGTCTTATCAGGTCTACAGGTTACAAACCTTGCCATTAAAAAGCCCGGATCGCggg                                                                                                                                        >  5:165842/1‑148 (MQ=255)
       cTTCACCAGAGGCGTATTGCCTATTGTTTGTTCTAATGTACTCACGATCTCTATCCATACGTGGTGTTGCCTGATGCGACGCTTGCGCGTCTTATCAGGTCTACAGGTTACAAACCTTGCCATTAAAAAGCCCGGATCGCGGGaa                                                                                                                                      <  7:162959/145‑1 (MQ=255)
       cTTCACCAGAGGCGTATTGCCTATTGTTTGTTCTAATGTACTCACGATCTCTATCCATACGTGGTGTTGCCTGATGCGACGCTTGCGCGTCTTATCAGGTCTACAGGTTACAAACCTTGCCATTAAAAAGCCCGGATCGCGGGaa                                                                                                                                      >  8:162959/1‑145 (MQ=255)
               gagGCGTATTGCCTATTGTTTGTTCTAATGTACTCACGATCTCTATCCATACGTGGTGTTGCCTGATGCGACGCTTGCGCGTCTTATCAGGTCTACAGGTTACAAACCTTGCCATTAAAAAGCCCGGATCGCGGGAAGCGCCTCCGGGc                                                                                                                          >  5:85674/1‑149 (MQ=255)
                          ccTATTGTTTGTTCTAATGTACTCACGATCTCTATCAATACGTGGTGTTGCCTGATGCGACGCTTGCGCGTCTTATCAGGTCTACAGGTTACAAACCTTGCCATTAAAAAGCCCGGATCGCGGGAAGCGCCTCCGGGCGTTTAACAtt                                                                                                                >  3:117663/1‑148 (MQ=255)
                            tATTGTTTGTTCTAATGTACTCACGATCTCTATCCATACGTGGTGTTGCCTGATGCGACGCTTGCGCGTCTTATCAGGTCTACAGGTTACAAACCTTGCCATTAAAAAGCCCGGATCGCGGGAAGCGCCTCCGGGCGTTTAACATTCAc                                                                                                             >  1:70149/1‑149 (MQ=255)
                            tATTGTTTGTTCTAATGTACTCACGATCTCTATCCATACGTGGTGTTGCCTGATGCGACGCTTGCGCGTCTTATCAGGTCTACAGGTTACAAACCTTGCCATTAAAAAGCCCGGATCGCGGGAAGCGCCTCCGGGCGTTTAACATTCAc                                                                                                             >  4:159789/1‑149 (MQ=255)
                              ttgtttgtTCTAATGTACTCACGATCTCTATCGATACGTGGTGTTGCCTGATGCGACGCTTGCGCGTCTTATCAGGTCTACAGGTTACAAACCTTGCCATTAAAAAGCCCGGATCGCGGGAAGCGCCTCCGGGCGTTTAACATTCAc                                                                                                             <  4:117663/147‑1 (MQ=255)
                                        tAATGTACTCACGATCTCTATCCATACGTGGTGTTGCCTGATGCGACGCTTGCGCGTCTTATCAGGTCTACAGGTTACAAACCTTGCCATTAAAAAGCCCGGATCGCGGGAAGCGCCTCCGGGCGTTTAACATTCACTCAACCTATCAg                                                                                                 >  8:21677/1‑149 (MQ=255)
                                        tAATGTACTAACGATCTCTAACCATACTTGGTGTTGCCTGATGCGACGCTTGCGCGTATTATCAGGTCTACAGGTTACAAACCTTTCCAATAAAAAGCCCGGATCGCGGTAAGCTCCTCCTGGAGTTTAACATTCACTCAACCTATAAg                                                                                                 >  6:134623/1‑149 (MQ=255)
                                             tACTCACGATCTCTATCCATACGTGGTGTTGCCTGATGCGACGCTTGCGCGTCTTATCAGGTCTACAGGTTACAAACCTTGCCATTAAAAAGCCCGGATCGCGGGAAGCGCCTCCGGGCGTTTAACATTCACTCAACCTATCAGGCGCt                                                                                            <  8:47187/149‑1 (MQ=255)
                                               cTCACGATCTCTATCCATACGTGGTGTTGCCTGATGCGACGCTTGCGCGTCTTATCAGGTCTACAGGTTACAAACCTTGCCATTAAAAAGCCCGGATCGCGGGAAGCGCCTCCGGGCGTTTAACATTCACTCAACCTATCAGGCGCttt                                                                                          <  7:21677/149‑1 (MQ=255)
                                                tCACGATCTCTATCCATACGTGGTGTTGCCTGATGCGACGCTTGCGCGTCTTATCAGGTCTACAGGTTACAAACCTTGCCATTAAAAAGCCCGGATCGCGGGAAGCGCCTCCGGGCGTTTAACATTCACTCAACCTATCAGGCGCTTtg                                                                                         >  6:44940/1‑149 (MQ=255)
                                                 cACGATCTCTATCCATACGTGGTGTTGCCTGATGCGACGCTTGCGCGTCTTATCAGGTCTACAGGTTACAAACCTTGCCATTAAAAAGCCCGGATCGCGGGAAGCGCCTCCGGGCGTTTAACATTCACTCAACCTATCAGGCGCTTtg                                                                                         >  8:130992/1‑148 (MQ=255)
                                                 cACGATCTCTATCCATACGTGGTGTTGCCTGATGCGACGCTTGCGCGTCTTATCAGGTCTACAGGTTACAAACCTTGCCATTAAAAAGCCCGGATCGCGGGAAGCGCCTCCGGGCGTTTAACATTCACTCAACCTATCAGGCGCTTtg                                                                                         <  7:130992/148‑1 (MQ=255)
                                                  aCGATCTCTATCCATACGTGGTGTTGCCTGATGCGACGCTTGCGCGTCTTATCAGGTCTACAGGTTACAAACCTTGCCATTAAAAAGCCCGGATCGCGGGAAGCGCCTCCGGGCGTTTAACATTCACTCAACCTATCAGGCGCTTtg                                                                                         <  3:159789/147‑1 (MQ=255)
                                                    gATCTCTATCCATACGTGGTGTTGCCTGATGCGACGCTTGCGCGTCTTATCAGGTCTACAGGTTACAAACCTTGCCATTAAAAAGCCCGGATCGCGGGAAGCGCCTCCGGGCGTTTAACAtt                                                                                                                >  3:117897/1‑122 (MQ=255)
                                                    gATCTCTATCCATACGTGGTGTTGCCTGATGCGACGCTTGCGCGTCTTATCAGGTCTACAGGTTACAAACCTTGCCATTAAAAAGCCCGGATCGCGGGAAGCGCCTCCGGGCGTTTAACAtt                                                                                                                <  4:117897/122‑1 (MQ=255)
                                                          tATCCATACGTGGTGTTGCCTGATGCGACGCTTGCGCGTCTTATCAGGTCTACAGGTTACAAACCTTGCCATTAAAAAGCCCGGATCGCGGGAAGCGCCTCCgg                                                                                                                            >  6:116749/1‑104 (MQ=255)
                                                          tATCCATACGTGGTGTTGCCTGATGCGACGCTTGCGCGTCTTATCAGGTCTACAGGTTACAAACCTTGCCATTAAAAAGCCCGGATCGCGGGAAGCGCCTCCgg                                                                                                                            <  5:116749/104‑1 (MQ=255)
                                                          tATCCATACGTGGTGTTGCCTGATGCGACGCTTGCGCGTCTTATCAGGTCTACAGGTTACAAACCTTGCCATTAAAAAGCCCGGATCGCGGGAAGCGCCTCCGGGCGTTTAACATTCACTCAACCTATCAGGCGCTTTGTGCGAGAGc                                                                                >  1:71328/1‑148 (MQ=255)
                                                              cATACGTGGTGTTGCCTGATGCGACGCTTGCGCGTCTTATCAGGTCTACAGGTTACAAACCTTGCCATTAAAAAGCCCGGATCGCGGGAAGCGCCTCCGGGCGTTTAACATTCACTCAAc                                                                                                        <  2:62995/120‑1 (MQ=255)
                                                               aTACGTGGTGTTGCCTGATGCGACGCTTGCGCGTCTTATCAGGTCTACAGGTTACAAACCTTGCCATTAAAAAGCCCGGATCGCGGGAAGCGCCTCCGGGCGTTTAACATTCACTCAACCTATCAGGCGCTTTGTGCGAGAGCAAGTTc                                                                          <  2:71328/149‑1 (MQ=255)
                                                                tACGTGGTGTTGCCTGATGCGACGCTTGCGCGTCTTATCAGGTCTACAGGTTACAAACCTTGCCATTAAAAAGCCCGGATCGCGGGAAGCGCCTCCGGGCGTTTAACATTCACTCAACCTATCAGGCGCTTTGTGCGAGAGCAAGTTcc                                                                         <  5:44940/149‑1 (MQ=255)
                                                                  cGTGGTGTTGCCTGATGCGACGCTTGCGCGTCTTATCAGGTCTACAGGTTACAAACCTTGCCATTAAAAAGCCCGGATCGCGGGAAGCGCCTCCGGGCGTTTAACATTCACTCAACCTATCAGGCGCTTTGTGCGAGAGCAAGTTCCt                                                                        >  8:88606/1‑148 (MQ=255)
                                                                         ttGCCTGATGCGACGCTTGCGCGTCTTATCAGGTCTACAGGTTACAAACCTTGCCATTAAAAAGCCCGGATCGCGGGAAGCGCCTCCGGGCGTTTAACATTCACTCAACCTATCAGGCGCTTTGTGCGAGAGCAAGTTCCTCATCGCgg                                                                <  6:165842/149‑1 (MQ=255)
                                                                          tGCCTGATGCGACGCTTGCGCGTCTTATCAGGTATACAGGTTACAAACCTTGCCATTAAAAAGCCCGGATCGCGGGAAGCGCCTCCGGGCGTTTAACATTCACTCAACCTATCAGGCGCTTTGTGCGAGAGCAAGTTCCTCATCGCgg                                                                <  5:134623/148‑1 (MQ=255)
                                                                                   cGACGCTTGCGCGTCTTATCAGGTCTACAGGTTACAAACCTTGCCATTAAAAAGCCCGGATCGCGGGAAGCGCCTCCGGGCGTTTAACATTCACTCAACCTATCAGGCGCTTTGTGCGAGAGCAAGTTCCTCATCGCGGGTTTCGATAc                                                      >  7:121883/1‑149 (MQ=255)
                                                                                    gACGCTTGCGCGTCTTATCAGGTCTACAGGTTACAAACCTTGCCATTAAAAAGCCCGGATCGCGGGAAGCGCCTCCGGGCGTTTAACATTCACTCAACCTATCAGGCGCTTTGTGCGAGAGc                                                                                >  6:116797/1‑122 (MQ=255)
                                                                                    gACGCTTGCGCGTCTTATCAGGTCTACAGGTTACAAACCTTGCCATTAAAAAGCCCGGATCGCGGGAAGCGCCTCCGGGCGTTTAACATTCACTCAACCTATCAGGCGCTTTGTGCGAGAGc                                                                                <  5:116797/122‑1 (MQ=255)
                                                                                     aCGCTTGCGCGTCTTATCAGGTCTACAGGTTACAAACCTTGCCATTAAAAAGCCCGGATCGCGGGAAGCGCCTCCGGGCGTTTAACATTCACTCAACCTATCAGGCGCTTTGTGCGAGAGCAAGTTCCTCATCGCGGGTTTCGATACGc                                                    <  6:85674/149‑1 (MQ=255)
                                                                                               gTCTTATCAGGTCTACAGGTTACAAACCTTGCCATTAAAAAGCCCGGATCGCGGGAAGCGCCTCCGGGCGTTTAACATTCACTCAACCTATCAGGCGCTTTGTGCGAGAGCAAGTTCCTCATCGCGGGTTTCGATACGCTCGTCGCCGt                                          <  2:70149/149‑1 (MQ=255)
                                                                                               gTCTTATCAGGTCTACAGGTTACAAACCTTGCCATTAAAAAGCCCGGATCGCGGGAAGCGCCTCCGGGCGTTTAACATTCACTCAACCTATCAGGCGCTTTGTGCGAGAGCAAGTTCCTCATCGCGGGTTTCGATACGCTCGTCGCCGt                                          >  4:76316/1‑149 (MQ=255)
                                                                                                  ttATCAGGTCTACAGGTTACAACCCTTGCCATTAAAAAGCCCGGATCGCGGGAAGCGCCTCCGGGCGTTTAACATTCACTCAACCTATCAGGCGCTTTGTGCGAGAGCAAGTTCCTCATCGCGGGTTTCGATACGCTCGTCGCCGTta                                        >  3:162455/1‑148 (MQ=255)
                                                                                                        ggTCTACAGGTTACAAACCTTGCCATTAAAAAGCCCGGATCGCGGGAAGCGCCTCCGGGCGTTTAACATTCACTCAACCTATCAGGCGCTTTGTGCGAGAGCAAGTTCCTCATCGCGGGTTTCGATACGCTCGTCGCCGTTATACAGcc                                 <  3:76316/149‑1 (MQ=255)
                                                                                                           cTACAGGTTACAACCCTTGCCATTAAAAAGCCCGGATCGCGGTACGCGCCTCCGGGCGTTTAGCATTCACTCAACCTATCAGGCGCTTTGTGCGCGAGCAAGTTCCTCATCGCGGGTTTCGATACGCTCGTCGCCGTTATACAGCcgcg                              >  3:44514/1‑149 (MQ=255)
                                                                                                            tACAGGTTACAAACCTTGCCATTAAAAAGCCCGGATCGCGGGAAGCGCCTCCGGGCGTTTAACATTCACTCAACCTATCAGGCGCTTTGTGCGAGAGCAAGTTCCTCATCGCGGGtt                                                             <  1:178982/117‑1 (MQ=255)
                                                                                                            tACAGGTTACAAACCTTGCCATTAAAAAGCCCGGATCGCGGGAAGCGCCTCCGGGCGTTTAACATTCACTCAACCTATCAGGCGCTTTGTGCGAGAGCAAGTTCCTCATCGCGGGtt                                                             >  2:178982/1‑117 (MQ=255)
                                                                                                              cAGGTTACAAACCTTGCCATTAAAAAGCCCGGATCGCGGGAAGCGCCTCCGGGCGTTTAACATTCACTCAACCTATCAGGCGCTTTGTGCGAGAGCAAGTTCCTCATCGCGGGTTTCGATACGCTCGTCGCCGTTATACAGCCGCGCa                            <  7:88606/148‑1 (MQ=255)
                                                                                                                 gTTACAAACCTTGCCATTAAAAAGCCCGGATCGCGGGAAGCGCCTCCGGGCGt                                                                                                                        <  6:75203/53‑1 (MQ=255)
                                                                                                                 gTTACAAACCTTGCCATTAAAAAGCCCGGATCGCGGGAAGCGCCTCCGGGCGt                                                                                                                        >  5:75203/1‑53 (MQ=255)
                                                                                                                 gTTACAAACCTTGCCATTAAAAAGCCCGGATCGCGGGAAGCGCCTCCGGGCGTTTAACATTCACTCAACCTATCAGGCGCTTTGTGCGAGAGCAAGTTCCTCATCGCGGGTTTCGATACGCTCGTCGCCGTTATACAGCCGCGCATGtt                        >  8:28398/1‑149 (MQ=255)
                                                                                                                 gTTACAAACCTTGCCATTAAAAAGCCCGGATCGCGGGAAGCGCCTCCGGGCGTTTAACATTCACTCAACCTATCAGGCGCTTTGTGCGAGAGCAAGTTCCTCATCGCGGGTTTCGATACGCTCGTCGCCGTTATACAGCCGCGCATGtt                        >  5:97415/1‑149 (MQ=255)
                                                                                                                    aCAAACCTTGCCATTAAAAAGCCCGGATCGCGGGAAGCGCCTCCGGGCGTTTAACATTCACTCAACCTATCAGGCGCTTTGTGCGAGAGCAAGTTCCTCATCGCGGGTTTCGATACGCTCGTCGCCGTTATACAGCCGCGCATGTTGc                      <  7:28398/148‑1 (MQ=255)
                                                                                                                      aaaCCTTGCCATTAAAAAGCCCGGATCGCGGGAAGCGCCTCCGGGCGTTTAACATTCACTCAACCTATCAGGCGCTTTGTGCGAGAGCAAGTTCCTCATCGCGGGTTTCGATACGCTCGTCGCCGTTATACAGCCGCGCATGTTGCAGa                   >  6:184926/1‑149 (MQ=255)
                                                                                                                       aaCCTTGCCATTAAAAAGCCCGGATCGCGGGAAGCGCCTCCGGGCGTTTAACATTCACTCAACCTATCAGGCGCTTTGTGCGAGAGCAAGTTCCTCATCGCGGGTTTCGATACGCTCGTCGCCGTTATACAGCCgcgc                             <  7:12915/138‑1 (MQ=255)
                                                                                                                       aaCCTTGCCATTAAAAAGCCCGGATCGCGGGAAGCGCCTCCGGGCGTTTAACATTCACTCAACCTATCAGGCGCTTTGTGCGAGAGCAAGTTCCTCATCGCGGGTTTCGATACGCTCGTCGCCGTTATACAGCCgcgc                             >  8:12915/1‑138 (MQ=255)
                                                                                                                       aaCCTTGCCATTAAAAAGCCCGGATCGCGGGAAGCGCCTCCGGGCGTTTAACATTCACTCAACCTATCAGGCGCTTTGTGCGAGAGCAAGTTCCTCATCGCGGGTTTCGATACGCGCGTCGCCGTTATACAGCCGCGCATGTTGCAGAc                  <  4:44514/149‑1 (MQ=255)
                                                                                                                             gCCATTAAAAAGCCCGGATCGCGGGAAGCGCCTCCGGGCGTTTAACATTCACTCAACCTATCAGGCGCTTTGTGCGAGAGCAAGTTCCTCATCGCGGGTTTCGATACGCTCGTCGCCGTTATACAGCCGCGCATGTTGCAGACCAAc              <  6:97415/147‑1 (MQ=255)
                                                                                                                               cATTAAAAAGCCCGGATCGCGGGAAGCGCCTCCGGGCGTTTAACATTCACTCAACCTATCAGGCGCTTTGTGCGAGAGCAAGTTCCTCATCGCGGGTTTCGATACGCTCGTCGCCGTTATACAGCCGCGCATGTTGCAGACCAACGAAt          <  7:21085/149‑1 (MQ=255)
                                                                                                                                       aGCCCGGATCGCGGGAAGCGCCTCCGGGCGTTTAACATTCACTCAACCTATCAGGCGCTTTGTGCGAGAGCAAGTTCCTCATCGCGGGTTTCGATACGCTCGTCGCCGTTATACAGCCGCGCATGTTGCAGACCAACGAATAAACGCTc  <  5:184926/149‑1 (MQ=255)
                                                                                                                                           cGGATCGCGGGAAGCGCCTCCGGGCGTTTAACATTCACTCAACCTATCAGGCGCTTTGTGCGAGAGCAAGTTCCTCATCGCGGGTTTCGATACGCTCGTCGCCGTTATACAGCCGCGCATGTTGCAGACCAACGAATAAACGCTc  <  6:88598/145‑1 (MQ=255)
                                                                                                                                           cGGATCGCGGGAAGCGCCTCCGGGCGTTTAACATTCACTCAACCTATCAGGCGCTTTGTGCGAGAGCAAGTTCCTCATCGCGGGTTTCGATACGCTCGTCGCCGTTATACAGCCGCGCATGTTGCAGACCAACGAATAAACGCTc  >  5:88598/1‑145 (MQ=255)
                                                                                                                                              |                                                                                                                                             
GCTGCAACTTCACCAGAGGCGTATTGCCTATTGTTTGTTCTAATGTACTCACGATCTCTATCCATACGTGGTGTTGCCTGATGCGACGCTTGCGCGTCTTATCAGGTCTACAGGTTACAAACCTTGCCATTAAAAAGCCCGGATCGCGGGAAGCGCCTCCGGGCGTTTAACATTCACTCAACCTATCAGGCGCTTTGTGCGAGAGCAAGTTCCTCATCGCGGGTTTCGATACGCTCGTCGCCGTTATACAGCCGCGCATGTTGCAGACCAACGAATAAACGCTC  >  NC_000913/2539532‑2539815

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