New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? REL606 3198500 =39 (0.710)10 (0.180) 5/62 NT 27.0% coding (563/990 nt) tdcB threonine dehydratase
?REL606 3424874 = 15 (0.240)intergenic (+92/‑34) nirD/nirC nitrite reductase small subunit/nitrite transporter

TAATAGATTTAATTGCCACCGCAATACCAGCAATTAAACCGCCACCACCAATTGGCACAATCACGTTATCG  >  REL606/3198465‑3198535
                                   |                                   
tAATAGATTTAATTGCCACCGCAATACCAGCAATTa                                     >  1:1899531/1‑36 (MQ=255)
tAATAGATTTAATTGCCACCGCAATACCAGCAATTa                                     >  1:3763404/1‑36 (MQ=255)
tAATAGATTTAATTGCCACCGCAATACCAGAAATTa                                     >  1:7315349/1‑36 (MQ=255)
 aaTAGATTTAATTGCCACCGCAATACCAGCAATTaa                                    >  1:7568191/1‑36 (MQ=255)
 aaTAGATTTAATTGCCACCGCAATACCAGCAATTaa                                    >  1:1322022/1‑36 (MQ=255)
  aTATATTTAAGTTCCATCGCAATACCAGCAATTaaa                                   <  1:1881475/36‑1 (MQ=16)
  aTAGATTTAATTGCCACCGCAATACCAGCAATTaaa                                   >  1:1744020/1‑36 (MQ=255)
  aTAGATTTAATTGCCACCGCAATACCAGCAATTaaa                                   <  1:6601788/36‑1 (MQ=255)
   tAGATTTAATTGCCACCGCAATACCAGCAATTAAAc                                  <  1:3394163/36‑1 (MQ=255)
   tAGATTTAATTGCCACCGCAATACCAGCAATTAAAc                                  >  1:8309014/1‑36 (MQ=255)
   tAGATTTAATTGCCACCGCAATACCAGCAATTAAAc                                  <  1:3227610/36‑1 (MQ=255)
   tAGATTTAATTGCCACCGCAATACCAGCAATTAAAc                                  >  1:4981492/1‑36 (MQ=255)
    aGATTTAATTGCCACCGCAATACCAGCAATTAAAcc                                 <  1:1088707/36‑1 (MQ=255)
    aGATTTAATTGCCACCGCAATACCAGCAATTAAAcc                                 <  1:6147502/36‑1 (MQ=255)
     gATTTAATTGCCACCGCAATACCAGCAATTAAACCg                                >  1:7546406/1‑36 (MQ=255)
      aTTTAATTGCCACCGCAATACCAGCAATTAAACCGc                               >  1:7401300/1‑36 (MQ=255)
       tttAATTGCCACCGCAATACCAGCAATTAAACCGcc                              >  1:7116062/1‑36 (MQ=255)
          aaTTGCCACCGCAATACCAGCAATTAAACCGccacc                           >  1:4750111/1‑36 (MQ=255)
          aaTTGCCACCGCAATACCAGCAATTAAACCGccacc                           >  1:911188/1‑36 (MQ=255)
          aaTTGCCACCGCAATACCAGCAATTAAACCGccacc                           >  1:1064478/1‑36 (MQ=255)
          aaTTGCCACCGCAATACCAGCAATTAAACCGccaac                           >  1:7567290/1‑34 (MQ=255)
           aTTGCCACCGCAATACCAGCAATTAAACCGccacct                          >  1:7411824/1‑35 (MQ=255)
           aTTGCCACCGCAATACCAGCAATTAAACCGccacca                          >  1:7983301/1‑36 (MQ=255)
             tGCCACCGCAATACCAGCAATTAAACCGccaccacc                        >  1:7718513/1‑36 (MQ=255)
             tGCCACCGCAATACCAGCAATTAAACCGccaccacc                        >  1:7706790/1‑36 (MQ=255)
              gCCACCGCAATACCAGCAATTAAACCGccaccacca                       >  1:4432574/1‑36 (MQ=255)
               ccACCGCAATACCAGCAATTAAACCGCCACCACCaa                      >  1:1983887/1‑36 (MQ=255)
                cACCGCAATACCAGCAATTAAACCGCCACCACCAAt                     >  1:2327766/1‑36 (MQ=255)
                cACCGCAATACCAGCAATTAAACCGCCACCACCAAt                     >  1:479402/1‑36 (MQ=255)
                cACCGCAATACCAGCAATTAAACCGCCACCACCAAt                     >  1:730579/1‑36 (MQ=255)
                 aCCGCAATACCAGCAATTAAACCGCCACCACCAAtt                    >  1:2259377/1‑36 (MQ=255)
                 aCCGCAATACCAGCAATTAAACCGCCACCACCAAtt                    >  1:5626561/1‑36 (MQ=255)
                 aCCGCAATACCAGCAATTAAACCGCCACCACCAAtt                    >  1:1829096/1‑36 (MQ=255)
                   caaaaTACAAGAAATAAAACCGCCACCACCAATTgg                  <  1:6865753/33‑1 (MQ=16)
                   cGCAATACCAGCAATTAAACCGCCACCACCAATTgg                  >  1:1222786/1‑36 (MQ=255)
                   cGCAATACCAGCAATTAAACCGCCACCACCAATTgg                  >  1:4114360/1‑36 (MQ=255)
                    gCAATACCAGCAATTAAACCGCCACCACCAATTGGc                 >  1:4437787/1‑36 (MQ=255)
                    gCAATACCAGCAATTAAACCGCCACCACCAATTGGc                 >  1:4056662/1‑36 (MQ=255)
                    gCAATACCAGCAATTAAACCGCCACCACAAATTGGc                 >  1:6760058/1‑36 (MQ=255)
                     aaaaacaaaaaataAAACCGCCACCACCAATTGGca                >  1:2665105‑M1/15‑36 (MQ=255)
                      aaTACCAGCAATTAAACCGCCACCACCAATTGGcac               >  1:4759113/1‑36 (MQ=255)
                      aaTACCAGCAATTAAACCGCCACCACCAATTGGcac               >  1:4837920/1‑36 (MQ=255)
                       aTACCAGAAATAAAACCGCCACCACCAATTGGcaca              <  1:5763206/36‑1 (MQ=37)
                        tACCAGCAATTAAACCGCCACCACCAATTGGCACaa             >  1:7864847/1‑36 (MQ=255)
                        tACCAGCAATTAAACCGCCACCACCAATTGGCACaa             >  1:1880121/1‑36 (MQ=255)
                        tACCAGCAATTAAACCGCCACCACCAATTGGCACaa             >  1:713692/1‑36 (MQ=255)
                         acaagaaataAAACCGCCACCACCAATTGGCACAAt            >  1:4145380‑M1/11‑36 (MQ=255)
                         acaaaaaattAAACCGCCACCACCAATTGGCACAAt            >  1:7820478‑M1/11‑36 (MQ=255)
                         aaaaaaaataAAACCGCCACCACCAATTGGCACAAt            >  1:3975229‑M1/11‑36 (MQ=255)
                         aCCAACAATAAAACCGCCACCACCAATTGGCACAAt            <  1:1422622/36‑1 (MQ=37)
                          ccAGAAATTAAACCGCCACCACCAATTGGCACAATc           <  1:774010/36‑1 (MQ=255)
                          aaagaaataAAACCGCCACCACCAATTGGCACAATc           >  1:1019580‑M1/10‑36 (MQ=255)
                           cAGCAATTAAACCGCCACCACCAATTGGCACAATCa          >  1:7536339/1‑36 (MQ=255)
                           cAGCAATTAAACCGCCACCACCAATTGGCACAATCa          >  1:2151836/1‑36 (MQ=255)
                              aaataAACCCGCCACCACCAATTGGCACAATCACGt       >  1:5694395‑M1/6‑36 (MQ=255)
                              aaaaaAAACCGCCACCACCAATTGGCACAATCACGt       >  1:5006785‑M1/6‑36 (MQ=255)
                               aaaaAAACCGCCACCACCAATTGGCACAATCACGtt      >  1:5965446‑M1/5‑36 (MQ=255)
                                 ttAAACCGCCACCACCAATTGGCACAATCACGTTAt    >  1:7028653/1‑36 (MQ=255)
                                  tAAACCGCCACCACCAATTGGCACAATCACGTTATc   >  1:4848602/1‑36 (MQ=255)
                                  aAATCCGCCACCACCAATTGGCACAATCACGTTATc   >  1:1777642‑M1/2‑36 (MQ=255)
                                  aAAACCGCCACCACCAATTGGCACAATCACGTTATc   >  1:7849497‑M1/2‑36 (MQ=255)
                                   AAACCGCCACCACCAATTGGCACAATCACGTTATCg  >  1:6964972‑M1/1‑36 (MQ=255)
                                   |                                   
TAATAGATTTAATTGCCACCGCAATACCAGCAATTAAACCGCCACCACCAATTGGCACAATCACGTTATCG  >  REL606/3198465‑3198535

Base quality scores:  ATCG  < 3 ≤  ATCG  < 11 ≤  ATCG  < 23 ≤  ATCG  < 30 ≤  ATCG  < 37 ≤  ATCG