New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? REL606 = 288616744 (0.890)6 (0.120) 4/56 NT 15.2% coding (68/2259 nt) xdhA xanthine dehydrogenase, molybdenum binding subunit
?REL606 3650766 = 23 (0.370)intergenic (+47/+140) cspA/hokA major cold shock protein/toxic polypeptide, small

AAGCACAGAATCTAAGATCCCTGCCATTTGGCGGGGATTTTTTTATTTGTTTTCAGGAAATAAATAATCGA  >  REL606/3650731‑3650801
                                   |                                   
cgggcacgatatactgacgattatgttatggcgggg                                     >  1:7560246‑M2/36‑36 (MQ=255)
aaGCACAGAATCTAAGATCCCTGCCATTTGGCgggg                                     >  1:2611542/1‑36 (MQ=255)
 gggcacgatatactgacgattatgttatggcggggg                                    >  1:448471‑M2/35‑35 (MQ=255)
 gggcacgatatactgacgattatgttatggcggggc                                    >  1:106900‑M2/35‑35 (MQ=255)
 atcacaGAATCTAAGATCCCTGCCATTTGGCGGGGa                                    <  1:4103835/34‑1 (MQ=255)
  gCACAGAATCTAAGATCCCTGCCATTTGGCGGGGAt                                   <  1:3541639/36‑1 (MQ=255)
    acaGAATCTAAGATCCCTGCCATTTGGCGGGGGttt                                 >  1:3038369/1‑36 (MQ=255)
    acaGAATCTAAGATCCCTGCCATTTGGCGGGGAttt                                 <  1:6646464/36‑1 (MQ=255)
    acaGAATCTAAGATCCCTGCCATTTGGCGGGGAttt                                 >  1:249523/1‑36 (MQ=255)
     caGAATCTAAGATCCCTGCCATTTGGCGGGGAtttt                                <  1:97748/36‑1 (MQ=255)
      cgatatactgacgattatgttatggcgggGGttttt                               >  1:6708099‑M2/30‑36 (MQ=255)
      aGAATCTAAGATCCCTGCCATTTGGCGGGGGttttt                               >  1:6118885/1‑36 (MQ=255)
      aGAATCTAAGATCCCTGCCATTTGGCGGGGAttttt                               <  1:1748063/36‑1 (MQ=255)
        atatactgacgattatgttatggcgggGAgtttttt                             >  1:6641227‑M2/28‑36 (MQ=255)
        aaTCTAAGATCCCTGCCATTTGGCGGGGAttttttt                             >  1:7571426/1‑36 (MQ=255)
        aaTCTAAGATCCCTGCCATTTGGCGGGGAttttttt                             <  1:579198/36‑1 (MQ=255)
        aaTCTAAGATCCCTGCCATTTGGCGGGGAttttttt                             <  1:1492897/36‑1 (MQ=255)
         aTCTAAGATCCCTGCCATTTGGCGGGGATTTTTTTa                            >  1:4845272/1‑36 (MQ=255)
          tCTAAGATCCCTGCCATTTGGCGGGGATTTTTTTAt                           >  1:62422/1‑36 (MQ=255)
           tactgacgattatgttatggcgggGATTTTTTTTtt                          >  1:2043834‑M2/25‑36 (MQ=255)
             aaGATCCCTGCCATTTGGCGGGGATTTTTTTAtttt                        >  1:577010/1‑35 (MQ=255)
             aaGATCCCTGCCATTTGGCGGGGATTTTTTTATTTg                        >  1:2428079/1‑36 (MQ=255)
               gATCCCTGCCATTTGGCGGGGATTTTTTTATTTTtt                      >  1:1340919/1‑36 (MQ=255)
               gATCCCTGCCATTTGGCGGGGATTTTTTTATTTGtt                      >  1:5133963/1‑36 (MQ=255)
               gATCCCTGCCATTTGGCGGGGATTTTTTTATTTGtt                      <  1:6324941/36‑1 (MQ=255)
                aTCCCTGCCATTTGGCGGGGATTTTTTTATTTGttt                     <  1:5009736/36‑1 (MQ=255)
                aTCCCTGCCATTGGGCGGGGATTTTTTTATTTGttt                     >  1:7234751/1‑36 (MQ=255)
                 acccTGCGATTTGGCGGGGATTTTTTTATTTGtttt                    <  1:3456975/35‑1 (MQ=255)
                   ccTGCCATTTGGCGGGGATTTTTTTATTTGTTTTCa                  >  1:887126/1‑36 (MQ=255)
                   ccTGCCATTTGGCGGGGATTTTTTTATTTGTTTTCa                  <  1:1375322/36‑1 (MQ=255)
                    cTGCCATTTGGCGGGGATTTTTTTATTTGTTTTCAg                 <  1:5328677/36‑1 (MQ=255)
                     tGCCATTTGGCGGGGATTTTTTTATTTGTTTTCAgg                >  1:7631654/1‑36 (MQ=255)
                      gCCATTTGGCGGGGATTTTTTTATTTGTTTTCatta               >  1:2541486/1‑33 (MQ=255)
                         aTTTGGCGGGGATTTTTTTATTTGTTTTCAGGaaat            <  1:1156393/36‑1 (MQ=255)
                         aTTTGGCGGGGATTTTTTTATTTGTTTTCAGGaaat            >  1:669980/1‑36 (MQ=255)
                          tttGGCGGGGATTTTTTTATTTGTTTTCAGGaaata           <  1:707006/36‑1 (MQ=255)
                           ttGGCGGGGATTTTTTTATTTGTTTTCAGGaaataa          <  1:2742746/36‑1 (MQ=255)
                           ttGGCGGGGATTTTTTTATTTGTTTTCAGGaaataa          <  1:2029750/36‑1 (MQ=255)
                            tggcgggGATTTTTTTATTTGTTTTCAGGaaataaa         <  1:3080534‑M2/29‑1 (MQ=255)
                             ggcgggGATTTTTTTATTTGTTTTCAGGaaataaat        <  1:7245473‑M2/30‑1 (MQ=255)
                               cgggGATTTTTTTATTTGTTTTCAGGAAataaataa      <  1:2649038‑M2/32‑1 (MQ=255)
                               cgggGATTTTTTTATTTGTTTTCAGGAAataaataa      <  1:7364016‑M2/32‑1 (MQ=255)
                                gggGATTTTTTTATTTGTTTTCAGGAAATAaataat     >  1:1682586‑M2/4‑36 (MQ=255)
                                gggGATTTTTTTATTTGTTTTCAGGAAATAaataat     <  1:7822298‑M2/33‑1 (MQ=255)
                                 ggGATTTTTTTATTTGTTTTCAGGAAATAAATAATc    <  1:4691756‑M2/34‑1 (MQ=255)
                                   gATTTTTTTATTTGTTTTCAGTAAATAAATAatcga  >  1:8170956/1‑36 (MQ=255)
                                   gATTTTTTTATTTGTTTTCAGGAAATAAATAatcga  >  1:4824947/1‑36 (MQ=255)
                                   |                                   
AAGCACAGAATCTAAGATCCCTGCCATTTGGCGGGGATTTTTTTATTTGTTTTCAGGAAATAAATAATCGA  >  REL606/3650731‑3650801

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