New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_014034 = 854430202 (0.920)2 (0.010) 2/180 28.3 0.01 coding (753/1134 nt) hypD hydrogenase expression/formation protein HypD
?NC_014034 1154525 = 137 (0.560)coding (889/1065 nt) lspL1 UDP‑glucuronate 5'‑epimerase
Rejected: Position hash score below cutoff.

Only 100 of 140 total aligned reads displayed.

TCGGCAATTCCGACAAGGTGCGGCTGCTGGATTTCGTCGAGGCGATCGAGGCCTGTCTGGGCAAGACGGCGCTCCGCAATTACCTGCCGATGCAAAAGGGTGATGTGCCCGCGACCTTTGCCGATGCGACGCTGCTGCACAGCCTGACCGG  >  NC_014034/1154475‑1154625
                                                  |                                                                                                    
tCGGCAATTCCGACAAGGTGCGGCTGCTGGATTTCGTCGAGGCGATCGAGGCCTGTCTGGGCAAGACGGCGCTCCGCAATTACCTGCCGATGCAAAAGGGt                                                    <  2:4585299/101‑1 (MQ=11)
tCGGCAATTCCGACAAGGTGCGGCTGCTGGATTTCGTCGAGGCGATCGAGGCCTGTCTGGGCAAGACGGCGCTCCGCAATTACCTGCCGATGCAAAAGGGt                                                    <  2:1925154/101‑1 (MQ=11)
 cGGCAATTCCGACAAGGTGCGGCTGCTGGATTTCGTCGAGGCGATCGAGGCCTGTCTGGGCAAGACGGCGCTCCGCAATTACCTGCCGATGCAAAAGGGTg                                                   >  1:437772/1‑101 (MQ=31)
 cGGCAATTCCGACAAGGTGCGGCTGCTGGATTTCGTCGAGGCGATCGAGGCCTGTCTGGGCAAGACGGCGCTCCGCAATTACCTGCCGATGCAAAAGGGTg                                                   <  2:4246781/101‑1 (MQ=31)
  ggCAATTCCGACAAGGTGCGGCTGCTGGATTTCGTCGAGGCGATCGAGGCCTGTCTGGGCAAGACGGCGCTCCGCAATTACCTGCCGATGCAAAAGGGTGa                                                  <  2:3106024/101‑1 (MQ=32)
   gCAATTCCGGCAAGGTGCGGCTGCTGGATTTCGTCGAGGCGATCGAGGCCTGTCTGGGCAAGACGGCGCTCCGCAATTACCTGCCGATGCAAAAGGGTGAt                                                 <  2:5001839/101‑1 (MQ=18)
   gCAATTCCGACAAGGTGCGGCTGCTGGATTTCGTCGAGGCGATCGAGGCCTGTCTGGGCAAGACGGCGCTCCGCAATTACCTGCCGATGCAAAAGGGTGAt                                                 <  1:3589807/101‑1 (MQ=33)
   gCAATTCCGACAAGGTGCGGCTGCTGGATTTCGTCGAGGCGATCGAGGCCTGTCTGGGCAAGACGGCGCTCCGCAATTACCTGCCGATGCAAAAGGGTGAt                                                 >  1:1493411/1‑101 (MQ=33)
   gCAATTCCGACAAGGTGCGGCTGCTGGATTTCGTCGAGGCGATCGAGGCCTGTCTGGGCAAGACGGCGCTCCGCAATTACCTGCCGATGCAAAAGGGTGAt                                                 <  2:4642622/101‑1 (MQ=33)
   gCAATTCCGACAAGGTGCGGCTGCTGGATTTCGTCGAGGCGATCGAGGCCTGTCTGGGCAAGACGGCGCTCCGCAATTACCCGCCGATGCAAAAGGGTGAt                                                 <  2:2030224/101‑1 (MQ=18)
    cAATTCCGACAAGGTGCGGCTGCTGGATTTCGTCGAGGCGATCGAGGCCTGTCTGGGCAAGACGGCGCTCCGCAATTACCTGCCGATGCAAAAGGGTGagg                                                >  2:3914137/1‑99 (MQ=17)
    cAATTCCGACAAGGTGCGGCTGCTGGATTTCGTCGAGGCGATCGAGGCCTGTCTGGGCAAGACGGCGCTCCGCAATTACCTGCCGATGCAAAAGGGTGAtg                                                <  2:4130397/101‑1 (MQ=34)
     aaTTCCGACAAGGTGCGGCTGCTGGATTTCGTCGAGGCGATCGAGGCCTGTCTGGGCAAGACGGCGCTCCGCAATTACCTGCCGATGCAAAAGGGTGAtgt                                               >  2:136976/1‑101 (MQ=35)
     aaTTCCGACAAGGTGCGGCTGCTGGATTTCGTCGAGGCGATCGAGGCCTGTCTGGGCAAGACGGCGCTCCGCAATTACCTGCCGATGCAAAAGGGTGAtgt                                               <  2:1586436/101‑1 (MQ=35)
      aTTCCGACAAGGTGCGTCTGCTGGATTTCGTCGAGGCGATCGAGGCCTGTCTGGGCAAGACGGCGCTCCGCAATTACCTGCCGATGCAAAAGGGTGAtgtg                                              <  2:1071163/101‑1 (MQ=25)
      aTTCCGACAAGGTGCGGCTGCTGGATTTCGTCGAGGCGATCGAGGCCTGTCTGGGCAAGACGGCGCTCCGCAATTACCTGCCGATGCAAAAGGGTGAtgtg                                              >  2:537144/1‑101 (MQ=35)
       ttCCGACAAGGTGCGGCTGCTGGATTTCGTCGAGGCGATCGAGGCCTGTCTGGGCAAGACGGCGCTCCGCAATTACCTGCCGATGCAAAAGGGTGATGTGc                                             <  1:4693594/101‑1 (MQ=35)
       ttCCGACAAGGTGCGGCTGCTGGATTTCGTCGAGGCGATCGAGGCCTGTCTGGGCAAGACGGCGCTCCGCAATTACCTGCCGATGCAAAAGGGTGATGTGc                                             <  2:3801181/101‑1 (MQ=35)
       ttCCGACAAGGTGCGGCTGCTGGATTTCGTCGAGGCGATCGAGGCCTGTCTGGGCAAGACGGCGCTCCGCAATTACCTGCCGATGCAAAAGGGTGATGTGc                                             <  1:2162479/101‑1 (MQ=35)
       ttCCGACAAGGTGCGGCTGCTGGATTTCGTCGAGGCGATCGAGGCCTGTCTGGGCAAGACGGCGCTCCGCAATTACCTGCCGATGCAAAAGGGTGATGTGc                                             <  2:4936812/101‑1 (MQ=35)
       ttCCGACAAGGTGCGGCTGCTGGATTTCGTCGAGGCGATCGAGGCCTGTCTGGGCAAGACGGCGCTCCGCAATTACCTGCCGATGCAAAAGGGTGAAGTGc                                             >  2:1796796/1‑101 (MQ=17)
        tCCGACAAGGTGCGGCTGCTGGATTTCGTCGAGGCGATCGAGGCCTGTCTGGGCAAGACGGCGCTCCGCAATTACCTGCCGATGCAAAAGGGGGATGTGcc                                            >  1:4232602/1‑101 (MQ=17)
         ccGACAAGGTGCGGCTGCTGGATTTCGTCGAGGCGATCGAGGCCTGTCTGGGCAAGACGGCGCTCCGCAATTACCTGCCGATGCAAAAGGGTGATGTGccc                                           <  2:2329390/101‑1 (MQ=35)
         ccGACAAGGTGCGGCTGCTGGATTTCGTCGAGGCGATCGAGGCCTGTCTGGGCAAGACGGCGCTCCGCAATTACCTGCCGATGCAAAAGGGTGATGTGccc                                           <  1:1363045/101‑1 (MQ=35)
          cGACAAGGTGCGGCTGCTGGATTTCGTCGAGGCGATCGAGGCCTGTCTGGGCAAGACGGCGCTCCGCAATTACCTGCCGATGCAAAAGGGTGATGTGCCcg                                          <  1:1971032/101‑1 (MQ=35)
           gccAAGGTGGGGCTGCTGGATTTCGTCGAGGCGATCGAGGCCTGTCTGGGCAAGACGGCGCTCCGCAATTTCCTGCCGATGCAAAAGGGTGATGTGCCcgc                                         <  1:3172000/99‑1 (MQ=255)
           gACAAGGTGCGGCTGCTGGATTTCGTCGAGGCGATCGAGGCCTGTCTGGGCAAGACGGCGCTCCGCAATTACCTGCCGATGCAAAAGGGTGATGTGCCcgc                                         <  1:3610968/101‑1 (MQ=35)
           gACAAGGTGCGGCTGCTGGATTTCGTCGAGGCGATCGAGGCCTGTCTGGGCAAGACGGCGCTCCGCAATTACCTGCCGATGCAAAAGGGTGATGTGCCcgc                                         >  1:361541/1‑101 (MQ=35)
           gACAAGGTGCGGCTGCTGGATTTCGTCGAGGCGATCGAGGCCTGTCTGGGCAAGACGGCGCTCCGCAATTACCTGCCGATGCAAAAGGGGGATGTGCCcgc                                         >  2:3839225/1‑101 (MQ=17)
           gACAAGGTGCGGCTGCTGGATTTCGTCGAGGCGATCGAGGCCTGGCTCGGCAAGACGGCGCTCCGCAATTACCTGCCGATGCAACAGGGTGATGTGCCcgc                                         >  2:827848/1‑101 (MQ=255)
             cAAGGTGCGGCTGCTGGATTTCGTCGAGGCGATCGAGGCCTGTCTGGGCAAGACGGCGCTCCGCAATTACCTGCCGATGCAAAAGGGTGATGTGCCCGCGa                                       <  2:4116105/101‑1 (MQ=37)
              aaGGTGCGGCTGCTGGATTTCGTCGAGGCGATCGAGGCCTGTCTGGGCAAGACGGCGCTCCGCAATTACCTGCCGATGCAAAAGGGTGATGTGCCCGCGAc                                      >  1:3285857/1‑101 (MQ=35)
              aaGGTGCGGCTGCTGGATTTCGTCGAGGCGATCGAGGCCTGTCTGGGCAAGACGGCGCTCCGCAATTACCTGCCGATGCAAAAGGGTGATGTGCCCGCGAc                                      <  1:844160/101‑1 (MQ=37)
               aGGTGCGGCTGCTGGATTTCGTCGAGGCGATCGAGGCCTGTCTGGGCAAGACGGCGCTCCGCAATTACCTGCCGATGCAAAAGGGTGATGTGCCCGCGAcc                                     <  2:265909/101‑1 (MQ=37)
               aGGTGCGGCTGCTGGATTTCGTCGAGGCGATCGAGGCCTGTCTGGGCAAGACGGCGCTCCGCAATTACCTGCCGATGCAAAAGGGTGATGTGCCCGCGAcc                                     >  1:5046128/1‑101 (MQ=35)
               aGGTGCGGCTGCGGGATTTCGTCGAGGCGATCGAGGCCTGTCTGGGCAAGACGGCGCTCCGCAATTACCTGCCGATGCAAAAGGGTGATGTGCCCGCGAcc                                     <  1:175532/101‑1 (MQ=37)
                 gTGCGGCTGCTGGATTTCGTCGAGGCGATCGAGGCCTGTCTGGGCAAGACGGCGCTCCGCAATTACCTGCCGATGCAAAAGGGTGATGTGCCCGCGACCtt                                   <  2:1081045/101‑1 (MQ=38)
                  tGCGGCTGCTGGATTTCGTCGAGGCGATCGAGGCCTGTCTGGGCAAGACGGCGCTCCGCAATTACCTGCCGATGCAAAAGGGTGATGTGCCCGCGACCttt                                  >  1:3745870/1‑101 (MQ=35)
                  ggCGGCTGCTGGATTTCGTCGAGGCGATCGAGGCCTGTCTGGGCAAGACGGCGCTCCGCAATTCCCTGCCGATGCAAAAGGGTGATGTGCCCGCGACCttt                                  <  1:2086251/100‑1 (MQ=37)
                   gCGGCTGCTGGATTTCGTCGAGGCGATCGAGGCCTGTCTGGGCAAGACGGCGCTCCGCAATTACCTGCCGATGCAAAAGGGTGATGTGCCCGCGACCTTTg                                 <  2:4643623/101‑1 (MQ=35)
                   gCGGCTGCTGGATTTCGTCGAGGCGATCGAGGCCTGTCTGGGCAAGACGGCGCTCCGCAATTACCTGCCGATGCAAAAGGGTGATGTGCCCGCGACCTTTg                                 <  2:1657687/101‑1 (MQ=38)
                   gCGGCTGCTGGATTTCGTCGAGGCGATCGAGGCCTGTCTGGGCAAGACGGCGCTCCGCAATTACCTGCCGATGCAAAAGGGTGATGTGCCCGCGACCTTTg                                 >  1:4209053/1‑101 (MQ=35)
                   gCGGCTGCTGGATTTCGTCGAGGCGATCGAGGCCTGTCTGGGCAAGACGGCGCTCCGCAATTACCTGCCGATGCAAAAGGGTGATGTGCCCGCGACCTTTg                                 <  2:3173087/101‑1 (MQ=35)
                    cGGCTGCTGGATTTCGTCGAGGCGATCGAGGCCTGTCTGGGCAAGACGGCGCTCCGCAATTACCTGCCGATGCAAAAGGGTGATGTGCCCGCGACCTTTGc                                >  2:502799/1‑101 (MQ=35)
                     ggCTGCTGGATTTCGTCGAGGCGATCGAGGCCTGTCTGGGCAAGACGGCGCTCCGCAATTACCTGCCGATGCAAAAGGGTGATGTGCCCGCGACCTTTGcc                               >  2:1071540/1‑101 (MQ=35)
                      gctgctGGATTTCGTCGAGGGGATCGAGGCCTGTCTGGGCAAGACGGCGCTCCGCAATTACCTGCCGATGCAAAAGGGTGATGTGCCCGCGACCTTTGCCg                              <  1:4091685/101‑1 (MQ=25)
                      gctgctGGATTTCGTCGAGGCGATCGAGGCCTGTCTGGGCAAGACGGCGCTCCGCAATTACCTGCCGATGCAAAAGGGTGATGTGCCCGCGACCTTTGCCg                              >  1:968807/1‑101 (MQ=35)
                       ctgctgGATTTCGTCGAGGCGATCGAGGCCTGTCTGGGCAAGACGGCGCTCCGCAATTACCTGCCGATGCAAAAGGGTGATGTGCCCGCGACCTTTGCCGa                             <  1:3784270/101‑1 (MQ=38)
                       ctgctgGATTTCGTCGAGGCGATCGAGGCCTGTCTGGGCAAGACGGCGCTCCGCAATTACCTGCCGATGCAAAAGGGTGATGTGCCCGCGACCTTTGCCGa                             >  1:4171923/1‑101 (MQ=35)
                       ctgctgGATTTCGTCGAGGCGATCGAGGCCTGTCTGGGCAAGACGGCGCTCCGCAATTACCTGCCGATGCAAAAGGGTGATGTGCCCGCGACCTTTGCCGa                             >  2:3355234/1‑101 (MQ=35)
                         gctgGATTTCGTCGAGGCGATCGAGGCCTGTCTGGGCAAGACGGCGCTCCGCAATTACCTGCCGATGCAAAAGGGTGATGTGCCCGCGACCTTTGCCGATg                           <  1:2491651/101‑1 (MQ=38)
                         gctgGATTTCGTCGAGGCGATCGAGGCCTGTCTGGGCAAGACGGCGCTCCGCAATTACCTGCCGATGCAAAAGGGTGATGTGCCCGCGACCTTTGCCGATg                           >  1:3256545/1‑101 (MQ=35)
                         gctgGATTTCGTCGAGGCGATCGAGGCCTGTCTGGGCAAGACGGCGCTCCGCAATTACCTGCCGATGCAAAAGGGTGATGTGCCCGCGACCTTTGCCGATg                           <  1:3406145/101‑1 (MQ=38)
                          ctgGATTTCGTCGAGGCGATCGAGGCCTGTCTGGGCAAGACGGCGCTCCGCAATTACCTGCCGATGCAAAAGGGTGATGTGCCCGCGACCTTTGCCGATGc                          >  1:3330574/1‑101 (MQ=35)
                          ctgGATTTCGTCGAGGCGATCGAGGCCTGTCTGGGCAAGACGGCGCTCCGCAATTACCTGCCGATGCAAAAGGGTGATGTGCCCGCGACCTTTGCCGATGc                          <  1:120133/101‑1 (MQ=38)
                          ctgGATTTCGTCGAGGCGATCGAGGCCTGTCTGGGCAAGACGGCGCTCCGCAATTACCTGCCGATGCAAAAGGGTGATGTGCCCGCGACCTTTGCCGATGc                          <  1:4107185/101‑1 (MQ=38)
                            ggATTTCGTCGAGGCGATCGAGGCCTGTCTGGGCAAGACGGCGCTCCGCAATTACCTGCCGATGCAAAAGGGTGATGTGCCCGCGACCTTTGCCGATGCGa                        >  2:2179069/1‑101 (MQ=35)
                            ggATTTCGTCGAGGCGATCGAGGCCTGTCTGGGCAAGACGGCGCTCCGCAATTACCTGCCGATGCAAAAGGGTGATGTGCCCGCGACCTTTGCCGATGCGa                        >  2:4731251/1‑101 (MQ=35)
                            ggATTTCGTCGAGGCGATCGAGGCCTGTCTGGGCAAGACGGCGCTCCGCAATTACCTGCCGATGCAAAAGGGTGATGTGCCCGCGACCTTTGCCGATGCGa                        <  1:1569307/101‑1 (MQ=38)
                            ggATTTCGTCGAGGCGATCGAGGCCTGTCTGGGCAAGACGGCGCTCCGCAATTACCTGCCGATGCAAAAGGGTGATGTGCCCGCGACCTTTGCCGATGCGa                        >  1:580083/1‑101 (MQ=35)
                             gATTTCGTCGAGGCGATCGAGGCCTGTCTGGGCAAGACGGCGCTCCGCAATTACCTGCCGATGCAAAAGGGTGATGTGCCCGCGACCTTTGCCGATGCGAc                       <  2:2860853/101‑1 (MQ=38)
                                ttCGTCGAGGCGATCGAGGCCTGTCTGGGCAAGACGGCGCTCCGCAATTACCTGCCGATGCAAAAGGGTGATGTGCCCGCGACCTTTGCCGATGCGACgct                    >  1:668447/1‑101 (MQ=35)
                                ttCGTCGAGGCGATCGAGGCCTGTCTGGGCAAGACGGCGCTCCGCAATTACCTGCCGATGCAAAAGGGTGATGTGCCCGCGACCTTTGCCGATGCGACgct                    >  1:4632394/1‑101 (MQ=35)
                                 tcgtggGGGCGTGCGGGGCCTGGCGGGTCAAGACGGCCCTCCGCAATTACCTGCCGATGCAAAAGGGTGATGTGCCCGCGGCCTTTGCCGATGCGACgctg                   <  2:3602044/101‑1 (MQ=11)
                                 tcgtcgAGGCGATCGAGGCCTGTCTGGGCAAGACGGCGCTCCGCAATTACCTGCCGATGCAAAAGGGTGATGTGCCCGCGACCTTTGCCGATGCGACgctg                   <  1:3539712/101‑1 (MQ=35)
                                  cgtcgAGGCGATCGAGGCCTGTCTGGGCAAGACGGCGCTCCGCAATTACCTGCCGATGCAAAAGGGTGATGTGCCCGCGACCTTTGCCGATGCGACgctgc                  >  2:4466026/1‑101 (MQ=35)
                                  cgtcgAGGCGATCGAGGCCTGTCTGGGCAAGACGGCGCTCCGCAATTACCTGCCGATGCAAAAGGGTGATGTGCCCGCGACCTTTGCCGATGCGACgctgc                  <  2:3374523/101‑1 (MQ=35)
                                  cgtcgAGGCGATCGAGGCCTGTCTGGGCAAGACGGCGCTCCGCAATTACCTGCCGATGCAAAAGGGTGATGTGCCCGCGACCTTTGCCGATGCGACgctgc                  >  2:3116245/1‑101 (MQ=35)
                                   gtcgAGGCGATCGAGGCCTGTCTGGGCAAGACGGCGCTCCGCAATTACCTGCCGATGCAAAAGGGTGATGTGCCCGCGACCTTTGCCGATGCGACgctgct                 >  1:1686759/1‑101 (MQ=35)
                                     cgAGGCGATCGAGGCCTGTCTGGGCAAGACGGCGCTCCGCAATTACCTGCCGATGCAAAAGGGTGATGTGCCCGCGACCTTTGCCGATGCGACGCtgctgc               >  1:4890548/1‑101 (MQ=35)
                                      gAGGCGGTCGAGGCCTGTCTGGGCAAGACGGCGCTCCGCAATTACCTGCCGATGCAAAAGGGTGATGTGCCCGCGACCTTTGCCGATGCGACGCTGCTGca              <  2:3465479/101‑1 (MQ=25)
                                      gAGGCGATCGAGGCCTGTCTGGGCAAGACGGCGCTCCGCAATTACCTGCCGATGCAAAAGGGTGATGTGCCCGCGACCTTTGCCGATGCGACGCTGCTGca              >  1:1032827/1‑101 (MQ=35)
                                       aGGCGATCGAGGCCTGTCTGGGCAAGACGGCGCTCCGCAATTACCTGCCGATGCAAAAGGGTGATGTGCCCGCGACCTTTGCCGATGCGACGCTGCTGcac             <  1:1071653/101‑1 (MQ=35)
                                       aGGCGATCGAGGCCTGTCTGGGCAAGACGGCGCTCCGCAATTACCTGCCGATGCAAAAGGGTGATGTGCCCGCGACCTTTGCCGATGCGACGCTGCTGcac             <  2:487506/101‑1 (MQ=35)
                                       aGGCGATCGAGGCCTGTCTGGGCAAGACGGCGCTCCGCAATTACCTGCCGATGCAAAAGGGTGATGTGCCCGCGACCTTTGCCGATGCGACGCTGCTGcac             <  1:2134922/101‑1 (MQ=35)
                                        ggCGATCGAGGCCTGTCTGGGCAAGGCGGCGCTCCGCAAGTACCTGCCGATGCAAAAGGGTGATGTGCCCGCGACCTTTGCCGATGCGACGCtgctgctcc            >  2:4189596/1‑98 (MQ=255)
                                        ggCGATCGAGGCCTGTCTGGGCAAGACGGCGCTCCGCAATTACCTGCCGATGCAAAAGGGTGATGTGCCCGCGACCTTTGCCGATGCGACGCTGCTGcaca            >  2:1595520/1‑101 (MQ=35)
                                         gCGATCGAGCCCTGTCTGGGCAAGACGGCGCTCCGCAATTACCTGCCGATGCAAAAGGGTGATGTGCCCGCGACCTTTGCCGATGCGACGCTGCTGCACAg           <  2:2020352/101‑1 (MQ=25)
                                          cGATCGAGGCCTGTCTGGGCAAGACGGCGCTCCGCAATTACCTGCCGATGCAAAAGGGTGATGTGCCCGCGACCTTTGCCGATGCGACGCTGCTGCACAGc          <  1:537223/101‑1 (MQ=35)
                                          cGAGCGAGGCCTGTCTGGGCAAGCCGGCGCTCTGCAATTACCTGCCGGTGCAAAAGCGTGATGTGCCCGCGACCTTTGCCGCTGCGACGCTGCTGCACAGc          >  2:820660/1‑101 (MQ=11)
                                           gagcGAGGCCTGTCTGGGCAAGACGGCGCTCCGCAATTACCTGCCGATGCAAAAGGGTGATGTGCCCGCGACCTTTGCCGATGCGACGCTGCTGCACAGcc         <  2:162496/98‑1 (MQ=25)
                                           gATCGAGGCCTGTCTGGGCAAGACGGCGCTCCGCAATTACCTGCCGATGCAAAAGGGTGATGTGCCCGCGACCTTTGCCGATGCGACGCTGCTGCACAGcc         >  1:5034751/1‑101 (MQ=35)
                                           gATCGAGGCCTGTCTGGGCAAGACGGCGCTCCGCAATTACCTGCCGATGCAAAAGGGTGATGTGCCCGCGACCTTTGCCGATGCGACGCTGCTGCACAGcc         >  1:3227070/1‑101 (MQ=35)
                                           gATCGAGGCCTGTCTGGGCAAGACGGCGCTCCGCAATTACCTGCCGATGCAAAAGGGTGATGTGCCCGCGACCTTTGCCGATGCGACGCTGCTGCACAGcc         <  2:3330311/101‑1 (MQ=35)
                                           gATCGAGGCCTGTCGGGGCAAGACGGCGCTCCGCAATTCCCTGCCGATGCAAAAGGGTGATGTGCCCGCGACCTTTGCCGATGCGACGCTGCTGCACAGcc         <  2:4689359/101‑1 (MQ=255)
                                            aTCGAGGCCTGTCTGGGCAAGGCGGCGCTCCGCAATTACCTGCCGATGCAAAAGGGTGATGTGCCCGCGACCTTTGCCGATGCGACGCTGCTGCACAGCCt        >  1:2870726/1‑101 (MQ=25)
                                            aTCGAGGCCTGTCTGGGCAAGACGGCGCCCCGCAATTACCTGCCGATGCAAAAGGGTGATGTGCCCGCGACCTTTGCCGATGCGACGCTGCTGCGCAGCCt        <  1:3914429/101‑1 (MQ=255)
                                             tggAGGCCTGTCTGGGCAAGACGGCGCTCCGCAATTACCTGCCGATGCAAAAGGGTGATGTGCCCGCGACCTTTGCCGATGCGACGCTGCTGCACAGCCTg       <  2:4905770/99‑1 (MQ=25)
                                             tCGAGGCCTGTCTGGGCAAGACGGCGCTCCGCAATTACCTGCCGATGCAAAAGGGTGATGTGCCCGCGACCTTTGCCGATGCGACGCTGCTGCACAGCCTg       >  1:4273587/1‑101 (MQ=35)
                                             tCGAGGCCTGTCTGGGCAAGACGGCGCTCCGCAATTACCTGCCGATGCAAAAGGGTGATGTGCCCGCGACCTTTGCCGATGCGACGCTGCTGCACAGCCTg       <  1:2548906/101‑1 (MQ=35)
                                             tCGAGGCCTGTCTGGGCAAGACGGCGCTCCGCAATTACCTGCCGATGCAAAAGGGTGATGTGCCCGCGACCTTTGCCGATGCGACGCTGCTGCACAGCCTg       <  2:3751701/101‑1 (MQ=35)
                                              cGAGGCCTGTCTGGGCAAGACGGCGCTCCGCAATTACCTGCCGATGCAAAAGGGTGATGTGCCCGCGACCTTTGCCGATGCGACGCTGCTGCACAGCCTGa      >  2:289413/1‑101 (MQ=35)
                                              cGAGGCCTGTCTGGGCAAGACGGCGCTCCGCAATTACCTGCCGATGCAAAAGGGTGATGTGCCCGCGACCTTTGCCGATGCGACGCTGCTGCACAGCCTGa      <  2:729998/101‑1 (MQ=35)
                                               gAGGCCTGTCTGGGCAAGACGGCGCTCCGCAATTACCTGCCGATGCAAAAGGGTGATGTGTCCGCGACCTTTGCCGATGCGACGCTGCTGCACAGCCTGAc     >  1:2586965/1‑101 (MQ=25)
                                               gAGGCCTGTCTGGGCAAGACGGCGCTCCGCAATTACCTGCCGATGCAAAAGGGTGATGTGCCCGCGACCTTTGCCGATGCGACGCTGCTGCACAGCCTGAc     >  2:1382480/1‑101 (MQ=35)
                                               gAGGCCTGTCGGGGCAAGACGGCTCTCCGCAATTACCTGCCGATGCAAAAGGGTGATGTGCCCGCGACCTTTGCCGATGCGACGCTGCTGCACAGCCTGAc     <  1:3233301/101‑1 (MQ=255)
                                               gAGGCCTGTCGGGGCAAGACGGCGCTCCGCAATTACCTGCCGATGCAAAAGGGTGATGTGCCCGCGACCTTTGCCGATGCGACGCTGCTGCACAGCCTGAc     <  2:739981/101‑1 (MQ=25)
                                                gggCCTGTCTGGGCAAGACGGCGCTCCGCAATTCCCTGCCGATGCAAAAGGGTGATGTGCCCGCGACCTTTGCCGATGCGACGCTGCTGCACAGCCTGAcc    <  1:2312488/100‑1 (MQ=25)
                                                  gCCTGTCTGGGCAAGACGGCGCTCCGCAATTACCTGCCGATGCAAAAGGGTGATGTGCCCGCGACCTTTGCCGATGCGACGCTGCTGCACAGCCTGACCgg  >  1:2549028/1‑101 (MQ=35)
                                                  gCCTGTCTGGGCAAGACGGCGCTCCGCAATTACCTGCCGATGCAAAAGGGTGATGTGCCCGCGACCTTTGCCGATGCGACGCTGCTGCACAGCCTGACCgg  <  1:2011807/101‑1 (MQ=35)
                                                  |                                                                                                    
TCGGCAATTCCGACAAGGTGCGGCTGCTGGATTTCGTCGAGGCGATCGAGGCCTGTCTGGGCAAGACGGCGCTCCGCAATTACCTGCCGATGCAAAAGGGTGATGTGCCCGCGACCTTTGCCGATGCGACGCTGCTGCACAGCCTGACCGG  >  NC_014034/1154475‑1154625

Base quality scores:  ATCG  < 3 ≤  ATCG  < 17 ≤  ATCG  < 25 ≤  ATCG  < 33 ≤  ATCG  < 40 ≤  ATCG