New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_014034 158228 =76 (0.330)59 (0.260) 28/194 16.1 0.51 intergenic (+31/‑199) parC/tuf1 DNA topoisomerase IV subunit A/translation elongation factor Tu
?NC_014034 = 2330172 39 (0.160)intergenic (+147/+165) RCAP_rcc02161/mrcA hypothetical protein/penicillin‑binding protein 1A
Rejected: Position hash score below cutoff.

CTGGGCAAGCGGGCGAGCGCCGGGCGGATGGCGCCGCGCGGTTTCCCGCGCGACAACCGGTTCCCCTGACGGGGTTTCGGGCGCTGCCGCGCCCGACCCGCGGGGGCTCCGCCCCCGCACCCCCGGCGTATTTTGCCAAGGAAAAGCAAGGTGCCCCCCGGACGTGTCGTTTCGGGGACGCGCGCTGCCCTGGCAGGCGG  >  NC_014034/158129‑158328
                                                                                                   |                                                                                                    
cTGGGCAAGCGGGCGAGCGCCGGGCGGATGGCGCCGCGCGGTTTCCCGCGCGACAACCGGTTCCCCTGACGGGGTTTCGGGCGCTGCCGCGCCCGACCcgc                                                                                                     <  1:659415/101‑1 (MQ=255)
   ggCAAGCGGGCGAGCGCCGGGCGGATGGCGCCGCGCGGTTTCCCGCGCGACAACCGGTTCCCCTGACGGGGTTTCGGGCGCTGCCGCGCCCGACCcgccgc                                                                                                  >  1:231352/1‑98 (MQ=255)
       aGCGGGCGAGCGCCGGGCGGATGGCGCCGCGCGGTTTCCCGCGCGACAACCGGTTCCCCTGACGCGGTTTCGGGCGCTGGCGCGCCCGACCcggcgtgcgc                                                                                              >  2:4615604/1‑93 (MQ=255)
         ggggCGAGCGCCGGGCGGATGGCGCCGCGCGGTTTCCCGCGCGACAACCGGTTCCCCTGACGGGGTTTCGGGCGCTGCCGCGCCCGACCCGCGGGGGCTcc                                                                                            <  1:2654591/100‑1 (MQ=255)
            gCGAGCGCCGGGCGGATGGCGCCGCGCGGTTTCCCGCGCGACAACCGGTTCCCCTGACGGGGTTTCGGGCGCTGCCGCGCCCGACCCGCGGGGGCTCCGcc                                                                                         <  1:1660004/101‑1 (MQ=255)
             cGAGCGCCGGGCGGATGGCGCCGCGCGGTTTCCCGCGCGACAACCGGTTCCCCTGACGGGGTTTCGGGCGCTGCCGCGCCCGACCCGCGGGGGCTCCGccc                                                                                        <  1:1802944/101‑1 (MQ=255)
              gAGCGCCGGGCGGATGGCGCCGCGCGTTTTCCCGCGCAACATCCGGTCACCCTGACGGGTTTTCGGGCGCTGCCCAGCCCGACCCGCGGGGGCTCCGcccc                                                                                       <  2:4223502/101‑1 (MQ=255)
               aGCGCCGGGCGGATGGCGCCGCGCGGTTTCCCGCGCGACAACCGGTTCCCCTGACGGGGTTTCGGGCGCTGCCGCGCCCGACCCGCGGGGGCTCCGccccc                                                                                      <  1:2627817/101‑1 (MQ=255)
               aGCGCCGGGCGGATGGCGCCGCGCGGTTTCCCGCGCGACAACCGGTTCCCCTGACGGGGTTTCGGGCGCTGCCGCGCCCGACCCGCGGGGGCTCCGccccc                                                                                      <  1:2462998/101‑1 (MQ=255)
               aGCGCCGGGCGGATGGCGCCGCGCGGTTTCCCGCGCGACAACCGGTTCCCCTGACGGGGTTTCGGGCGCTGCCGCGCCCGACCCGCGGGGGCTCCGccccc                                                                                      <  2:2934903/101‑1 (MQ=255)
                  gcCGGGCGGATGGCGCCGCGCGGTTTCCCGCGCGACAACCGGTTCCCCTGACGGGGTTTCGGGCGCTGCCGCGCCCGACCCGCGGGGGCTCCGCCCCCGCa                                                                                   <  1:463631/101‑1 (MQ=255)
                   ccGGGCGGATGGCGCCGCGCGGTTTCCCGCGCGACAACCGGTTCCCCTGACGGGGTTTCGGGCGCTGCCGCGCCCGACCCGCGGGGGCTCCGCCCCCGCAc                                                                                  <  1:4026073/101‑1 (MQ=255)
                   ccGGGCGGATGGCGCCGCGCGGTTTCCCGCGCGACAACCGGTTCCCCTGACGGGGTTTCGGGCGCTGCCGCGCCCGACCCGCGGGGGCTCCGCCCCCGCAc                                                                                  <  2:4156326/101‑1 (MQ=255)
                   ccGGGCGGATGGCGCCGCGCGGTTTCCCGCGCGACAACCGGTTCCCCTGACGGGGTTTCGGGCGCTGCCGCGCCCGACCCGCGGGGGCTCCGCCCCCGCAc                                                                                  <  2:4586447/101‑1 (MQ=255)
                     gggCGGGTGGCGCCGCGCGGTTTCCCGCGCGACAACCGGTTCCCCTGACGGGGTTTCGGGCGCTGCCGCGCCCGACCCGCGGGGGCTCCGCCCCCGCAccc                                                                                <  1:3528419/101‑1 (MQ=255)
                     gggCGGATGGCGCCGCGCGGTTTCCCGCGCGACAACCGGTTCCCCTGACGGGGTTTCGGGCGCTGCCGCGCCCGACCCGCGGGGGCTCCGCCCCCGCAccc                                                                                <  2:2654883/101‑1 (MQ=255)
                      ggCGGATGGCGCCGCGCGGTTTCCCGCGCGACAACCGGTTCCCCTGACGGGGTTTCGGGCGCTGCCGCGCCCGACCCGCGGGGGCTCCGCCCCCGCAcccc                                                                               <  2:4477783/101‑1 (MQ=255)
                      ggCGGATGGCGCCGCGCGGTTTCCCGCGCGACAACCGGTTCCCCTGACGGGGTTTCGGGCGCTGCCGCGCCCGACCCGCGGGGGCTCCGCCCCCGCAcccc                                                                               <  2:1997999/101‑1 (MQ=255)
                       gCGGATGGCGCCGCGCGGTTTCCCGCGCGACAACCGGTTCCCCTGACGGGGTTTCGGGCGCTGCCGCGCCCGACCCGCGGGGGCTCCGCCCCCGCAccccc                                                                              <  2:765397/101‑1 (MQ=255)
                        cGGATGGCGCCGCGCGGTTTCCCGCGCGACAACCGGTTCCCCTGACGGGGTTTCGGGCGCTGCCGCGCCCGACCCGCGGGGGCTCCGCCCCCGCACCCCCg                                                                             <  2:2654471/101‑1 (MQ=255)
                        cGGATGGCGCCGCGCGGTTTCCCGCGCGACAACCGGTTCCCCTGACGGGGTTTCGGGCGCTGCCGCGCCCGACCCGCGGGGGCTCCGCCCCCGCACCCCCg                                                                             <  1:2256001/101‑1 (MQ=255)
                        cGGATGGCGCCGCGCGGTTTCCCGCGCGACAACCGGTTCCCCTGACGGGGTTTCGGGCGCTGCCGCGCCCGACCCGCGGGGGCTCCGCCCCCGCACCCCCg                                                                             <  1:1051429/101‑1 (MQ=255)
                         ggATGGCGGCGGGCCGTTTCCCGCGCGGCAACCGGTTCCCCGGACGGGGTTTCGGGCGCTGCCGCGCCCGGCCCGCGGGGGCTCCGCCCCCGCACCCCCgg                                                                            <  1:547552/101‑1 (MQ=255)
                         ggATGGCGCCGCGCGGTTTCCCGCGCGACAACCGGTTCCCCTGACGGGGTTTCGGGCGCTGCCGCGCCCGACCCGCGGGGGCTCCGCCCCCGCACCCCCgg                                                                            <  1:1365608/101‑1 (MQ=255)
                          gATGGCGCCGCGCGGTTTCCCGCGCGACAACCGGTTCCCCTGACGGGGTTTCGGGCGCTGCCGCGCCCGACCCGCGGGGGCTCCGCCCCCGCACCCCCGGc                                                                           <  2:1552987/101‑1 (MQ=255)
                          gATGGCGCCGCGCGGTTTCCCGCGCGACAACCGGTTCCCCTGACGGGGTTTCGGGCGCTGCCGCGCCCGACCCGCGGGGGCTCCGCCCCCGCACCCCCGGc                                                                           <  2:2754956/101‑1 (MQ=255)
                          gATGGCGCCGCGCGGTTTCCCGCGCGACAACCGGTTCCCCTGACGGGGTTTCGGGCGCTGCCGCGCCCGACCCGCGGGGGCTCCGCCCCCGCACCCCCGGc                                                                           <  2:3347737/101‑1 (MQ=255)
                            tGGCGCCGCGCGGTTTCCCGCGCGACAACCGGTTCCCCTGACGGGGTTTCGGGCGCTGCCGCGCCCGACCCGCGGGGGCTCCGCCCCCGCACCCCCGGCGt                                                                         <  1:3181840/101‑1 (MQ=255)
                            tGGCGCCGCGCGGTTTCCCGCGCGACAACCGGTTCCCCTGACGGGGTTTCGGGCGCTGCCGCGCCCGACCCGCGGGGGCTCCGCCCCCCCACCCCCGGCGt                                                                         <  1:2576432/101‑1 (MQ=255)
                              gcgcCGCGCGGTTTCCCGCGCGACAACCGGTTCCCCTGACGGGGTTTCGGGCGCTGCCGCGCCCGACCCGCGGGGGCTCCGCCCCCGCACCCCCGGCGTAt                                                                       <  1:950213/101‑1 (MQ=255)
                              gcgcCGCGCGGTTTCCCGCGCGACAACCGGTTCCCCTGACGGGGTTTCGGGCGCTGCCGCGCCCGACCCGCGGGGGCTCCGCCCCCGCACCCCCGGCGTAt                                                                       <  2:3433452/101‑1 (MQ=255)
                               cgccgcGCGGTTTCCCGCGCGACAACCGGTTCCCCTGACGGGGTTTCGGGCGCTGCCGCGCCCGACCCGCGGGGGCTCCGCCCCCGCACCCCCGGCGTAtt                                                                      <  1:4238081/101‑1 (MQ=255)
                                gccgcGCGGTTTCCCGCGCGACAACCGGTTCCCCTGACGGGGTTTCGGGCGCTGCCGCGCCCGACCCGCGGGGGCTCCGCCCCCGCACCCCCGGCGTAttt                                                                     <  1:1459084/101‑1 (MQ=255)
                                gccgcGCGGTTTCCCGCGCGACAACCGGTTCCCCTGACGGGGTTTCGGGCGCTGCCGCGCCCGACCCGCGGGGGCTCCGCCCCCGCACCCCCGGCGTAttt                                                                     <  1:2798091/101‑1 (MQ=255)
                                 ccgcGCGGTTTCCCGCGCGACAACCGGTTCCCCTGACGGGGTTTCGGGCGCTGCCGCGCCCGACCCGCGGGGGCTCCGCCCCCGCACCCCCGGCGTAtttt                                                                    <  1:3702962/101‑1 (MQ=255)
                                 ccgcGCGGTTTCCCGCGCGACAACCGGTTCCCCTGACGGGGTTTCGGGCGCTGCCGCGCCCGACCCGCGGGGGCTCCGCCCCCGCACCCCCGGCGTAtttt                                                                    <  1:503150/101‑1 (MQ=255)
                                  cgcgcgGTTTCCCGCGCGACAACCGGTTCCCCTGACGGGGTTTCGGGCGCTGCCGCGCCCGACCCGCGGGGGCTCCGCCCCCGCACCCCCGGCGTATTTTg                                                                   <  1:379331/101‑1 (MQ=255)
                                      cgGTTTCCCGCGCGACAACCGGTTCCCCTGACGGGGTTTCGGGCGCTGCCGCGCCCGACCCGCGGGGGCTCCGCCCCCGCACCCCCGGCGTATTTTGCCaa                                                               <  2:4224515/101‑1 (MQ=255)
                                          ttCCCGCGCGACAACCGGTTCCCCTGACGGGGTTTCGGGCGCTGCCGCGCCCGACCCGCGGGGGCTCCGCCCCCGCACCCCCGGCGTATTTTGCCAAGGaa                                                           <  1:4251563/101‑1 (MQ=255)
                                          ttCCCGCGCGACAACCGGTTCCCCTGACGGGGTTTCGGGCGCTGCCGCGCCCGACCCGCGGGGGCTCCGCCCCCGCACCCCCGGCGTATTTTGCCAAGGaa                                                           <  1:700260/101‑1 (MQ=255)
                                            cccGCGCGACAACCGGTTCCCCTGACGGGGTTTCGGGCGCTGCCGCGCCCGACCCGCGGGGGCTCCGCCCCCGCACCCCCGGCGTATTTTGCCAAGGaaaa                                                         <  1:22254/101‑1 (MQ=255)
                                             ccGCGCGACAGCCGGTTNCCCTGACGGGGTTTCGGGCGCTGCCGCGCCCGACCCGCGGGGGCTCCGCCCCCGCACCCCCGGCGTATTTTGCCAAGGAAAAg                                                        <  2:509984/101‑1 (MQ=255)
                                              cgcgcgACAACCGGTTCCCCTGACGGGGTTTCGGGCGCTGCCGCGCCCGACCCGCGGGGGCTCCGCCCCCGCACCCCCGGCGTATTTTGCCAAGGAAAAGc                                                       <  2:948045/101‑1 (MQ=255)
                                               gcgcgACAACCGGTTCCCCTGACGGGGTTTCGGGCGCTGCCGCGCCCGACCCGCGGGGGCTCCGCCCCCGCACCCCCGGCGTATTTTGCCAAGGAAAAGCa                                                      <  1:809882/101‑1 (MQ=255)
                                               gcgcgACAACCGGTTCCCCTGACGGGGTTTCGGGCGCTGCCGCGCCCGACCCGCGGGGGCTCCGCCCCCGCACCCCCGGCGTATTTTGCCAAGGAAAAGCa                                                      <  1:33822/101‑1 (MQ=255)
                                                cgcgACAACCGGTTCCCCTGACGGGGTTTCGGGCGCTGCCGCGCCCGACCCGCGGGGGCTCCGCCCCCGCACCCCCGGCGTATTTTGCCAAGGAAAAGCaa                                                     <  2:1921979/101‑1 (MQ=255)
                                                cgcgACAACCGGTTCCCCTGACGGGGTTTCGGGCGCTGCCGCGCCCGACCCGCGGGGGCTCCGCCCCCGCACCCCCGGCGTATTTTGCCAAGGAAAAGCaa                                                     <  1:4344989/101‑1 (MQ=255)
                                                 gcgACAACCGGTTCCCCTGACGGGGTTTCGGGCGCTGCCGCGCCCGACCCGCGGGGGCTCCGCCCCCGCACCCCCGGCGTATTTTGCCAAGGAAAAGCAAg                                                    <  2:1742443/101‑1 (MQ=255)
                                                  cgACAACCGGTTCCCCTGACGGGGTTTCGGGCGCTGCCGCGCCCGACCCGCGGGGGCTCCGCCCCCGCACCCCCGGCGTATTTTGCCAAGGAAAAGCAAgg                                                   <  2:2643563/101‑1 (MQ=255)
                                                  cgACAACCGGTTCCCCTGACGGGGTTTCGGGCGCTGCCGCGCCCGACCCGCGGGGGCTCCGCCCCCGCACCCCCGGCGTATTTTGCCAAGGAAAAGCAAgg                                                   <  1:1726743/101‑1 (MQ=255)
                                                       aCCGGTTCCCCTGACGGGGTTTCGGGCGCTGCCGCGCCCGACCCGCGGGGGCTCCGCCCCCGCACCCCCGGCGTATTTTGCCAAGGAAAAGCAAGGTGccc                                              <  1:1920774/101‑1 (MQ=255)
                                                        ccGGTTCCCCTGACGGGGTTTCGGGCGCTGCCGCGCCCGACCCGCGGGGGCTCCGCCCCCGCACCCCCGGCGTATTTTGCCAAGGAAAAGCAAGGTGcccc                                             <  1:2764652/101‑1 (MQ=255)
                                                         cGGTTCCCCTGGCGGGTTTTCGGGCGCTCCCGCGCCCGACCCGCGGGGGCTCCGCCCCCGCCCCCCCGGCGTATTTTGCCAAGGAAAAGCAAGGTGccccc                                            <  1:1876604/101‑1 (MQ=37)
                                                         cGGTTCCCCTGACGGGGTTTCGGGCGCTGCCGCGCCCGACCCGCGGGGGCTCCGCCCCCGCACCCCCGGCGTATTTTGCCAAGGAAAAGCAAGGTGccccc                                            <  1:1793819/101‑1 (MQ=255)
                                                         cGGTTCCCCTGACGGGGTTTCGGGCGCTGCCGCGCCCGACCCGCGGGGGCTCCGCCCCCGCACCCCCGGCGTATTTTGCCAAGGAAAAGCAAGGTGccccc                                            <  1:945910/101‑1 (MQ=255)
                                                         cGGTTCCCCTGACGGGGTTTCGGGCGCTGCCGCGCCCGACCCGCGGGGGCTCCGCCCCCGCACCCCCGGCGTATTTTGCCAAGGAAAAGCAAGGTGccccc                                            <  2:3158768/101‑1 (MQ=255)
                                                          gggtCCCCTGACGGGGTTTCGGGCGCTGCCGCGCCCGACCCGCGGGGGCTCCGCCCCCGCACCCCCGGCGTATTTTGCCAAGGAAAAGCAAGGTGcccccc                                           <  2:624910/98‑1 (MQ=255)
                                                            ttCCCCTGACGGGGTTTCGGGCGCTGCCGCGCCCGACCCGCGGGGGCTCCGCCCCCGCACCCCCGGCGTATTTTGCCAAGGAAAAGCAAGGTGCCCCCCgg                                         <  1:576466/101‑1 (MQ=255)
                                                            ttCCCCTGACGGGGTTTCGGGCGCTGCCGCGCCCGACCCGCGGGGGCTCCGCCCCCGCACCCCCGGCGTATTTTGCCAAGGAAAAGCAAGGTGCCCCCCgg                                         <  1:509977/101‑1 (MQ=255)
                                                             tCCCCTGACGGGGTTTCGGGCGCTGCCGCGCCCGACCCGCGGGGGCTCCGCCCCCGCACCCCCGGCGTATTTTGCCAAGGAAAAGCAAGGTGCCCCCCGGa                                        <  1:4350608/101‑1 (MQ=255)
                                                              ccccTGACGGGGTTTCGGGCGCTGCCGCGCCCGACCCGCGGGGGCTCCGCCCCCGCACCCCCGGCGTATTTTGCCAAGGAAAAGCAAGGTGCCCCCCGGAc                                       <  2:4032393/101‑1 (MQ=255)
                                                              ccccTGACGGGGTTTCGGGCGCTGCCGCGCCCGACCCGCGGGGGCTCCGCCCCCGCACCCCCGGCGTATTTTGCCAAGGAAAAGCAAGGTGCCCCCCGGAc                                       <  1:2812329/101‑1 (MQ=255)
                                                                 cTGACGGGGTTTCGGGCCCTGCCGCGCCCGACCCGCGGGGGCTCCGCCCCCGCACCCCCGGCGTATTTTGCCAAGGAAAAGCAAGGTGCCCCCCGGACgtg                                    <  2:1066601/101‑1 (MQ=255)
                                                                  tGACGGGGTTTCGGGCGCTGCCGCGCCCGCCCCGCGGGGGCTCCGCCCCCGCACCCCCGGCGTATTTTGCCAAGGAAAAGCAAGGTGCCCCCCGGACgtgt                                   <  1:3613726/101‑1 (MQ=255)
                                                                      ggggTTTCGGGCGCTGCCGCGCCCGACCCGCGGGGGCTCCGCCCCCGCACCCCCGGCGTATTTTGCCAAGGAAAAGCAAGGTGCCCCCCGGACGTGTCGtt                               <  2:1740449/101‑1 (MQ=255)
                                                                       gggTTTCGGGCGCTGCCGCGCCCGACCCGCGGGGGCTCCGCCCCCGCACCCCCGGCGTATTTTGCCAAGGAAAAGCAAGGTGCCCCCCGGACGTGTCGttt                              <  2:955504/101‑1 (MQ=255)
                                                                       gggTTTCGGGCGCTGCCGCGCCCGACCCGCGGGGGCTCCGCCCCCGCACCCCCGGCGTATTTTGCCAAGGAAAAGCAAGGTGCCCCCCGGACGTGTCGttt                              <  2:2386421/101‑1 (MQ=255)
                                                                        ggTTTCGGGCGCTGCCGCGCCCGACCCGCGGGGGCTCCGCCCCCGCACCCCCGGCGTATTTTGCCAAGGAAAAGCAAGGTGCCCCCCGGACGTGTCGTTTc                             <  2:106191/101‑1 (MQ=255)
                                                                            tCGGGCGCTGCCGCGCCCGACCCGCGGGGGCTCCGCCCCCGCACCCCCGGCGTATTTTGCCAAGGAAAAGCAAGGTGCCCCCCGGACGTGTCGTTTCgggg                         <  2:754308/101‑1 (MQ=255)
                                                                              gggCGCTGCCGCGCCCGACCCGCGGGGGCTCCGCCCCCGCACCCCCGGCGTATTTTGCCAAGGAAAAGCAAGGTGCCCCCCGGACGTGTCGTTTCGGGGAc                       <  1:3484078/101‑1 (MQ=255)
                                                                                  gcTGCCGCGCCCGACCCGCGGGGGCTCCGCCCCCGCACCCCCGGCGTATTTTGCCAAGGAAAAGCAAGGTGCCCCCCGGACGTGTCGTTTCGGGGAcgcgc                   <  1:4171648/101‑1 (MQ=255)
                                                                                     ccgccgcCCTGCCCGCGGGTGCTCCGCCCCCGCACCCCCGGCGTATTTTGCCAAGGAAAAGCAAGGTGCCCCCCTGACGTGTCGTTTCGGGGACGCGCGCt                <  2:4339759/97‑1 (MQ=37)
                                                                                         cgcCCGACCCGCGGGGGCTCCGCCCCCGCACCCCCGGCGTATTTTGCCAAGGAAAAGCAAGGTGCCCCCCGGACGTGTCGTTTCGGGGACGCGCGCTGccc            <  1:4421610/101‑1 (MQ=255)
                                                                                           cccGACCCGCGGGGGCTCCGCCCCCGCACCCCCGGCGTATTTTGCCAAGGAAAAGCAAGGTGCCCCCCGGACGTGTCGTTTCGGGGACGCGCGCTGCCCTg          <  2:871659/101‑1 (MQ=255)
                                                                                            ccGACCCGCGGGGGCTCCGCCCCCGCACCCCCGGCGTATTTTGCCAAGGAAAAGCAAGGTGCCCCCCGGACGTGTCGTTTCGGGGACGCGCGCTGCCCTgg         <  1:2851218/101‑1 (MQ=255)
                                                                                             cGACCCGCGGGGGCTCCGCCCCCGCACCCCCGGCGTATTTTGCCAAGGAAAAGCAAGGTGCCCCCCGGACGTGTCGTTTCGGGGACGCGCGCTGCCCTGGc        <  2:2829242/101‑1 (MQ=255)
                                                                                              gccccGCGGGGCCCCCCCCCCCCCCCCCCCGGCGTATTTTGCCAAGGAAAAGCAAGGTCCCCCCCGGACGTGTCGTTTCGGGGACGCGCGCTGCCCTGGCa       <  2:1961780/99‑1 (MQ=25)
                                                                                                 ccGCGGGGGCTCCGCCCCCGCACCCCCGGCGTATTTTGCCAAGGAAAAGCAAGGTGCCCCCCGGACGTGTCGTTTCGGGGACGCGCGCTGCCCTGGCAggc    <  1:3262342/101‑1 (MQ=255)
                                                                                                  cgcgGGGGCTCCGCCCCCGCACCCCCGGCGTATTTTGCCAAGGAAAAGCAAGGTGCCCCCCGGACGTGTCGTTTCGGGGACGCGCGCTGCCCTGGCAggcg   <  1:153196/101‑1 (MQ=255)
                                                                                                   gcgGGGGCTCCGCCCCCGCACCCCCGGCGTATTTTGCCAAGGAAAAGCAAGGTGCCCCCCGGACGTGTCGTTTCGGGGACGCGCGCTGCCCTGGCAggcgg  <  1:4615857/101‑1 (MQ=255)
                                                                                                   |                                                                                                    
CTGGGCAAGCGGGCGAGCGCCGGGCGGATGGCGCCGCGCGGTTTCCCGCGCGACAACCGGTTCCCCTGACGGGGTTTCGGGCGCTGCCGCGCCCGACCCGCGGGGGCTCCGCCCCCGCACCCCCGGCGTATTTTGCCAAGGAAAAGCAAGGTGCCCCCCGGACGTGTCGTTTCGGGGACGCGCGCTGCCCTGGCAGGCGG  >  NC_014034/158129‑158328

Base quality scores:  ATCG  < 3 ≤  ATCG  < 18 ≤  ATCG  < 24 ≤  ATCG  < 33 ≤  ATCG  < 39 ≤  ATCG