Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_014034 702,786 G→C 100% R812R (CGC→CGG RCAP_rcc00643 ← diguanylate cyclase/phosphodiesterase

Read alignment evidence...
  seq id position change freq score reads annotation genes product
*NC_014034702,7860G→C100.0% 815.2 214R812R (CGC→CGGRCAP_rcc00643diguanylate cyclase/phosphodiesterase
Reads supporting (aligned to +/- strand):  new base (114/100):  ref base (0/0):  total (114/100)

Only 100 of 226 total aligned reads displayed.

ATCTTCAGCGTCGTCACCGGCATCCGGCGCAGATAGGTCAGGCTGGAATAGCCGGTGCCGAAATCATCGATCGCCAGCGCGATTCCCATCTTGCGCAGGCGCTGCAGCGCCTCCATCGCGCCCGGGATGTCGACGGCCAGCTGATCCTCGGTCAGTTCGATCTGCAGCCCGAAGGGCGGCGCGACATCGTCCGAAG  >  NC_014034/702688‑702883
                                                                                                  |                                                                                                 
aTCTTCAGCGTCGTCACCGGCATCCGGCGCAGATAGGTCAGGCTGGAATAGCCGGTGCCGAAATCATCGATCGCCAGCGCGATTCCCATCTTGCGCAGccg                                                                                                 <  1:3952377/101‑1 (MQ=255)
aTCTTCAGCGTCGTCACCGGCATCCGGCGCAGATAGGTCAGGCTGGAATAGCCGGTGCCGAAATCATCGATCGCCAGCGCGATTCCCATCTTGCGCAGccg                                                                                                 >  2:3830079/1‑101 (MQ=255)
  cTTCAGCGTCGTCACCGGCATCCGGCGCAGATAGGTCAGGCTGGAATAGCCGGTGCCGAAATCATCGATCGCCAGCGCGATTCCCATCTTGCGCAGCCGCt                                                                                               >  2:4561555/1‑101 (MQ=255)
  cTTCAGCGTCGTCACCGGCATCCGGCGCAGATAGGTCAGGCTGGAATAGCCGGTGCCGAAATCATCGATCGCCAGCGCGATTCCCATCTTGCGCAGCCGCt                                                                                               >  1:597290/1‑101 (MQ=255)
   ttCAGCGTCGTCACCGGCATCCGGCGCAGATAGGTCAGGCTGGAATAGCCGGTGCCGAAATCATCGATCGCCAGCGCGATTCCCATCTTGCGCAGCCGCTg                                                                                              >  2:2871822/1‑101 (MQ=255)
    tCAGCGTCGTCACCGGCATCCGGCGCAGATAGGTCAGGCTGGAATAGCCGGTGCCGAAATCATCGATCGCCAGCGCGATTCCCATCTTGCGCAGCCGCTGc                                                                                             <  2:2778274/101‑1 (MQ=255)
    tCAGCGTCGTCACCGGCATCCGGCGCAGATAGGTCAGGCTGGAATAGCCGGTGCCGAAATCATCGATCGCCAGCGCGATTCCCATCTTGCGCAGCCGCTGc                                                                                             >  2:742095/1‑101 (MQ=255)
       gCGTCGTCACCGGCATCCGGCGCAGATAGGTCAGGCTGGAATAGCCGGTGCCGAAATCATCGATCGCCAGCGCGATTCCCATCTTGCGCAGCCGCTGCAgc                                                                                          <  1:3155660/101‑1 (MQ=255)
       gCGTCGTCACCGGCATCCGGCGCAGATAGGTCAGGCTGGAATAGCCGGTGCCGAAATCATCGATCGCCAGCGCGATTCCCATCTTGCGCAGCCGCTGCAgc                                                                                          <  1:1810223/101‑1 (MQ=255)
        cgtcgtCACCGGCATCCGGCGCAGATAGGTCAGGCTGGAATAGCCGGTGCCGAAATCATCGATCGCCAGCGCGATTCCCATCTTGCGCAGCCGCTGCAgcg                                                                                         >  1:3730958/1‑101 (MQ=255)
          tcgtcACCGGCATCCGGCGCAGATAGGTCAGGCTGGAATAGCCGGTGCCGAAATCATCGATCGCCAGCGCGATTCCCATCTTGCGCAGCCGCTGCAGCGcc                                                                                       <  2:3005025/101‑1 (MQ=255)
              cACCGGCATCCGGCGCAGATAGGTCAGGCTGGAATAGCCGGTGCCGAAATCATCGATCGCCAGCGCGATTCCCATCTTGCGCAGCCGCTGCAGCGCCTCCa                                                                                   <  2:901610/101‑1 (MQ=255)
              cACCGGCATCCGGCGCAGATAGGTCAGGCTGGAATAGCCGGTGCCGAAATCATCGATCGCCAGCGCGATTCCCATCTTGCGCAGCCGCTGCAGCGCCTCCa                                                                                   <  2:413709/101‑1 (MQ=255)
               aCCGGCATCCGGCGCAGATAGGTCAGGCTGGAATAGCCGGTGCCGAAATCATCGATCGCCAGCGCGATTCCCATCTTGCGCAGCCGCTGCAGCGCCTCCAt                                                                                  <  1:4426905/101‑1 (MQ=255)
                 cGGCATCCGGCGCAGATAGGTCAGGCTGGAATAGCCGGTGCCGAAATCATCGATCGCCAGCGCGATTCCCATCTTGCGCAGCCGCTGCAGCGCCTCCATcg                                                                                >  1:3344160/1‑101 (MQ=255)
                   gCATCCGGCGCAGATAGGTCAGGCTGGAATAGCCGGTGCCGAAATCATCGATCGCCAGCGCGATTCCCATCTTGCGCAGCCGCTGCAGCGCCTCCATcgcg                                                                              <  2:2753921/101‑1 (MQ=255)
                    cATCCGGCGCAGATAGGTCAGGCTGGAATAGCCGGTGCCGAAATCATCGATCGCCAGCGCGATTCCCATCTTGCGCAGCCGCTGCAGCGCCTCCATCgcgc                                                                             >  1:2327052/1‑101 (MQ=255)
                     aTCCGGCGCAGATAGGTCAGGCTGGAATAGCCGGTGCCGAAATCATCGATCGCCAGCGCGATTCCCATCTTGCGCAGCCGCTGCAGCGCCTCCATCGCGcc                                                                            >  1:2215959/1‑101 (MQ=255)
                      tCCGGCGCAGATAGGTCAGGCTGGAATAGCCGGTGCCGAAATCATCGATCGCCAGCGCGATTCCCATCTTGCGCAGCCGCTGCAGCGCCTCCATCGCGccc                                                                           <  1:3277992/101‑1 (MQ=255)
                        cGGCGCAGATAGGTCAGGCTGGAATAGCCGGTGCCGAAATCATCGATCGCCAGCGCGATTCCCATCTTGCGCAGCCGCTGCAGCGCCTCCATCGCGCCCgg                                                                         <  2:91472/101‑1 (MQ=255)
                        cGGCGCAGATAGGTCAGGCTGGAATAGCCGGTGCCGAAATCATCGATCGCCAGCGCGATTCCCATCTTGCGCAGCCGCTGCAGCGCCTCCATCGCGCCCgg                                                                         >  2:972236/1‑101 (MQ=255)
                           cgcAGATAGGTCAGGCTGGAATAGCCGGTGCCGAAATCATCGATCGCCAGCGCGATTCCCATCTTGCGCAGCCGCTGCAGCGCCTCCATCGCGCCCGGGAt                                                                      <  2:3039871/101‑1 (MQ=255)
                            gcAGATAGGTCAGGCTGGAATAGCCGGTGCCGAAATCATCGATCGCCAGCGCGATTCCCATCTTGCGCAGCCGCTGCAGCGCCTCCATCGCGCCCGGGATg                                                                     >  2:4360939/1‑101 (MQ=255)
                            gcAGATAGGTCAGGCTGGAATAGCCGGTGCCGAAATCATCGATCGCCAGCGCGATTCCCATCTTGCGCAGCCGCTGCAGCGCCTCCATCGCGCCCGGGATg                                                                     <  1:1148210/101‑1 (MQ=255)
                            gcAGATAGGTCAGGCTGGAATAGCCGGTGCCGAAATCATCGATCGCCAGCGCGATTCCCATCTTGCGCAGCCGCTGCAGCGCCTCCATCGCGCCCGGGATg                                                                     >  1:3490150/1‑101 (MQ=255)
                            gcAGATAGGTCAGGCTGGAATAGCCGGTGCCGAAATCATCCATCGCCAGCGCGATTCCCATCTTGCGCAGCCGCTGCAGCGCCTCCATCGCGCCCGGGATg                                                                     <  1:1558154/101‑1 (MQ=255)
                              aGATAGGTCAGGCTGGAATAGCCGGTGCCGAAATCATCGATCGCCAGCGCGATTCCCATCTTGCGCAGCCGCTGCAGCGCCTCCATCGCGCCCGGGATGTc                                                                   >  2:1665843/1‑101 (MQ=255)
                              aGATAGGTCAGGCTGGAATAGCCGGTGCCGAAATCATCGATCGCCAGCGCGATTCCCATCTTGCGCAGCCGCTGCAGCGCCTCCATCGCGCCCGGGATGTc                                                                   >  1:1626088/1‑101 (MQ=255)
                              aGATAGGTCAGGCTGGAATAGCCGGTGCCGAAATCATCGATCGCCAGCGCGATTCCCATCTTGCGCAGCCGCTGCAGCGCCTCCATCGCGCCCGGGATGTc                                                                   <  1:4151394/101‑1 (MQ=255)
                              aGATAGGTCAGGCTGGAATAGCCGGTGCCGAAATCATCGATCGCCAGCGCGATTCCCATCTTGCGCAGCCGCTGCAGCGCCTCCATCGCGCCCGGGATGTc                                                                   >  2:4288683/1‑101 (MQ=255)
                               gATAGGTCAGCCTGGAATAGCCGGTGCCGAAATCATCGATCGCCAGCGCGATTCCCATCTTGCGCAGCCGCTGCAGCGCCTCCATCGCGCCCGGGATGTCg                                                                  >  1:834948/1‑101 (MQ=255)
                                aTAGGTCAGGCTGGAATAGCCGGTGCCGAAATCATCGATCGCCAGCGCGATTCCCATCTTGCGCAGCCGCTGCAGCGCCTCCATCGCGCCCGGGATGTCGa                                                                 >  1:4298036/1‑101 (MQ=255)
                                aTAGGTCAGGCTGGAATAGCCGGTGCCGAAATCATCGATCGCCAGCGCGATTCCCATCTTGCGCAGCCGCTGCAGCGCCTCCATCGCGCCCGGGATGTCGa                                                                 <  1:4215347/101‑1 (MQ=255)
                                 tAGGTCAGGCTGGAATAGCCGGTGCCGAAATCATCGATCGCCAGCGCGATTCCCATCTTGCGCAGCCGCTGCAGCGCCTCCATCGCGCCCGGGATGTCGAc                                                                >  1:1041243/1‑101 (MQ=255)
                                  acgTCAGGCTGGAATAGCCGGTGCCGAAATCATCGATCCCCAGCGCGATTCCCATCTTCCGCAGCCGCTGCAGCGCCTCCATCGCGCCCGGGATGTCGACg                                                               <  1:1103237/99‑1 (MQ=255)
                                   ggTCAGGCTGGAATAGCCGGTGCCGAAATCATCGATCGCCAGCGCGATTCCCATCTTGCGCAGCCGCTGCAGCGCCTCCATCGCGCCCGGGATGTCGACgg                                                              >  1:932163/1‑101 (MQ=255)
                                     tCAGGCTGGAATAGCCGGTGCCGAAATCATCGATCGCCAGCGCGATTCCCATCTTGCGCAGCCGCTGCAGCGCCTCCATCGCGCCCGGGATGTCGACGGcc                                                            >  2:4176421/1‑101 (MQ=255)
                                     tCAGGCTGGAATAGCCGGTGCCGAAATCATCGATCGCCAGCGCGATTCCCATCTTGCGCAGCCGCTGCAGCGCCTCCATCGCGCCCGGGATGTCGACGGcc                                                            <  1:304982/101‑1 (MQ=255)
                                       aGGCTGGAATAGCCGGTGCCGAAATCATCGATCGCCAGCGCGATTCCCATCTTGCGCAGCCGCTGCAGCGCCTCCATCGCGCCCGGGATGTCGACGGCCAg                                                          <  1:871804/101‑1 (MQ=255)
                                         gcgggAATAGCCGGTGCCGAAATCATCGATCGCCAGCGCGATTCCCATCTTGCGCAGCCGCTGCAGCGCCTCCATCGCGCCCGGGATGTCGACGGCCAGCt                                                        <  1:320985/98‑1 (MQ=255)
                                         gCTGGAATAGCCGGTGCCGAAATCATCGATCGCCAGCGCGATTCCCATCTTGCGCAGCCGCTGCAGCGCCTCCATCGCGCCCGGGATGTCGACGGCCAGCt                                                        <  2:245824/101‑1 (MQ=255)
                                           tGGAATAGCCGGTGCCGAAATCATCGATCGCCAGCGCGATTCCCATCTTGCGCAGCCGCTGCAGCGCCTCCATCGCGCCCGGGATGTCGACGGCCAGCTGa                                                      >  1:2120215/1‑101 (MQ=255)
                                           tGGAATAGCCGGTGCCGAAATCATCGATCGCCAGCGCGATTCCCATCTTGCGCAGCCGCTGCAGCGCCTCCATCGCGCCCGGGATGTCGACGGCCAGCTGa                                                      >  1:511374/1‑101 (MQ=255)
                                           tGGAATAGCCGGTGCCGAAATCATCGATCGCCAGCGCGATTCCCATCTTGCGCAGCCGCTGCAGCGCCTCCATCGCGCCCGGGATGTCGACGGCCAGCTGa                                                      <  2:2909611/101‑1 (MQ=255)
                                           tGGAATAGCCGGTGCCGAAATCATCGATCGCCAGCGCGATTCCCATCTTGCGCAGCCGCTGCAGCGCCTCCATCGCGCCCGGGATGTCGACGGCCAGCTGa                                                      >  2:3269351/1‑101 (MQ=255)
                                            ggAATAGCCGGTGCCGAAATCATCGATCGCCAGCGCGATTCCCATCTTGCGCAGCCGCTGCAGCGCCTCCATCGCGCCCGGGATGTCGACGGCCAGCTGAt                                                     >  2:2795642/1‑101 (MQ=255)
                                            ggAATAGCCGGTGCCGAAATCATCGATCGCCAGCGCGATTCCCATCTTGCGCAGCCGCTGCAGCGCCTCCATCGCGCCCGGGATGTCGACGGCCAGCTGAt                                                     >  2:2858353/1‑101 (MQ=255)
                                                   ccGGTGCCGAAATCATCGATCGCCAGCGCGATTCCCATCTTGCGCAGCCGCTGCAGCGCCTCCATCGCGCCCGGGATGTCGACGGCCAGCTGATCCTCGGt                                              >  1:1425878/1‑101 (MQ=255)
                                                      gTGCCGAAATCATCGATCGCCAGCGCGATTCCCATCTTGCGCAGCCGCTGCAGCGCCTCCATCGCGCCCGGGATGTCGACGGCCAGCTGATCCTCGGTCAg                                           >  1:1762018/1‑101 (MQ=255)
                                                      gTGCCGAAATCATCGATCGCCAGCGCGAGTCCCATCTTGCGCAGCCGCTGCAACTCCCCCATCCCGCCCGGGATGTCGACGGCCAGCTGATCCACGGACAg                                           >  2:3943168/1‑101 (MQ=255)
                                                       tGCCGAAATCATCGATCGCCAGCGCGATTCTCATTTTGCGCAGCCGCTGCAGCGCCTCCATCGCGCCCGGGATGTCGACGGCCAGCTGATCCTCGGTCAGt                                          <  1:332714/101‑1 (MQ=255)
                                                        gCCGAAATCATCGATTTCAAGAGCGATTCCCATCTTGCACAGCCGCTGCAGCTCCTCATTCGCGCAAGGGATGCCGATGCGCAGATGAACCTAGGCCAGtt                                         >  2:678189/1‑101 (MQ=255)
                                                          cGAAATCATCGATCGCCAGCGCGATTCCCATCTTGCGCAGCCGCTGCAGCGCCTCCATCGCGCCCGGGATGTCGACGGCCAGCTGATCCTCGGTCAGTTCg                                       <  1:2426769/101‑1 (MQ=255)
                                                          cGAAATCATCGATCGCCAGCGCGATTCCCATCTTGCGCAGCCGCTGCAGCGCCTCCATCGCGCCCGGGATGTCGACGGCCAGCTGATCCTCGGTCAGTTCg                                       >  2:358952/1‑101 (MQ=255)
                                                          cGAAATCATCGATCGCCAGCGCGATTCCCATCTTGCGCAGCCGCTGCAGCGCCTCCATCGCGCCCGGGATGTCGACGGCCAGCTGATCCTCGGTCAGTTCg                                       <  1:4481419/101‑1 (MQ=255)
                                                          cGAAATCATCGATCGCCAGCGCGATTCCCATCTTGCGCAGCCGCTGCAGCGCCTCCATCGCGCCCGGGATGTCGACGGCCAGCTGATCCTCGGTCAGTTCg                                       <  1:1235824/101‑1 (MQ=255)
                                                           gAAATCATCGATCGCCAGCGCGATTCCCATCTTGCGCAGCCGCTGCAGCGCCTCCATCGCGCCCGGGATGTCGACGGCCAGCTGATCCTCGGTCAGTTCGa                                      >  1:1659915/1‑101 (MQ=255)
                                                           gAAATCATCGATCGCCAGCGCGATTCCCATCTTGCGCAGCCGCTGCAGCGCCTCCATCGCGCCCGGGATGTCGACGGCCAGCTGATCCTCGGTCAGTTCGa                                      <  1:4383043/101‑1 (MQ=255)
                                                            aaaTCATCGATCGCCAGCGCGATTCCCATCTTGCGCAGCCGCTGCAGCGCCTCCATCGCGCCCGGGATGTCGACGGCCAGCTGATCCTCGGTCAGGTCGAt                                     >  2:4683415/1‑101 (MQ=255)
                                                             aaTCATCGATCGCCAGCGCGATTCCCATCTTGCGCAGCCGCTGCAGCGCCTCCATCGCGCCCGGGATGTCGACGGCCAGCTGATCCTCGGTCAGTTCGATc                                    >  2:2940861/1‑101 (MQ=255)
                                                             aaTCATCGATCGCCAGCGCGATTCCCATCTTGCGCAGCCGCTGCAGCGCCTCCATCGCGCCCGGGATGTCGACGGCCAGCTGATCCTCGGTCAGTTCGATc                                    >  1:4333175/1‑101 (MQ=255)
                                                              atcatcGATCGCCAGCGCGATTCCCATCTTGCGCAGCCGCTGCAGCGCCTCCATCGCGCCCGGGATGTCGACGGCCAGCTGATCCTCGGTCAGTTCGATCt                                   >  1:2743940/1‑101 (MQ=255)
                                                              ataatcgatcgCTCACTCAATTCCCATCTTGGCCCGCCGCTGCAGCACCTCCATCGCTCCCGGGATGTCGACGGCCAGCTGAGACTCGGTCAGTTCGATCt                                   <  2:3333113/98‑1 (MQ=255)
                                                               tcatcGATCGCCAGCGCGATTCCCATCTTGCGCAGCCGCTGCAGCGCCTCCATCGCGCCCGGGATGTCGACGGCCAGCTGATCCTCGGTCAGTTCGATCTg                                  <  2:1384140/101‑1 (MQ=255)
                                                                catcGATCGCCAGCGCGATTCCCATCTTGCGCAGCCGCTGCAGCGCCTCCATCGCGCCCGGGATGTCGACGGCCAGCTGATCCTCGGTCAGTTCGATCTGc                                 >  2:2641878/1‑101 (MQ=255)
                                                                  tcgatcgCCAGCGCGATTCCCATCTTGCGCAGCCGCTGCAGCGCCTCCATCGCGCCCGGGATGTCGACGGCCAGCTGATCCTCGGTCAGTTCGATCTGCAg                               <  2:1213142/101‑1 (MQ=255)
                                                                  tcgatcgCCAGCGCGATTCCCATCTTGCGCAGCCGCTGCAGCGCCTCCATCGCGCCCGGGATGTCGACGGCCAGCTGATCCTCGGTCAGTTCGATCTGCAg                               <  2:1744861/101‑1 (MQ=255)
                                                                   cgatcgCCAGCGCGATTCCCATCTTGCGCAGCCGCTGCAGCGCCTCCATCGCGCCCGGGATGTCGACGGCCAGCTGATCCTCGGTCAGTTCGATCTCCAGc                              >  1:3439409/1‑101 (MQ=255)
                                                                     atcgCCAGCGCGATTCCCATCTTGCGCAGCCGCTGCAGCGCCTCCATCGCGCCCGGGATGTCGACGGCCAGCTGATCCTCGGTCAGTTCGATCTGCAGccc                            <  2:3008822/101‑1 (MQ=255)
                                                                      tcgCCCGCGCGATTCCCATCTTGCGCAGCCGCTGCAGCGCCTCCATCGCGCCCGGGATGTCGACGGCCAGCTGATCCTCGGTCAGTTCGATCTGCAGCCCg                           <  2:1329156/101‑1 (MQ=255)
                                                                      tcgCCAGCGCGATTCCCATCTTGCGCAGCCGCTGCAGCGCCTCCATCGCGCCCGGGATGTCGACGGCCAGCTGATCCTCGGTCAGTTCGATCTGCAGCCCg                           >  1:1925224/1‑101 (MQ=255)
                                                                      tcgCCAGCGCGATTCCCATCTTGCGCAGCCGCTGCAGCGCCTCCATCGCGCCCGGGATGTCGACGGCCAGCTGATCCTCGGTCAGTTCGATCTGCAGCCCg                           <  2:806027/101‑1 (MQ=255)
                                                                      tcgCCAGCGCGATTCCCATCTTGCGCAGCCGCTGCAGCGCCTCCATCGCGCCCGGGATGTCGACGGCCAGCTGATCCTCGGTCAGTTCGATCTGCAGCCCg                           <  2:366649/101‑1 (MQ=255)
                                                                        gCCAGCGCGATTCCCATCTTGCGCAGCCGCTGCAGCGCCTCCATCGCGCCCGGGATGTCGACGGCCAGCTGATCCTCGGTCAGTTCGATCTGCAGCCCGaa                         <  1:4020127/101‑1 (MQ=255)
                                                                        gCCAGCGCGATTCCCATCTTGCGCAGCCGCTGCAGCGCCTCCATCGCGCCCGGGATGTCGACGGCCAGCTGATCCTCGGTCAGTTCGATCTGCAGCCCGaa                         >  2:3444656/1‑101 (MQ=255)
                                                                          cAGCGCGATTCCCATCTTGCGCAGCCGCTGCAGCGCCTCCATCGCGCCCGGGATGTCGACGGCCAGCTGATCCTCGGTCAGTTCGATCTGCAGCCCGAAgg                       <  1:2920562/101‑1 (MQ=255)
                                                                               cgATTCCCATCTTGCGCAGCCGCTGCAGCGCCTCCATCGCGCCCGGGATGTCGACGGCCAGCTGATCCTCGGTCAGTTCGATCTGCAGCCCGAAGggcggc                  <  2:4469629/101‑1 (MQ=255)
                                                                                gATTCCCATCTTGCGCAGCCGCTGCAGCGCCTCCATCGCGCCCGGGATGTCGACGGCCAGCTGATCCTCGGTCAGTTCGATCTGCAGCCCGAAGGgcggcg                 <  1:1005316/101‑1 (MQ=255)
                                                                                gATTCCCATCTTGCGCAGCCGCTGCAGCGCCTCCATCGCGCCCGGGATGTCGACGGCCAGCTGATCCTCGGTCAGTTCGATCTGCAGCCCGAAGGgcggcg                 <  2:2235701/101‑1 (MQ=255)
                                                                                gATTCCCATCTTGCGCAGCCGCTGCAGCGCCTCCATCGCGCCCGGGATGTCGACGGCCAGCTGATCCTCGGTCAGTTCGATCTGCAGCCCGAAGGgcggcg                 >  1:1270076/1‑101 (MQ=255)
                                                                                 aTTCCCATCTTGCGCAGCCGCTGCAGCGCCTCCATCGCGCCCGGGATGTCGACGGCCAGCTGATCCTCGGTCAGTTCGATCTGCAGCCCGAAGGGCGgcgc                >  2:2340896/1‑101 (MQ=255)
                                                                                  ttCCCATCTTGCGCAGCCGCTGCAGCGCCTCCATCGCGCCCGGGATGTCGACGGCCAGCTGATCCTCGGTCAGTTCGATCTGCAGCCCGAAGGGCGGcgcg               <  2:4695829/101‑1 (MQ=255)
                                                                                  ttCCCATCTTGCGCAGCCGCTGCAGCGCCTCCATCGCGCCCGGGATGTCGACGGCCAGCTGATCCTCGGTCAGTTCGATCTGCAGCCCGAAGGGCGGcgcg               >  1:3888264/1‑101 (MQ=255)
                                                                                  ttCCCATCTTGCGCAGCCGCTGCAGCGCCTCCATCGCGCCCGGGATGTCGACGGCCAGCTGATCCTCGGTCAGTTCGATCTGCAGCCCGAAGGGCGGcgcg               >  2:2733753/1‑101 (MQ=255)
                                                                                   tCCCATCTTGCGCAGCCGCTGCAGCGCCTCCATCGCGCCCGGGATGTCGACGGCCAGCTGATCCTCGGTCAGTTCGATCTGCAGCCCGAAGGGCGGCGCGa              >  2:1707842/1‑101 (MQ=255)
                                                                                    cccATCTTGCGCAGCCGCTGCAGCGCCTCCATCGCGCCCGGGATGTCGACGGCCAGCTGATCCTCGGTCAGTTCGATCTGCAGCCCGAAGGGCGGCGCGAc             >  1:680000/1‑101 (MQ=255)
                                                                                     ccATCTTGCGCAGCCGCTGCAGCGCCTCCATCGCGCCCGGGATGTCGACGGCCAGCTGATCCTCGGTCAGTTCGATCTGCAGCCCGAAGGGCGGCGCGACa            <  1:4066016/101‑1 (MQ=255)
                                                                                     ccATCTTGCGCAGCCGCTGCAGCGCCTCCATCGCGCCCGGGATGTCGACGGCCAGCTGATCCTCGGTCAGTTCGATCTGCAGCCCGAAGGGCGGCGCGACa            >  2:3356064/1‑101 (MQ=255)
                                                                                      cATCTTGCGCAGCCGCTGCAGCGCCTCCATCGCGCCCGGGATGTCGACGGCCAGCTGATCCTCGGTCAGTTCGATCTGCAGCCCGAAGGGCGGCGCGACAt           >  2:1809720/1‑101 (MQ=255)
                                                                                       aTCTTGCGCAGCCGCTGCAGCGCCTCCATCGCGCCCGGGATGTCGACGGCCAGCTGATCCTCGGTCAGTTCGATCTGCAGCCCGAAGGGCGGCGCGACATc          <  1:1485467/101‑1 (MQ=255)
                                                                                       aTCTTGCGCAGCCGCTGCAGCGCCTCCATCGCGCCCGGGATGTCGACGGCCAGCTGATCCTCGGTCAGTTCGATCTGCAGCCCGAAGGGCGGCGCGACATc          <  1:3597186/101‑1 (MQ=255)
                                                                                          ttGCGCAGCCGCTGCAGCGCCTCCATCGCGCCCGGGATGTCGACGGCCAGCTGATCCTCGGTCAGTTCGATCTGCAGCCCGAAGGGCGGCGCGACATCGTc       >  1:3534306/1‑101 (MQ=255)
                                                                                           tGCGCAGCCGCTGCAGCGCCTCCATCGCGCCCGGGATGTCGACGGCCAGCTGATCCTCGGTCAGTTCGATCTGCAGCCCGAAGGGCGGCGCGACATCGTcc      <  1:2547027/101‑1 (MQ=255)
                                                                                           tGCGCAGCCGCTGCAGCGCCTCCATCGCGCCCGGGATGTCGACGGCCAGCTGATCCTCGGTCAGTTCGATCTGCAGCCCGAAGGGCGGCGCGACATCGTcc      >  1:1194169/1‑101 (MQ=255)
                                                                                            gcgcAGCCGCTGCAGCGCCTCCATCGCGCCCGGGATGTCGACGGCCAGCTGATCCTCGGTCAGTTCGATCTGCAGCCCGAAGGGCGGCGCGACATCGTCCg     <  2:1007225/101‑1 (MQ=255)
                                                                                              gcAGCCGCTGCAGCGCCTCCATCGCGCCCGGGATGTCGACGGCCAGCTGATCCTCGGTCAGTTCGATCTGCAGCCCGAAGGGCGGCGCGACATCGTCCGaa   <  1:2784862/101‑1 (MQ=255)
                                                                                               cAGCCGCTGCAGCGCCTCCATCGCGCCCGGGATGTCGACGGCCAGCTGATCCTCGGTCAGTTCGATCTGCAGCCCGAAGGGCGGCGCGACATCGTCCGAAg  <  2:597226/101‑1 (MQ=255)
                                                                                               cAGCCGCTGCAGCGCCTCCATCGCGCCCGGGATGTCGACGGCCAGCTGATCCTCGGTCAGTTCGATCTGCAGCCCGAAGGGCGGCGCGACATCGTCCGAAg  >  1:1889630/1‑101 (MQ=255)
                                                                                               cAGCCGCTGCAGCGCCTCCATCGCGCCCGGGATGTCGACGGCCAGCTGATCCTCGGTCAGTTCGATCTGCAGCCCGAAGGGCGGCGCGACATCGTCCGAAg  <  2:1068690/101‑1 (MQ=255)
                                                                                               cAGCCGCTGCAGCGCCTCCATCGCGCCCGGGATGTCGACGGCCAGCTGATCCTCGGTCAGTTCGATCTGCAGCCCGAAGGGCGGCGCGACATCGTCCGAAg  >  1:4681646/1‑101 (MQ=255)
                                                                                                  |                                                                                                 
ATCTTCAGCGTCGTCACCGGCATCCGGCGCAGATAGGTCAGGCTGGAATAGCCGGTGCCGAAATCATCGATCGCCAGCGCGATTCCCATCTTGCGCAGGCGCTGCAGCGCCTCCATCGCGCCCGGGATGTCGACGGCCAGCTGATCCTCGGTCAGTTCGATCTGCAGCCCGAAGGGCGGCGCGACATCGTCCGAAG  >  NC_014034/702688‑702883

Base quality scores:  ATCG  < 3 ≤  ATCG  < 20 ≤  ATCG  < 29 ≤  ATCG  < 34 ≤  ATCG  < 41 ≤  ATCG