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% 746.3 198R812R (CGC→CGGRCAP_rcc00643diguanylate cyclase/phosphodiesterase
Reads supporting (aligned to +/- strand):  new base (100/98):  ref base (0/0):  total (100/98)

Only 100 of 210 total aligned reads displayed.

ATCTTCAGCGTCGTCACCGGCATCCGGCGCAGATAGGTCAGGCTGGAATAGCCGGTGCCGAAATCATCGATCGCCAGCGCGATTCCCATCTTGCGCAGGCGCTGCAGCGCCTCCATCGCGCCCGGGATGTCGACGGCCAGCTGATCCTCGGTCAGTTCGATCTGCAGCCCGAAGGGCGGCGCGACATCGTCCGAA  >  NC_014034/702688‑702882
                                                                                                  |                                                                                                
aTCTTCAGCGTCGTCACCGGCATCCGGCGCAGATAGGTCAGGCTGGAATAGCCGGTGCCGAAATCATCGATCGCCAGCGCGATTCCCATCTTGCGCAGccg                                                                                                >  1:2806707/1‑101 (MQ=255)
aTCTTCAGCGTCGTCACCGGCATCCGGCGCAGATAGGTCAGGCTGGAATAGCCGGTGCCGAAATCATCGATCGCCAGCGCGATTCCCATCTTGCGCAGccg                                                                                                >  2:348648/1‑101 (MQ=255)
  cTTCAGCGTCGTCACCGGCATCCGGCGCAGATAGGTCAGGCTGGAATAGCCGGTGCCGAAATCATCGACCGCCAGCGCGATTCCCATCTTGCGCAGCCGCt                                                                                              <  2:792769/101‑1 (MQ=255)
   ttCAGCGTCGTCACCGGCATCCGGCGCAGATAGGTCAGGCTGGAATAGCCGGTGCCGAAATCATCGATCGCCAGCGCGATTCCCATCTTGCGCAGCCGCTg                                                                                             >  1:925583/1‑101 (MQ=255)
   ttCAGCGTCGTCACCGGCATCCGGCGCAGATAGGTCAGGCTGGAATAGCCGGTGCCGAAATCATCGATCGCCAGCGCGATTCCCATCTTGCGCAGCCGCTg                                                                                             <  2:3081836/101‑1 (MQ=255)
    tCAGCGTCGTCACCGGCATCCGGCGCAGATAGGTCAGGCTGGAATAGCCGGTGCCGAAATCATCGATCGCCAGCGCGATTCCCATCTTGCGCAGCCGCTGc                                                                                            >  1:852479/1‑101 (MQ=255)
    tCAGCGTCGTCACCGGCATCCGGCGCAGATAGGTCAGGCTGGAATAGCCGGTGCCGAAATCATCGATCGCCAGCGCGATTCCCATCTTGCGCAGCCGCTGc                                                                                            >  1:3026568/1‑101 (MQ=255)
      aGCGTCGTCACCGGCATCCGGCGCAGATAGGTCAGGCTGGAATAGCCGGTGCCGAAATCATCGATCGCCAGCGCGATTCCCATCTTGCGCAGCCGCTGCAg                                                                                          >  1:3296834/1‑101 (MQ=255)
         gtcgtcACCGGCATCCGGCGCAGATAGGTCAGGCTGGAATAGCCGGTGCCGAAATCATCGATCGCCAGCGCGATTCCCATCTTGCGCAGCCGCTGCAgcgc                                                                                       >  2:3423454/1‑101 (MQ=255)
         gtcgtcACCGGCATCCGGCGCAGATAGGTCAGGCTGGAATAGCCGGTGCCGAAATCATCGATCGCCAGCGCGATTCCCATCTTGCGCAGCCGCTGCAgcgc                                                                                       <  2:2391460/101‑1 (MQ=255)
            gtcACCGGCATCCGGCGCAGATAGGTCAGGCTGGAATAGCCGGTGCCGAAATCATCGATCGCCAGCGCGATTCCCATCTTGCGCAGCCGCTGCAGCGCCTc                                                                                    <  1:1087291/101‑1 (MQ=255)
             tcACCGGCATCCGGCGCAGATAGGTCAGGCTGGAATAGCCGGTGCCGAAATCATCGATCGCCAGCGCGATTCCCATCTTGCGCAGCCGCTGCAGCGCCTcc                                                                                   <  1:4022753/101‑1 (MQ=255)
             tcACCGGCATCCGGCGCAGATAGGTCAGGCTGGAATAGCCGGTGCCGAAATCATCGATCGCCAGCGCGATTCCCATCTTGCGCAGCCGCTGCAGCGCCTcc                                                                                   <  1:2494507/101‑1 (MQ=255)
              cACCGGCATCCGGCGCAGATAGGTCAGGCTGGAATAGCCGGTGCCGAAATCATCGATCGCCAGCGCGATTCCCATCTTGCGCAGCCGCTGCAGCGCCTCCa                                                                                  >  2:2169309/1‑101 (MQ=255)
              cACCGGCATCCGGCGCAGATAGGTCAGGCTGGAATAGCCGGTGCCGAAATCATCGATCGCCAGCGCGATTCCCATCTTGCGCAGCCGCTGCAGCGCCTCCa                                                                                  <  1:98384/101‑1 (MQ=255)
               cccGGCATCCGGCGCAGATAGGTCAGGCTGGAATAGCCGGTGCCGAAATCATCGATCGCCAGCGCGATTCCCATCTTGCGCAGCCGCTGCAGCGCCTCCAt                                                                                 <  1:3723139/100‑1 (MQ=255)
                ccGGCATCCGGCGCAGATAGGTCAGGCTGGAATAGCCGGTGCCGAAATCATCGATCGCCAGCGCGATTCCCATCTTGCGCAGCCGCTGCAGCGCCTCCATc                                                                                >  2:3571706/1‑101 (MQ=255)
                 cGGCATCCGGCGCAGATAGGTCAGGCTGGAATAGCCGGTGCCGAAATCATCGATCGCCAGCGCGATTCCCATCTTGCGCAGCCGCTGCAGCGCCTCCATcg                                                                               <  2:3026343/101‑1 (MQ=255)
                   gCATCCGGCGCAGATAGGTAAGGCTGGAATAGCCGGTGCCGAAATCATCGATCGCCAGCGCGATTCCCATCTTGCGCAGCCGCTGCAGCGCCTCCATcgcg                                                                             <  1:1992950/101‑1 (MQ=255)
                     aTCCGGCGCAGATAGGTCAGGCTGGAATAGCCGGTGCCGAAATCATCGATCGCCAGCGCGATTCCCATCTTGCGCAGCCGCTGCAGCGCCTCCATCGCGcc                                                                           <  1:3881533/101‑1 (MQ=255)
                     aTCCGGCGCAGATAGGTCAGGCTGGAATAGCCGGTGCCGAAATCATCGATCGCCAGCGCGATTCCCATCTTGCGCAGCCGCTGCAGCGCCTCCATCGCGcc                                                                           <  2:3842845/101‑1 (MQ=255)
                     aTCCGGCGCAGATAGGTCAGGCTGGAATAGCCGGTGCCGAAATCATCGATCGCCAGCGCGATTCCCATCTTGCGCAGCCGCTGCAGCGCCTCCATCGCGcc                                                                           >  1:1041267/1‑101 (MQ=255)
                     aTCCGGCGCAGATAGGTCAGGCTGGAATAGCCGGTGCCGAAATCATCGATCGCCAGCGCGATTCCCATCTTGCGCAGCCGCTGCAGCGCCTCCATCGCGcc                                                                           <  1:1844183/101‑1 (MQ=255)
                      tCCGGCGCAGATAGGTCAGGCTGGAATAGCCGGTGCCGAAATCATCGATCGCCAGCGCGATTCCCATCTTGCGCAGCCGCTGCAGCGCCTCCATCGCGccc                                                                          <  1:4099868/101‑1 (MQ=255)
                        cGGCGCAGATAGGTCAGGCTGGAATAGCCGGTGCCGAAATCATCGATCGCCAGCGCGATTCCCATCTTGCGCAGCCGCTGCAGCGCCTCCATCGCGCCCgg                                                                        >  1:3171637/1‑101 (MQ=255)
                         ggCGCAGTTAGGTCAGGCGGGAATAGCCGGTCCCGAAATCATCGATCGCCAGCGCGATTCCCATCTTGCGCAGCCGCTGCAGCGCCTCCACCGCGCCCggg                                                                       >  2:1773347/1‑101 (MQ=255)
                         ggCGCAGATAGGTCAGGCTGGAATAGCCGGTGCCGAAATCATCGATCGCCAGCGCGATTCCCATCTTGCGCAGCCGCTGCAGCGCCTCCATCGCGCCCggg                                                                       >  1:3153269/1‑101 (MQ=255)
                           cgcAGATAGGTCAGGCTGGAATAGCCGGTGCCGAAATCATCGATCGCCAGCGCGATTCCCATCTTGCGCAGCCGCTGCAGCGCCTCCATCGCGCCCGGGAt                                                                     <  2:1432928/101‑1 (MQ=255)
                             cAGATAGGTCAGGCTGGAATAGCCGGTGCCGAAATCATCGATCGCCAGCGCGATTCCCATCTTGCGCAGCCGCTGCAGCGCCTCCATCGCGCCCGGGATGt                                                                   <  2:2323828/101‑1 (MQ=255)
                               gATAGGTCAGGCTGGAATAGCCGGTGCCGAAATCATCGATCGCCAGCGCGATTCCCATCTTGCGCAGCCGCTGCAGCGCCTCCATCGCGCCCGGGATGTCg                                                                 <  2:3894575/101‑1 (MQ=255)
                                       aGGCTGGAATAGCCGGTGCCGAAATCATCGATCGCCAGCGCGATTCCCATCTTGCGCAGCCGCTGCAGCGCCTCCATCGCGCCCGGGATGTCGACGGCCAg                                                         >  1:420703/1‑101 (MQ=255)
                                       aGGCTGGAATAGCCGGTGCCGAAATCATCGATCGCCAGCGCGATTCCCATCTTGCGCAGCCGCTGCAGCGCCTCCATCGCGCCCGGGATGTCGACGGCCAg                                                         >  2:3681952/1‑101 (MQ=255)
                                       aGGCTGGAATAGCCGGTGCCGAAATCATCGATCGCCAGCGCGATTCCCATCTTGCGCAGCCGCTGCAGCGCCTCCATCGCGCCCGGGATGTCGACGGCCAg                                                         >  2:772331/1‑101 (MQ=255)
                                            ggAATAGCCGGTGCCGAAATCATCGATCGCCAGCGCGATTCCCATCTTGCGCAGCCGCTGCAGCGCCTCCATCGCGCCCGGGATGTCGACGGCCAGCTGAt                                                    >  2:1590378/1‑101 (MQ=255)
                                              aaTAGCCGGTGCCGAAATCATCGATCGCCAGCGCGATTCCCATCTTGCGCAGCCGCTGCAGCGCCTCCATCGCGCCCGGGATGTCGACGGCCAGCTGATct                                                  <  2:1856408/101‑2 (MQ=255)
                                              aaTAGCCGGTGCCGAAATCATCGATCGCCAGCGCGATTCCCATCTTGCGCAGCCGCTGCAGCGCCTCCATCGCGCCCGGGATGTCGACGGCCAGCTGATcc                                                  >  2:1277264/1‑101 (MQ=255)
                                               aTAGCCGGTGCCGAAATCATCGATCGCCAGCGCGATTCCCATCTTGCGCAGCCGCTGCAGCGCCTCCATCGCGCCCGGGATGTCGACGGCCAGCTGATCCt                                                 >  2:1341276/1‑101 (MQ=255)
                                               aTAGCCGGTGCCGAAATCATCGATCGCCAGCGCGATTCCCATCTTGCGCAGCCGCTGCAGCGCCTCCATCGCGCCCGGGATGTCGACGGCCAGCTGATCCt                                                 >  2:1145189/1‑101 (MQ=255)
                                                taaccGGTGCCGAAATCATCGATCGCCAGCGCGATTCCCATCTTGCGCAGCCGCTGCAGCGCCTCCATCGCGCCCGGGATGTCGACGGCCAGCTGATCCTc                                                >  1:2420367/4‑101 (MQ=255)
                                                tAGCCGGTGCCGAAATCATCGATCGCCAGCGCGATTCCCATCTTGCGCAGCCGCTGCAGCGCCTCCATCGCGCCCGGGATGTCGACGGCCAGCTGATCCTc                                                <  1:1255152/101‑1 (MQ=255)
                                                tAGCCGGTGCCGAAATCATCGATCGCCAGCGCGATTCCCATCTTGCGCAGCCGCTGCAGCGCCTCCATCGCGCCCGGGATGTCGACGGCCAGCTGATCCTc                                                <  2:2482198/101‑1 (MQ=255)
                                                 aGCCGGTGCCGAAATCATCGATCGCCAGCGCGATTCCCATCTTGCGCAGCCGCTGCAGCGCCTCCATCGCGCCCGGGATGTCGACGGCCAGCTGATCCTCg                                               >  2:1472784/1‑101 (MQ=255)
                                                 aGCCGGTGCCGAAATCATCGATCGCCAGCGCGATTCCCATCTTGCGCAGCCGCTGCAGCGCCTCCATCGCGCCCGGGATGTCGACGGCCAGCTGATCCTCg                                               >  1:1865421/1‑101 (MQ=255)
                                                  gCCGGTGCCGAAACCACCGATCGCCAGCGCGATTCCCACCTTGCGCAGCCGCTGCAGCGCCTCCATCGCGCCCGGGATGTCGACGGCCAGCTGATCCTCgg                                              <  2:2525193/101‑1 (MQ=255)
                                                  gCCGGAGCCGAAATCATCGATCGCCAGCGCGATTCCCATCTTGCGCAGCCGCTGCAGCGCCTCCATCGCGCCCGGGATGTCGACGGCCAGCTGATCCTCgg                                              >  1:285438/1‑101 (MQ=255)
                                                    ctgTGCCGAAATCATCGATCGCCAGCGCGATTCCCATCTTGCGCAGCCGCTGCAGCGCCTCCATCGCGCCCGGGATGTCGACGGCCAGCTGATCCTCggcc                                            >  2:1489228/3‑99 (MQ=255)
                                                     ggTGCCGAAATCATCGATCGCCAGCGCGATTCCCATCTTGCGCAGCCGCTGCAGCGCCTCCATCGCGCCCGGGATGTCGACGGCCAGCTGATCCTCGGTCa                                           >  1:4133020/1‑101 (MQ=255)
                                                     ggTGCCGAAATCATCGATCGCCAGCGCGATTCCCATCTTGCGCAGCCGCTGCAGCGCCTCCATCGCGCCCGGGATGTCGACGGCCAGCTGATCCTCGGTCa                                           >  1:3696837/1‑101 (MQ=255)
                                                     ggTGCCGAAATCATCGATCGCCAGCGCGATTCCCATCTTGCGCAGCCGCTGCAGCGCCTCCATCGCGCCCGGGATGTCGACGGCCAGCTGATCCTCGGTCa                                           >  1:3479495/1‑101 (MQ=255)
                                                     ggTGCCGAAATCATCGATCGCCAGCGCGATTCCCATCTTGCGCAGCCGCTGCAGCGCCTCCATCGCGCCCGGGATGTCGACGGCCAGCTGATCCTCGGTCa                                           >  1:3063757/1‑101 (MQ=255)
                                                      gtcccGAAATCATCGATCGCCAGCGCGATTCCCATCTTGCGCAGCCGCTGCAGCGCCTCCATCGCGCCCGGGATGTCGACGGCCAGCTGATCCTCGGTCAg                                          <  2:2685370/98‑1 (MQ=255)
                                                      gTGCCGAAATCATCGATCGCCAGCGCGATTCCCATCTTGCGCAGCCGCTGCAGCGCCTCCATCGCGCCCGGGATGTCGACGGCCAGCTGATCCTCGGTCAg                                          <  2:1455159/101‑1 (MQ=255)
                                                        gCCGAAATCATCGAACGCCGGCGCGATTCCCATCTTGCGCAGCCGCTGCAGCGCCTCCATCGCGCCCGGGATGTCGACGGCCAGCTGATCCTCGGTCAGtt                                        <  2:4043816/101‑1 (MQ=255)
                                                         ccGAAATCATCGATCGCCAGCGCGATTCCCATCTTGCGCAGCCGCTGCAGCGCCTCCATCGCGCCCGGGATGTCGACGGCCAGCTGATCCTCGGTCAGTTc                                       <  1:3042704/101‑1 (MQ=255)
                                                         ccGAAATCATCGATCGCCAGCGCGATTCCCATCTTGCGCAGCCGCTGCAGCGCCTCCATCGCGCCCGGGATGTCGACGGCCAGCTGATCCTCGGTCAGTTc                                       <  1:1341376/101‑1 (MQ=255)
                                                          cGAAATCATCGATCGCCAGCGCGATTCCCATCTTGCGCAGCCGCTGCAGCGCCTCCATCGCGCCCGGGATGTCGACGGCCAGCTGATCCTCGGTCAGTTCg                                      <  2:1326713/101‑1 (MQ=255)
                                                           gAAATCATCGATCGCCAGCGCGATTCCCATCTTGCGCAGCCGCTGCAGCGCCTCCATCGCGCCCGGGATGTCGACGGCCAGCTGATCCTCGGCCAGTTCGa                                     >  1:123675/1‑101 (MQ=255)
                                                            aaaTCATCGATCGCCAGCGCGATTCCCATCTTGCGCAGCCGCTGCAGCGCCTCCATCGCGCCCGGGATGTCGACGGCCAGCTGATCCTCGGTCAGTTCGAt                                    >  2:3984549/1‑101 (MQ=255)
                                                                catcGATCGCCAGCGCGATTCCCATCTTGCGCAGCCGCTGCAGCGCCTCCATCGCGCCCGGGATGTCGACGGCCAGCTGATCCTCGGTCAGTTCGATCTGc                                >  1:3233734/1‑101 (MQ=255)
                                                                 atcgagcgCCGGCGCGATTCCCATCTTGCGCAGCCGCTGCAGCGCCTCCATCGCGCCCGGGATGTCGACGGCCAGCTGATCCTCGGTAAGTTCGATCTGCa                               <  1:627379/101‑1 (MQ=38)
                                                                 atcgaccgCCAGCGCGATTCCCATCTTGCGCAGCCGCTGCAGCGCCTCCATCGCGCCCGGGATGTCGACGGCCAGCTGATCCTCGGTCAGTTCGATCTGCa                               <  2:3063532/101‑1 (MQ=255)
                                                                    gatcgCCAGCGCGATTCCCATCTTGCGCAGCCGCTGCAGCGCCTCCATCGCGCCCGGGATGTCGACGGCCAGCTGATCCTCGGTCAGTTCGATCTGCAGcc                            >  2:3390/1‑101 (MQ=255)
                                                                    gatcgCCAGCGCGATTCCCATCTTGCGCAGCCGCTGCAGCGCCTCCATCGCGCCCGGGATGTCGACGGCCAGCTGATCCTCGGTCAGTTCGATCTGCAGcc                            <  1:2171089/101‑1 (MQ=255)
                                                                    gatcgCCAGCGCGATTCCCATCTTGCGCAGCCGCTGCAGCGCCTCCATCGCGCCCGGGATGTCGACGGCCAGCTGATCCTCGGTCAGTTCGATCTGCAGcc                            <  1:2758825/101‑1 (MQ=255)
                                                                      tcgCCAGCGCGATTCCCATCTTGCGCAGCCGCTGCAGCGCCTCCATCGCGCCCGGGATGTCGACGGCCAGCTGATCCTCGGTCAGTTCGATCTGCAGCCCg                          <  2:1892014/101‑1 (MQ=255)
                                                                      tcgCCAGCGCGATTCCCATCTTGCGCAGCCGCTCCAGCGCCTCCATCGCGCCCGGGATGTCGACGGCCAGCTGATCCTCGGTCAGTTCGATCTGCAGCCCg                          <  1:3658077/101‑1 (MQ=255)
                                                                      ccgatctcgcgATTCCCATCTTGCGCAGCCGCTGCAGCGCCTCCATCGCGCCCGGGATGTCGACGGCCAGCTGATCCTCGGTCAGTTCGATCTGCAGCCCg                          <  1:775803/94‑1 (MQ=255)
                                                                       cgCCAGCGCGATTCCCATCTTGCGCAGCCGCTGCAGCGCCTCCATCGCGCCCGGGATGTCGACGGCCAGCTGATCCTCGGTCAGTTCGATCTGCAGCCCGa                         >  2:477874/1‑101 (MQ=255)
                                                                        gCCAGCGCGATTCCCATCTTGCGCAGCCGCTGCAGCGCCTCCATCGCGCCCGGGATGTCGACGGCCAGCTGATCCTCGGTCAGTTCGATCTGCAGCCCGaa                        <  2:293364/101‑1 (MQ=255)
                                                                          cAGCGCGATTCCCATCTTGCGCAGCCGCTGCAGCGCCTCCATCGCGCCCGGGATGTCGACGGCCAGCTGATCCTCGGTCAGTTCGATCTGCAGCCCGAAgg                      <  2:2806488/101‑1 (MQ=255)
                                                                          cAGCGCGATTCCCATCTTGCGCAGCCGCTGCAGCGCCTCCATCGCGCCCGGGATGTCGACGGCCAGCTGATCCTCGGTCAGTTCGATCTGCAGCCCGAAgg                      >  1:3959834/1‑101 (MQ=255)
                                                                          cAGCGCGATTCCCATCTTGCGCAGCCGCTGCAGCGCCTCCATCGCGCCCGGGATGTCGACGGCCAGCTGATCCTCGGTCAGTTCGATCTGCAGCCCGAAgg                      <  1:1795408/101‑1 (MQ=255)
                                                                           aGCGCGATTCCCATCTTGCGCAGCCGCTGCAGCGCCTCCATCGCGCCCGGGATGTCGACGGCCAGCTGATCCTCGGTCAGTTCGATCTGCAGCCCGAAggg                     <  2:4140377/101‑1 (MQ=255)
                                                                           aGCGCGATTCCCATCTGGCGCAGCCGCTGCAGCGCCTCCATCGCGCCCGGGATGTCGACGGCCAGCTGATCCTCGGTCAGTTCGATCTGCAGCCCGAAggg                     <  1:1698065/101‑1 (MQ=255)
                                                                            gcgcGATTCCCATCTTGCGCAGCCGCTGCAGCGCCTCCATCGCGCCCGGGATGTCGACGGCCAGCTGATCCTCGGTCAGTTCGATCTGCAGCCCGAAGggc                    >  2:2117813/1‑101 (MQ=255)
                                                                              gcgATTCCCATCTTGCGCAGCCGCTGCAGCGCCTCCATCGCGCCCGGGATGTCGACGGCCAGCTGATCCTCGGTCAGTTCGATCTGCAGCCCGAAGggcgg                  >  1:2610082/1‑101 (MQ=255)
                                                                              gcgATTCCCATCTTGCGCAGCCGCTGCAGCGCCTCCATCGCGCCCGGGATGTCGACGGCCAGCTGATCCTCGGTCAGTTCGATCTGCAGCCCGAAGggcgg                  >  1:2215373/1‑101 (MQ=255)
                                                                              gcgATTCCCATCTTGCGCAGCCGCTGCAGCGCCTCCATCGCGCCCGGGATGTCGACGGCCAGCTGATCCTCGGTCAGTTCGATCTGCAGCCCGAAGggcgg                  >  2:1755674/1‑101 (MQ=255)
                                                                               ggtttCCCATCTTGCGCAGCCGCTGCAGCGCCTCCATCGCGCCCGGGATGTCGACGGCCAGCTGATCCTCGGTCAGTTCGATCTGCAGCCCGAAGggcggc                 <  2:925498/98‑1 (MQ=255)
                                                                               cgATTCCCATCTTGCGCAGCCGCTGCAGCGCCTCCATCGCGCCCGGGATGTCGACGGCCAGCTGATCCTCGGTCAGTTCGATCTGCAGCCCGAAGggcggc                 >  1:3382058/1‑101 (MQ=255)
                                                                                gATTCCCATCTTGCGCAGCCGCTGCAGCGCCTCCATCGCGCCCGGGATGTCGACGGCCAGCTGATCCTCGGTCAGTTCGATCTGCAGCCCGAAGGgcggcg                <  2:3171398/101‑1 (MQ=255)
                                                                                gATTCCCATCTTGCGCAGCCGCTGCAGCGCCTCCATCGCGCCCGGGATGTCGACGGCCAGCTGATCCTCGGTCAGTTCGATCTGCAGCCCGAAGGgcggcg                >  2:3643758/1‑101 (MQ=255)
                                                                                   tCCCATCTTGCGCAGCCGCTGCAGCGCCTCCATCGCGCCCGGGATGTCGACGGCCAGCTGATCCTCGGTCAGTTCGATCTGCAGCCCGAAGGGCGGCGCGa             >  1:3549142/1‑101 (MQ=255)
                                                                                   tCCCATCTTGCCCAGCCGCTGCAGCGCCTCCATCGCGCCCGGGATGTCGACGGCCAGCTGATCCTCGGTCAGTTCGATCTGCAGCCCGAAGGGCGGCGCGa             <  1:4041303/101‑1 (MQ=255)
                                                                                    cccATCTTGCGCAGCCGCTGCAGCGCCTCCATCGCGCCCGGGATGTCGACGGCCAGCTGATCCTCGGTCAGTTCGATCTGCAGCCCGAAGGGCGGCGCGAc            <  2:2753547/101‑1 (MQ=255)
                                                                                    cccATCTTGCGCAGCCGCTGCAGCGCCTCCATCGCGCCCGGGATGTCGACGGCCAGCTGATCCTCGGTCAGTTCGATCTGCAGCCCGAAGGGCGGCGCGAc            >  2:3278379/1‑101 (MQ=255)
                                                                                     ccATCTTGCGCAGCCGCTGCAGCGCCTCCATCGCGCCCGGGATGTCGACGGCCAGCTGATCCTCGGTCAGTTCGATCTGCAGCCCGAAGGGCGGCGCGACa           >  2:670413/1‑101 (MQ=255)
                                                                                     ccATCTTGCGCAGCCGCTGCAGCGCCTCCATCGCGCCCGGGATGTCGACGGCCAGCTGATCCTCGGTCAGTTCGATCTGCAGCCCGAAGGGCGGCGCGACa           <  2:1963150/101‑1 (MQ=255)
                                                                                     ccATCTTGCGCAGCCGCTGCAGCGCCTCCATCGCGCCCGGGATGTCGACGGCCAGCTGATCCTCGGTCAGTTCGATCTGCAGCCCGAAGGGCGGCGCGACa           >  2:2921652/1‑101 (MQ=255)
                                                                                      cATCTTGCGCAGCCGCTGCAGCGCCTCCATCGCGCCCGGGATGTCGACGGCCAGCTGATCCTCGGTCAGTTCGATCTGCAGCCCGAAGGGCGGCGCGACAt          >  2:2010512/1‑101 (MQ=255)
                                                                                      cATCTTGCGCAGCCGCTGCAGCGCCTCCATCGCGCCCGGGATGTCGACGGCCAGCTGATCCTCGGTCAGTTCGATCTGCAGCCCGAAGGGCGGCGCGACAt          >  1:159213/1‑101 (MQ=255)
                                                                                       aTCTTGCGCAGCCGCTGCAGCGCCTCCATCGCGCCCGGGATGTCGACGGCCAGCTGATCCTCGGTCAGTTCGATCTGCAGCCCGAAGGGCGGCGCGACATc         <  1:1455995/101‑1 (MQ=255)
                                                                                        tCTTGCGCAGCCGCTGCAGCGCCTCCATCGCGCCCGGGATGTCGACGGCCAGCTGATCCTCGGTCAGTTCGATCTGCAGCCCGAAGGGCGGCGCGACATCg        <  1:1487472/101‑1 (MQ=255)
                                                                                        ccTTGCGCAGCCGCTGCAGCGCCTCCATCGCGCCCGGGATGTCGACGGCCAGCTGATCCTCGGTCAGTTCGATCTGCAGCCCGAAGGGCGGCGCGACATCg        <  2:1849347/100‑1 (MQ=255)
                                                                                         cTTGCGCAGCCGCTGCAGCGCCTCCATCGCGCCCGGGATGTCGACGGCCAGCTGATCCTCGGTCAGTTCGATCTGCAGCCTGAAGGGCGGCGCGACATCGt       >  1:2951699/1‑101 (MQ=255)
                                                                                         cTTGCGCAGCCGCTGCAGCGCCTCCATCGCGCCCGGGATGTCGACGGCCAGCTGATCCTCGGTCAGTTCGATCTGCAGCCCGAAGGGCGGCGCGACATCGt       >  1:2500149/1‑101 (MQ=255)
                                                                                          ttGCGCAGCCGCTGCAGCGCCTCCATCGCGCCCGGGATGTCGACGGCCAGCTGATCCTCGGTCAGTTCGATCTGCAGCCCGAAGGGCGGCGCGACATCGTc      <  1:1738625/101‑1 (MQ=255)
                                                                                          ttGCGCAGCCGCTGCAGCGCCTCCATCGCGCCCGGGATGTCGACGGCCAGCTGATCCTCGGTCAGTTCGATCTGCAGCCCGAAGGGCGGCGCGACATCGTc      <  1:2580147/101‑1 (MQ=255)
                                                                                           tGCGCAGCCGCTGCAGCGCCTCCATCGCGCCCGGGATGTCGACGGCCAGCTGATCCTCGGTCAGTTCGATCTGCAGCCCGAAGGGCGGCGCGACATCGTcc     <  1:4189171/101‑1 (MQ=255)
                                                                                             cgcAGCCGCTGCAGCGCCTCCATCGCGCCCGGGATGTCGACGGCCAGCTGATCCTCGGTCAGTTCGATCTGCAGCCCGAAGGGCGGCGCGACATCGTCCGa   <  2:2197757/101‑1 (MQ=255)
                                                                                                  |                                                                                                
ATCTTCAGCGTCGTCACCGGCATCCGGCGCAGATAGGTCAGGCTGGAATAGCCGGTGCCGAAATCATCGATCGCCAGCGCGATTCCCATCTTGCGCAGGCGCTGCAGCGCCTCCATCGCGCCCGGGATGTCGACGGCCAGCTGATCCTCGGTCAGTTCGATCTGCAGCCCGAAGGGCGGCGCGACATCGTCCGAA  >  NC_014034/702688‑702882

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