New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? AP009048 39915 =49 (0.980)4 (0.090) 3/92 NT 8.3% coding (472/1143 nt) caiA crotonobetaine reductase subunit II, FAD‑binding
?AP009048 39938 = 44 (0.980)coding (449/1143 nt) caiA crotonobetaine reductase subunit II, FAD‑binding

AGCCCCGTCGCGCGCCATCACCACGATGTACGGGGTGTAGGCGCTGCTGGTAATAAAACACTTACTACCATTAAGATAAATCTTACCATTTCTACGGGTATAAGTCGTTTTCAGGCTACCCACGTCGGAGCCCGCGCCC  >  AP009048/39868‑40006
                                                                      |                                                                    
aGCCCCGTCGCGCGCCATCACCACGATGTACGGGGTGTAGGCGCTGCTGGTAATAAAACACTTACTACCAt                                                                      <  1:2831217/71‑1 (MQ=255)
  cctgaaaacgacttatacccgtagaaatggtaagatttatcttaaTGGTAATAAAACACTTACTACCAt                                                                      >  1:115925‑M1/46‑69 (MQ=255)
    ccGTCGCGCGCCATCACCACGATGTACGGGGTGTAGGCGCTGCTGGTAATAAAACACTTACTACCAt                                                                      >  1:4336027/1‑67 (MQ=255)
     cGTCGCGCGCCATCACCACGATGTACGGGGTGTAGGCGCTGCTGGTAATAAAACACTTACTACCAt                                                                      >  1:2052212/1‑66 (MQ=255)
             gcCATCACCACGATGTACGGGGTGTAGGCGCTGCTGGTAATAAAACACTTACTACCATTAAGATAAATCt                                                          >  1:3634928/1‑70 (MQ=255)
              ccATCACCACGATGTACGGGGTGTAGGCGCTGCTGGTAATAAAACACTTACTACCATTAAGATAAATCTTa                                                        <  1:1314130/71‑1 (MQ=255)
              ccATCACCACGATGTACGGGGTGTAGGCGCTGCTGGTAATAAAACACTTACTACCATTAAGATAAATCTTa                                                        <  1:3018168/71‑1 (MQ=255)
               cATCACCACGATGTACGGGGTGTAGGCGCTGCTGGTAATAAAACACTTACTACCATTaa                                                                   >  1:608417/1‑59 (MQ=255)
               cATCACCACGATGTACGGGGTGTAGGCGCTGCTGGTAATAAAACACTTACTACCATTAAGATAAATCTTAc                                                       >  1:304979/1‑71 (MQ=255)
                aTCACCACGATGTACGGGGTGTAGGCGCTGCTGGTAATAAAACACTTACTACCATTAAGATAAATCTTAc                                                       >  1:4612820/1‑70 (MQ=255)
                aTCACCACGATGTACGGGGTGTAGGCGCTGCTGGTAATAAAACACTTACTACCATTAAGATAAATCTTAc                                                       >  1:4339495/1‑70 (MQ=255)
                 tCACCACGATGTACGGGGTGTAGGCGCTGCTGGTAATAAAACACTTACTACCATTAAGATAAATCTTAc                                                       >  1:2234883/1‑69 (MQ=255)
                 tCACCACGATGTACGGGGTGTAGGCGCTGCTGGTAATAAAACACTTACTACCATTAAGATAAATCTTACCa                                                     >  1:1453891/1‑71 (MQ=255)
                 tCACCACGATGTACGGGGTGTAGGCGCTGCTGGTAATAAAACACTTACTACCATTAAGATAAATCTTACCa                                                     >  1:4105214/1‑71 (MQ=255)
                 tCACCACGATGTACGGGGTGTAGGCGCTGCTGGTAATAAAACACTTACTACCATTAAGATAAATCTTACCa                                                     >  1:3589858/1‑71 (MQ=255)
                   accacGATGTACGGGGTGTAGGCGCTGCTGGTAATAAAACACTTACTACCATTAAGATAAATCTTACCAtt                                                   >  1:4481212/1‑71 (MQ=255)
                    ccacGATGTACGGGGTGTAGGCGCTGCTGGTAATAAAACACTTACTACCATTAAGATaa                                                              <  1:3672003/59‑1 (MQ=255)
                     cacGATGTACGGGGTGTAGGCGCTGCTGGTAATAAAACACTTACTACCATTAAGATAAATCTTACCATTTc                                                 >  1:4312013/1‑71 (MQ=255)
                      acGATGTACGGGGTGTAGGCGCTGCTGGTAATAAAACACTTACTACCATTAAGATAAATCTTACCATTTc                                                 >  1:2053991/1‑70 (MQ=255)
                       cGATGTACGGGGTGTAGGCGCTGCTGGTAATAAAACACTTACTACCATTAAGATAAATCTTACCAt                                                    >  1:3415085/1‑66 (MQ=255)
                       cGATGTACGGGGTGTAGGCGCTGCTGGTAATAAAACACTTACTACCATTAAGATAAATCTTACCATTTCTa                                               >  1:96179/1‑71 (MQ=255)
                       cGATGTACGGGGTGTAGGCGCTGCTGGTAATAAAACACTTACTACCATTAAGATAAATCTTACCATTTCTa                                               >  1:3384904/1‑71 (MQ=255)
                            tACGGGGTGTAGGCGCTGCTGGTAATAAAACACTTACTACCATTAAGATAAATCTTAc                                                       >  1:3534057/1‑58 (MQ=255)
                            tACGGGGTGTAGGCGCTGCTGGTAATAAAACACTTACTACCATTAAGATAAATCTTACCATTTCTACGGGt                                          <  1:3194909/71‑1 (MQ=255)
                              cGGGGTGTAGGCGCTGCTGGTAATAAAACACTTACTACCATTAAGATAAATCTTACCATTTCTACGGGtat                                        <  1:4386230/71‑1 (MQ=255)
                              cGGGGTGTAGGCGCTGCTGGTAATAAAACACTTACTACCATTAAGATAAATCTTACCATTTCTACGGGtat                                        <  1:581549/71‑1 (MQ=255)
                              cGGGGTGTAGGCGCTGCTGGTAATAAAACACTTACTACCATTAAGATAAATCTTACCATTTCTACGGGtat                                        <  1:2238544/71‑1 (MQ=255)
                                 ggTGTAGGCGCTGCTGGTAATAAAACACTTACTACCATTAAGATAAATCTTACCATTTCTACGGGTATAAg                                     <  1:3953924/71‑1 (MQ=255)
                                  gtgtAGGCGCTGCTGGTAATAAAACACTTACTACCATTAAGATAAATCTTACCATTTCTACGGGTATaa                                      <  1:554921/69‑1 (MQ=255)
                                     tAGGCGCTGCTGGTAATAAAACACTTAATACCATTAAGATAAATCTTACCATTTCTACGGGTATAAGTCGt                                 <  1:2532777/71‑1 (MQ=255)
                                       ggCGCTGCTGGTAATAAAACACTTACTACCATTAAGATAAATCTTACCATTTCTACGGGTATAAGTCGtt                                <  1:4317787/70‑1 (MQ=255)
                                        gcgcTGCTGGTGATAAAACACTTACTACCATTAAGATAAATCTTACCATTTCTACGGGTATAAGTCGtt                                <  1:2796667/69‑1 (MQ=255)
                                        gcgcTGCTGGTAATAAAACACTTACTACCATTAAGATAAATCTTACCATTTCTACGGGTATAAGTCGtttt                              >  1:2320017/1‑71 (MQ=255)
                                        gcgcTGCTGGTAATAAAACACTTACTACCATTAAGATAAATCTTACCATTTCTACGGGTATAAGTCGtttt                              <  1:183275/71‑1 (MQ=255)
                                        gcgcTGCTGGTAATAAAACACTTACTACCATTAAGATAAATCTTACCATTTCTACGGGTATAAGTCGtttt                              <  1:825165/71‑1 (MQ=255)
                                        gcgcTGCTGGTAATAAAACACTTACTACCATTAAGATAAATCTTACCATTTCTACGGGTATAAGTCGtttt                              <  1:3726406/71‑1 (MQ=255)
                                        gcgcTGCTGGTAATAAAACACTTACTACCATTAAGATAAATCTTACCATTTCTACGGGTATAAGTCGtttt                              <  1:4348925/71‑1 (MQ=255)
                                             gctgGTAATAAAACACTTACTACCATTAAGATAAATCTTACCATTTc                                                 >  1:1402226/1‑47 (MQ=255)
                                             gctgGTAATAAAACACTTACTACCATTAAGATAAATCTTACCATTTc                                                 >  1:2958107/1‑47 (MQ=255)
                                             gctgGTAATAAAACACTTACTACCATTAAGATAAATCTTACCATTTCTACGGGTATAAGTCGTTTTCAGGc                         >  1:1632159/1‑71 (MQ=255)
                                              ctgGTAATAAAACACTTACTACCATTAAGATAAATCTTACCAttt                                                  >  1:527988/1‑45 (MQ=255)
                                              atggtagtaagtgttttattaccaTTAAGATAAATCTTACCATTTCTACGGGTATAAGTCGTTTTCAgg                          >  1:115925‑M2/25‑69 (MQ=255)
                                               tggtagtaagtgttttattaccaTTAAGATAAATCTTACCATTTCTACGGGTATaa                                      >  1:1129623‑M2/24‑56 (MQ=255)
                                               tggtagtaagtgttttattaccaTTAAGATAAATCTTACCATTTCTACGGGTATaa                                      >  1:3829255‑M2/24‑56 (MQ=255)
                                               TGGTAATAAAACACTTACTACCATTAAGATAAATc                                                           >  1:3986228‑M1/1‑35 (MQ=255)
                                               TGGTAATAAAACACTTACTACCATTAAGATAAATCTTACCATTTCTACGGGTATAAGTCGTTTTc                             <  1:93343‑M1/65‑1 (MQ=255)
                                                 GTAATAAAACACTTACTACCATTAAGATAAATCTTACCATTTCTACGGGTATAAGTCGTTTTCAGGCTAcc                     <  1:2366161‑M1/71‑1 (MQ=255)
                                                 GTAATAAAACACTTACTACCATTAAGATAAATCTTACCATTTCTACGGGTATAAGTCGTTTTCAGGCTAcc                     <  1:1960043‑M1/71‑1 (MQ=255)
                                                 GTAATAAAACACTTACTACCATTAAGATAAATCTTACCATTTCTACGGGTATAAGTCGTTTTCAGGCTAc                      >  1:145801‑M1/1‑70 (MQ=255)
                                                  TAATAAAACACTTACTACCATTAAGATAAATCTTACCATTTCTACGGGTATAAGTCGTTTTCAGGCTAcc                     <  1:3326324‑M1/70‑1 (MQ=255)
                                                   AATAAAACACTTACTACCATTAAGATAAATCTTACCAtt                                                   <  1:973587‑M1/39‑1 (MQ=255)
                                                   AATAAAACACTTACTACCATTAAGATAAATCTTACCATTTc                                                 >  1:4395858‑M1/1‑41 (MQ=255)
                                                   AATAAAACACTTACTACCATTAAGATAAATCTTACCATTTCTACGGGTATAAGTCGTTTTCAGGCTACCCa                   <  1:816800‑M1/71‑1 (MQ=255)
                                                    gtaagtgttttattaccaTTAAGATAAATCTTACCATTTCTACGGGTATAAGTCGtttt                              <  1:4001500‑M2/41‑1 (MQ=255)
                                                      aaaaCACTTACTACCATTAAGATAAATCTTACCATTTCTACgg                                            >  1:3559103/1‑43 (MQ=255)
                                                       aaaCACTTACTACCATTAAGATAAATCTTACCATTTCTACgg                                            >  1:2133334/1‑42 (MQ=255)
                                                       aaaCACTTACTACCATTAAGATAAATCTTACCATTTCTACGGGTATAAGTCGTTTTCAGGCTAcc                     >  1:2728974/1‑65 (MQ=255)
                                                       aaaCACTTACTACCATTAAGATAAATCTTACCATTTCTACGGGTATAAGTCGTTTTCAGGCTACCCACGTc               >  1:266500/1‑71 (MQ=255)
                                                        aaCACTTACTACCATTAAGATAAATCTTACCATTTCTACGGGTATAAGTCGTTTTCAGGCTACCCACGTCg              <  1:1989757/71‑1 (MQ=255)
                                                         acacTTACTACCATTAAGATAAATCTTACCATTTc                                                 <  1:2968222/35‑1 (MQ=255)
                                                         acacTTACTACCATTAAGATAAATCTTACCATTTCTACGGGTATAAGTCGTTTTCa                            >  1:315657/1‑56 (MQ=255)
                                                         acacTTACTACCATTAAGATAAATCTTACCATTTCTACGGGTATAAGTCGTTTTCAGGCTACCCACGTCg              >  1:1657648/1‑70 (MQ=255)
                                                            cTTACTACCATTAAGATAAATCTTACCATTTCTACGGGTATAAGTCGTTTTCAGGCTACCCACGTCGGAGc          <  1:4092156/71‑1 (MQ=255)
                                                             ttACTACCATTAAGATAAATCTTACCATTTCTACg                                             >  1:144962/1‑35 (MQ=255)
                                                             ttACTACCATTAAGATAAATCTTACCATTTCTACGGGTATAAGTCGTTTTCa                            >  1:4324853/1‑52 (MQ=255)
                                                             ttACTACCATTAAGATAAATCTTACCATTTCTACGGGTATAAGTCGTTTTCAGGCTACCCAc                  <  1:1750586/62‑1 (MQ=255)
                                                                 taccaTTAAGATAAACCTTACCATTTCTACGGGTATAAg                                     <  1:4418050‑M2/34‑1 (MQ=255)
                                                                   ccaTTAAGATAAATCTTACCATTTCTACGGGTATAAGTCGtttt                              <  1:1863055‑M2/41‑1 (MQ=255)
                                                                   ccaTTAAGATAAATCTTACCATTTCTACGGGTATAAGTCGttt                               >  1:4155411‑M2/4‑43 (MQ=255)
                                                                   ccaTTAAGATAAATCTTACCATTTCTACGGGTATAAGTCGTTTTCa                            >  1:4401924‑M2/4‑46 (MQ=255)
                                                                    caTTAAGATAAATCTTACCATTTCTACGGGTATAAGTCGTTTTCAGGCTACCCACGTCGGAGCCCGCGccc  <  1:1692508‑M2/69‑1 (MQ=255)
                                                                      ttAAGATAAATCTTACCATTTCTACGGGTATAAGTCGTTTTc                             <  1:3867736/42‑1 (MQ=255)
                                                                      ttAAGATAAATCTTACCATTTCTACGGGTATAAGTCGTTTTCa                            >  1:3141157/1‑43 (MQ=255)
                                                                      |                                                                    
AGCCCCGTCGCGCGCCATCACCACGATGTACGGGGTGTAGGCGCTGCTGGTAATAAAACACTTACTACCATTAAGATAAATCTTACCATTTCTACGGGTATAAGTCGTTTTCAGGCTACCCACGTCGGAGCCCGCGCCC  >  AP009048/39868‑40006

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