New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? REL606 632818 =34 (0.660)7 (0.140) 4/60 NT 18.9% coding (17/297 nt) citD citrate lyase, acyl carrier (gamma) subunit
?REL606 1241885 = 26 (0.430)coding (337/1737 nt) ycgO potassium/proton antiporter

TAATAAGCCTTCAATCAAATCAAGATTAAACAGCCACGCCGCCATCATGCCGGTTAAACCAGAGGTGAT  >  REL606/1241850‑1241918
                                   |                                 
taataaGCCTTAAATCAAATCAAGATTAAACAGCCa                                   <  1:405648/36‑1 (MQ=255)
 aataaGCCTTCAATCAAATCAAGATTAAACAGCCAc                                  >  1:6108061/1‑36 (MQ=255)
 aataaGCCTTCAATCAAATCAAGATTAAACAGCCAc                                  >  1:5217116/1‑36 (MQ=255)
  ataaGCCTCAAATAAAATCAAAATAAAACAGCCAcg                                 <  1:7708364/36‑1 (MQ=16)
   taaGCCTTCAATCAAATCAAGATTAAACAGCCAcgc                                >  1:5082394/1‑36 (MQ=255)
    aaGCCTTCAATCAAATCAAGATTAAACAGCCAcgcc                               >  1:2862/1‑36 (MQ=255)
     aGCCTTCAATCAAATCAAGATTAAACAGCCAcgccg                              >  1:1442875/1‑36 (MQ=255)
     aGCCTTAAAAAAAATAAAGATTAAACAGCCAcgccg                              <  1:6181097/36‑1 (MQ=19)
      gcattCAAACAAATAAAGATAAAACAGCCAcgccgc                             <  1:5072613/33‑1 (MQ=16)
        cTTCAAACAAATCAAGATAAAACAGCCACGCCGCca                           <  1:799964/36‑1 (MQ=37)
           caatcaaATCAAGATTAAACAGCCACGCCGCcatca                        >  1:6261181/1‑36 (MQ=255)
            aatcaaatcaAGATTAAACAGCCACGCCGCtatcat                       >  1:6779148/1‑36 (MQ=255)
            aatcaaatcaAGATTAAACAGCCACGCCGCcatcat                       >  1:5123563/1‑36 (MQ=255)
            aatcaaatcaAGATTAAACAGCCACGCCGCcatcat                       >  1:4458825/1‑36 (MQ=255)
            aatcaaatcaAGATTAAACAGCCACGCCGCcatcat                       >  1:163031/1‑36 (MQ=255)
            aatcaaagaaAGACTAAACAGCCACGCCGCcatcat                       <  1:2935209/36‑1 (MQ=25)
              aaaaataaaAATAAAACAGCCCCGCCGCCATCATGc                     <  1:1156205/34‑1 (MQ=12)
               caaatcaaGATTAAACAGCCACGCCGCCATCATGcc                    >  1:2129594/1‑36 (MQ=255)
                aaataaaAATAAAACAGCCACGCCGCCATCATGCCg                   <  1:2037549/36‑1 (MQ=21)
                  aaaaaGATTAAACAGCCACGCCGCCATCATGCCGGt                 <  1:4280612/33‑1 (MQ=255)
                   taaaGATAAACCAGCCACGCCGCCATCATGCCGGtt                <  1:3158675/34‑1 (MQ=25)
                   acaaGATTAAACAGCCACGCCGCCATCATGCCGGtt                <  1:5279818/35‑1 (MQ=255)
                   acaaGATAAAACAGCCACGCCGCCATCATGCCGGtt                <  1:2467572/35‑1 (MQ=255)
                    caaGATTAAACAGCCACGCCGCCATCATGCCGGTTa               >  1:1254592/1‑36 (MQ=255)
                     aaaataaaccagccACGCCGCCATCATGCCGGTTaa              <  1:3733130‑M2/22‑1 (MQ=255)
                     aaaataaaacagccACGCCGCCATCATGCCGGTTaa              <  1:5831416‑M2/22‑1 (MQ=255)
                     aaGATTAAACAGCCACCCCGCCATCATGCCGGTTaa              >  1:744536/1‑36 (MQ=255)
                      aaataaaccagccACGCCGCCATCATGCCGGTTaaa             <  1:5173318‑M2/23‑1 (MQ=255)
                      aGATTAAACAGCCACGCCGCCATCATGCCGGTTaaa             >  1:5293257/1‑36 (MQ=255)
                       aataaaacagccACGCCGCCATCATGCCGGTTAAAc            <  1:6889132‑M2/24‑1 (MQ=255)
                       aataaaacagccACGCCGCCATCATGCCGGTTAAAc            <  1:3493249‑M2/24‑1 (MQ=255)
                        aataaccagccACGCCGCCATCATGCCGGTTAAAcc           <  1:1563016‑M2/25‑1 (MQ=255)
                        aatAAACAGCCACGCCGCCATCATGCCGGTTAAAcc           <  1:24120/34‑1 (MQ=255)
                         ttAAACAGCCACGCCGCCATCATGCCGGTTAAACCa          >  1:4160340/1‑36 (MQ=255)
                          tAAACAGCCACGCCGCCATCATGCCGGTTAAACAag         >  1:280160/1‑36 (MQ=255)
                          aaaaCAGCCACGCCGCCATCATGCCGGTTAAACCag         <  1:3770772/35‑1 (MQ=255)
                            accagccACGCCGCCATCATGCCGGTTAAACCagag       <  1:1264106‑M2/29‑1 (MQ=255)
                            aaCAGCCACGCCGCCATCATGCCGGTTAAACCatag       >  1:6871150/1‑36 (MQ=255)
                            aaCAGCCACGCCGCCATCATGCCGGTTAAACCagag       >  1:7030429/1‑36 (MQ=255)
                              cagccACGCCGCCATCATGCCGGTTAAACCAGAGGt     <  1:5608804‑M2/31‑1 (MQ=255)
                               agccACGCCGCCATCATGCCGGTTAAACCAGAGGTg    >  1:5781260‑M2/5‑36 (MQ=255)
                               agccACGCCGCCATCATGCCGGTTAAACCAGAGGTg    >  1:4447655‑M2/5‑36 (MQ=255)
                                 ccACGCCGCCATCATGCCGGTTAAACCAGAGGTGAt  >  1:212481‑M2/3‑36 (MQ=255)
                                 ccACGCCGCCATCATGCCGGTTAAACCAGAGGTGAt  >  1:4764480‑M2/3‑36 (MQ=255)
                                   |                                 
TAATAAGCCTTCAATCAAATCAAGATTAAACAGCCACGCCGCCATCATGCCGGTTAAACCAGAGGTGAT  >  REL606/1241850‑1241918

Base quality scores:  ATCG  < 3 ≤  ATCG  < 6 ≤  ATCG  < 19 ≤  ATCG  < 30 ≤  ATCG  < 37 ≤  ATCG