New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? REL606 1946059 =30 (0.580)5 (0.100) 4/60 NT 12.7% coding (149/390 nt) cheY chemotaxis regulator transmitting signal to flagellar motor component
?REL606 = 3240660 39 (0.650)coding (392/1245 nt) mtr tryptophan transporter of high affinity

GGGCATGTTCCAGTCGGAGATAACAAATCCATAACCGCCTGCCTGCAACTTATTGAGAGCGTCGACGCC  >  REL606/1946025‑1946093
                                  |                                  
gggCATGTTCCAGTCGGAGATAACAAATCCATAAcc                                   <  1:2079555/36‑1 (MQ=255)
gggCATGTTCCAGTCGGAGATAACAAATCCATAAcc                                   <  1:4483589/36‑1 (MQ=255)
accacaaacgctaccagcagtgcaaaaccaaaaccc                                   <  1:83551‑M1/2‑1 (MQ=255)
accacaaacgctaccagcagtgcaaaaccaaaaccc                                   <  1:5861371‑M1/2‑1 (MQ=255)
 ggCATGTTCCAGTCGGAGATAACAAATCCATAACCg                                  <  1:6343455/36‑1 (MQ=255)
 ccacaaacgctaccagcagtgcaaaaccaaaacCCg                                  <  1:4245601‑M1/3‑1 (MQ=255)
  gaaTGTTCCAGTCGGAGATAACAAATCCATAACCgc                                 <  1:7139170/34‑1 (MQ=255)
  gCATGTTCCAGTCGTATATAACAAATCCATAACCgc                                 <  1:765438/36‑1 (MQ=37)
  gCATGTTCCAGTCGGAGATAACAAATCCATAACCgc                                 <  1:3178493/36‑1 (MQ=255)
  cacaaacgctaccagcagtgcaaaaccaaaacCCgc                                 >  1:7663078‑M1/33‑36 (MQ=255)
  cacaaacgctaccagcagtgcaaaaccaaaacCCgc                                 <  1:729355‑M1/4‑1 (MQ=255)
  cacaaacgctaccagcagtgcaaaaccaaaacCCgc                                 <  1:4935313‑M1/4‑1 (MQ=255)
   cATGTTCCAGTCGGAGATAACAAATCCATAACCgcc                                >  1:3685236/1‑36 (MQ=255)
   acaaacgctaccagcagtgcaaaaccaaaacCCgcc                                <  1:3245079‑M1/5‑1 (MQ=255)
   acaaacgctaccagcagtgcaaaaccaaaacCCgcc                                >  1:6227802‑M1/32‑36 (MQ=255)
   acaaacgctaccagcagtgcaaaaccaaaacCCgcc                                >  1:7837913‑M1/32‑36 (MQ=255)
     tGTTCAAGTCGGAGATAACAAATCCATAACCgcctg                              <  1:4174170/36‑1 (MQ=255)
      gTTCCAGTCGGAGATAACAAATCCATAACCgcctgc                             >  1:7759210/1‑36 (MQ=255)
       ttCCAGTCGGAGATAACAAATCCATAACCgcctgcc                            >  1:128015/1‑36 (MQ=255)
        tCCAGTCGGAGATAACAAATCCATAACCgcctgcct                           >  1:399802/1‑36 (MQ=255)
          cAGTCGGAGATAACAAATCCATAACCGCctgcctgc                         >  1:3778232/1‑36 (MQ=255)
          cAGTCGGAGATAACAAATCCATAACCGCctgcctgc                         >  1:2507201/1‑36 (MQ=255)
           aGTCGGAGATAACAAATCCATAACCGCCTGCCTGCa                        >  1:4280661/1‑36 (MQ=255)
           aGTCGGAGATAACAAATCCATAACCGCCTGCCTGCa                        >  1:5704521/1‑36 (MQ=255)
            gTCGGAGATAACAAATCCATAACCGCCTGCCTGCaa                       >  1:7452550/1‑36 (MQ=255)
            gTCGGAGATAACAAATCCATAACCGCCTGCCTGCaa                       >  1:1618146/1‑36 (MQ=255)
            gTCGGAGATAACAAATCCATAACCGCCTGCCTGCaa                       >  1:4970479/1‑36 (MQ=255)
               ggAGATAACAAATCCATAACCGCCTGCCTGCAACtt                    >  1:6565430/1‑36 (MQ=255)
               ggAGATAACAAATCCATAACCGCCTGCCTGCAACtt                    >  1:1701088/1‑36 (MQ=255)
               ggAGATAACAAATCCATAACCGCCTGCCTGCAACtt                    >  1:609874/1‑36 (MQ=255)
               ggAAATAAAAAATCAATAACCGCCTGCCTGCAACtt                    <  1:2621946/36‑1 (MQ=21)
                gagaTAACAAATCCATAACCGCCTGCCTGCAACTTa                   >  1:751320/1‑36 (MQ=255)
                 agaTAACAAATCCATAACCGCCTGCCTGCAACTTAt                  >  1:7417598/1‑36 (MQ=255)
                 agaTAACAAATCCATAACCGCCTGCCTGCAACTTAt                  >  1:6125988/1‑36 (MQ=255)
                  gaTAACAAATCCATAACCGCCTGCCTGCAACTTAtt                 >  1:2864631/1‑36 (MQ=255)
                    tAACAAATCCATAACCGCCTGCCTGCAACTTATTga               >  1:6989504/1‑36 (MQ=255)
                     aaCAAATCCATAACCGCCTGCCTGCAACTTATTgag              >  1:7494753/1‑36 (MQ=255)
                      aCAAATCCATAACAGCCTGCCTGCTTATTATTgagc             >  1:153377/1‑35 (MQ=16)
                       aaaaccaaaaaCCGCCTGCCTGCAACTTATTGagag            >  1:3930422‑M1/12‑36 (MQ=255)
                       aaaaccaaaaaCCGCCTGCCTGCAACTTATTGagag            >  1:938316‑M1/12‑36 (MQ=255)
                        aaaccaaaacCCGCCTGCCTGCAACTTATTGAGAGc           >  1:2937989‑M1/11‑36 (MQ=255)
                        aaaTCCATAACCGCCTGCCTGCAACTTATTGatagt           >  1:3872329/1‑35 (MQ=255)
                        aaaTCCATAACCGCCTGCCTGCAACTTATTGAGAGc           >  1:2611308/1‑36 (MQ=255)
                         aaTCCATAACCGCCTGCCTGCAACTTATTGAGAGCg          >  1:5466923/1‑36 (MQ=255)
                          accaaaacCCGCCTGCCTGCAACTTATTGAGAGCGt         >  1:3165523‑M1/9‑36 (MQ=255)
                              aaacCCGCCTGCCTGCAACTTATTGAGAGCGTCGAc     >  1:989506‑M1/5‑36 (MQ=255)
                              aTAACCGCCTGCCTGCAACTTATTGAGAGCGTCGAc     >  1:2838383/1‑36 (MQ=255)
                              aTAACCGCCTGCCTGCAACTTATTGAGAGCGTCGAc     >  1:4784309/1‑36 (MQ=255)
                                aaCCGCCTGCCTGCAACTTATTGAGAGCGTCGACGc   >  1:2473148/1‑36 (MQ=255)
                                 aCCGCCTGCCTGCAACTTATTGAGAGCGTCGACGcc  >  1:1511778/1‑36 (MQ=255)
                                 aCCGCCTGCCTGCAACTTATTGAGAGCGTCGACGcc  <  1:401069/36‑1 (MQ=255)
                                  |                                  
GGGCATGTTCCAGTCGGAGATAACAAATCCATAACCGCCTGCCTGCAACTTATTGAGAGCGTCGACGCC  >  REL606/1946025‑1946093

Base quality scores:  ATCG  < 3 ≤  ATCG  < 11 ≤  ATCG  < 22 ≤  ATCG  < 32 ≤  ATCG  < 37 ≤  ATCG