Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_014034 834,164 A→G 100% intergenic (‑830/‑236) mcpA1 ← / → hoxH methyl‑accepting chemotaxis protein McpA/NAD‑reducing hydrogenase HoxS subunit beta

Read alignment evidence...
  seq id position change freq score reads annotation genes product
*NC_014034834,1640A→G100.0% 1009.2 263intergenic (‑830/‑236)mcpA1/hoxHmethyl‑accepting chemotaxis protein McpA/NAD‑reducing hydrogenase HoxS subunit beta
Reads supporting (aligned to +/- strand):  new base (130/133):  ref base (0/0):  total (130/133)

Only 100 of 270 total aligned reads displayed.

TCGGACGCAGCGGGCAATTCGTTTGACCCGCCCGGCGTGCTGGCGCAATTTTGATCCATCCCAAGGGAAGGGTGGAGGATGGAGAGGTTCGGGCCGCGAAGGGCCGCAGATCTTGGAAGATGACGACGGGGTCCGGGTCGGACACAGGGCCCCCGAACCATCGCGTCGACAGCGGGAGAGATCCCGCGCCGACCAC  >  NC_014034/834066‑834261
                                                                                                  |                                                                                                 
 cGGACGCAGCGGGCAATTCGTTTGACCCGCCCGGCGTGCTGGCGCAATTTTGATCCATCCCAAGGGAAGGGTGGAGGATGGAGAGGTTCGGGCCGCGGAgg                                                                                                <  1:4875361/101‑1 (MQ=255)
 cGGACGCAGCGGGCAATTCGTTTGACCCGCCCGGCGTGCTGGCGCAATTTTGATCCATCCCAAGGGAAGGGTGGAGGATGGAGAGGTTCGGGCCGCGGAgg                                                                                                >  1:4317301/1‑101 (MQ=255)
  ggACGCAGCGGGCAATTCGTTTGACCCGCCCGGCGTGCTGGCGCAATTTTGATCCATCCCACGGGAAGGGTGGAGGTTGGAGGGGTTGGGGCCGCGGAggg                                                                                               >  2:1352191/1‑101 (MQ=255)
  ggACGCAGCGGGCAATTCGTTTGACCCGCCCGGCGTGCTGGCGCAATTTTGATCCATCCCAAGGGAAGGGTGGAGGATGGAGAGGTTCGGGCCGCGGAggg                                                                                               <  2:512296/101‑1 (MQ=255)
   gACGCAGCGGGCAATTCGTTTGACCCGCCCGGCGTGCTGGCGCAATTTTGATCCATCCCAAGGGAAGGGTGGAGGATGGAGAGGTTCGGGCCGCGGAGGGc                                                                                              <  1:4809818/101‑1 (MQ=255)
     cGCAGCGGGCAATTCGTTTGACCCGCCCGGCGTGCTGGCGCAATTTTGATCCATCCCAAGGGAAGGGTGGAGGATGGAGAGGTTCGGGCCGCGGAGGGCCg                                                                                            <  2:3115695/101‑1 (MQ=255)
     cGCAGCGGGCAATTCGTTTGACCCGCCCGGCGTGCTGGCGCAATTTTGATCCATCCCAAGGGAAGGGTGGAGGATGGAGAGGTTCGGGCCGCGGAGGGCCg                                                                                            >  1:4639432/1‑101 (MQ=255)
       cAGCGGGCAATTCGTTTGACCCGCCCGGCGTGCTGGCGCAATTTTGATCCATCCCAAGGGAAGGGTGGAGGATGGAGAGGTTCGGGCCGCGGAGGGCCGCa                                                                                          >  2:822327/1‑101 (MQ=255)
       cAGCGGGCAATTCGTTTGACCCGCCCGGCGTGCTGGCGCAATTTTGATCCATCCCAAGGGAAGGGTGGAGGATGGAGAGGTTCGGGCCGCGGAGGGCCGCa                                                                                          <  1:3898017/101‑1 (MQ=255)
       cAGCGGGCAATTCGTTTGACCCGCCCGGCGTGCTGGCGCAATTTTGATCCATCCCAAGGGAAGGGTGGAGGATGGAGAGGTTCGGGCCGCGGAGGGCCGCa                                                                                          >  1:3599104/1‑101 (MQ=255)
        aGCGGGCAATTCGTTTGACCCGCCCGGCGTGCTGGCGCAATTTTGATCCATCCCAAGGGAAGGGTGGAGGATGGAGAGGTTCGGGCCGCGGAGGGCCGCAg                                                                                         <  1:1551385/101‑1 (MQ=255)
         gCGGGCAATTCGTTTGACCCGCCCGGCGTGCTGGCGCAATTTTGATCCATCCCAAGGGAAGGGTGGAGGATGGAGAGGTTCGGGCCGCGGAGGGCCGCAGa                                                                                        <  2:5446755/101‑1 (MQ=255)
         gCGGGCAATTCGTTTGACCCGCCCGGCGTGCTGGCGCAATTTTGATCCATCCCAAGGGAAGGGTGGAGGATGGAGAGGTTCGGGCCGCGGAGGGCCGCAGa                                                                                        <  2:3834809/101‑1 (MQ=255)
          cGGGCAATTCGTTTGACCCGCCCGGCGTGCTGGCGCAATTTTGATCCATCCCAAGGGAAGGGTGGAGGATGGAGAGGTTCGGGCCGCGGAGGGCCGCAGAt                                                                                       >  1:4592767/1‑101 (MQ=255)
           gggCAATTCGTTTGACCCGCCCGGCGTGCTGGCGCAATTTTGATCCATCCCAAGGGAAGGGTGGAGGATGGAGAGGTTCGGGCCGCGGAGGGCCGCAGATc                                                                                      >  1:2866953/1‑101 (MQ=255)
           gggCAATTCGTTTGACCCGCCCGGCGTGCTGGCGCAATTTTGATCCATCCCAAGGGAAGGGTGGAGGATGGAGAGGTTCGGGCCGCGGAGGGCCGCAGATc                                                                                      <  1:10084/101‑1 (MQ=255)
            ggCAATTCGTTTGACCCGCCCGGCGTGCTGGCGCAATTTTGATCCATCCCAAGGGAAGGGTGGAGGATGGAGAGGTTCGGGCCGCGGAGGGCCGCAGATCt                                                                                     <  1:2967882/101‑1 (MQ=255)
              cAATTCGTTTGACCCGCCCGGCGTGCTGGCGCAATTTTGATCCATCCCAAGGGAAGGGTGGAGGATGGAGAGGTTCGGGCCGCGGAGGGCCGCAGATCtag                                                                                   >  1:3029219/1‑99 (MQ=255)
                aTTCGTTTGACCCGCCCGGCGTGCTGGCGCAATTTTGATCCATCCCAAGGGAAGGGTGGAGGATGGAGAGGTTCGGGCCGCGGAGGGCCGCAGATCTTGGa                                                                                 <  1:749579/101‑1 (MQ=255)
                  tCGTTTGACCCGCCCGGCGTGCTGGCGCAATTTTGATCCATCCCAAGGGAAGGGTGGAGGATGGAGAGGTTCGGGCCGCGGAGGGCCGCAGATCTTGGAAg                                                                               <  1:401495/101‑1 (MQ=255)
                  tCGTTTGACCCGCCCGGCGTGCTGGCGCAATTTTGATCCATCCCAAGGGAAGGGTGGAGGATGGAGAGGTTCGGGCCGCGGAGGGCCGCAGATCTTGGAAg                                                                               <  1:3738014/101‑1 (MQ=255)
                    gTTTGACCCGCCCGGCGTGCTGGCGCAATTTTGATCCATCCCAAGGGAAGGGTGGAGGATGGAGAGGTTCGGGCCGCGGAGGGCCGCAGATCTTGGAAGAt                                                                             <  2:5399626/101‑1 (MQ=255)
                    gTTTGACCCGCCCGGCGTGCTGGCGCAATTTTGATCCATCCCAAGGGAAGGGTGGAGGATGGAGAGGTTCGGGCCGCGGAGGGCCGCAGATCTTGGAAGAt                                                                             >  1:2327830/1‑101 (MQ=255)
                    gTTTGACCCGCCCGGCGTGCTGGCGCAATTTTGATCCATCCCAAGGGAAGGGTGGAGGATGGAGAGGTTCGGGCCGCGGAGGGCCGCAGATCTTGGAAGAt                                                                             <  2:3598887/101‑1 (MQ=255)
                     tttGACCCGCCCGGCGTGCTGGCGCAATTTTGATCCATCCCAAGGGAAGGGTGGAGGATGGAGAGGTTCGGGCCGCGGAGGGCCGCAGATCTTGGAAGATg                                                                            <  1:5067680/101‑1 (MQ=255)
                     tttGACCCGCCCGGCGTGCTGGCGCAATTTTGATCCATCCCAAGGGAAGGGTGGAGGATGGAGAGGTTCGGGCCGCGGAGGGCCGCAGATCTTGGAAGATg                                                                            <  2:2973771/101‑1 (MQ=255)
                     tttGACCCGCCCGGCGTGCTGGCGCAATTTTGATCCATCCCAAGGGAAGGGTGGAGGATGGAGAGGTTCGGGCCGCGGAGGGCCGCAGATCTTGGAAGATg                                                                            >  1:1840423/1‑101 (MQ=255)
                     tttGACCCGCCCGGCGTGCTGGCGCAATTTTGATCCATCCCAAGGGAAGGGTGGAGGATGGAGAGGTTCGGGCCGCGGAGGGCCGCAGATCTTGGAAGATg                                                                            >  1:3773867/1‑101 (MQ=255)
                       tGACCCGCCCGGCGTGCTGGCGCAATTTTGATCCATCCCAAGGGAAGGGTGGAGGATGGAGAGGTTCGGGCCGCGGAGGGCCGCAGATCTTGGAAGATgac                                                                          <  1:1123666/101‑1 (MQ=255)
                         aCCCGCCCGGCGTGCTGGCGCAATTTTGATCCATCCCAAGGGAAGGGTGGAGGATGGAGAGGTTCGGGCCGCGGAGGGCCGCAGATCTTGGAAGATgacga                                                                        >  2:2183980/1‑101 (MQ=255)
                          cccgcccgGCGTGCTGGCGCAATTTTGATCCATCCCAAGGGAAGGGTGGAGGATGGAGAGGTTCGGGCCGCGGAGGGCCGCAGATCTTGGAAGATgacgac                                                                       <  1:5449637/101‑1 (MQ=255)
                             gcccgGCGTGCTGGCGCAATTTTGATCCATCCCAAGGGAAGGGTGGAGGATGGAGAGGTTCGGGCCGCGGAGGGCCGCAGATCTTGGAAGATGACGACggg                                                                    >  1:5728112/1‑101 (MQ=255)
                              cccgGCGTGCTGGCGCAATTTTGATCCATCCCAAGGGAAGGGTGGAGGATGGAGAGGTTCGGGCCGCGGAGGGCCGCAGATCTTGGAAGATGACGACgggg                                                                   >  2:4303554/1‑101 (MQ=255)
                               ccgGCGTGCTGGCGCAATTTTGATCCATCCCAAGGGAAGGGTGGAGGATGGAGAGGTTCGGGCCGCGGAGGGCCGCAGATCTTGGAAGATGACGACggggg                                                                  >  2:1866305/1‑100 (MQ=255)
                                cgGCGTGCTGGCGCAATTTTGATCCATCCCAAGGGAAGGGTGGAGGATGGAGAGGTTCGGGCCGCGGAGGGCCGCAGATCTTGGAAGATGACGACGGGGTc                                                                 <  2:1541544/101‑1 (MQ=255)
                                cgGCGTGCTGGCGCAATTTTGATCCATCCCAAGGGAAGGGTGGAGGATGGAGAGGTTCGGGCCGCGGAGGGCCGCAGATCTTGGAAGATGACGACGGGGTc                                                                 >  2:2520681/1‑101 (MQ=255)
                                cgGCGTGCTGGCGCAATTTTGATCCATCCCAAGGGAAGGGCGGAGGATGGAGAGGTTCGGGCCGCGGAGGGCCGCAGATCTTGGAAGATGACGACggaggc                                                                 >  1:1866501/1‑97 (MQ=255)
                                 ggCGTGCTGGCGCAATTTTGATCCATCCCAAGGGAAGGGTGGAGGATGGAGAGGTTCGGGCCGCGGAGGGCCGCAGATCTTGGAAGATGACGACGGGGTcc                                                                <  1:3505536/101‑1 (MQ=255)
                                  gCGTGCTGGCGCAATTTTGATCCATCCCAAGGGAAGGGTGGAGGATGGAGAGGTTCGGGCCGCGGAGGGCCGCAGATCTTGGAAGATGACGACGGGGTCCg                                                               <  1:2674320/101‑1 (MQ=255)
                                   cGTGCTGGCGCAATTTTGATCCATCCCAAGGGAAGGGTGGAGGATGGAGAGGTTCGGGCCGCGGAGGGCCGCAGATCTTGGAAGATGACGACGGGGTCCgg                                                              <  1:5593320/101‑1 (MQ=255)
                                      gCTGGCGCAATTTTGATCCATCCCAAGGGAAGGGTGGAGGATGGAGAGGTTCGGGCCGCGGAGGGCCGCAGATCTTGGAAGATGACGACGGGGTCCGGGTc                                                           <  1:4746192/101‑1 (MQ=255)
                                       cTGGCGCAATTTTGATCCATCCCAAGGGAAGGGTGGAGGATGGAGAGGTTCGGGCCGCGGAGGGCCGCAGATCTTGGAAGATGACGACGGGGTCCGGGTCg                                                          >  1:4235273/1‑101 (MQ=255)
                                        tGGCGCAATTTTGATCCATCCCAAGGGAAGGGTGGAGGATGGAGAGGTTCGGGCCGCGGAGGGCCGCAGATCTTGGAAGATGACGACGGGGTCCGGGTCgg                                                         >  1:1692793/1‑101 (MQ=255)
                                        tGGCGCAATTTTGATCCATCCCAAGGGAAGGGTGGAGGATGGAGAGGTTCGGGCCGCGGAGGGCCGCAGATCTTGGAAGATGACGACGGGGTCCGGGTCgg                                                         >  2:3406841/1‑101 (MQ=255)
                                          gcgcAATTTTGATCCATCCCAAGGGAAGGGTGGAGGATGGAGAGGTTCGGGCCGCGGAGGGCCGCAGATCTTGGAAGATGACGACGGGGTCCGGGTCGGac                                                       <  1:1661177/101‑1 (MQ=255)
                                          gcgcAATTTTGATCCATCCCAAGGGAAGGGTGGAGGATGGAGAGGTTCGGGCCGCGGAGGGCCGCAGATCTTGGAAGATGACGACGGGGTCCGGGTCGGac                                                       >  2:5471505/1‑101 (MQ=255)
                                           cgcAATTTTGATCCATCCCAAGGGAAGGGTGGAGGATGGAGAGGTTCGGGCCGCGGAGGGCCGCAGATCTTGGAAGATGACGACGGGGTCCGGGTCGGaca                                                      <  2:4235013/101‑1 (MQ=255)
                                           cgcAATTTTGATCCATCCCAAGGGAAGGGTGGAGGATGGAGAGGTTCGGGCCGCGGAGGGCCGCAGATCTTGGAAGATGACGACGGGGTCCGGGTCGGaca                                                      <  2:2995531/101‑1 (MQ=255)
                                             cAATTTTGATCCATCCCAAGGGAAGGGTGGAGGATGGAGAGGTTCGGGCCGCGGAGGGCCGCAGATCTTGGAAGATGACGACGGGGTCCGGGTCGGAcaca                                                    <  2:693736/101‑1 (MQ=255)
                                              aaTTTTGATCCATCCCAAGGGAAGGGTGGAGGATGGAGAGGTTCGGGCCGCGGAGGGCCGCAGATCTTGGAAGATGACGACGGGGTCCGGGTCGGACACAg                                                   >  2:593493/1‑101 (MQ=255)
                                               aTTTTGATCCATCCCAAGGGAAGGGTGGAGGATGGAGAGGTTCGGGCCGCGGAGGGCCGCAGATCTTGGAAGATGACGACGGGGTCCGGGTCGGACACAgg                                                  <  2:2122985/101‑1 (MQ=255)
                                               aTTTTGATCCATCCCAAGGGAAGGGTGGAGGATGGAGAGGTTCGGGCCGCGGAGGGCCGCAGATCTTGGAAGATGACGACGGGGTCCGGGTCGGACACAgg                                                  <  1:4647124/101‑1 (MQ=255)
                                                ttttGATCCATCCCAAGGGAAGGGTGGAGGATGGAGAGGTTCGGGCCGCGGAGGGCCGCAGATCTTGGAAGATGACGACGGGGTCCGGGTCGGACACAggg                                                 >  1:1417321/1‑101 (MQ=255)
                                                 tttGATCCATCCCAAGGGAAGGGTGGAGGATGGAGAGGTTCGGGCCGCGGAGGGCCGCAGATCTTGGAAGATGACGACGGGGTCCGGGTCGGACACAGGGc                                                <  1:1988769/101‑1 (MQ=255)
                                                  ttGATCCATCCCAAGGGAAGGGTGGAGGATGGAGAGGTTCGGGCCGCGGAGGGCCGCAGATCTTGGAAGATGACGACGGGGTCCGGGTCGGACACAGGGcc                                               >  2:1510638/1‑101 (MQ=255)
                                                      tccatccCAAGGGAAGGGTGGAGGATGGAGAGGTTCGGGCCGCGGAGGGCCGCAGATCTTGGAAGATGACGACGGGGTCCGGGTCGGACACAGGGCCCCCg                                           >  2:491021/1‑101 (MQ=255)
                                                        catccCAAGGGAAGGGTGGAGGATGGAGAGGTTCGGGCCGCGGAGGGCCGCAGATCTTGGAAGATGACGACGGGGTCCGGGTCGGACACAGGGCCCCCGaa                                         >  1:5332475/1‑101 (MQ=255)
                                                         atccCAAGGGAAGGGTGGAGGATGGAGAGGTTCGGGCCGCGGAGGGCCGCAGATCTTGGAAGATGACGACGGGGTCCGGGTCGGACACAGGGCCCCCGAAc                                        >  2:1972314/1‑101 (MQ=255)
                                                          tccCAAGGGAAGGGTGGAGGATGGAGAGGTTCGGGCCGCGGAGGGCCGCAGATCTTGGAAGATGACGACGGGGTCCGGGTCGGACACAGGGCCCCCGAAcc                                       <  1:3407046/101‑1 (MQ=255)
                                                          tccCAAGGGAAGGGTGGAGGATGGAGAGGTTCGGGCCGCGGAGGGCCGCAGATCTTGGAAGATGACGACGGGGTCCGGGTCGGACACAGGGCCCCCGAAcc                                       <  2:2310890/101‑1 (MQ=255)
                                                           cccAAGGGAAGGGTGGAGTATGGAGAGGTTCGGGCCGCGGAGGGCCGCAGATCTTGAAAGATGACGACGGGGTCCGGGTCGGACACAGGGCCCCCGAACCa                                      <  2:3900085/101‑1 (MQ=255)
                                                           cccAAGGGAAGGGTGGAGGATGGAGAGGTTCGGGCCGCGGAGGGCCGCAGATCTTGGAAGATGACGACGGGGTCCGGGTCGGACACAGGGCCCCCGAACCa                                      <  2:448437/101‑1 (MQ=255)
                                                             cAAGGGAAGGGTGGAGGATGGAGAGGTTCGGGCCGCGGAGGGCCGCAGATCTTGGAAGATGACGACGGGGTCCGGGTCGGACACAGGGCCCCCGAACCATc                                    <  1:4274818/101‑1 (MQ=255)
                                                             cAAGGGAAGGGTGGAGGATGGAGAGGTTCGGGCCGCGGAGGGCCGCAGATCTTGGAAGATGACGACGGGGTCCGGGTCGGACACAGGGCCCCCGAACCATc                                    >  2:2476250/1‑101 (MQ=255)
                                                              aagggaagggTGGAGGATGGAGAGGTTCGGGCCGCGGAGGGCCGCAGATCTTGGAAGATGACGACGGGGTCCGGGTCGGACACAGGGCCCCCGAACCATcg                                   >  1:2156786/1‑101 (MQ=255)
                                                              aagggaagggTGGAGGATGGAGAGGTTCGGGCCGCGGAGGGCCGCAGATCTTGGAAGATGACGACGGGGTCCGGGTCGGACACAGGGCCCCCGAACCATcg                                   >  1:1052959/1‑101 (MQ=255)
                                                              aagggaagggTGGAGGATGGAGAGGTTCGGGCCGCGGAGGGCCGCAGATCTTGGAAGATGACGACGGGGTCCGGGTCGGACACAGGGCCCCCGAACCATcg                                   <  1:2393464/101‑1 (MQ=255)
                                                                gggaagggTGGAGGATGGAGAGGTTCGGGCCGCGGAGGGCCGCAGATCTTGGAAGATGACGACGGGGTCCGGGTCGGACACAGGGCCCCCGAACCATcgcg                                 >  2:4293415/1‑101 (MQ=255)
                                                                 ggaagggTGGAGGATGGAGAGGTTCGGGCCGCGGAGGGCCGCAGATCTTGGAAGATGACGACGGGGTACGGGTCGGACACAGGGCCCCCGAACCATCGCGt                                >  1:4987685/1‑101 (MQ=255)
                                                                   aagggTGGAGGATGGAGAGGTTCGGGCCGCGGAGGGCCGCAGATCTTGGAAGATGACGACGGGGTCCGGGTCGGACACAGGGCCCCCGAACCATCGCGTCg                              >  1:1467688/1‑101 (MQ=255)
                                                                   aagggTGGAGGATGGAGAGGTTCGGGCCGCGGAGGGCCGCAGATCTTGGAAGATGACGACGGGGTCCGGGTCGGACACAGGGCCCCCGAACCATCGCGTCg                              <  1:3188342/101‑1 (MQ=255)
                                                                   aagggTGGAGGATGGAGAGGTTCGGGCCGCGGAGGGCCGCAGATCTTGGAAGATGACGACGGGGTCCGGGTCGGACACAGGGCCCCCGAACCATCGCGTCg                              <  2:2619206/101‑1 (MQ=255)
                                                                   aagggTGGAGGATGGAGAGGTTCGGGCCGCGGAGGGCCGCAGATCTTGGAAGATGACGACGGGGTCCGGGTCGGACACAGGGCCCCCGAACCATCGCGTCg                              >  2:412010/1‑101 (MQ=255)
                                                                      ggTGGAGGATGGAGAGGTTCGGGCCGCGGAGGGCCGCAGATCTTGGAAGATGACGACGGGGTCCGGGTCGGACACAGGGCCCCCGAACCATCGCGTCGACa                           <  2:3754912/101‑1 (MQ=255)
                                                                        tGGAGGATGGAGAGGTTCGGGCCGCGGAGGGCCGCAGATCTTGGAAGATGACGACGGGGTCCGGGTCGGACACAGGGCCCCCGAACCATCGCGTCGACAGc                         <  2:3700346/101‑1 (MQ=255)
                                                                             gaTGGAGAGGTTCGGGCCGCGGAGGGCCGCAGATCTTGGAAGATGACGACGGGGTCCGGGTCGGACACAGGGCCCCCGAACCATCGCGTCGACAGCGGgag                    >  1:2749729/1‑101 (MQ=255)
                                                                             gaTGGAGAGGTTCGGGCCGCGGAGGGCCGCAGATCTTGGAAGATGACGACGGGGTCCGGGTCGGACACAGGGCCCCCGAACCATCGCGTCGACAGCGGgag                    >  1:2035113/1‑101 (MQ=255)
                                                                              aTGGAGAGGTTCGGGCCGCGGAGGGCCGCAGATCTTGGAAGATGACGACGGGGTCCGGGTCGGACACAGGGCCCCCGAACCATCGCGTCGACAGCGGgaga                   >  2:3056763/1‑101 (MQ=255)
                                                                              aTGGAGAGGTTCGGGCCGCGGAGGGCCGCAGATCTTGGAAGATGACGACGGGGTCCGGGTCGGACACAGGGCCCCCGAACCATCGCGTCGACAGCGGgaga                   >  2:4820058/1‑101 (MQ=255)
                                                                              aTGGAGAGGTTCGGGCCGCGGAGGGCCGCAGATCTTGGAAGATGACGACGGGGTCCGGGTCGGACACAGGGCCCCCGAACCATCGCGTCGACAGCGGgaga                   <  2:2365964/101‑1 (MQ=255)
                                                                              aTGGAGAGGTTCGGGCCGCGGAGGGCCGCAGATCTTGGAAGATGACGACGGGGTCCGGGTCGGACACAGGGCCCCCGAACCATCGCGTCGACAGCGGgaga                   <  1:954382/101‑1 (MQ=255)
                                                                              aTGGAGAGGTTCGGGCCGCGGAGGGCCGCAGATCTTGGAAGATGACGACGGGGTCCGGGTCGGACACAGCGCCCCCGAACCATCGCGTCGACAGCGGgaga                   <  1:5127090/101‑1 (MQ=255)
                                                                               tGGAGAGGTTCGGGCCGCGGAGGGCCGCAGATCTTGGAAGATGACGACGGGGTCCGGGTCGGACACAGGGCCCCCGAACCATCGCGTCGACAGCgattgcg                  <  2:1692694/101‑7 (MQ=255)
                                                                                ggAGAGGTTCGGGCCGCGGAGGGCCGCAGATCTTGGAAGATGACGACGGGGTCCGGGTCGGACACAGGGCCCCCGAACCATCGCGTCGACAGCGGgagaga                 >  1:5502890/1‑101 (MQ=255)
                                                                                ggAGAGGTTCGGGCCGCGGAGGGCCGCAGATCTTGGAAGATGACGACGGGGTCCGGGTCGGACACAGGGCCCCCGAACCATCGCGTCGACAGCGGgagaga                 >  1:4130427/1‑101 (MQ=255)
                                                                                ggAGAGGTTCGGGCCGCGGAGGGCCGCAGATCTTGGAAGATAACGACGGGGTCCGGGTCGGACACAGGGCCCCCGAACCATCGCGTCGACAGCGGgaggga                 >  1:4371248/1‑101 (MQ=255)
                                                                                  agagGTTCGGGCCGCGGAGGGCCGCAGATCTTGGAAGATGACGACGGGGTCCGGGTCGGACACAGGGCCCCCGAACCATCGCGTCGACAGCGGGAGAGATc               <  2:1838958/101‑1 (MQ=255)
                                                                                  agagGTTCGGGCCGCGGAGGGCCGCAGATCTTGGAAGATGACGACGGGGTCCGGGTCGGACACAGGGCCCCCGAACCATCGCGTCGACAGCGGGAGAGATc               >  1:2508630/1‑101 (MQ=255)
                                                                                   gagGTTCGGGCCGCGGAGGGCCGCAGATCTTGGAAGATGACGACGGGGTCCGGGTCGGACACAGGGCCCCCGAACCATCGCGTCGACAGCGGGAGAGATcc              <  2:1135250/101‑1 (MQ=255)
                                                                                   gagGTTCGGGCCGCGGAGGGCCGCAGATCTTGGAAGATGACGACGGGGTCCGGGTCGGACACAGGGCCCCCGAACCATCGCGTCGACAGCGGGAGAGATcc              <  2:5776622/101‑1 (MQ=255)
                                                                                     ggTTCGGGCCGCGGAGGGCCGCAGATCTTGGAAGATGACGACGGGGTCCGGGTCGGACACAGGGCCCCCGAACCATCGCGTCGACAGCGGGAGAGATCCcg            <  2:1291673/101‑1 (MQ=255)
                                                                                     ggTTCGGGCCGCGGAGGGCCGCAGATCTTGGAAGATGACGACGGGGTCCGGGTCGGACACAGGGCCCCCGAACCATCGCGTCGACAGCGGGAGAGATCCcg            <  1:330008/101‑1 (MQ=255)
                                                                                      gTTCGGGCCGCGGAGGGCCGCAGATCTTGGAAGATGACGACGGGGTCCGGGTCGGACACAGGGCCCCCGAACCATCGCGTCGACAGCGGGAGAGATCCcgc           >  2:4656214/1‑101 (MQ=255)
                                                                                          gggCCGCGGAGGGCCGCAGATCTTGGAAGATGACGACGGGGTCCGGGTCGGACACAGGGCCCCCGAACCATCGCGTCGACAGCGGGAGGGATCCCGCGCCg       >  1:4344614/1‑101 (MQ=255)
                                                                                           ggCCGCGGAGGGCCGCAGATCTTGGAAGATGACGACGGGGTCCGGGTCGGACACAGGGCCCCCGAACCATCGCGTCGACAGCGGGAGAGATCCCGCGCCGa      >  2:3226669/1‑101 (MQ=255)
                                                                                            gCCGCGGAGGGCCGCAGATCTTGGAAGATGACGACGGGGTCCGGGTCGGACACAGGGCCCCCGAACCATCGCGTCGACAGCGGGAGAGATCCCGCGCCGAc     <  1:2650396/101‑1 (MQ=255)
                                                                                             ccGCGGAGGGCCGCAGATCTTGGAAGATGACGACGGGGTCCGGGTCGGACACAGGGCCCCCGAACCATCGCGTCGACAGCGGGAGAGATCCCGCGCCGAcc    <  1:3437218/101‑1 (MQ=255)
                                                                                              cgcgGAGGGCCGCAGATCTTGGAAGATGACGACGGGGTCCGGGTCGGACACAGGGCCCCCGAACCATCGCGTCGACAGCGGGAGAGATCCCGCGCCGAcc  >  2:1439813/1‑100 (MQ=255)
                                                                                               gcgGAGGGCCGCAGATCTTGGAAGATGACGACGGGGTCCGGGTCGGACACAGGGCCCCCGAACCATCGCGTCGACAGCGGGAGAGATCCCGCGCCGACcac  >  1:828442/1‑101 (MQ=255)
                                                                                               gcgGAGGGCCGCAGATCTTGGAAGATGACGACGGGGTCCGGGTCGGACACAGGGCCCCCGAACCATCGCGTCGACAGCGGGAGAGATCCCGCGCCGACcac  <  1:4453202/101‑1 (MQ=255)
                                                                                                  |                                                                                                 
TCGGACGCAGCGGGCAATTCGTTTGACCCGCCCGGCGTGCTGGCGCAATTTTGATCCATCCCAAGGGAAGGGTGGAGGATGGAGAGGTTCGGGCCGCGAAGGGCCGCAGATCTTGGAAGATGACGACGGGGTCCGGGTCGGACACAGGGCCCCCGAACCATCGCGTCGACAGCGGGAGAGATCCCGCGCCGACCAC  >  NC_014034/834066‑834261

Base quality scores:  ATCG  < 3 ≤  ATCG  < 23 ≤  ATCG  < 30 ≤  ATCG  < 34 ≤  ATCG  < 41 ≤  ATCG