New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_014034 = 770327214 (0.970)2 (0.010) 2/186 37.3 0.01 coding (159/552 nt) bsaA1 glutathione peroxidase
?NC_014034 = 2811345 224 (0.940)coding (1180/1650 nt) sulP sulphate transporter
Rejected: Position hash score below cutoff.

Only 100 of 226 total aligned reads displayed.

CCGCCGCCCGCGCGACCGAGGCGGCGGCGCCGAGCTTCACCTTCGCCTCGATCGAGGGCGGGGTGATCGACACCGCCGCCTATCGCGGTCATCCGGTGCTGGTCGTGAACACCGCGTCGCTCTGCGGCTTCACCCCGCAGCTCGAAGGGCTGCAGGCGCTGCACAAGGCCTGGGGGCCGAAGGGGCTGATCGTGCTGGCCG  >  NC_014034/770227‑770427
                                                                                                    |                                                                                                    
ccggcgCCCGCACGACCCAGGCAGCTGCAAGAAGAGTCAACTTCGCCTGGATCGTGGCTGCGGTGATCGACACCGCCGCCTATCGCGATCAAGCGGTTCTg                                                                                                      <  2:813696/101‑1 (MQ=255)
ccgccgCCCGCGCGACCGAGGCGGCGGCGCCGAGCTTCACCTTCGCCTCGATCGAGGGCGGGGTGATCGACACCGCCGCCTATCGCGGTCATCCGGTGCTg                                                                                                      <  1:4205421/101‑1 (MQ=255)
ccgccgCCCGCGCGACCGAGGCGGCGGCGCCGAGCTTCACCTTCGCCTCGATCGAGGGCGGGGTGATCGACACCGCCGCCTATCGCGGTCATCCGGTGCTg                                                                                                      <  2:1619726/101‑1 (MQ=255)
 cgccgcCCGCGCGACCGAGGCGGCGGCGCCGAGCTTCACCTTCGCCTCGATCGAGGGCGGGGTGATCGACACCGCCGCCTATCGCGGTCATCCGGTGCTgg                                                                                                     <  1:1172053/101‑1 (MQ=255)
 cgccgcCCGCGCGACCGAGGCGGCGGCGCCGAGCTTCACCTTCGCCTCGATCGAGGGCGGGGTGATCGACACCGCCGCCTATCGCGGTCATCCGGCGCTgg                                                                                                     >  1:3559626/1‑101 (MQ=255)
   ccgcccgcGCGACCGAGGCGGCGGCGCCGAGCTTCACCTTCGCCTCGATCGAGGGCGGGGTGATCGACACCGCCGCCTATCGCGGTCATCCGGTGCTGGTc                                                                                                   <  2:1281235/101‑1 (MQ=255)
    cgcccgcGCGCCCGAGGCGGCGGCGCCGAGCTTCACCTTCGCCTCGATCGAGGGCGGGGTGATCGACACCGCCGCCTATCGCGGTCATCCGGTGCTGGTCg                                                                                                  <  2:2106689/101‑1 (MQ=255)
    cgcccgcGCGACCGAGGCGGCGGCGCCGAGCTTCACCTTCGCCTCGATCGAGGGCGGGGTGATCGACACCGCCGCCTATCGCGGTCATCCGGTGCTGGTCg                                                                                                  <  2:1248034/101‑1 (MQ=255)
         gcgcgACCGAGGCGGCGGCGCCGAGCTTCACCTTCGCCTCGATCGAGGGCGGGGTGATCGACACCGCCGCCTATCGCGGTCATCCGGTGCTGGTCGTGAac                                                                                             <  2:4345205/101‑1 (MQ=255)
          cgcgACCGAGGCGGCGGCGCCGAGCTTCACCTTCGCCTCGATCGAGGGCGGGGTGATCGACACCGCCGCCTATCGCGGTCATCCGGTGCTGGTCGTGAaca                                                                                            <  2:104458/101‑1 (MQ=255)
          cgcgACCGAGGCGGCGGCGCCGAGCTTCACCTTCGCCTCGATCGAGGGCGGGGTGATCGACACCGCCGCCTATCGCGGTCATCCGGTGCTGGTCGTGAaca                                                                                            >  2:2232144/1‑101 (MQ=255)
            cgACCGAGGCGGCGGCGCCGAGCTTCACCTTCGCCTCGATCGAGGGCGGGGTGATCGACACCGCCGCCTATCGCGGTCATCCGGTGCTGGTCGTGAACAcc                                                                                          >  2:1876285/1‑101 (MQ=255)
            cgACCGAGGCGGCGGCGCCGAGCTTCACCTTCGCCTCGATCGAGGGCGGGGTGATCGACACCGCCGCCTATCGCGGTCATCCGGTGCTGGTCGTGAACAcc                                                                                          >  2:451420/1‑101 (MQ=255)
             gACCGAGGCGGCGGCGCCGAGCTTCACCTTCGCCTCGATCGAGGGCGGGGTGATCGACACCGCCGCCTATCGCGGTCATCCGGTGCTGGTCGTGAACACcg                                                                                         >  1:4454346/1‑101 (MQ=255)
              gccGAGGCGGCGGCGCCGAGCTTCACCTTCGCCTCGATCGAGGGCGGGGTGATCGACACCGCCGCCTATCGCGGTCATCCGGTGCTGGTCGTGAACACcgc                                                                                        <  1:2184735/100‑1 (MQ=255)
               ccGAGGCGGCGGCGCCGAGCTTCACCTTCGCCTCGATCGAGGGCGGGGTGATCGACACCGCCGCCTATCGCGGTCATCCGGTGCTGGTCGTGAACACcgcg                                                                                       <  2:1492851/101‑1 (MQ=255)
               ccGAGGCGGCGGCGCCGAGCTTCACCTTCGCCTCGATCGAGGGCGGGGTGATCGACACCGCCGCCTATCGCGGTCATCCGGTGCTGGTCGTGAACACcgcg                                                                                       <  1:2050699/101‑1 (MQ=255)
                 gAGGCGGCGGCGCCGAGCTTCACCTTCGCCTCGATCGAGGGCGGGGTGATCGACACCGCCGCCTATCGCGGTCATCCGGTGCTGGTCGTGAACACCGCGTc                                                                                     >  1:1430295/1‑101 (MQ=255)
                  aGGCGGCGGCGCCGAGCTTCACCTTCGCCTCGATCGAGGGCGGGGTGATCGACACCGCCGCCTATCGCGGTCATCCGGTGCTGGTCGTGAACACCGCGTCg                                                                                    >  1:1954004/1‑101 (MQ=255)
                  aGGCGGCGGCGCCGAGCTTCACCTTCGCCTCGATCGAGGGCGGGGTGATCGACACCGCCGCCTATCGCGGTCATCCGGTGCTGGTCGTGAACACCGCGTCg                                                                                    >  1:1567388/1‑101 (MQ=255)
                  aGGCGGCGGCGCCGAGCTTCACCTTCGCCTCGATCGAGGGCGGGGTGATCGACACCGCCGCCTATCGCGGTCATCCGGTGCTGGTCGTGAACACCGCGTCg                                                                                    >  2:3033590/1‑101 (MQ=255)
                  aGGCGGCGGCGCCGAGCTTCACCTTCGCCTCGATCGAGGGCGGGGTGATCGACACCGCCGCCTATCGCGGTCATCCGGTGCTGGTCGTGAACACCGCGCCg                                                                                    >  1:1625308/1‑101 (MQ=255)
                   ggcggcggcGCCGAGCTTCACCTTCGCCTCGATCGAGGGCGGGGTGATCGACACCGCCGCCTATCGCGGTCATCCGGTGCTGGTCGTGAACACCGCGTCGc                                                                                   >  2:4218898/1‑101 (MQ=255)
                   ggcggcggcGCCGAGCTTCACCTTCGCCTCGATCGAGGGCGGGGTGATCGACACCGCCGCCTATCGCGGTCATCCGGTGCTGGTCGTGAACACCGCGTCGc                                                                                   >  1:3769599/1‑101 (MQ=255)
                   ggcggcggcGCCGAGCTTCACCTTCGCCTCGATCGAGGGCGGGGTGATCGACACCGCCGCCTATCGCGGTCATCCGGTGCTGGTCGTGAACACCGCGTCGc                                                                                   <  2:2991854/101‑1 (MQ=255)
                      ggcggcgCCGAGCTTCACCTTCGCCTCGATCGAGGGCGGGGTGATCGACACCGCCGCCTATCGCGGTCATCCGGTGCTGGTCGTGAATACCGCGTCGctct                                                                                >  2:3911123/1‑101 (MQ=255)
                      ggcggcgCCGAGCTTCACCTTCGCCTCGATCGAGGGCGGGGTGATCGACACCGCCGCCTATCGCGGTCATCCGGTGCTGGTCGTGAACACCGCGTCGctct                                                                                >  2:2029639/1‑101 (MQ=255)
                       gcggcgTCGAGCTTCACCTTCGCCTCGATCGAGGGCGGGGTGATCGACACCGCCGCCTATCGCGGTCATCCGGTGCTGGTCGTGAACACCGCGTCGCTCTg                                                                               <  1:740326/101‑1 (MQ=255)
                       gcggcgCCGAGCTTCACCTTCGCCTCGATCGAGGGCGGGGTGATCGACACCGCCGCCTATCGCGGTCATCCGGTGCTGGTCGTGAACACCGCGTCGCTCTg                                                                               <  1:1469706/101‑1 (MQ=255)
                       gcggcgCCGAGCTTCACCTTCGCCTCGATCGAGGGCGGGGTGATCGACACCGCCGCCTATCGCGGTCATCCGGTGCTGGTCGTGAACACCGCGTCGCTCTg                                                                               >  2:3497074/1‑101 (MQ=255)
                         ggcgCCGAGCTTCACCTTCGCCTCGATCGAGGGCGGGGTGATCGACACCGCCGCCTATCGCGGTCATCCGGTGCTGGTCGTGAACACCGCGTCGCTCTGCg                                                                             <  2:327126/101‑1 (MQ=255)
                         ggcgCCGAGCTTCACCTTCGCCTCGATCGAGGGCGGGGTGATCGACACCGCCGCCTATCGCGGTCATCCGGTGCTGGTCGTGAACACCGCGTCGCTCTGCg                                                                             >  1:2505410/1‑101 (MQ=255)
                         ggcgCCGAGCTTCACCTTCGCCTCGATCGAGGGCGGGGTGATCGACACCGCCGCCTATCGCGGTCATCCGGTGCTGGTCGTGAACACCGCGTCGCTCTGCg                                                                             >  2:273590/1‑101 (MQ=255)
                         ggcgCCGAGCTTCACCTTCGCCTCGATCGAGGGCGGGGTGATCGACACCGCCGCCTATCGCGGTCATCCGGTGCTGGTCGTGAACACCGCGTCGCTCTGCg                                                                             >  1:2941932/1‑101 (MQ=255)
                            gcCGAGCTTCACCTTCGCCTCGATCGAGGGCGGGGTGATCGACACCGCCGCCTATCGCGGTCATCCGGTGCTGGTCGTGAACACCGCGTCGCTCTGCGGCt                                                                          >  1:1246517/1‑101 (MQ=255)
                              cGATCTTCACCTTCGCCTGGATCGAGAGCGTGGTGATCGACACCGCCGTCTATCGCGGTCATCCGGTACTGGTCGTGAACACCGCGTCGCTCTGCGGCGTc                                                                        <  2:4666986/101‑1 (MQ=255)
                               gAGCTTCACCTTCGCCTCGATCGAGGGCGGGGTGATCGACACCGCCGCCTATCGCGGTCATCCGGTGCTGGTCGTGAACACCGCGTCGCTCTGCGGCTTCa                                                                       >  1:2574571/1‑101 (MQ=255)
                               gAGCTTCACCTTCGCCTCGATCGAGGGCGGGGTGATCGACACCGCCGCCTATCGCGGTCATCCGGTGCTGGTCGTGAACACCGCGTCGCTCTGCGGCTTCa                                                                       >  1:4421421/1‑101 (MQ=255)
                               gAGCTTCACCTTCGCCTCGATCGAGGGCGGGGTGATCGACACCGCCGCCTATCGCGGTCATCCGGTGCTGGTCGTGAACACCGCGTCGCTCTGCGGCTTCa                                                                       >  2:3357472/1‑101 (MQ=255)
                               gAGCTTCACCTTCGCCTCGATCGAGGGCGGGGTGATCGACACCGCCGCCTATCGCGGTCATCCGGTGCTGGTCGTGAACACCGCGTCGCTCTGCGGCTTCa                                                                       <  2:1668586/101‑1 (MQ=255)
                                  cTTCACCTTCGCCTCGATCGAGGGCGGGGTGATCGACACGGCCGCCTATCGCGGTCATCCGGTGCTGGTCGTGAACACCGCGTCGCTCTGCGGCTTCAccc                                                                    <  2:2619510/101‑1 (MQ=255)
                                  cTTCACCTTCGCCTCGATCGAGGGCGGGGTGATCGACACCGCCGCCTATCGCGGTCATCCGGTGCTGGTCGTGAACACCGCGTCGCTCTGCGGCTTCAccc                                                                    >  1:3867417/1‑101 (MQ=255)
                                    tCACCTTCGCCTCGATCGAGGGCGGGGTGATCGACACCGCCGCCTATCGCGGTCATCCGGTGCTGGTCGTGAACACCGCGTCGCTCTGCGGCTTCACCCCg                                                                  <  2:1068455/101‑1 (MQ=255)
                                      cccTTCGCCTCGATCGAGGGCGGGGTGATCGACACCGCCGCCTATCGCGGTCATCCGGTGCTGGTCGTGAACACCGCGTCGCTCTGCGGCTTCACCCCGCa                                                                <  2:3601922/100‑1 (MQ=255)
                                      aCCTTCGCCTCGATCGAGGGCGGGGTGATCGACCCCGCCGCCTATCGCGGTCATCCGGTGCTGGTCGTGAACACCGCGTCGCTCTGCGGCTTCACCCCGCa                                                                <  1:3976756/101‑1 (MQ=255)
                                         ttCGCCTCGATCGAGGGCGGGGTGATCGACACCGCCGCCTATCGCGGTCATCCGGTGCTGGTCGTGAACACCGCGTCGCTCTGCGGCTTCACCCCGCAGCt                                                             <  1:2232306/101‑1 (MQ=255)
                                         ttCGCCTCGATCGAGGGCGGGGTGATCGACACCGCCGCCTATCGCGGTCATCCGGTGCTGGTCGTGAACACCGCGTCGCTCTGCGGCTTCACCCCGCAGCt                                                             <  1:4608158/101‑1 (MQ=255)
                                           cGCCTCGATCGAGGGCGGGGTGATCGACACCGCCGCCTATCGCGGTCATCCGGTGCTGGTCGTGAACACCGCGTCGCTCTGCGGCTTCACCCCGCAGCTCg                                                           >  1:685463/1‑101 (MQ=255)
                                           cGCCTCGATCGAGGGCGGGGTGATCGACACCGCCGCCTATCGCGGTCATCCGGTGCTGGTCGTGAACACCGCGTCGCTCTGCGGCTTCACCCCGCAGCTCg                                                           >  1:3546408/1‑101 (MQ=255)
                                           cGCCTCGATCGAGGGCGGGGTGATCGACACCGCCGCCTATCGCGGTCATCCGGTGCTGGTCGTGAACACCGCGTCGCTCTGCGGCTTCACCCCGCAGCTCg                                                           <  1:1762676/101‑1 (MQ=255)
                                             ccTCGATCGAGGGCGGGGTGATCGACACCGCCGCCTATCGCGGTCATCCGGTGCTGGTCGTGAACACCGCGTCGCTCTGCGGCTTCACCCCGCAGCTCGaa                                                         >  1:3738548/1‑101 (MQ=255)
                                              cTCGATCGAGGGCGGGGTGATCGACACCGCCGCCTATCGCGGTCATCCGGTGCTGGTCGTGAACACCGCGTCGCTCTGCGGCTTCACCCCGCAGCTCGAAg                                                        >  2:4451387/1‑101 (MQ=255)
                                               tcgatcgaGGGCGGGGTGATCGACACCGCCGCCTATCACGGTCATCCGGCGCTGGTCGTGAACACCGCGTCGCTCCGCGGCTTCCCCCCGCAGCCCGCAgg                                                       >  2:3411248/1‑101 (MQ=255)
                                                  atcgaGGGCGGGGTGATCGACACCGCCGCCTATCGCGGTCATCCGGTGCTGGTCGTGAACACCGCGTCGCTCTGCGGCTTCACCCCGCAGCTCGAAGGGCt                                                    <  1:963678/101‑1 (MQ=255)
                                                   tcgaGGGCGGGGTGATCGACACCGCCGCCTATCGCGGTCATCCGGTGCTGGTCGTGAACACCGCGTCGCTCTGCGGCTTCACCCCGCAGCTCGAAGGGCTg                                                   <  1:3509857/101‑1 (MQ=255)
                                                    cgaGGGCGGGGTGATCGACACCGCCGCCTATCGCGGTCATCCGGTGCTGGTCGTGAACACCGCGTCGCTCTGCGGCTTCACCCCGCAGCTCGAAGGGCTGc                                                  >  2:1339389/1‑101 (MQ=255)
                                                     gaGGGCGGGGTGATCGACACCGCCGCCTATCGCGGTCATCCGGTGCTGGTCGTGAACACCGCGTCGCTCTGCGGCTTCACCCCGCAGCTCGAAGGGCTGCa                                                 >  2:3561185/1‑101 (MQ=255)
                                                      aGGGCGGGGTGATCGACACCGCCGCCTATCGCGGTCATCCGGTGCTGGTCGTGAACACCGCGTCGCTCTGCGGCTTCACCCCGCAGCTCGAAGGGCTGCAg                                                >  1:571759/1‑101 (MQ=255)
                                                       gggCGGGGTGATCGACACCGCCGCCTATCGCGGTCATCCGGTGCTGGTCGTGAACACCGCGTCGCTCTGCGGCTTCACCCCGCAGCTCGAAGGGCTGCAgg                                               >  1:3518031/1‑101 (MQ=255)
                                                           ggggTGATCGACACCGCCGCCTATCGCGGTCATCCGGTGCTGGTCGTGAACACCGCGTCGCTCTGCGGCTTCACCCCGCAGCTCGACGGGCTGCAGGCGCt                                           >  1:2653988/1‑101 (MQ=255)
                                                           ggggTGATCGACACCGCCGCCTATCGCGGTCATCCGGTGCTGGTCGTGAACACCGCGTCGCTCTGCGGCTTCACCCCGCAGCTCGAAGGGCTGCAGGCGCt                                           >  1:3003062/1‑101 (MQ=255)
                                                            gggTGATCGACACCGCCGCCTATCGCGGTCATCCGGTGCTGGTCGTGAACACCGCGTCGCTCTGCGGCTTCACCCCGCAGCTCGAAGGGCTGCAGGCGCTg                                          <  1:1837614/101‑1 (MQ=255)
                                                              gTGATCGACACCGCCGCCTATCGCGGTCATCCGGTGCTGGTCGTGAACACCGCGTCGCTCTGCGGCTTCACCCCGCAGCTCGAAGGGCTGCAGGCGCTGca                                        >  1:2254001/1‑101 (MQ=255)
                                                               tGATCGACACCGCCGCCTATCGCGGTCATCCGGTGCTGGTCGTGAACACCGCGTCGCTCTGCGGCTTCACCCCGCAGCTCGAAGGGCTGCAGGCGCTGcac                                       >  1:661515/1‑101 (MQ=255)
                                                                 aTCGACACCGCCGCCTATCGCGGTCATCCGGTGCTGGTCGTGAACACCGCGTCGCTCTGCGGCTTCACCCCGCAGCTCGAAGGGCTGCAGGCGCTGCACaa                                     >  1:1812995/1‑101 (MQ=255)
                                                                  tCGACACCGCCGCCTATCGCGGTCATCCGGTGCTGGTCGTGAACACCGCGTCGCTCTGCGGCTTCACCCCGCAGCTCGAAGGGCTGCAGGCGCTGCACAAg                                    <  1:3290240/101‑1 (MQ=255)
                                                                     acacCGCCGCCTATCGCGGTCATCCGGTGCTGGTCGTGAACACCGCGCCGCTCTGCGGCTTCACCCCGCAGATCGAAGGGCTGGCGGCGCTGCACCAGGcc                                 >  2:773707/1‑101 (MQ=255)
                                                                      cacCGCCGCCTATCGCGGTCATCCGGTGCTGGTCGTGAACACCGCGTCGCTCTGCGGCTTCACCCCGCAGCTCGAAGGGCTGCAGGCGCTGCACAAGGCCt                                <  2:3886939/101‑1 (MQ=255)
                                                                        ccgccgCCTATCGCGGTCATCCGGTGCTGGTCGTGAACACCGCGTCGCTCTGCGGCTTCACCCCGCAGCTCGAAGGGCTGCAGGCGCTGCACAAGGCCTgg                              <  1:3185319/101‑1 (MQ=255)
                                                                        ccgccgCCTATCGCGGTCATCCGGTGCTGGTCGTGAACACCGCGTCGCTCTGCGGCTTCACCCCGCAGCTCGAAGGGCTGCAGGCGCTGCACAAGGCCTgg                              <  1:3357721/101‑1 (MQ=255)
                                                                         tctcgcaTATCGCGGTCATCCGGTGCTGGTCGTGAACACCGCGTCGCTCTGCGGCTTCACCCCGCAGCTCGAAGGGCTGCAGGCGCTGCACAAGGCCTggg                             <  1:762611/98‑1 (MQ=255)
                                                                          gcggccTATCGCGCTCATCCGGTGCTGGTCGTGACCACCGGGTCCCTCTGCGGCTTCCCCCCGCACCTCG.AGGGCTGCAGGCGCTGCACAAGGCCTggggg                           <  2:2495029/98‑1 (MQ=255)
                                                                           ccgccTATCGCGGTCATCCGGTGCTGGTCGTGAACACCGCGTCGCTCTGCGGCTTCACCCCGCAGCTCGAAGGGCTGCAGGCGCTGCACAAGGCCTggggg                           <  1:4383897/101‑1 (MQ=255)
                                                                           ccgccTATCGCGGTCATCCGGTGCTGGTCGTGAACACCGCGTCGCTCTGCGGCTTCACCCCGCAGCTCGAAGGGCTGCAGGCGCTGCACAAGGCCTggggg                           <  1:1315300/101‑1 (MQ=255)
                                                                            cgccTATCGCGGTCATCCGGTGCTGGTCGTGAACACCGCGTCGCTCTGCGGCTTCACCCCGCAGCTCGAAGGGCTGCAGGCGCTGCACAAGGCCTGGGGGc                          <  2:3858721/101‑1 (MQ=255)
                                                                             gccTATCGCGGTCATCCGGTGCTGGTCGTGAACACCGCGTCGCTCTGCGGCTTCCCCCCGCAGCTCGAAGGGCTGCAGGCGCTGCACAAGGCCTGGGGGcc                         <  1:910084/101‑1 (MQ=255)
                                                                             gccTATCGCGGTCATCCGGTGCTGGTCGTGAACACCGCGTCGCTCTGCGGCTTCACCCCGCAGCTCGAAGGGCTGCAGGCGCTGCACAAGGCCTGGGGGcc                         <  2:3236724/101‑1 (MQ=255)
                                                                              ccTATCGCGGTCATCCGGTGCTGGTCGTGAACACCGCGTCGCTCTGCGGCTTCACCCCGCAGCTCGAAGGGCTGCAGGCGCTGCACAAGGCCTGGGGGCCg                        <  1:4220413/101‑1 (MQ=255)
                                                                               cTATCGCGGTCATCCGGTGCTGGTCGTGAACACCGCGTCGCTCTGCGGCTTCACCCCGCAGCTCGAAGGGCTGCAGGCGCTGCACAAGGCCTGGGGGCCGa                       >  2:2881971/1‑101 (MQ=255)
                                                                                tATCGCGGTCATCCGGTGCTGGTCGTGAACACCGCGTCGCTCTGCGGCTTCACCCCGCAGCTCGAAGGGCTGCAGGCGCTGCACAAGGCCTGGGGGCCGaa                      >  2:673715/1‑101 (MQ=255)
                                                                                tATCGCGGTCATCCGGTGCTGGTCGTGAACACCGCGTCGCTCTGCGGCTTCACCCCGCAGCTCGAAGGGCTGCAGGCGCTGCACAAGGCCTGGGGGCCGaa                      >  2:1426968/1‑101 (MQ=255)
                                                                                  tCGCGGTCATCCGGTGCTGGTCGTGAACACCGCGTCGCTCTGCGGCTTCACCCCGCAGCTCGAAGGGCTGCAGGCGCTGCACAAGGCCTGGGGGCCGAAgg                    <  2:2777685/101‑1 (MQ=255)
                                                                                  tCGCGGTCATCCGGTGCTGGTCGTGAACACCGCGTCGCTCTGCGGCTTCACCCCGCAGCTCGAAGGGATGCAGGCGCTGCACAAGGACTGGGGGCCGAAgg                    >  1:4023752/1‑101 (MQ=255)
                                                                                   cgcgGTCATCCGGTGCTGGTCGTGAACACCGCGTCGCTCTGCGGCTTCACCCCGCAGCTCGAAGGGCTGCAGGCGCTGCACAAGGCCTGGGGGCCGAAggg                   >  1:4130447/1‑101 (MQ=255)
                                                                                   cgcgGTCATCCGGTGCTGGTCGTGAACACCGCGTCGCTCTGCGGCTTCACCCCGCAGCTCGAAGGGCTGCAGGCGCTGCACAAGGCCTGGGGGCCGAAggg                   <  2:1225308/101‑1 (MQ=255)
                                                                                    gcgGTCATCCGGTGCTGGTCGTGAACACCGCGTCGCTCTGCGGCTTCACCCCGCAGCTCGAAGGGCTGCAGGCGCTGCACAAGGCCTGGGGGCCGAAgggg                  >  1:2708168/1‑101 (MQ=255)
                                                                                     cgGTCATCCGGTGCTGGTCGTGAACACCGCGTCGCTCTGCGGCTTCACCCCGCAGCTCGAAGGGCTGCAGGCGCTGCACAAGGCCTGGGGGCCGAAGGGGc                 <  2:4140689/101‑1 (MQ=255)
                                                                                      ggTCATCCGGTGCTGGTCGTGAACACCGCGTCGCTCTGCGGCTTCACCCCGCAGCTCGAAGGGCTGCAGGCGCTGCACAAGGCCTGGGGGCCGAAGGGGCt                >  2:2839388/1‑101 (MQ=255)
                                                                                        tCATCCGGTGCTGGTCGTGAACACCGCGTCGCTCTGCGGCTTCACCCCGCAGCTCGAAGGGCTGCAGGCGCTGCACAAGGCCTGGGGGCCGAAGGGGCTGa              >  1:3691673/1‑101 (MQ=255)
                                                                                        tCATCCGGAGCTGGTCGTGAACACCGCGTCGCTCTGCGGCTTCACCCCGCAGCTCGAAGGGCTGCAGGCGCTGCACAAGGCCTGGGGGCCGAAGGGGCTGa              >  2:1738681/1‑101 (MQ=255)
                                                                                         cATCCGGTGCTGGTCGTGAACACCGCGTCGCTCTGCGGCTTCACCCCGCAGCTCGAAGGGCTGCAGGCGCTGCACAAGGCCTGGGGGCCGAAGGGGCTGAt             >  2:1458717/1‑101 (MQ=255)
                                                                                         cATCCGGTGCTGGTCGTGAACACCGCGTCGCTCTGCGGCTTCACCCCGCAGCTCGAAGGGCTGCAGGCGCTGCACAAGGCCTGGGGGCCGAAGGGGCTGAt             <  1:2882187/101‑1 (MQ=255)
                                                                                          aTCCGGTGCTGGTCGTGAACACCGCGTCGCTCTGCGGCTTCACCCCGCAGCTCGACGGGCTGCAGGCGCTGCACAAGGCCTGGGGGCCGCAGGGGCTGCTc            >  2:4246864/1‑101 (MQ=255)
                                                                                          aTCCGGTGCTGGTCGTGAACACCGCGTCGCTCTGCGGCTTCACCCCGCAGCTCGAAGGGCTGCAGGCGCTGCACAAGGCCTGGGGGCCGAAGGGGCTGATc            <  1:1876411/101‑1 (MQ=255)
                                                                                            ccGGTGCTGGTCGTGAACACCGCGTCGCTCTGCGGCTTCACCCCGCAGCTCGAAGGGCTGCAGGCGCTGCACAAGGCCTGGGGGCCGAAGGGGCTGATCGt          <  2:4023467/101‑1 (MQ=255)
                                                                                                tGCTGGTCGTGAACACCGCGTCGCTCTGCGGCTTCACCCCGCAGCTCGAAGGGCTGCAGGCGCTGCACAAGGCCTGGGGGCCGAAGGGGCTGATCGTGCTg      <  2:373759/101‑1 (MQ=255)
                                                                                                tGCTGGTCGTGAACACCGCGTCGCTCTGCGGCTTCACCCCGCAGCTCGAAGGGCTGCAGGCGCTGCACAAGGCCTGGGGGCCGAAGGGGCTGATCGTGCTg      <  2:4001913/101‑1 (MQ=255)
                                                                                                tGCTGGTCGTGAACACCGCGTCGCTCTGCGGCTTCACCCCGCAGCTCGAAGGGCTGCAGGCGCTGCACAAGGCCTGGGGGCCGAAGGGGCTGATCGTGCTg      <  1:404624/101‑1 (MQ=255)
                                                                                                tGCTGGTCGTGAACACCGCGTCGCTCTGCGGCTTCACCCCGCAGCTCGAAGGGCTGCAGGCGCTGCACAAGGCCTGGGGGCCGAAGGGGCTGATCGTGCTg      <  1:1058021/101‑1 (MQ=255)
                                                                                                  cTGGTCGTGAACACCGCGTCGCTCTGCGGCTTCACCCCGCAGCTCGAAGGGCTGCAGGCGCTGCACAAGGCCTGGGGGCCGAAGGGGCTGATCGTGCTGGc    >  2:2536670/1‑101 (MQ=255)
                                                                                                    |                                                                                                    
CCGCCGCCCGCGCGACCGAGGCGGCGGCGCCGAGCTTCACCTTCGCCTCGATCGAGGGCGGGGTGATCGACACCGCCGCCTATCGCGGTCATCCGGTGCTGGTCGTGAACACCGCGTCGCTCTGCGGCTTCACCCCGCAGCTCGAAGGGCTGCAGGCGCTGCACAAGGCCTGGGGGCCGAAGGGGCTGATCGTGCTGGCCG  >  NC_014034/770227‑770427

Base quality scores:  ATCG  < 3 ≤  ATCG  < 21 ≤  ATCG  < 29 ≤  ATCG  < 33 ≤  ATCG  < 40 ≤  ATCG