New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NZ_CP009273 322260 =34 (1.250)3 (0.130) 3/150 2.1 10.0% coding (699/1671 nt) betA choline dehydrogenase
?NZ_CP009273 322363 = 24 (1.000)coding (596/1671 nt) betA choline dehydrogenase
Rejected: Frequency below/above cutoff threshold.

GAGGCAATCGCGCCTGCACATAACAGCACTTCTTTGTTGGCCGTTGCGCGGGTTGGGATGGTGCTGTCGCCTTCCAGCCATTCGACGCCCACCGCGCGTTTGCCGTCAAAAATGATGTGATCGGTCATAGCGTGAGTACGAATGGTCAGGTTAGGACGCGATTTGGCCTGATCG  >  NZ_CP009273/322171‑322344
                                                                                         |                                                                                    
gAGGCAATCGCGCCTGCACATAACAGCACTTCTTTGTTGGCCGTTGCGCGGGTTGGGATGGTGCTGTCGCCTTCCAGCCATTCGACGccc                                                                                      >  2:573296/1‑90 (MQ=255)
gAGGCAATCGCGCCTGCACATAACAGCACTTCTTTGTTGGCCGTTGCGCGGGTTGGGATGGTGCTGTCGCCTTCCAGCCATTCGACGccc                                                                                      >  2:204798/1‑90 (MQ=255)
gAGGCAATCGCGCCTGCACATAACAGCACTACTTTGTTGGCCGTTGCGCGGGTTGGGATGGTGCTGTCGCCTTCCAGCCATTCGACGccc                                                                                      <  2:770840/90‑1 (MQ=255)
       gcacggacgatctcaacggttatcagcaggaaggttttggtccgatggatcgcaccgtcacgccgcagggccgtcgcgccagCACcgcgc                                                                               <  1:437323‑M1/8‑1 (MQ=255)
                    tAACAGCACTTCTTTGTTGGCCGTTGCGCGGGTTGGGATGGTGCTGTCGCCTTCCAGCCATTCGACGCCCACCGCGCGTTTGCCGTCaaa                                                                  >  1:259372/1‑90 (MQ=255)
                    tAACAGCACTTCTTTGTTGGCCGTTGCGCGGGTTGGGATGGTGCTGTCGCCTTCCAGCCATTCGACGCCCACCGCGCGTTTGCCGTCaaa                                                                  >  2:3635/1‑90 (MQ=255)
                    tAACAGCACTTCTTTGTTGGCCGTTGCGCGGGTTGGGATGGTGCTGTCGCCTTCCAGCCATTCGACGCCCACCGCGCGTTTGCCGTCaaa                                                                  >  1:87482/1‑90 (MQ=255)
                        aGCACTTCTTTGTTGGCCGTTGCGCGGGTTGGGATGGTGCTGTCGCCTTCCAGCCATTCGACGCCCACCGCGCGTTTGCCGTCAAAAatg                                                              <  1:334747/90‑1 (MQ=255)
                          ttatcagcaggaaggttttggtccgatggatcgcaccgtcacgccgcagggccggcgcgccagGACCGCGCGTTTGCCGCCAAAAatgat                                                            <  1:211768‑M1/27‑1 (MQ=255)
                            cttcttTGTTGGCCGTTGCGCGGGTTGGGATGGTGCTGTCGCCTTCCAGCCATTCGACGCCCACCGCGCGTTTGCCGTCAAAAATGAtgt                                                          >  1:86339/1‑90 (MQ=255)
                            cttcttTGTTGGCCGTTGCGCGGGTTGGGATGGTGCTGTCGCCTTCCAGCCATTCGACGCCCACCGCGCGTTTGCCGTCAAAAATGAtgt                                                          <  2:322630/90‑1 (MQ=255)
                              tcttTGTTGGCCGTTGCGCGGGTTGGGATGGTGCTGTCGCCTTCCAGCCATTCGACGCCCACCGCGCGTTTGCCGTCAAAAatgat                                                            >  2:334138/1‑86 (MQ=255)
                              tcttTGTTGGCCGTTGCGCGGGTTGGGATGGTGCTGTCGCCTTCCAGCCATTCGACGCCCACCGCGCGTTTGCCGTCAAAAatgat                                                            <  1:334138/86‑1 (MQ=255)
                                tttGTTGGCCGTTGCGCGGGTTGGGATGGTGCTGTCGCCTTCCAGCCATTCGACGCCCACCGCGCGTTTGCCGTCAAAAATGATGTGATc                                                      <  1:573296/90‑1 (MQ=255)
                                tttGTTGGCCGTTGCGCGGGTTGGGATGGTGCTGTCGCCTTCCAGCCATTCGACGCCCACCGCGCGTTTGCCGTCAAAAATGATGTGATc                                                      <  2:758430/90‑1 (MQ=255)
                                tttGTTGGCCGTTGCGCGGGTTGGGATGGTGCTGTCGCCTTCCAGCCATTCGACGCCCACCGCGCGTTTGCCGTCAAAAATGATGTGATc                                                      <  2:727702/90‑1 (MQ=255)
                                tttGTTGGCCGTTGCGCGGGTTGGGATGGTGCTGTCGCCTTCCAGCCATTCGACGCCCACCGCGCGTTTGCCGTCAAAAATGATGTGATc                                                      <  2:207384/90‑1 (MQ=255)
                                tttGTTGGCCGTTGCGCGGGTTGGGATGGTGCTGTCGCCTTCCAGCCATTCGACGCCCACCGCGCGTTTGCCGTCAAAAATGATGTGATc                                                      <  2:259372/90‑1 (MQ=255)
                                  tgttgGCCGGTGCGCGGGTTGGGATGGTGCTGTCGCCTTCCAGCCATTCGACGCCCACCGCGCGTTTGCCGTCAAAAATGATGTGATCgg                                                    <  1:13900/90‑1 (MQ=255)
                                      ggCCGTTGCGCGGGTTGGGATGGTGCTGTCGCCTTCCAGCCATTCGACGCCCACCGCGCGTTTGCCGTCAAAAATGATGTGATCGGTCAt                                                <  1:442791/90‑1 (MQ=255)
                                      ggCCGTTGCGCGGGTTGGGATGGTGCTGTCGCCTTCCAGCCATTCGACGCCCACCGCGCGTTTGCCGTCAAAAATGATGTGATCGGTCAt                                                >  1:442061/1‑90 (MQ=255)
                                      ggCCGTTGCGCGGGTTGGGATGGTGCTGTCGCCTTCCAGCCATTCGACGCCCACCGCGCGTTTGCCGTCAAAAATGATGTGATCGGTCAt                                                >  1:251493/1‑90 (MQ=255)
                                        ccGTTGCGCGGGTTGGGATGGTGCTGTCGCCTTCCAGCCATTCGACGCCCACCGCGCGTTTGCCGTCAAAAATGATGTGATCGGTc                                                  >  2:96926/1‑86 (MQ=255)
                                        ccGTTGCGCGGGTTGGGATGGTGCTGTCGCCTTCCAGCCATTCGACGCCCACCGCGCGTTTGCCGTCAAAAATGATGTGATCGGTc                                                  <  1:96926/86‑1 (MQ=255)
                                           ttGCGCGGGTTGGGATGGTGCTGTCGCCTTCCAGCCATTCGACGCCCACCGCGCGTTTGCCGTCAAAAATGATGTGATCGGTCATAGCGt                                           >  2:100010/1‑90 (MQ=255)
                                             gcgcGGGTTGGGATGGTGCTGTCGCCTTCCAGCCATTCGACGCCCACCGCGCGTTTGCCGTCAAAAATGATGTGATCGGTCATAGCGTGa                                         <  1:634764/90‑1 (MQ=255)
                                                    ttGGGATGGTGCTGTCGCCTTCCAGCCATTCGACGCCCACCGCGCGTTTGCCGTCAAAAATGATGTGATCGGTCATAGCGTGAGTACGaa                                  <  1:3635/90‑1 (MQ=255)
                                                         aTGGTGCTGTCGCCTTCCAGCCATTCGACGCCCACCGCGCGTTTGCCGTCAAAAATGATGTGATCGGTCATAGCGTGAGTACGAATGGTc                             <  2:86339/90‑1 (MQ=255)
                                                            gTGCTGTCGCCTTCCAGCCATTCGACGCCCACCGCGCGTTTGCCGTCAAAAATGATGTGATCGGTc                                                  <  2:452621/66‑1 (MQ=255)
                                                            gTGCTGTCGCCTTCCAGCCATTCGACGCCCACCGCGCGTTTGCCGTCAAAAATGATGTGATCGGTc                                                  >  1:452621/1‑66 (MQ=255)
                                                            gTGCTGTCGCCTTCCAGCCATTCGACGCCCACCGCGCGTTTGCCGTCAAAAATGATGTGATCGGTCATAGCGTGAGTACGAATGGTCAga                          >  2:324217/1‑89 (MQ=255)
                                                              gCTGTCGCCTTCCAGCCATTCGACGCCCACCGCGCGTTTGCCGTCAAAAATGATGTGATCGGTCATAGCGTGAGTACGAATGGTCAGGtt                        <  2:442061/90‑1 (MQ=255)
                                                                   cgccgcagggccgtcgcgccagCACCGCGCGTTTGCCGTCAAAAATGATGTGATCGGTCATAGCGTGAGTACGAATGGTCAGGTTAGGAc                   <  2:675060‑M1/68‑1 (MQ=255)
                                                                      cTTCCAGCCATTCGACGCCCACCGCGCGTTTGCCGTCAAAAATGATGTGATCGGTCATAGCGTGAGTACGAATGATCAGGTTAGGAcgcg                >  2:742557/1‑90 (MQ=255)
                                                                           aGCCATTCGACGCCCACCGCGCGTTTGCCGTCAAAAATGATGTGATCGGTCATAGCGTGAGTACGAATGGTCAGGTTAGGACGCGATTTg           >  2:471557/1‑90 (MQ=255)
                                                                                ttCGACGCCCACCGCGCGTTTGCCGTCAAAAATGATGTGATCGGTCATAGCGTGAGTACGAATGGTCAGGTTAGGACGCGATTTGGCCTg      >  1:333884/1‑90 (MQ=255)
                                                                                   gACGCCCACCGCGCGTTTGCCGTCAAAAATGATGTGATCGGTCATAGCGTGAGTACGAATGGTCAGGTTAGGACGCGATTTGGCCTGATc   <  1:606741/90‑1 (MQ=255)
                                                                                    gccagCACCGCGCGTTTGCCGTCAAAAATGATGTGATCGGTCATAGCGTGAGTACGAATGGTCAGGTTAGGACGCGATTTGGCCTGATCg  <  2:415431‑M1/85‑1 (MQ=255)
                                                                                         |                                                                                    
GAGGCAATCGCGCCTGCACATAACAGCACTTCTTTGTTGGCCGTTGCGCGGGTTGGGATGGTGCTGTCGCCTTCCAGCCATTCGACGCCCACCGCGCGTTTGCCGTCAAAAATGATGTGATCGGTCATAGCGTGAGTACGAATGGTCAGGTTAGGACGCGATTTGGCCTGATCG  >  NZ_CP009273/322171‑322344

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