New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_014034 158479 =NA (NA)3 (0.010) 3/194 28.2 0.03 coding (53/1176 nt) tuf1 translation elongation factor Tu
?NC_014034 3557414 = 191 (0.780)coding (380/984 nt) RCAP_rcc03360 iron siderophore/cobalamin ABC transporter permease
Rejected: Position hash score below cutoff.

Only 100 of 125 total aligned reads displayed.

GCCATATTCACTGGGCCGACGGCGGCCCCCGAACCAGAGGCGCCCGAAATGGCAAAGGCAAAGTTTGAACGGAACAAGCCGCACGTGAACATCGGCACGATCGGCCACGTGGACCACGGCAAGACGACGCTGACGGCGGCGATCACGA  >  NC_014034/158379‑158526
                                                                                                    |                                               
gcCATATTCACTGGGCCGACGGCGGCCCCCGAACCAGAGGCGCCCGAAATGGCAAAGGCAAAGTTTGAACGGAACAAGCCGCACGTGAACATCGGCACGAt                                                 >  1:3586845/1‑101 (MQ=255)
gcCATATTCACTGGGCCGACGGCGGCCCCCGAACCAGAGGCGCCCGAAATGGCAAAGGCAAAGTTTGAACGGAACAAGCCGCACGTGAACATCGGCACGAt                                                 >  2:24770/1‑101 (MQ=255)
 ccATATTCACTGGGCCGACGGCGGCCCCCGAACCAGAGGCGCCCGAAATGGCAAAGGCAAAGTTTGAACGGAACAAGCCGCACGTGAACATCGGCACGATc                                                >  1:3928303/1‑101 (MQ=255)
 ccATATTCACTGGGCCGACGGCGGCCCCCGAACCAGAGGCGCCCGAAATGGCAAAGGCAAAGTTTGAACGGAACAAGCCGCACGTGAACATCGGCACGATc                                                <  1:76917/101‑1 (MQ=255)
  cATATTCACTGGGCCGACGGCGGCCCCCGAACCAGAGGCGCCCGAAATGGCAAAGGCAAAGTTTGAACGGAACAAGCCGCACGTGAACATCGGCACGATCg                                               >  1:338213/1‑101 (MQ=255)
  cATATTCACTGGGCCGACGGCGGCCCCCGAACCAGAGGCGCCCGAAATGGCAAAGGCAAAGTTTGAACGGAACAAGCCGCACGTGAACATCGGCACGATCg                                               >  1:1726365/1‑101 (MQ=255)
  cATATTCACTGGGCCGACGGCGGCCCCCGAACCAGAGGCGCCCGAAATGGCAAAGGCAAAGTTTGAACGGAACAAGCCGCACGTGAACATCGGCACGATCg                                               >  2:1835183/1‑101 (MQ=255)
  cATATTCACTGGGCCGACGGCGGCCCCCGAACCAGAGGCGCCCGAAATGGCAAAGGCAAAGTTTGAACGGAACAAGCCGCACGTGAACATCGGCACGATCg                                               >  1:4642936/1‑101 (MQ=255)
   atatTCACTGGGCCGACGGCGGCCCCCGAACCAGAGGCGCCCGAAATGGCAAAGGCAAAGTTTGAACGGAACAAGCCGCACGTGAACATCGGCACGATCgg                                              >  1:660548/1‑101 (MQ=255)
   atatTCACTGGGCCGACGGCGGCCCCCGAACCAGAGGCGCCCGAAATGGCAAAGGCAAAGTTTGAACGGAACAAGCCGCACGTGAACATCGGCACCATCgg                                              >  1:2274300/1‑101 (MQ=255)
    tatTCACTGGGCCGACGGCGGCCCCCGAACCAGAGGCGCCCGAAATGGCAAAGGCAAAGTTTGAACGGAACAAGCCGCACGTGAACATCGGCACGATCGGc                                             >  1:754395/1‑101 (MQ=255)
     atTCACTGGGCCGACGGCGGCCCCCGAACCAGAGGCGCCCGAAATGGCAAAGGCAAAGTTTGAACGGAACAAGCCGCACGTGAACATCGGCACGATCGGcc                                            >  2:258221/1‑101 (MQ=255)
     atTCACTGGGCCGACGGCGGCCCCCGAACCAGAGGCGCCCGAAATGGCAAAGGCAAAGTTTGAACGGAACAAGCCGCACGTGAACATCGGCACGATCGGcc                                            >  1:733786/1‑101 (MQ=255)
      ttCACTGGGCCGACGGCGGCCCCCGAACCAGAGGCGCCCGAAATGGCAAAGGCAAAGTTTGAACGGAACAAGCCGCACGTGAACATCGGCACGATCGGCCa                                           <  1:3213246/101‑1 (MQ=255)
      ttCACTGGGCCGACGGCGGCCCCCGAACCAGAGGCGCCCGAAATGGCAAAGGCAAAGTTTGAACGGAACAAGCCGCACGTGAACATCGGCACGATCGGCCa                                           >  1:4247845/1‑101 (MQ=255)
      ttCACTGGGCCGACGGCGGCCCCCGAACCAGAGGCGCCCGAAATGGCAAAGGCAAAGTTTGAACGGAACAAGCCGCACGTGAACATCGGCACGATCGGCCa                                           <  1:2516345/101‑1 (MQ=255)
      ttCACTGGGCCGACGGCGGCCCCCGAACCAGAGGCGCCCGAAATGGCAAAGGCAAAGTTTGAACGGAACAAGCCGCACGTGAACATCGGCACGATCGGCCa                                           >  1:211496/1‑101 (MQ=255)
      ttCACTGGGCCGACGGCGGCCCCCGAACCAGAGGCGCCCGAAATGGCAAAGGCAAAGTTTGAACGGAACAAGCCGCACGTGAACATCGGCACGATCGGCCa                                           >  2:1921279/1‑101 (MQ=255)
      ttCACGGGGCCGACGGCGGCCCCCGAACCAGAGGCGCCCGAAATGGCAAAGGCAAAGTTTGAACGGAACAAGCCGCACGTGAACATCGGCACGATCGGCCa                                           <  1:1080579/101‑1 (MQ=255)
       tCACTGGGCCGACGGCGGCCCCCGAACCAGAGGCGCCCGAAATGGCAAAGGCAAAGTTTGAACGGAACAAGCCGCACGTGAACATCGGCACGATCGGCCAc                                          <  1:3050011/101‑1 (MQ=255)
       tCACTGGGCCGACGGCGGCCCCCGAACCAGAGGCGCCCGAAATGGCAAAGGCAAAGTTTGAACGGAACAAGCCGCACGTGAACATCGGCACGATCGGCCAc                                          <  2:4994535/101‑1 (MQ=255)
       tCACGGGGCCGACGGCGGCCCCCGAACCAGAGGCGCCCGAAATGGCAAAGGCAAAGTTTGAACGGAACAAGCCGCACGTGAACATCGGCACGATCGGCCAc                                          <  2:2607134/101‑1 (MQ=255)
        cACTGGGCCGACGGCGGCCCCCGAACCAGAGGCGCCCGAAATGGCAAAGGCAAAGTTTGAACGGAACAAGCCGCACGTGAACATCGGCACGATCGGCCACg                                         >  1:4557846/1‑101 (MQ=255)
        cACTGGGCCGACGGCGGCCCCCGAACCAGAGGCGCCCGAAATGGCAAAGGCAAAGTTTGAACGGAACAAGCCGCACGTGAACATCGGCACGATCGGCCACg                                         >  1:1020833/1‑101 (MQ=255)
        cACTGGGCCGACGGCGGCCCCCGAACCAGAGGCGCCCGAAATGGCAAAGGCAAAGTTTGAACGGAACAAGCCGCACGTGAACATCGGCACGATCGGCCACg                                         >  1:1233677/1‑101 (MQ=255)
         aCTGGGCCGACGGCGGCCCCCGAACCAGAGGCGCCCGAAATGGCAAAGGCAAAGTTTGAACGGAACAAGCCGCACGTGAACATCGGCACGATCGGCCACgg                                        >  2:5032045/1‑100 (MQ=255)
           tGGGCCGACGGCGGCCCCCGAACCAGAGGCGCCCGAAATGGCAAAGGCAAAGTTTGAACGGAACAAGCCGCACGTGAAGATCGGCACGATCGGCCACGTgg                                      >  2:832023/1‑101 (MQ=255)
           tGGGCCGACGGCGGCCCCCGAACCAGAGGCGCCCGAAATGGCAAAGGCAAAGTTTGAACGGAACAAGCCGCACGTGAACATCGGCACGATCggatcggaag                                      >  1:3462192/1‑93 (MQ=255)
            gggCCGACGGCGGCCCCCGAACCAGAGGCGCCCGAAATGGCAAAGGCAAAGTTTGAACGGAACAAGCCGCACGTGAACATCGGCACGATCGGCCACGTGGa                                     <  1:2374267/101‑1 (MQ=255)
            gggCCGACGGCGGCCCCCGAACCAGAGGCGCCCGAAATGGCAAAGGCAAAGTTTGAACGGAACAAGCCGCACGTGAACATCGGCACGATCGGCCACGTGGa                                     >  2:3183647/1‑101 (MQ=255)
            gggCCGACGGCGGCCCCCGAACCAGAGGCGCCCGAAATGGCAAAGGCAAAGTTTGAACGGAACAAGCCGCACGTGAACATCGGCACGATCGGCCACGTGGa                                     >  1:2754732/1‑101 (MQ=255)
            gggCCGACGGCGGCCCCCGAACCAGAGGCGCCCGAAATGGCAAAGGCAAAGTTTGAACGGAACAAGCCGCACGTGAACATCGGCACGATCGGCCACGTGGa                                     >  1:4108349/1‑101 (MQ=255)
            ggcccGACGGCGGCCCCCGAACCAGAGGCGCCCGAAATGGCAAAGGCAAAGTTTGAACGGAACAAGCCGCACGTGAACATCGGCACGATCGGCCACGTGGa                                     <  2:1233560/98‑1 (MQ=255)
              gCCGGCGGCGGCCCCCGAACCAGGGGCGCCCGAAATGGCAAAGGCAAAGTTTGAACGGAACAAGCCGCACGTGAACATCGGCACGATCGGCCACGTGGAcc                                   <  2:2331754/101‑1 (MQ=255)
              gCCGACGGCGGCCCCCGAACCAGAGGCGCCCGAAATGGCAAGGGCAAAGTTTGAACGGAACAAGCCGCACGTGAACATCGGCACGATCGGCCACGTGGAcc                                   <  2:660459/101‑1 (MQ=255)
               ccGAGGGCGGCCCCCGATCAAGGGGCGCCCGAATGGGCAAATGAAAAGTTTGAACGGAACAAGCCGCACGTGAACATCGGCACGATCGGCCACGTGGACCa                                  <  2:471610/101‑1 (MQ=17)
                cGACGGCGGCCCCCGAACCAGAGGCGCCCGAAATGGCAAAGGCAAAGTTTGAACGGAACAAGCCGCACGTGAACATCGGCACGATCGGCCACGTGGACCAc                                 <  2:339543/101‑1 (MQ=255)
                 gACGGCGGCCCCCGAACCAGAGGCGCCCGAAATGGCAAAGGCAAAGTTTGAACGGAACAAGCCGCACGTGAACATCGGCACGATCGGCCACGTGGACCACg                                >  1:429016/1‑101 (MQ=255)
                  aCGGCGGCCCCCGAACCAGAGGCGCCCGAAATGGCAAAGGCAAAGTTTGAACGGAACAAGCCGCACGTGAACATCGGCACGATCGGCCACGTGGACCACgg                               >  1:2174684/1‑101 (MQ=255)
                  aCGGCGGCCCCCGAACCAGAGGCGCCCGAAATGGCAAAGGCAAAGTTTGAACGGAACAAGCCGCACGTGAACATCGGCACGATCGGCCACGTGGACCACgg                               >  2:1687174/1‑101 (MQ=255)
                   cggcggCCCCCGAACCAGAGGCGCCCGAAATGGCAAAGGCAAAGTTTGAACGGAACAAGCCGCACGTGAACATCGGCACGATCGGCCACGTGGACCACGGc                              <  1:4127026/101‑1 (MQ=255)
                   cggcggCCCCCGAACCAGAGGCGCCCGAAATGGCAAAGGCAAAGTTTGAACGGAACAAGCCGCACGTGAACATCGGCACGATCGGCCACGTGGACCACGGc                              >  2:2183354/1‑101 (MQ=255)
                   cggcggCCCCCGAACCAGAGGCGCCCGAAATGGCAAAGGCAAAGTTTGAACGGAACAAGCCGCACGTGAACATCGGCACGATCGGCCACGTGGACCACGGc                              <  1:2055852/101‑1 (MQ=255)
                   cggcggCCCCCGAACCAGAGGCGCCCGAAATGGCAAAGGCAAAGTTTGAACGGAACAAGCCGCACGTGAACATCGGCACGATCGGCCACGTGGACCACGGc                              <  1:255466/101‑1 (MQ=255)
                   cggcggCCCCCGAACCAGAGGCGCCCGAAATGGCAAAGGCAAAGTTTGAACGGAACAAGCCGCACGTGAACATCGGCACGATCGGCCACGTGGACCACGGc                              <  1:4648465/101‑1 (MQ=255)
                    tctggcCCCCGAACCAGAGGCCCCCGAAATGGCAAAGGCAAAGTTTGAACGGAACAAGCCGCACGTGAACATCGGCACGATCGGCCACGTGGACCACGGCa                             <  1:135222/98‑1 (MQ=255)
                    ggcggcCCCCGAACCAGAGGCGCCCGAAATGGCAAAGGCAAAGTTTGAACGGAACAAGCCGCACGTGAACATCGGCACGATCGGCCACGTGGACCACGGCa                             >  1:410187/1‑101 (MQ=255)
                    ggcggcCCCCGAACCAGAGGCGCCCGAAATGGCAAAGGCAAAGTTTGAACGGAACAAGCCGCACGTGAACATCGGCACGATCGGCCACGTGGACCACGGCa                             <  1:3306375/101‑1 (MQ=255)
                      cggcCCCCGAACCAGAGGCGCCCGAAATGGCAAAGGCAAAGTTTGAACGGAACAAGCCGCACGTGAACATCGGCACGATCGGCCACGTGGACCACGGCAAg                           <  1:873234/101‑1 (MQ=255)
                      cggcCCCCGAACCAGAGGCGCCCGAAATGGCAAAGGCAAAGTTTGAACGGAACAAGCCGCACGTGAACATCGGCACGATCGGCCACGTGGACCACGGCAAg                           <  1:2810956/101‑1 (MQ=255)
                      cggcCCCCGAACCAGAGGCGCCCGAAATGGCAAAGGCAAAGTTTGAACGGAACAAGCCGCACGTGAACATCGGCACGATCGGCCACGTGGACCACGGCAAg                           >  1:1999117/1‑101 (MQ=255)
                       ggcCCCCGAACCAGAGGCGCCCGAAATGGCAAAGGCAAAGTTTGAACGGAACAAGCCGCACGTGAACATCGGCACGATCGGCCACGTGGACCACGGCAAga                          >  2:135206/1‑101 (MQ=255)
                        gcCCCCGAACCAGAGGCGCCCGAAATGGCAAAGGCAAAGTTTGAACGGAACAAGCCGCACGTGAACATCGGCACGATCGGCCACGTGGACCACGGCAAgac                         <  1:689332/101‑1 (MQ=255)
                        gcCCCCGAACCAGAGGCGCCCGAAATGGCAAAGGCAAAGTTTGAACGGAACAAGCCGCACGTGAACATCGGCACGATCGGCCACGTGGACCACGGCAAgac                         <  1:4767483/101‑1 (MQ=255)
                        gcCCCCGAACCAGAGGCGCCCGAAATGGCAAAGGCAAAGTTTGAACGGAACAAGCCGCACGTGAACATCGGCACGATCGGCCACGTGGACCACGGCAAgac                         >  2:1084128/1‑101 (MQ=255)
                        gcCCCCGAACCAGAGGCGCCCGAAATGGCAAAGGCAAAGTTTGAACGGAACAAGCCGCACGTGAACATCGGCACGATCGGCCACGTGGACCACGGCAAgac                         >  2:4091013/1‑101 (MQ=255)
                         cccccGAACCAGAGGCGCCCGAAATGGCAAAGGCAAAGTTTGAACGGAACAAGCCGCACGTGAACATCGGCACGATCGGCCACGTGGACCACGGCAAgacg                        <  1:4421688/101‑1 (MQ=255)
                          ccccGAACCAGAGGCGCCCGAAATGGCAAAGGCAAAGTTTGAACGGAACAAGCCGCACGTGAACATCGGCACGATCGGCCACGTGGACCACGGCAAgacga                       <  1:486960/101‑1 (MQ=255)
                          ccccGAACCAGAGGCGCCCGAAATGGCAAAGGCAAAGTTTGAACGGAACAAGCCGCACGTGAACATCGGCACGATCGGCCACGTGGACCACGGCAAgacga                       >  1:3529790/1‑101 (MQ=255)
                           cccGACCCAGGGGCGCCCGAAATGGAAAAGGAAAAGTTTGAACGAAACAAGCCGCACGTGAACATCGGCACGATCGGCCACGTGGACCACGGCAAgccgat                      <  2:2303237/101‑2 (MQ=14)
                           cccGAACCAGAGGCGCCCGAAATGGCAAAGGCAAAGTTTGAACGGAACAAGCCGCACGTGAACATCGGCACGATCGGCCACGTGGACCACGGCAAgacgac                      <  2:4108032/101‑1 (MQ=255)
                           cccGAACCAGAGGCGCCCGAAATGGCAAAGGCAAAGTTTGAACGGAACAAGCCGCACGTGAACATCGGCACGATCGGCCACGTGGACCACGGCAAgacgac                      <  2:1613198/101‑1 (MQ=255)
                           cccGAACCAGAGGCGCCCGAAATGGCAAAGGAAAAGTTTGAACGGAACAAGCCGCACGTGAACATCGGCACGATCGGCCACGTGGACCACGGCAAgacgac                      <  2:120197/101‑1 (MQ=255)
                            ccGAACCAGAGGCGCCCGAAATGGCAAAGGCAAAGTTTGAACGGAACAAGCCGCACGTGAACATCGGCACGATCGGCCACGTGGACCACGGCAAGacgacg                     >  1:2029400/1‑101 (MQ=255)
                             cGAACCAGAGGCGCCCGAAATGGCAAAGGCAAAGTTTGAACGGAACAAGCCGCACGTGAACATCGGCACGATCGGCCACGTGGACCACGGCAAGACGACGc                    <  2:1405824/101‑1 (MQ=255)
                              gAACCAGAGGCGCCCGAAATGGCAAAGGCAAAGTTTGAACGGAACAAGCCGCACGTGAACATCGGCACGATCGGCCACGTGGACCACGGCAAGACGACGCt                   >  1:4396149/1‑101 (MQ=255)
                              gAACCAGAGGCGCCCGAAATGGCAAAGGCAAAGTTTGAACGGAACAAGCCGCACGTGAACATCGGCACGATCGGCCACGTGGACCACGGCAAGACGACGCt                   >  1:2433950/1‑101 (MQ=255)
                              gAACCAGAGGCGCCCGAAATGGCAAAGGCAAAGTTTGAACGGAACAAGCCGCACGTGAACATCGGCACGATCGGCCACGTGGACCACGGCAAGACGACGCt                   >  2:4817879/1‑101 (MQ=255)
                              gAACCAGAGGCGCCCGAAATGGCAAAGGCAAAGTTTGAACGGAACAAGCCGCACGTGAACATCGGCACGATCGGCCACGTGGACCACGGCAAGACGACGCt                   <  2:718084/101‑1 (MQ=255)
                              gAACCAGAGGCGCCCGAAATGGCAAAGGCAAAGTTTGAACGGAACAAGCCGCACGTGAACATCGGCACGATCGGCCACGTGGACCACGGCAAGACGACGCt                   >  2:4882769/1‑101 (MQ=255)
                               aaCCGGAGGGGCCCGAAATGGCAAAGGCAAAGTTTGAACGGAACAAGCCGCACGTGAACATCGGCACGATCGGCCACGTGGACCACGGCAAGACGACGCTg                  <  2:2274145/101‑1 (MQ=255)
                               aaCCAGAGGCGCCCGAAATGGCAAAGGCAAAGTTTGAACGGAACAAGCCGCACGTGAACATCGGCACGATCGGCCACGTGGACCACGGCAAGACGACGCTg                  >  2:2177932/1‑101 (MQ=255)
                                aCCAGAGGCGCCCGAAATGGCAAAGGCAAAGTTTGAACGGAACAAGCCGCACGTGAACATCGGCACGATCGGCCACGTGGACCACGGCAAGACGACGCTGa                 >  1:3484145/1‑101 (MQ=255)
                                 ccAGAGGCGCCCGAAATGGCAAAGGCAAAGTTTGAACGGAACAAGCCGCACGTGAACATCGGCACGATCGGCCACGTGGACCACGGCAAGACGACGCTGAc                >  2:2153530/1‑101 (MQ=255)
                                 ccAGAGGCGCCCGAAATGGCAAAGGCAAAGTTTGAACGGAACAAGCCGCACGTGAACATCGGCACGATCGGCCACGTGGACCACGGCAAGACGACGCTGAc                <  1:3436301/101‑1 (MQ=255)
                                  cAGAGGCGCCCGAAATGGCAAAGGCAAAGTTTGAACGGAACAAGCCGCACGTGAACATCGGCACGATCGGCCACGTGGACCACGGCAAGACGACGCTGAcg               <  2:2043209/101‑1 (MQ=255)
                                  cAGAGGCGCCCGAAATGGCAAAGGCAAAGTTTGAACGGAACAAGCCGCACGTGAACATCGGCACGATCGGCCACGTGGACCACGGCAAGACGACGCTGAcg               >  1:554388/1‑101 (MQ=255)
                                  cAGAGGCGCCCGAAATGGCAAAGGCAAAGTTTGAACGGAACAAGCCGCACGTGAACATCGGCACGATCGGCCACGTGGACCACGGCAAGACGACGCTGAcg               <  2:310356/101‑1 (MQ=255)
                                   agagGCGCCCGAAATGGCAAAGGCAAAGTTTGAACGGAACAAGCCGCACGTGAACATCGGCACGATCGGCCACGTGGACCACGGCAAGACGACGCTGAcgg              >  2:4793914/1‑101 (MQ=255)
                                   agagGCGCCCGAAATGGCAAAGGCAAAGTTTGAACGGAACAAGCCGCACGTGAACATCGGCACGATCGGCCACGTGGACCACGGCAAGACGACGCTGAcgg              <  2:2427705/101‑1 (MQ=255)
                                   agagGCGCACGAAATGGCAAAGGCAAAGTTTGAACGGAACAAGCCGCACGTGAACATCGGCACGATCGGCCACGTGGACCACGGCAAGACGACGCTGAcgg              >  2:4419792/1‑101 (MQ=255)
                                    gagGCGCCCGAAATGGCAAAGGCAAAGTTTGAACGGAACAAGCCGCACGTGACCATCGGCACGATCGGCCACGTGGACCACGGGAAgacgacacgaccggg             >  2:18878/1‑92 (MQ=14)
                                    gagGCGCCCGAAATGGCAAAGGCAAAGTTTGAACGGAACAAGCCGCACGTGAACATCGGCACGATCGGCCACGTGGACCACGGCAAGACGACGCTGAcggc             >  1:1585434/1‑101 (MQ=255)
                                     agGCGCCCGAAATGGCAAAGGCAAAGTTTGAACGGAACAAGCCGCACGTGAACATCGGCACGATCGGCCACGTGGACCACGGCAAGACGACGCTGAcggcg            >  2:4863760/1‑101 (MQ=40)
                                         gcCCGAAATGGCAAAGGCAAAGTTTGAACGGAACAAGCCGCACGTGAACATCGGCACGATCGGCCACGTGGACCACGGCAAGACGACGCTGACGGCGGCGa        >  1:961991/1‑101 (MQ=37)
                                           ccGAAAGGGCAAAGGCAAAGTTTGAACGGAACAAGCCGCACGTGAACATCGGCACGATCGGCCACGTGGACCACGGCAAGACGACGCTGACGGCGGCGATc      <  2:4670102/101‑1 (MQ=18)
                                            cGAAATGGCAAAGGCAAAGTTTGAACGGAACAAGCCGCACGTGGACATCGGCACGATCGGCCAGGTGGACCACGGCAAGACGACGCTGACGGCGGCGATCa     >  1:4321175/1‑101 (MQ=16)
                                            cGAAATGGCAAAGGCAAAGTTTGAACGGAACAAGCCGCACGTGAACATCGGCACGATCGGCCACGTGGACCACGGCAAGACGACGCTGACGGCGGCGATCa     <  1:3137967/101‑1 (MQ=33)
                                             gAAATGGCAAAGGCAAAGTTTGAACGGAACAAGCCGCACGTGAACATCGGCACGATCGGCCACGTGGACCACGGCAAGACGACGCTGACGGCGGCGATCAc    >  2:839280/1‑101 (MQ=32)
                                             gAAATGGCAAAGGCAAAGTTTGAACGGAACAAGCCGCACGTGAACATCGGCACGATCGGCCACGTGGACCACGGCAAGACGACGCTGACGGCGGCGATCAc    >  2:3768820/1‑101 (MQ=32)
                                             gAAATGGCAAAGGCAAAGTTTGAACGGAACAAGCCGCACGTGAACATCGGCACGATCGGCCACGTGGACCACGGCAAGACGACGCTGACGGCGGCGATCAc    <  2:1517138/101‑1 (MQ=32)
                                             gAAATGGCAAAGGCAAAGTTTGAACGGAACAAGCCGCACGTGAACATCGGCACGATCGGCCACGTGGACCACGGCAAGACGACGCTGACGGCGGCGATCAc    >  1:368351/1‑101 (MQ=32)
                                              aaaTGGCAAAGGCAAAGTTTGCACGGAACAAGCCGCACGTGAACATCGGCACGATCGGCCACGTGGACCACGGCAAGACGACGCTGACGGCGGCGATCACg   >  1:2874982/1‑101 (MQ=14)
                                              aaaTGGCAAAGGCAAAGTTTGAACGGAACAAGCCGCACGTGAACATCGGCACGATCGGCCACGTGGACCACGGCAAGACGACGCTGACGGCGGCGATCACg   >  2:1665427/1‑101 (MQ=31)
                                              aaaTGGCAAAGGCAAAGTTTGAACGGAACAAGCCGCACGTGAACATCGGCACGATCGGCCACGTGGACCACGGCAAGACGACGCTGACGGCGGCGATCACg   >  1:912320/1‑101 (MQ=31)
                                              aaaTGGCAAAGGCAAAGTTTGAACGGAACAAGCCGCACGTGAACATCGGCACGATCGGCCACGTGGACCACGGCAAGACGACGCTGACGGCGGCGATCACg   <  1:2178087/101‑1 (MQ=31)
                                              aaaTGGCAAAGGCAAAGTTTGAACGGAACAAGCCGCACGTGAACATCGGCACGATCGGCCAAGTGGACCACGGCAAGACGACGCTGACGGCGGCCATCACg   >  2:3981968/1‑101 (MQ=9)
                                               aaTGGCAAAGGCAAAGTTTGAACGGAACAAGCCGCACGTGAACATCGGCACGATCGGCCACGTGGACCACGGCAAGACGACGCTGACGGCGGCGATCACGa  <  2:338160/101‑1 (MQ=11)
                                               aaTGGCAAAGGCAAAGTTTGAACGGAACAAGCCGCACGTGAACATCGGCACGATCGGCCACGTGGACCACGGCAAGACGACGCTGACGGCGGCGATCACGa  <  1:790630/101‑1 (MQ=11)
                                               aaTGCCAAAGGAAAAGGTTGAACGGAACAAGCCGCACGTGAACATCGGCACGATCGGCCACGTGGAACACGGCAATACGACGCTGACGGCGGTGATCACGa  <  2:151950/101‑1 (MQ=11)
                                                                                                    |                                               
GCCATATTCACTGGGCCGACGGCGGCCCCCGAACCAGAGGCGCCCGAAATGGCAAAGGCAAAGTTTGAACGGAACAAGCCGCACGTGAACATCGGCACGATCGGCCACGTGGACCACGGCAAGACGACGCTGACGGCGGCGATCACGA  >  NC_014034/158379‑158526

Base quality scores:  ATCG  < 3 ≤  ATCG  < 20 ≤  ATCG  < 29 ≤  ATCG  < 33 ≤  ATCG  < 41 ≤  ATCG