Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_014034 2,394,542 2 bp→CC 100% coding (603‑604/708 nt) modC2 → molybdate ABC transporter ATP‑binding protein ModC

Read alignment evidence...
  seq id position change freq score reads annotation genes product
*NC_0140342,394,5420G→C100.0% 629.6 168E201D (GAG→GAC) modC2molybdate ABC transporter ATP‑binding protein ModC
Reads supporting (aligned to +/- strand):  new base (70/98):  ref base (0/0):  total (70/98)
*NC_0140342,394,5430A→C100.0% 628.6 168T202P (ACC→CCC) modC2molybdate ABC transporter ATP‑binding protein ModC
Reads supporting (aligned to +/- strand):  new base (70/98):  ref base (0/0):  total (70/98)

Only 100 of 173 total aligned reads displayed.

TCGGCGCTTGATCCCGACACCCGGGCCCGGCTGCAAGATCAGCTGCGCCGCCTGCACGATCTCTTTGGCGCCACGACGCTGATGACAAGCCACGAGACCGCCGAGGCGCGGCGCCTGGCCGACCGCGTGCTGCGGCTGGAGGCCGGGCGGCTGGGCGCCGCGCCGCCGCCGCCGCCATCCCTGCGCCTC  >  NC_014034/2394447‑2394635
                                                                                               ||                                                                                            
 cGGCGCTTGATCCCGCCACCCGGGCCCGGCTGCAAGATCAGCTGCGCCGCCTGCACGATCTCTTTGGCGCCACGACGCTGATGACAAGCCACGACCccgcc                                                                                         <  2:463241/101‑1 (MQ=255)
 cGGCGCTTGATCCCGACCCCCGGGCCCGGCTGCAAGATCAGCTGCGCCGCCTGCACGATCTCTTTGGCGCCACGACGCTGATGACAAGCCACGACCccgcc                                                                                         <  1:2917200/101‑1 (MQ=255)
 cGGCGCTTGATCCCGACACCCGGGCCCGGCTGCAAGATCAGCTGCGCCGCCTGCACGATCTCTTTGGCGCCACGACGCTGATGACAAGCCACGACCccgcc                                                                                         >  1:1012589/1‑101 (MQ=255)
 cGGCGCTTGATCCCGACACCCGGGCCCGGCTGCAAGATCAGCTGCGCCGCCTGCACGATCTCTTTGGCGCCACGACGCTGATGACAAGCCACGACCccgcc                                                                                         >  2:6308766/1‑101 (MQ=255)
  ggCGCTTGATCCCGACACCCGGGCCCGGCTGCAAGATCAGCTGCGCCGCCTGCACGATCTCTTTGGCGCCACGACGCTGATGACAAGCCACGACCccgccg                                                                                        >  2:6736784/1‑101 (MQ=255)
   gcgcTTGATCCCGACACCCGGGCCCGGCTGCAAGATCAGCTGCGCCGCCTGCACGATCTCTTTGGCGCCACGACGCTGATGACAAGCCACGACCCCGCCGa                                                                                       >  2:6512515/1‑101 (MQ=255)
   gcgcTTGATCCCGACACCCGGGCCCGGCTGCAAGATCAGCTGCGCCGCCTGCACGATCTCTTTGGCGCCACGACGCTGATGACAAGCCACGACCCCGCCGa                                                                                       >  1:2203039/1‑101 (MQ=255)
    cgcTTGATCCCGACCCCCGGGCCCGGCTGCAAGATCAGCTGCGCCGCCTGCACGATCTCTTTGGCGCCACGACGCTGATGACAAGCCACGACCCCGCCGAg                                                                                      <  2:3013727/101‑1 (MQ=255)
    cgcTGGATCCCGCCACCGGGGCCGGGCTGCAAGATCAGCTGCGCCGCCTGCACGATCTCTTTGGCGCCACGCCGCTGATGACAAGCCACGACCCCGCCGAg                                                                                      <  2:3568572/101‑1 (MQ=255)
      cTTGATCCCGACACCCGGGCCCGGCTGCAAGATCAGCTGCGCCGCCTGCACGATCTCTTTGGCGCCACGACGCTGATGACAAGCCACGACCCCGCCGAggc                                                                                    >  2:1112819/1‑101 (MQ=255)
      cTTGATCCCGACACCCGGGCCCGGCTGCAAGATCAGCTGCGCCGCCTGCACGATCTCTTTGGCGCCACGACGCTGATGACAAGCCACGACCCCGCCGAggc                                                                                    >  2:487272/1‑101 (MQ=255)
       ttGATCCCGACCCCCGGGCCCGGCTGCAAGATCAGCTGCGCCGCCTGCACGATCTCTTTGGCGCCACGACGCTGATGACAAGCCACGACCCCGCCGAggcg                                                                                   <  1:2903162/101‑1 (MQ=255)
        tGATCCCGACCCCCGGGCCCGGCTGCAAGATCAGCTGCGCCGCCTGCACGATCTCTTTGGCGCCACGACGCTGATGACAAGCCACGACCCCGCCGAggcgc                                                                                  <  1:6209122/101‑1 (MQ=255)
         taccccGCCCCCCGGGCCCGGCTGCAAGATCATCTGCGCCGCCTGCACGATCTCTTTGGCGCCACGACGCTGATGACAAGCCACGACCCCGCCGAggcgcg                                                                                 <  1:4207625/98‑1 (MQ=255)
         gATCCCGACACCCGGGCCCGGCTGCAAGATCAGCTGCGCCGCCTGCACGATCTCTTTGGCGCCACGACGCTGATGACAAGCCACGACCCCGCCGAggcgcg                                                                                 <  2:592133/101‑1 (MQ=255)
         gATCCCGACACCCGGGCCCGGCTGCAAGATCAGCTGCGCCGCCTGCACGATCTCTTTGGCGCCACGACGCTGATGACAAGCCACGACCCCGCCGAggcgcg                                                                                 >  2:2973513/1‑101 (MQ=255)
          aTCCCGACACCCGGGCCCGGCTGCAAGATCAGCTGCGCCGCCTGCACGATCTCTTTGGCGCCACGACGCTGATGACAAGCCACGACCCCGCCGAggcgcgg                                                                                >  1:2229439/1‑101 (MQ=255)
              cGACACCCGGGCCCGGCTGCAAGATCAGCTGCGCCGCCTGCACGATCTCTTTGGCGCCACGACGCTGATGACAAGCCACGACCCCGCCGAGGCGCGGCGcc                                                                            <  1:5797768/101‑1 (MQ=255)
              cGACACCCGGGCCCGGCTGCAAGATCAGCTGCGCCGCCTGCACGATCTCTTTGGCGCCACGACGCTGATGACAAGCCACGACCCCGCCCAGGCGCGGCGcc                                                                            <  1:5118112/101‑1 (MQ=255)
                  ccccGGGCCCGGCTGCAAGATCAGCTGCGCCGCCTGCACGACCTCTTTGGCGCCACGACGCTGATGACAAGCCACGCCCCCGCCGAGGCGCGGCGCCTGGc                                                                        <  1:6352152/100‑1 (MQ=255)
                   cccGGGCCCGGCTGCAAGATCAGCTGCGCCGCCTGCACGATCTCTTTGGCGCCACGACGCTGATGACAAGCCACGACCCCGCCGAGGCGCGGCGCCTGGcc                                                                       <  2:383072/101‑1 (MQ=255)
                   cccGCGCCCGGCTGCAAGATCACCTGCGCCGCCTGCACGATCTCTTTGGCGCCACGACGCTGATGACAAGCCACGACCCCGCCGAGGCGCGGCGCCTGGcc                                                                       <  2:2114196/101‑1 (MQ=255)
                    ccGGGCCCGGCTGCAAGATCAGCTGCGCCGCCTGCACGATCTCTTTGGCGCCACGACGCTGATGACAAGCCACGACCCCGCCGAGGCGCGGCGCCTGGCCg                                                                      <  2:946530/101‑1 (MQ=255)
                      gggCCCGGCTGCAAGATCAGCTGCGCCGCCTGCACGATCTCTTTGGCGCCACGACGCTGATGACAAGCCACGACCCCGCCGAGGCGCGGCGCCTGGCCgtc                                                                    <  2:2552666/101‑3 (MQ=255)
                      gggCCCGGCTGCAAGATCAGCTGCGCCGCCTGCACGATCTCTTTGGCGCCACGACGCTGATGACAAGCCACGACCCCGCCGAGGCGCGGCGCCTGGCCGAc                                                                    <  2:5504374/101‑1 (MQ=255)
                       ggCCCGGCTGCAAGATCAGCTGCGCCGCCTGCACGATCTCTTTGGCGCCACGACGCTGATGACAAGCCACGACCCCGCCGAGGCGCGGCGCCTGGCCGAcc                                                                   <  2:6115349/101‑1 (MQ=255)
                        gCCCGGCTGCAAGATCAGCTGCGCCGCCTGCACGATCTCTTTGGCGCCACGACGCTGATGACAAGCCACGACCCCGCCGAGGCGCGGCGCCTGGCCGACcg                                                                  >  1:6068978/1‑101 (MQ=255)
                        gCCCGGCTGCAAGATCAGCTGCGCCGCCTGCACGATCTCTTTGGCGCCACGACGCTGATGACAAGCCACGACCCCGCCGAGGCGCGGCGCCTGGCCGACcg                                                                  <  2:4773117/101‑1 (MQ=255)
                         cccGGCTGCAAGATCAGCTGCGCCGCCTGCACGATCTCTTTGGCGCCACGACGCTGATGACAAGCCACGACCCCGCCGAGGCGCGGCGCCTGGCCGACcgc                                                                 >  2:5557677/1‑101 (MQ=255)
                            ggCTGCAAGATCAGCTGCGCCGCCTGCACGATCTCTTTGGCGCCACGACGCTGATGACAAGCCACGACCCCGCCGAGGCGCGGCGCCTGGCCGACCGCGtg                                                              <  1:1486982/101‑1 (MQ=255)
                            ggCTGCAAGATCAGCTGCGCCGCCTGCACGATCTCTTTGGCGCCACGACGCTGATGACAAGCCACGACCCCGCCGAGGCGCGGCGCCTGGCCGACCGCGtg                                                              >  1:5356397/1‑101 (MQ=255)
                              cTGCAAGATCAGCTGCGCCGCCTGCACGATCTCTTTGGCGCCACGACGCTGATGACAAGCCACGACCCCGCCGAGGCGCGGCGCCTGGCCGACCGCGtgct                                                            >  1:931978/1‑101 (MQ=255)
                              cTGCAAGACCAGCTGCGCCGCCTGCACGATCTCTTTGGCGCCACGACGCTGATGACAAGCCACGACCCCGCCGAGGCGCGGCGCCTGGCCGACCGCGtgct                                                            <  1:6249563/101‑1 (MQ=255)
                               tgtaatgtCAGCTGCGCCGCCTGCACGTTCTCTTTGGCGCCCCGGCGCTGATGACAAGCCACGACCCCGCCGAGGCGCGGCGCCTGGCCGACCGCGtgctg                                                           <  1:4839905/94‑1 (MQ=255)
                               tGCAAGATCAGCTGCGCCGCCTGCACGATCTCTTTGGCGCCACGACGCTGATGACAAGCCACGCCCCCGCCGAGGCGCGGCGCCTGGCCGACCGCGtgctg                                                           <  1:4497975/101‑1 (MQ=255)
                               tGCAAGACCGGCGGCGCCGCCTGCACGATCTCTTTGGCGCCACGACGCTGATGACAAGCCACGACCCCGCCGAGGCGCGGCGCCTGGCCGACCGCGtgctg                                                           <  1:487284/101‑1 (MQ=255)
                               ggCAAGATCAGCTGCGCCGCCTGCACGATCTCTTTGGCGCCACGACGCTGATGACAAGCCACGACCCCGCCGAGGCGCGGCGCCTGGCCGACCGCGtgctg                                                           <  2:3964698/100‑1 (MQ=255)
                                 cAAGATCAGCTGCGCCGCCTGCACGATCTCTTTGGCGCCACGACGCTGATGACAAGCCACGACCCCGCCGAGGCGCGGCGCCTGGCCGACCGCGTGCTGCg                                                         >  2:3333040/1‑101 (MQ=255)
                                 cAAGATCAGCTGCGCCGCCTGCACGATCTCTTTGGCGCCACGACGCTGATGACAAGCCACGACCCCGCCGAGGCGCGGCGCCTGGCCGACCGCGTGCTGCg                                                         >  1:836411/1‑101 (MQ=255)
                                  aaGATCAGCTGCGCCGCCTGCACGATCTCTTTGGCGCCACGACGCTGATGACAAGCCACGACCCCGCCGAGGCGCGGCGCCTGGCCGACCGCGTGCTGCgg                                                        <  1:517850/101‑1 (MQ=255)
                                  aaGATCAGCTGCGCCGCCTGCACGATCTCTTTGGCGCCACGACGCTGATGACAAGCCACGACCCCGCCGAGGCGCGGCGCCTGGCCGACCGCGTGCTGCgg                                                        <  1:5562469/101‑1 (MQ=255)
                                  aaGATCAGCTGCCCCGCCTGCACGATCTCTTTGGCGCCACGACGCTGATGACAAGCCACGACCCCGCCGAGGCGCGGCGCCTGGCCGACCGCGTGCTGCgg                                                        <  2:735049/101‑1 (MQ=255)
                                   aGATCAGCTGCGCCGCCTGCACGATCTCTTTGGCGCCACGACGCTGATGACAAGCCACGACCCCGCCGAGGCGCGGCGCCTGGCCGACCGCGTGCTGCGGc                                                       >  2:52204/1‑101 (MQ=255)
                                   aGATCAGCTGCGCCGCCTGCACGATCTCTTTGGCGCCACGACGCTGATGACAAGCAACGACCCCGCCGAGGCGCGGCGCCTGGCCGACCGCGTGCTGCGGc                                                       <  1:5515503/101‑1 (MQ=255)
                                    gATCAGCTGCCCCGCCTGCACGATCTCTTTGGCGCCACGACGCTGATGACAAGCCACGACCCCGCCGAGGCGCGGCGCCTGGCCGACCGCGTGCTGCGGCt                                                      <  1:3525553/101‑1 (MQ=255)
                                     aTCAGCTGCGCCGCCTGCACGATCTCTTTGGCGCCACGACGCTGATGACAAGCCACGACCCCGCCGAGGCGCGGCGCCTGGCCGACCGCGTGCTGCGGCTg                                                     >  1:1600952/1‑101 (MQ=255)
                                      tcggCTGCGCCGCCTGCACGATCTCTTTGGCGCCACGACGCTGATGACAAGCCACGACCCCGCCGAGGCGCGGCGCCTGGCCGACCGCGTGCTGCGGCTgg                                                    <  2:5982934/98‑1 (MQ=255)
                                      tCAGCTGCGCCGCCTGCACGATCTCTTTGGCGCCACGACGCTGATGACAAGCCACGACCCCGCCGAGGCGCGGCGCCTGGCCGACCGCGTGCTGCGGCTgg                                                    <  1:3617057/101‑1 (MQ=255)
                                       cAGCTGCGCCGCCTGCAGGATCTCTTTGGCGCCAGGACGCTGATGACAAGCCACGACCCCGCCGAGGCGCGGCGCCTGGCCGACCGCGTGCTGCGGCTGGa                                                   <  2:3685260/101‑1 (MQ=255)
                                       cAGCTGCGCCGCCTGCACGATCTCTTTGGCGCCCCGACGCTGATGACAAGCCACGACCCCGCCGAGGCGCGGCGCCTGGCCGACCGCGTGCTGCGGCTGGa                                                   <  2:4971627/101‑1 (MQ=255)
                                        aGCTGCGCCGCCTGCACGATCTCTTTGGCGCCACGCCGCTGATGACAAGCCACGACCCCGCCGAGGCGCGGCGCCTGGCCGACCGCGTGCTGCGGCTGGAg                                                  <  2:6160627/101‑1 (MQ=255)
                                        aGCTGCGCCGCCTGCACGATCTCTTTGGCGCCACGACGCTGATGACAAGCCACGACCCCGCCGAGGCGCGGCGCCTGGCCGACCGCGTGCTGCGGCTGGAg                                                  >  1:4993134/1‑101 (MQ=255)
                                          cTGCGCCGCCTGCACGATCTCTTTGGCGCCACGACGCTGATGACAAGCCACGACCCCGCCGAGGCGCGGCGCCTGGCCGACCGCGTGCTGCGGCTGGCGGc                                                >  1:3039563/1‑101 (MQ=255)
                                           tGCGCCGCCTGCACGATCTCTTTGGCGCCACGACGCTGATGACAAGCCACGACCCCGCCGAGGCGCGGCGCCTGGCCGACCGCGTGCTGCGGCTGGAGGcc                                               >  1:4038969/1‑101 (MQ=255)
                                           tGCGCCGCCTGCACGATCTCTTTGGCGCCACGACGCTGATGACAAGCCACGACCCCGCCGAGGCGCGGCGCCTGGCCGACCGCGTGCTGCGGCTGGAGGcc                                               >  1:4205650/1‑101 (MQ=255)
                                           gccgccgcCTGCACGATCTCTTTGGCGCCACGACGCTGATGACAAGCCACGACCCCGCCGAGGCGCGGCGCCTGGCCGACCGCGTGCTGCGGCTGGAGGcc                                               <  1:3622130/99‑1 (MQ=255)
                                           gccgccgcCTGCACGATCTCTTTGGCGCCACGACGCTGATGACAAGCCACGACCCCGCCGAGGCGCGGCGCCTGGCCGACCGCGTGCTGCGGCTGGAGGcc                                               <  2:429316/99‑1 (MQ=255)
                                            gcgcCGCCTGCACGATCTCTTTGGCGCCACGACGCTGATGACAAGCCACGACCCCGCCGAGGCGCGGCGCCTGGCCGACCGCGTGCTGCGGCTGGAGGCCg                                              <  1:10835/101‑1 (MQ=255)
                                            gccccgccTGCACGATCTCTTTGGCGCCACGACGCTGATGACAAGCCACGACCCCGCCGAGGCGCGGCGCCTGGCCGACCGCGTGCTGCGGCTGGAGGCCg                                              <  1:241378/98‑1 (MQ=255)
                                                   cTGCACGATCTCTTTGGCGCCACGACGCTGATGACAAGCCACGACCCCGCCGAGGCGCGGCGCCTGGCCGACCGCGTGCTGCGGCTGGAGGCCGgggggcc                                       >  1:5980198/1‑100 (MQ=255)
                                                    tGCACGATCTCTTGGGGGCCACGACGCTGATGACAAGCCACGACCCCGCCGAGGCGCGGCGCCTGGCCGACCGCGTGCTGCGGCTGGAGGCCGGGCGGCTg                                      <  2:3428520/101‑1 (MQ=255)
                                                     gCACGATCTCTTTGGCGCCACGACGCTGATGACAAGCCACGACCCCGCCGAGGCGCGGCGCCTGGCCGACCGCGTGCTGCGGCTGGCGGCCGGGCGGCTgg                                     >  1:6778526/1‑101 (MQ=255)
                                                     gCACGATCTCTTTGGCGCCACGACGCTGATGACAAGCCACGACCCCGCCGAGGCGCGGCGCCTGGCCGACCGCGTGCTGCGGCTGGAGGCCGGGCGGCTgg                                     >  1:4554737/1‑101 (MQ=255)
                                                      cACGATCTCTTGGGCGCCACGACGCTGATGACAAGCCACGACCCCGCCGAGGCGCGGCGCCTGGCCGACCGCGTGCTGCGGCTGGAGGCCGGGCGGCTggg                                    <  1:1112854/101‑1 (MQ=255)
                                                        cGATCTCTTTGGCGCCACGACGCTGATGACAAGCCACGACCCCGCCGAGGCGCGGCGCCTGGCCGACCGCGTGCTGCGGCTGGAGGCCGGGCGGCTGGgcg                                  >  2:5559521/1‑101 (MQ=255)
                                                          acctgttgggCACCACGCCGCGGATACCAAGCCACGACCCCGCCGAGGCGCGGCGCCTGGCCGACCGCGTGCTGCGGCTGGAGGCCGGGCGGCTGGgcgcc                                <  2:4205455/93‑1 (MQ=255)
                                                          aTCTCTTTGGCGCCACGACGCTGATGACAAGCCACGACCCCGCCGAGGCGCGGCGCCTGGCCGACCGCGTGCTGCGGCTGGCGGCCGGGCGGCGGGgcgcc                                >  1:3231152/1‑101 (MQ=255)
                                                          aTCTCTTTGGCGCCACGACGCTGATGACAAGCCACGACCCCGCCGAGGCGCGGCGCCTGGCCGACCGCGTGCTGCGGCTGGAGGCCGGGCGGCTGGgcgcc                                <  1:2712296/101‑1 (MQ=255)
                                                          aTCTCTGTGGCGCCACGCCGCTGATGACAAGCCACGACCCCGCCGAGGCGCGGCGCCTGGCCGACCGCGTGCTGCGGCTGGAGGCCGGGCGGCTGGgcgcc                                <  2:6063232/101‑1 (MQ=255)
                                                           tctcTTTGGCCCCACGACGCTGATGACAAGCCACGACCCCGCCGAGGCGCGGCGCCTGGCCGACCGCGTGCTGCGGCTGGAGGCCGGGCGGCTGGgcgccg                               <  1:2130584/101‑1 (MQ=255)
                                                             tctTTGGCGCCACGACGCTGATGACAAGCCACGACCCCGCCGAGGCGCGGCGCCTGGCCGACCGCGTGCTGCGGCTGGAGGCCGGGCGGCTGGgcgccgcg                             <  2:4554498/101‑1 (MQ=255)
                                                             tctTTGGCGCCACGACGCTGATGACAAGCCACGACCCCGCCGAGGCGCGGCGCCTGGCCGACCGCGTGCTGCGGCTGGAGGCCGGGCGGCTGGgcgccgcg                             >  1:4020614/1‑101 (MQ=255)
                                                               tttGGCGCCACGACGCTGATGACAAGCCACGACCCCGCCGAGGCGCGGCGCCTGGCCGACCGCGTGCTGCGGCTGGAGGCCGGGCGGCTGGgcgccgcgcc                           >  1:4278753/1‑101 (MQ=255)
                                                               tttGGCGCCACGACGCTGATGACAAGCCACGACCCCGCCGAGGCGCGGCGCCTGGCCGACCGCGTGCTGCGGCTGGAGGCCGGGCGGCTGGgcgccgcgcc                           <  1:785257/101‑1 (MQ=255)
                                                               tttGGCGCCACGACGCTGATGACAAGCCACGACCCCGCCGAGGCGCGGCGCCTGGCCGACCGCGTGCTGCGGCTGGAGGCCGGGCGGCTGGgcgccgcgcc                           <  2:3959291/101‑1 (MQ=255)
                                                                ttGGCGCCACGACGCTGATGACAAGCCACGACCCCGCCGAGGCGCGGCGCCTGGCCGACCGCGTGCTGCGGCTGGAGGCCGGGCGGCTGGGcgccgcgccg                          <  2:1883884/101‑1 (MQ=255)
                                                                tgggCGCCACGACGCTGATGACAAGCCACGACCCCGCCGAGGCGCGGCGCCTGGCCGACCGCGTGCTGCGGCTGGAGGCCGGGCGGCTGGGcgccgcgccg                          <  1:3472109/99‑1 (MQ=255)
                                                                 tGGCGCCACGACGCTGATGATAAGCCACGACCCCGCCGAGGCGCGGCGCCTGGCCGACCGCGTGCTGCGGCTGGAGGCCGGGCGGTTGGGCgccgcgtcgc                         <  1:6856383/101‑1 (MQ=255)
                                                                 tGGCGCCACGACGCTGATGACAAGCCACGACCCCGCCGAGGCGCGGCGCCTGGCCGACCGCGTGCTGCGGCTGGAGGCCGGGCGGCTGGGCgccgcgccgc                         >  2:1521851/1‑101 (MQ=255)
                                                                  ggCGCCACGCCGCTGATGACAAGCCACGACCCCGCCGAGGCGCGGCGCCTGGCCGACCGCGTGCTGCGGCTGGAGGCCGGGCGGCTGGGCGCCGcgccgcc                        <  1:5528040/101‑1 (MQ=255)
                                                                  ggCGCCACGACGCTGATGACAAGCCACGACCCCGCCGAGGCGCGGCGCCTGGCCGACCGCGTGCTGCGGCTGGAGGCCGGGCGGCTGGGCGCCGcgccgcc                        <  2:783289/101‑1 (MQ=255)
                                                                  ggCCCCACGACGCTGATGACAAGCCACGACCCCGCCGAGGCGCGGCGCCTGGCCGACCGCGTGCTGCGGCTGGAGGCCGGGCGGCTGGGCGCCGcgccgcc                        <  1:3686439/101‑1 (MQ=255)
                                                                  gccgcCACGACGCGGATGACACGCCACGACCCCGCCGAGGCGCGGCGCCTGGCCGACCGCGTGCTGCGGCTGGAGGCCGGGCGGCTGGGCGCCGcgccgcc                        <  2:2865263/99‑1 (MQ=255)
                                                                   gcgcCACGACGCTGATGACAAGCCACGACCCCGCCGAGGCGCGGCGCCTGGCCGACCGCGTGCTGCGGCTGGAGGCCGGGAGGCTGGGCGCCGcgccgccg                       >  1:5603450/1‑101 (MQ=255)
                                                                   gccccACGCCGCTGATGACAGGCCACGACCCCGCCGAGGCGCGGCGCCTGGCCGACCGCGTGCTGCGGCTGGAGGCCGGGCGGCTGGGCGCCGcgtctccg                       <  1:5621886/98‑1 (MQ=255)
                                                                     cccccgccgCTGATGACAAGCCACGCCCCCGCCGAGGCTCGGCGCCTGGCCGACCGCGTGCTGCGGCTGGAGGCCGGGCGGCTGGGCGCCGcgccgccgcc                     <  1:5097099/94‑1 (MQ=255)
                                                                      cccggacgCTGAGGACAGGCCCCGCCCCCGCCGAGGCGCGGCGCCTGGCCGACCGCGTGCTGCGGCTGGAGGCCGGGCGGCTGGGCGCCGcgccgccgccg                    <  1:3565351/97‑1 (MQ=255)
                                                                          gacgCTGATGACAAGCCACGACCCCGCCGAGGCGCGGCGCCTGGCCGACCGCGTGCTGCGGCTGGAGGCCGGGCGGCTGGGCGCCGcgccgccgccgccgc                >  1:397167/1‑101 (MQ=255)
                                                                           ccgCTGATGCCAAGCCACGCCCCCGCCGAGGCGCGGCGCCTGGCCGACCGCGTGCTGCGGCTGGAGGCCGGGCGGCTGGGCGCCGCgccgccgccgccgcc               <  2:2133888/100‑1 (MQ=255)
                                                                           acgCTGATGACAAGCCACGACCCCGCCGAGGCGCGGCGCCTGGCCGACCGCGTGCTGCGGCTGGAGGCCGGGCGGCTGGGCGCCGCgccgccgccgccgcc               >  2:50802/1‑101 (MQ=255)
                                                                           acgCTGAGGACAAGCCACGACCCCGCCGAGGCGCGGCGCCTGGCCGACCGCGTGCTGCGGCTGGAGGCCGGGCGGCTGGGCGCCGCgccgccgccgccgcc               <  1:3811730/101‑1 (MQ=255)
                                                                            cgCTGATGACAAGCCACGACCCCGCCGAGGCGCGGCGCCTGGCCGACCGCGTGCTGCGGCTGGAGGCCGGGCGGGTGGGCGGCGgcccgccgcagcccgac              >  1:2287968/1‑97 (MQ=255)
                                                                             gCTGATGCAAAGCCACGACCCCGCCGAGGCGCGGCGCCTGGCCGACCGCGTGCTGCGGCTGGAGGCCGGGCGGCTGGGCGCCGCGCCGCCGCCGCCGCCAt             <  2:2404334/101‑1 (MQ=255)
                                                                               gagggaCAAGCCACGCCCCCGCCGAGGCGCGGCGCCTGGCCGCCCGCGTGCTGCGGCTGGAGGCCGGGCGGCTGGGCGCCGCGCCGCCGCCGCCGCCATcc           <  2:4591829/97‑1 (MQ=255)
                                                                               gaatcaCAAGCCACGACCCCGCCGAGGCGCGGCGCCTGGCCGACCGCGTGCTGCGGCTGGAGGCCGGGCGGCTGGGCGCCGCGCCGCCGCCGCCGCCATcc           <  2:548859/96‑1 (MQ=255)
                                                                                gtgaccAAGCCACGACCCCGCCGGGGCGCGGCGGCTGGCCGACCGCGTGCTGCGCCTGGAGGCCGGGCGGCTGGGCGCCGCGCCGCCGCCGCCGCCATccc          <  2:7058531/96‑1 (MQ=255)
                                                                                    aaaaGCCACGCCCCCGCCGAGGCGCGGCGCCTGGCCGACCGCGTGCTGCGGCTGGAGGCCGGGCGGCTGGGCGCCGCGCCGCCGCCGCCGCCATccccatg      <  2:2202944/99‑5 (MQ=255)
                                                                                      aaGCCACGCCCCCGCCGAGGCGCGGCGCCTGGCCGACCGCGTGCTGCGGCTGGAGGCCGGGCGGCTGGGCGCCGCGCCGCCGCCGCCGCCATCCCTGCGcc    <  2:4389859/101‑1 (MQ=255)
                                                                                       aGCCACGACCCCGCCGAGGCGCGGCGCCTGGCCGACCGCGTGGTGCGGCTGGAGGCCGGGCGGCTGGGCGCCGCGCCGCCGCCGCCGCCATCCCTGCGCCt   >  2:3441990/1‑101 (MQ=255)
                                                                                        gCCACGACCCCGCCGAGGCGCGGCGCCTGGCCGACCGCGTGCTGCGGCTGGAGGCCGGGCGGCTGGGCGCCGCGCCGCCGCCGCCGCCATCCCTGCGCCTc  >  2:6605058/1‑101 (MQ=255)
                                                                                               ||                                                                                            
TCGGCGCTTGATCCCGACACCCGGGCCCGGCTGCAAGATCAGCTGCGCCGCCTGCACGATCTCTTTGGCGCCACGACGCTGATGACAAGCCACGAGACCGCCGAGGCGCGGCGCCTGGCCGACCGCGTGCTGCGGCTGGAGGCCGGGCGGCTGGGCGCCGCGCCGCCGCCGCCGCCATCCCTGCGCCTC  >  NC_014034/2394447‑2394635

Base quality scores:  ATCG  < 3 ≤  ATCG  < 15 ≤  ATCG  < 23 ≤  ATCG  < 32 ≤  ATCG  < 40 ≤  ATCG