New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_002947 = 519644027 (1.440)4 (0.230) 3/226 NT 14.3% coding (123/690 nt) ybgJ allophanate hydrolase subunit 1
?NC_002947 = 5196465 23 (1.330)coding (98/690 nt) ybgJ allophanate hydrolase subunit 1

CTGTGTCAGGCTGCAAGCCGTCCAGCGCCTGAGCGATTCGCGCCCTGGCCTCGCCTGGCGCCAGGGCAGACTGCACCATCAGCGTAGTATAAGACGGCACCAGGTCCAGCAAATGCTCGCCGAACGCCGCGCTCAAGCGCCGGCTGGCGGCGAGCATCCATGGCATGTTGACCTCATCGATACGGTCGAACAGGCGTACCATCAGGCTGTCGATGGCGACCACTTCGATTCGCGGCTTCATCACACCTCCAGCCCATCGAGGGCCTGGCGAATCTGCC  >  NC_002947/5196319‑5196596
                                                                                                                                                  |                                                                                                                                   
cTGTGTCAGGCTGCAAGCCGTCCAGCGCCTGAGCGATTCGCGCCCTGGCCTCGCCTGGCGCCAGGGCAGACTGCACCATCAGCGTAGTATAAGACGGCACCAGGTCCAGCAAATGCTCGCCGAACGCCGCGCTCAAGCGCCGGCTGGcgg                                                                                                                                  >  2:405173‑M2/1‑147 (MQ=255)
            gCAAGCCGTCCAGCGCCTGAGCGATTCGCGCCCTGGCCTCGCCTGGCGCCAGGGCAGACTGCACCATCAGCGTAGTATAAGACGGCACCAGGTCCAGCAAATGCTCGCCGAACGCCGCGCTCAAGCGCCGGCTGGCGGCGAGCATCCATg                                                                                                                      <  2:221571/150‑1 (MQ=255)
                               aGCGATTCGCGCCCTGGCCTCGCCTGGCGCCAGGGCAGACTGCACCATCAGCGTAGTATAAGACGGCACCAGGTCCAGCAAATGCTCGCCGAACGCCGCGCTCAAGCGCCGGCTGGCGGCGAGCATCCATGGCATGtt                                                                                                               >  2:480471/1‑138 (MQ=255)
                               aGCGATTCGCGCCCTGGCCTCGCCTGGCGCCAGGGCAGACTGCACCATCAGCGTAGTATAAGACGGCACCAGGTCCAGCAAATGCTCGCCGAACGCCGCGCTCAAGCGCCGGCTGGCGGCGAGCATCCATGGCATGtt                                                                                                               <  1:480471/138‑1 (MQ=255)
                                    ttCGCGCCCTGGCCTCGCCTGGCGCCAGGGCAGACTGCACCATCAGCGTAGTATAAGACTGCACCAGGTCCAGCAAATGCTCGCCGAACGCCGCGCTCAAGCGCCGGCTGGCGGCGAGCATCCATGGCATGTTGACCTCATCGATACGGt                                                                                              >  1:216365/1‑150 (MQ=255)
                                     tCGCGCCCTGGCCTCGCCTGGCGCCAGGGCAGACTGCACCATCAGCGTAGTATAAGACGGCACCAGGTCCAGCAAATGCTCGCCGAACGCCGCGCTCAAGCGCCGGCTGGCGGCGAGCATCCATGGCATGTTGACCTCATCGATACGGTc                                                                                             <  2:388273/150‑1 (MQ=255)
                                                    gCCTGGCGCCAGGGCAGACTGCACCATCAGCGTAGTATAAGACGGCACCAGGTCCAGCAAATGCTCGCCGAACGCCGCGCTCAAGCGCCGGCTGGCGGCGAGCATCCATGGCATGTTGACCTCATCGATACGGTCGAACAGGCGTACCAt                                                                              >  1:186434/1‑150 (MQ=255)
                                                              aGGGCAGACTGCACCATCAGCGTAGTATAAGACGGCACCAGGTCCAGCAAATGCTCGCCGAACGCCGCGCTCAAGCGCCGGCTGGCGGCGAGCATCCATGGCATGTTGACCTCATCGATACGGTCGAACAGGCGTACCATCAGGCt                                                                        >  2:383669/1‑146 (MQ=255)
                                                              aGGGCAGACTGCACCATCAGCGTAGTATAAGACGGCACCAGGTCCAGCAAATGCTCGCCGAACGCCGCGCTCAAGCGCCGGCTGGCGGCGAGCATCCATGGCATGTTGACCTCATCGATACGGTCGAACAGGCGTACCATCAGGCt                                                                        <  1:383669/146‑1 (MQ=255)
                                                                      cTGCACCATCAGCGTAGTATAAGACGGCACCAGGTCCAGCAAATGCTCGCCGAACGCCGCGCTCAAGCGCCGGCTGGCGGCGAGCATCCATGGCATGTTGACCTCATCGATACGGTCGAACAGGCGTACCATCAGGCTGTCGATGGCGAc                                                            >  2:339689/1‑150 (MQ=255)
                                                                       tGCACAATCAGCGTAGTATAAGACGGCACCAGGTCCAGCAAATGCTCGCCGAACGCCGCGCTCAAGCGCCGGCTGGCGGCGAGCATCCATGGCATGTTGACCTCATCGATACGGTCGAACAGGCGTACCATCAGGCTGTCGATGGCGAcc                                                           >  2:506982/1‑150 (MQ=255)
                                                                                    tagtaTAAGACGGCACCAGGTCCAGCAAATGCTCGCCGAACGCCGCGCTCAAGCGCCGGCTGGCGGCGAGCATCCATGGCATGTTg                                                                                                              <  2:65755/86‑1 (MQ=255)
                                                                                    tagtaTAAGACGGCACCAGGTCCAGCAAATGCTCGCCGAACGCCGCGCTCAAGCGCCGGCTGGCGGCGAGCATCCATGGCATGTTg                                                                                                              >  1:65755/1‑86 (MQ=255)
                                                                                     agtaTAAGACGGCACCAGGTCCAGCAAATGCTCGCCGAACGCCGCGCTCAAGCGCCGGCTGGCGGCGAGCATCCATGGCATGTTGACCTCATCGATACGGTCGAACAGGCGTACCATCAGGCTGTCGATGGCGACCACTTCGATTCGCgg                                             <  1:333790/150‑1 (MQ=255)
                                                                                      gtaTAAGACGGCACCAGGTCCAGCAAATGCTCGCCGAACGCCGCGCTCAAGCGCCGGCTGGCGGCGAGCATCCATGGCATGTTGACCTCATCGATACGGTCGAACAGGCGTACCATCAGGCTGTCGATGGCGACCACTTCGATTCGCGGc                                            <  2:126425/150‑1 (MQ=255)
                                                                                       aataAGACGGCACCAGGTCCAGCAAATGCTCGCCGAACGCCGCGCTCAAGCGCCGGCTGGCGGCGAGCATCCATGGCATGTTGACCTCATCGATACGGTCGAACAGGCGTACCATCAGGCTGTCGATGGCGACCACTTCGATTCGCGGCt                                           <  2:216365/149‑1 (MQ=255)
                                                                                                gCACCAGGTCCAGCAAATGCTCGCCGAACGCCGCGCTCAAGCGCCGGCTGGCGGCGAGCATCCATGGCATGTTGACCTCATCGATACGGTCGAACAGGCGTACCATCAGGCTGTCGATGGCGACCACTTCGATTCGCGGCTTCATCacac                                  >  1:22557/1‑150 (MQ=255)
                                                                                                     aGGTCCAGCAAATGCTCGCCGAACGCCGCGCTCAAGCGCCGGCTGGCGGCGAGCATCCATGGCATGTTGACCTCATCGATACGGTCGAACAGGCGTACCATCAGGCTGTCGATGGCGa                                                             >  1:214646/1‑118 (MQ=255)
                                                                                                     aGGTCCAGCAAATGCTCGCCGAACGCCGCGCTCAAGCGCCGGCTGGCGGCGAGCATCCATGGCATGTTGACCTCATCGATACGGTCGAACAGGCGTACCATCAGGCTGTCGATGGCGa                                                             <  2:214646/118‑1 (MQ=255)
                                                                                                          cAGCAAATGCTCGCGGAACGCCGCGCTCAAGCGCCGGCTGGcggcgagcatttgctggacctggtgccgtcttatactacgctgatggtgcagtctgccctggcgccaggcgaggcaagggcgcgaatcgctcaggcgctggacggcttg                        <  1:17513‑M2/150‑110 (MQ=255)
                                                                                                           aGCAAATGCTCGCCGAACGCCGCGCTCAAGCGCCGGCTGGCGGCGAGCATCCATGGCATGTTGACCTCATCGATACGGTCGAACAGGCGTACCATCAg                                                                           >  1:495750/1‑98 (MQ=255)
                                                                                                           aGCAAATGCTCGCCGAACGCCGCGCTCAAGCGCCGGCTGGCGGCGAGCATCCATGGCATGTTGACCTCATCGATACGGTCGAACAGGCGTACCATCAg                                                                           <  2:495750/98‑1 (MQ=255)
                                                                                                            gCAAATGCTCGCCGAACGCCGCGCTCAAGCGCCGGCTGGCGGCGAGCATCCATGGCATGTTGACCTCATCGATACGGTCGAACAGGCGTACCATCAGGCTGTCGATGGCGACCACTTCGATTCGCGGCTTCATCACACCTCCAGCCCATc                      <  2:216830/150‑1 (MQ=255)
                                                                                                               aaTGCTCGCCGAACGCCGCGCTCAAGCGCTGGCTGGCGGCGAGCATCCATGGCATGTTGACCTCATCGATACGGTCGAACAGGCGTACCATCAGGCTGTCGATGGCGACCACTTCGATTCGCGGCTTCATCACACCTCCAGCCCATCGAg                   <  1:206147/150‑1 (MQ=255)
                                                                                                               aaTGCTCGCCGAACGCCGCGCTCAAGCGCCGGCTGGCGGCGAGCATCCATGGCATGTTGACCTCATCGATACGGTCGAACAGGCGTACCATCAGGCTGTCGATGGCGACCACTTCGATTCGCGGCTTCATCACACCTCCAGCCCATCGAg                   <  1:339689/150‑1 (MQ=255)
                                                                                                                    tCGCCGAACGCCGCGCTCAAGCGCCGGCTGGcggcgagcatttgctggacctggtgccgtcttatactacgcttatggtgcagtctgccctggcgccaggcgaggccagggcgcgaatcgct                                          <  1:190456‑M2/122‑92 (MQ=255)
                                                                                                                    tCGCCGAACGCCGCGCTCAAGCGCCGGCTGGcggcgagcatttgctggacctggtgccgtcttatactacgctgatggtgcagtctgccctggcgccaggcgaggccagggcgcgaatcgct                                          >  2:190456‑M2/1‑31 (MQ=255)
                                                                                                                    tCGCCGAACGCCGCGCTCAAGCGCCGGCTGGcggcgagcatttgctggacctggtgccgtcttataataggctgatggtgcagtctgccctggcgccaggcgaggccagggcgcgagtcgctcaggcgctggacggcttgcagcctaaca              <  2:198563‑M2/150‑120 (MQ=255)
                                                                                                                                gcgcTCAAGCGCCGGCTGGCGGCGAGCATCCATGGCATGTTGACCTCATCGATACGGTCGAACAGGCGTACCATCAGGCTGTCGATGGCGACCACTTCGATTCGCGGCTTCATCACACCTCCAGCCCATCGAGGGCCTGGCGAATCTGcc  <  1:94289/150‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                   
CTGTGTCAGGCTGCAAGCCGTCCAGCGCCTGAGCGATTCGCGCCCTGGCCTCGCCTGGCGCCAGGGCAGACTGCACCATCAGCGTAGTATAAGACGGCACCAGGTCCAGCAAATGCTCGCCGAACGCCGCGCTCAAGCGCCGGCTGGCGGCGAGCATCCATGGCATGTTGACCTCATCGATACGGTCGAACAGGCGTACCATCAGGCTGTCGATGGCGACCACTTCGATTCGCGGCTTCATCACACCTCCAGCCCATCGAGGGCCTGGCGAATCTGCC  >  NC_002947/5196319‑5196596

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