New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913_3_pae_pgi = 374914715 (0.930)4 (0.260) 4/240 1.6 20.9% coding (571/663 nt) sgbH 3‑keto‑L‑gulonate 6‑phosphate decarboxylase
?NC_000913_3_pae_pgi = 3749165 16 (1.040)coding (589/663 nt) sgbH 3‑keto‑L‑gulonate 6‑phosphate decarboxylase

ACAGGCCAGCGGGCAACAGTGGGGCGAAGCCGATCTGGCACGCATGAAGGCGCTTTCAGATATCGGCCTTGAGCTTTCCATTACTGGCGGGATTACCCCTGCTGACCTGCCGCTGTTTAAAGATATCCGCGTGAAAGCGTTTATTGCCGGGCGGGCACTGGCAGGCGCGGCAAATCCGGCGCAAGTGGCTGGCGATTTCCATGCGCAAATCGACGCTATCTGGGGAGGTGCGCGTGCGTAATCATCAGTTAGGGATTTATGAAAAAGCGCTGGCGAAAGAT  >  NC_000913_3_pae_pgi/3748999‑3749279
                                                                                                                                                    |                                                                                                                                    
acaGGCCAGCGGGCAACAGTGGGGCGAAGCCGATCTGGCACGCATGAAGGCGCTTTCAGATATCGGCCTTGAGCTTTCCATTACTGGCGGGATTACCCCTGCTGACCTGCCGCTGTTTAAAGATATCCGCGTGAAAGCGTTTATTGCcgg                                                                                                                                     >  4:22672/1‑150 (MQ=255)
                          aaGCCGATCTGGCACGCATGAAGGCGCTTTCAGATATCGGCCTTGAGCTTTCCATTACTGGCGGGATTACCCCTGCTGACCTGCCGCTGTTTAAAGATATCCGCGTGAAAGCGTTTATTGCCGGGCGGGCACTGGCAGGCGCGGCAAATc                                                                                                           <  3:22672/150‑1 (MQ=255)
                               gATCTGGCACGCATGAAGGCGCTTTCAGATATCGGCCTTGAGCTTTCCATTACTGGCGGGATTACCCCTGCTGACCTGCCGCTGTTTAAAGATATCCGCGTGAAAGCGTTTATTGCCGGGCGGGCACTGGCAGGCGCGGCAAATCCGgcg                                                                                                      >  2:89957/1‑150 (MQ=255)
                                       aCGCATGAAGGCGCTTTCAGATATCGGCCTTGAGCTTTCCATTACTGGCGGGATTACCCCTGCTGACCTGCCGCTGTTTAAAGATATCCGCGTGAAAGCGTTTATTGCCGcgcctgccagtgcccgcccggcaataaa                                                                                                          <  3:14033‑M1/138‑29 (MQ=255)
                                       aCGCATGAAGGCGCTTTCAGATATCGGCCTTGAGCTTTCCATTACTGGCGGGATTACCCCTGCTGACCTGCCGCTGTTTAAAGATATCCGCGTGAAAGCGTTTATTGCCGcgcctgccagtgcccgcccggcaataaa                                                                                                          >  4:14033‑M1/1‑110 (MQ=255)
                                               aGGCGCTTTCAGATATCGGCCTTGAGCTTTCCATTACTGGCGGGATTACCCCTGCTGACCTGCCGCTGTTTAAAGATATCCGCGTGAAAGCGTTTATTGCCGGGCGGGCACTGGCAGGCGCGGCAAATCCGGCGCAAGTGGCTGGCGAtt                                                                                      <  1:200512/150‑1 (MQ=255)
                                                     tttCAGATATCGGCCTTGAGCTTTCCATTACTGGCGGGATTACCCCTGCTGACCTGCCGCTGTTTAAAGATATCCGCGTGAAAGCGTTTATTGCCGGGCGGGCACTGGCAGGCGCGGCAAATCCGGCGCAAGTGGCTGGCGATTTCCATg                                                                                >  1:6937/1‑150 (MQ=255)
                                                        cAGATATCGGCCTTGAGCTTTCCATTACTGGCGGGATTACCCCTGCTGACCTGCCGCTGTTTAAAGATATCCGCGTGAAAGCGTTTATTGCCGcgcctgccagtgcccgcccggcaa                                                                                                              <  1:42613‑M1/117‑25 (MQ=255)
                                                        cAGATATCGGCCTTGAGCTTTCCATTACTGGCGGGATTACCCCTGCTGACCTGCCGCTGTTTAAAGATATCCGCGTGAAAGCGTTTATTGCCGcgcctgccagtgcccgcccggcaa                                                                                                              >  2:42613‑M1/1‑93 (MQ=255)
                                                                 gCCTTGAGCTTTCCATTACTGGCGGGATTACCCCTGCTGACCTGCCGCTGTTTAAAGATATCCGCGTGAAAGCGTTTATTGCCGGGCGGGCACTGGCAGGCGCGGCAAATCCGGCGCAAGTGGCTGGCGATTTCCATGCGCAAATCGACg                                                                    >  1:134675/1‑150 (MQ=255)
                                                                          tttCCATTACTGGCGGGATTACCCCTGCTGACCTGCCGCTGTTTAAAGATATCCGCGTGAAAGCGTTTATTGCCGGGCGGGCACTGGCAGGCGCGGCAAATCCGGCGCAAGTGGCTGGCGATTTCCATGCGCAAATCGACGCTATCTggg                                                           >  4:111502/1‑150 (MQ=255)
                                                                            tCCATTACTGGCGGGATTACCCCTGCTGACCTGCCGCTGTTTAAAGATATCCGCGTGAAAGCGTTTATTGCCGGGCGGGCACTGGCAGGCGc                                                                                                                   <  2:157344‑M2/92‑2 (MQ=255)
                                                                            tCCATTACTGGCGGGATTACCCCTGCTGACCTGCCGCTGTTTAAAGATATCCGCGTGAAAGCGTTTATTGCCGGGCGGGCACTGGCAGGCGc                                                                                                                   >  1:157344‑M2/1‑91 (MQ=255)
                                                                                       cGGGATTACCCCTGCTGACCTGCCGCTGTTTAAAGATATCCGCGTGAAAGCGTTTATTGCCGGGCGGGCACTGGCAGGCGCGGCAAATCCGGCGCAAGTGGCTGGCGATTTCCATGCGCAAATCGACGCTATCTGGGGAGGTGCgcgtgc                                              <  2:134675/150‑1 (MQ=255)
                                                                                        gggATTACCCCTGCTGACCTGCCGCTGTTTAAAGATATCCGCGTGAAAGCGTTTATTGCCGGGCGGGCACTGGCAGGCGCGGCAAATCCGGCGCAAGTGGCTGGCGATTTCCATGCGCAAATCGACGCTATCTGGGGAGGCGCgcgtgcg                                             >  2:28820/1‑150 (MQ=255)
                                                                                               ccccTGCTGACCTGCCGCTGTTTAAAGATATCCGCGTGAAAGCGTTTATTGCCGGGCGGGCACTGGCAGGCGCGGCAAATCCGGCGCAAGTGGCTGGCGATTTCCATGCGCAAATCGACGCTATCTGGGGAGGTgcgc                                                  <  1:191117/138‑1 (MQ=255)
                                                                                               ccccTGCTGACCTGCCGCTGTTTAAAGATATCCGCGTGAAAGCGTTTATTGCCGGGCGGGCACTGGCAGGCGCGGCAAATCCGGCGCAAGTGGCTGGCGATTTCCATGCGCAAATCGACGCTATCTGGGGAGGTgcgc                                                  >  2:191117/1‑138 (MQ=255)
                                                                                                 ccTGCTGACCTGCCGCTGTTTAAAGATATCCGCGTGAAAGCGTTTATTGCCGGGCGGGCACTGGCAGGCGCGGCAAATCCGGCGCAAGTGGCTGGCGATTTCCATGCGCAAATCGACGCTATCTGGGGAGGTGCGCGTGCGTAAtcatca                                    >  4:99827/1‑150 (MQ=255)
                                                                                                   tgctgACCTGCCGCTGTTTAAAGATATCCGCGTGAAAGCGTTTATTGCCGGGCGGGCACTGGCAGGCGCGGCAAATCCGGCGCAAGTGGCTGGCGATTTCCATGCGCAAATNGACNCTATCTGGNGAGGTGCGCGTGCGTAATCATCAGt                                  <  3:99827/150‑1 (MQ=255)
                                                                                                                   tttAAAGATATCCGCGTGAAAGCGTTTATTGCCGGGCGGGCACTGGCAGGCGCGGCAAATCCGGCGCAAGTGGCTGGCGATTTCCATGCGCAAATCGACGCTATCTGGGGAGGTGCGCGTGCGTAATCATCAGTTAGGGATTTATGaaaa                  >  1:39039/1‑150 (MQ=255)
                                                                                                                                   tGAAAGCGTTTATTGCCGGGCGGGCACTGGCAGGCGCGGCAAATCCGGCGCAAGTGGCTGGCGATTTCCATGCGCAAATCGACGCTATCTGGGGAGGCGCGCGTGCGTAATCATCAGTTAGGGATTTATGAAAAAGCGCTGGCGAAAGAt  <  1:28820/150‑1 (MQ=255)
                                                                                                                                           tttATTGCCGGGCGGGCACTGGCAGGCGcggcaataaacgctttcacgcggatatctttaaacagcggcaggtcagcaggggtaatcccgccagtaatggaaagctcaaggccgatatctgaaagcgccttcatgcgt      <  3:14033‑M2/138‑111 (MQ=255)
                                                                                                                                           tttATTGCCGGGCGGGCACTGGCAGGCGcggcaataaacgctttcacgcggatatctttaaacagcggcaggtcagcaggggtaatcccgccagtaatggaaagctcaaggccgatatctgaaagcgccttcatgcgt      >  4:14033‑M2/1‑28 (MQ=255)
                                                                                                                                               ttGCCGGGCGGGCACTGGCAGGCGcggcaataaacgctttcacgcggatatctttaaacagcggcaggtcagcaggggtaatcccgccagtaatggaaagctcaaggccgatatctg                       >  2:42613‑M2/1‑24 (MQ=255)
                                                                                                                                               ttGCCGGGCGGGCACTGGCAGGCGcggcaataaacgctttcacgcggatatctttaaacagcggcaggtcagcaggggtaatcccgccagtaatggaaagctcaaggccgatatctg                       <  1:42613‑M2/117‑94 (MQ=255)
                                                                                                                                                    gcgcctgccagtgcccgcccggcaataaacgctttcacgcggatatctttaaacagcggcaggtcagcaggggtaatcccgccagtaatgga                                           <  2:157344‑M1/92‑92 (MQ=255)
                                                                                                                                                    gcgcctgccagtgcccgcccggcaataaacgctttcacgcggatatctttaaacagcggcaggtcagcaggggtaatcccgccagtaatgga                                           >  1:157344‑M1/1‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                    
ACAGGCCAGCGGGCAACAGTGGGGCGAAGCCGATCTGGCACGCATGAAGGCGCTTTCAGATATCGGCCTTGAGCTTTCCATTACTGGCGGGATTACCCCTGCTGACCTGCCGCTGTTTAAAGATATCCGCGTGAAAGCGTTTATTGCCGGGCGGGCACTGGCAGGCGCGGCAAATCCGGCGCAAGTGGCTGGCGATTTCCATGCGCAAATCGACGCTATCTGGGGAGGTGCGCGTGCGTAATCATCAGTTAGGGATTTATGAAAAAGCGCTGGCGAAAGAT  >  NC_000913_3_pae_pgi/3748999‑3749279

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