New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NZ_CP009273 3022788 =10 (0.460)3 (0.150) 3/156 1.6 23.8% coding (1503/1920 nt) ygfT formate‑dependent uric acid utilization protein YgfT
?NZ_CP009273 3022819 = 10 (0.500)coding (1472/1920 nt) ygfT formate‑dependent uric acid utilization protein YgfT

GTACGAATCAGGCCCACCGCAGTTAAACGCCCATCTTCGTCACAAGCGATATATTGCGGTTGAACATTGAACTGAAACTCGACACCTTCCTCGCGCGCATTAACCACCTCTTTGCGCGAGCCCGGCATACTGACTTCATCACGACGATACGCGCAGGTCACGCTGGCGGCATTGA  >  NZ_CP009273/3022699‑3022873
                                                                                         |                                                                                     
acgacaatggattgtttgcggacttccatccgcctcaatgccgccagcgtgacctgcgcgtatcgtcgtgatgaagtcagtatgccgggc                                                                                       >  1:296887‑M1/90‑90 (MQ=255)
 tACGAATCAGGCCCACCGCAGTTAAACGCCCATCTTCGTCACAAGCGTTATATTGCGGTTGAACATTGAACTGAACCTCGACACCTTCCt                                                                                      <  1:239430/90‑1 (MQ=255)
   acaatggattgtttgcggacttccatccgcctcaatgccgccagcgtgacctgcgcgtatcgtcgtgatggagtcagtatgccgggCTcg                                                                                    >  1:231613‑M1/87‑90 (MQ=255)
                             cccATCTTCGTCACAAGCGATATATTGCGGTTGAACATTGAACTGAAACTCGACACCTTCCTCGCGCGCATTAACCACCTCTTTGCGCGa                                                          >  1:194166/1‑90 (MQ=255)
                              cgcctcaatgccgccagcgtgacctgcgcgtatcgtcgtgatgaagtcagtatgccgggCCCGCGCGCATTAACCACCTCTTTGCGCGAg                                                         <  1:375413‑M1/31‑1 (MQ=255)
                               cATCTTCGTCACAAGCGATATATTGCGGTTGAACATTGAACTGAAACTCGACACCTTCCTCGCGCGCATTAACCACCTCTTTGCGCGAGc                                                        >  2:17143/1‑90 (MQ=255)
                                             gCGATATATTGCGGTTGAACATTGAACTGAAACTCGACACCTTCCTCGCGCGCATTAACCACCTCTTTGCGCGAGCCCGGCATACTGACt                                          <  2:194166/90‑1 (MQ=255)
                                             gCGATATATTGCGGTTGAACATTGAACTGAAACTCGACACCTTCCTCGCGCGCATTAACCACCTCTTTGCGCGAGCCCGGCATACTGACt                                          <  2:73345/90‑1 (MQ=255)
                                                        cGGTTGAACATTGAACTGAAACTCGACACCTTCCTCGCGCGCATTAACCACCTCTTTGCGCGAGCCCGGCATACTGACTTCATCacgacg                               >  1:210683/1‑90 (MQ=255)
                                                            tatcgtcgtgatgaagtcagtatgccgggCTCGCGCGCATTAACCACCTCTTTGCGAGAGCCCGGCATACTGACTTCATCACGACGATAc                           <  2:332508‑M1/61‑1 (MQ=255)
                                                            tGAACATTGAACTGAAACTCGACACCTTCCTCGCGCGCATTAACCACCTCTTTGCGCGAGCCCGGCATACTGACTTCATCACGACGATAc                           >  2:435839/1‑90 (MQ=255)
                                                              tcgtcgtgatgaagtcagtatgccgggCTCGCGCGCATTAACCACCTCTTTGCGCGAGCCCGGCATACTGACTTCATCACGACGATAcgc                         <  1:241563‑M1/63‑1 (MQ=255)
                                                                cATTGACCTGAAACTCGACACCTTCCTCGCGCGCATTAACCTCCTCTTTGCGCGAGCCCGGCATACTGACTTCATCACGACGATACgcgc                       <  1:600485/90‑1 (MQ=255)
                                                                         gAAACTCGACACCTTCCTCGCGCGCATTAACCACCTCTTTGCGCGAGCCCGGCATACTGACTTCATCACGACGATACGCGCAGGTCACGc              >  2:407648/1‑90 (MQ=255)
                                                                                     cTTCCTCGCGCGCATTAACCACCTCTTTGCGCGAGCCCGGCATACTGACTTCATCACGACGATACGCGCAGGTCACGCTGGCGGCATTGa  <  1:573659/90‑1 (MQ=255)
                                                                                         |                                                                                     
GTACGAATCAGGCCCACCGCAGTTAAACGCCCATCTTCGTCACAAGCGATATATTGCGGTTGAACATTGAACTGAAACTCGACACCTTCCTCGCGCGCATTAACCACCTCTTTGCGCGAGCCCGGCATACTGACTTCATCACGACGATACGCGCAGGTCACGCTGGCGGCATTGA  >  NZ_CP009273/3022699‑3022873

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