Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_014034 2,821,448 +G 2.7% coding (197/1062 nt) RCAP_rcc02629 ← diguanylate cyclase/phosphodiesterase

Read alignment evidence...
  seq id position change freq score reads annotation genes product
*NC_0140342,821,4481.→G 2.7% 5.1 ‑3.8 ‑0.0147coding (197/1062 nt)RCAP_rcc02629diguanylate cyclase/phosphodiesterase
Reads supporting (aligned to +/- strand):  new base (4/0):  ref base (54/89):  total (58/89)
Fisher's exact test strand distribution p-value = 2.27e-02
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 9.76e-01

Only 100 of 179 total aligned reads displayed.

GAACATCGGATAGCGCACCAGCGTCATCACCAGCCGCCGCCCGGTCAGCGCGACCTCGCGCCCGGG.CTCGACGCCCCGGGCCCGGCGCAGGGTGAAATGGC.GGGCG..AAATCCGCCCCCTCCGCCCGCTCCGGCCCCTGACCGACGATCCGCGCCAGCGTCGGCCCCATGGCGCGGATCCGCCCGCGCGCATCCAGCCACAG  >  NC_014034/2821348‑2821548
                                                                                                      |                                                                                                      
gAACATCGGATAGCGCACCAGCGTCATCACCAGCCGCCGCCCGGTCAGCGCGACCTCGCGCCCGGG.CTCGACGCCCCGGGCCCGGCGCAGGGTGAAATggc                                                                                                         >  1:1596789/1‑101 (MQ=255)
gAACATCGGATAGCGCACCAGCGTCATCACCAGCCGCCGCCCGGTCAGCGCGACCTCGCGCCCGGG.CTCGACGCCCCGGGCCCGGCGCAGGGTGAAATggc                                                                                                         >  1:3302281/1‑101 (MQ=255)
 aaCATCGGATAGCGCACCAGCGTCATCACCAGCCGCCGCCCGGTCAGCGCGACCTCGCGCCCGGG.CTCGACGCCCCGGGCCCGGCGCAGGGTGAAATggc.g                                                                                                       >  1:2522922/1‑101 (MQ=255)
  aCATCGGATAGCGCACCAGCGTCATCACCAGCCGCCGCCCGGTCAGCGCGACCTCGCGCCCGGG.CTCGACGCCCCGGGCCCGGCGCAGGGTGAAATggc.gg                                                                                                      >  2:1151208/1‑101 (MQ=255)
   cATCGGATAGCGCACCAGCGTCATCACCAGCCGCCGCCCGGTCAGCGCGACCTCGCGCCCGGG.CTCGACGCCCCGGGCCCGGCGCAGGGTGAAATggc.ggg                                                                                                     >  2:2579509/1‑101 (MQ=255)
   aaTCGGATAGCGCCCCAGCGTCCTCACCAGCCGCCGCCCGGTCGGCGCGACCCCGCGCCCGGGCCTCGACGCCCCGGGCCCGGCGCAGGGTGAAATggc.gg                                                                                                      <  2:1670258/100‑1 (MQ=255)
     tCGGATAGCGCACCAGCGTCATCACCAGCCGCCGCCCGGTCAGCGCGACCTCGCGCCCGGG.CTCGACGCCCCGGGCCCGGCGCAGGGTGAAATggc.gggcg                                                                                                   <  1:3081310/101‑1 (MQ=255)
      cGGATAGCGCACCAGCGTCATCACCAGCCGCCGCCCGGTCAGCGCGACCTCGCGCCCGGG.CTCGACGCCCCGGGCCCGGCGCAGGGTGAAATGGC.GGGCG..a                                                                                                >  2:4662820/1‑101 (MQ=255)
       ggATAGCGCACCAGCGTCATCACCAGCCGCCGCCCGGTCAGCGCGACCTCGCGCCCGGG.CTCGACGCCCCGGGCCCGGCGCAGGGTGAAATGGC.GGGCG..aa                                                                                               <  2:354154/101‑1 (MQ=255)
        gATAGCGCACCAGCGTCATCACCAGCCGCCGCCCGGTCAGCGCGACCTCGCGCCCGGG.CTCGACGCCCCGGGCCCGGCGCAGGGTGAAATGGC.GGGCG..gaa                                                                                              >  1:2885372/1‑101 (MQ=255)
         aTAGCGCACCAGCGTCATCACCAGCCGCCGCCCGGTCAGCGCGACCTCGCGCCCGGG.CTCGACGCCCCGGGCCCGGCGCAGGGTGAAATggc.ggggggaaa                                                                                               >  2:2087417/1‑95 (MQ=255)
          tAGCGCACCAGCGTCATCACCAGCCGCCGCCCGGTCAGCGCGACCTCGCGCCCGGG.CTCGACGCCCCGGGCCCGGCGCAGGGTGAAATGGC.GGGCG..AAATc                                                                                            <  2:4326705/101‑1 (MQ=255)
             cgcACCAGCGTCATCACCAGCCGCCGCCCGGTCAGCGCGACCTCGCGCCCGGG.CTCGACGCCCCGGGCCCGGCGCAGGGTGAAATGGC.GGGCG..GAAACCGc                                                                                         >  2:4627997/1‑101 (MQ=255)
             cgcACCAGCGTCATCACCAGCCGCCGCCCGGTCAGCGCGACCTCGCGACCGGG.CTCGACGCCCCGCGCCCGGCGCAGGGTGAAATggc.ggggggagacgac                                                                                           >  1:2347689/1‑91 (MQ=255)
              gcACCAGCGTCATCACCAGCCGCCGCCCGGTCAGCGCGACCTCGCGCCCGGG.CTCGACGCCCCGGGCCCGGCGCAGGGTGAAATGGC.GGGCG..AAATCCGcc                                                                                        <  2:2155560/101‑1 (MQ=255)
              gcACCAGCGTCATCACCAGCCGCCGCCCGGTCAGCGCGACCTCGCGCCCGGG.CTCGACGCCCCGGGCCCGGCGCAGGGTGAAATGGC.GGGCG..AAATCCGcc                                                                                        <  1:3849139/101‑1 (MQ=255)
                aCCAGCGTCATCACCAGCCGCCGCCCGGTCAGCGCGACCTCGCGCCCGGG.CTCGACGCCCCGGGCCCGGCGCAGGGTGAAATGGC.GGGCG..AAATCCGcccc                                                                                      <  1:4340364/101‑1 (MQ=255)
                aCCAGCGTCATCACCAGCCGCCGCCCGGTCAGCGCGACCTCGCGCCCGGG.CTCGACGCCCCGGGCCCGGCGCAGGGTGAAATGGC.GGGCG..AAATCCGcccc                                                                                      <  2:4327462/101‑1 (MQ=255)
                  cAGCGTCATCACCAGCCGCCGCCCGGTCAGCGCGACCTCGCGCCCGGG.CTCGACGCCCCGGGCCCGGCGGAGGGTGAAATGGC.GGGCG..AAATCCGCCCCCt                                                                                    <  1:4295947/101‑1 (MQ=255)
                  cAGCGTCATCACCAGCCGCCGCCCGGTCAGCGCGACCTCGCGCCCGGG.CTCGACGCCCCGGGCCCGGCGCAGGGTGAAATGGC.GGGCG..AAATCCGCCCCCt                                                                                    <  1:955151/101‑1 (MQ=255)
                    gCGTCATCACCAGCCGCCGCCCGGTCAGCGCGACCTCGCGCCCGGG.CTCGACGCCCCGGGCCCGGCGCAGGGTGAAATGGC.GGGCG..AAATCCGCCCCCTcc                                                                                  <  2:2633435/101‑1 (MQ=255)
                    gCGTCATCACCAGCCGCCGCCCGGTCAGCGCGACCTCGCGCCCGGG.CTCGACGCCCCGGGCCCGGCGCAGGGTGAAATGGC.GGGCG..AAATCCGCCCCCTcc                                                                                  <  1:4274158/101‑1 (MQ=255)
                    gCGTCATCACCAGCCGCCGCCCGGTCAGCGCGACCTCGCGCCCGGG.CTCGACGCCCCGGGCCCGGCGCAGGGTGAAATGGC.GGGCG..AAATCCGCCCCCTcc                                                                                  <  1:4109627/101‑1 (MQ=255)
                    gCGTCATCACCAGCCGCCGCCCGGTCAGCGCGACCTCGCGCCCGGG.CTCGACGCCCCGGGCCCGGCGCAGGGTGAAATGGC.GGGCG..AAATCCGCCCCCTcc                                                                                  <  2:615715/101‑1 (MQ=255)
                     cGTCATCACCAGCCGCCGCCCGGTCAGCGCGACCTCGCGCCCGGG.CTCGACGCCCCGGGCCCGGCGCAGGGTGAAATGGC.GGGCG..AAATCCGCCCCCTccg                                                                                 <  2:4006184/101‑1 (MQ=255)
                      gTCATCACCAGCCGCCGCCCGGTCAGCGCGACCTCGCGCCCGGG.CTCGACGCCCCGGGCCCGGCGCAGGGTGAAATGGC.GGGGG..GAAACCGCCCCCCcccc                                                                                >  1:1000552/1‑99 (MQ=255)
                      gTCATCACCAGCCGCCGCCCGGTCAGCGCGACCTCGCGCCCGGG.CTCGACGCCCCGGGCCCGGCGCAGGGGGAAATGGCGGGGGG..AAAACCCCCCCCCccc                                                                                 >  1:2057326/1‑100 (MQ=255)
                       tcatcaCCAGCCGCCGCCCGGTCAGCGCGACCTCGCGCCCGGG.CTCGACGCCCCGGGCCCGGCGCAGGGTGAAATGGC.GGGCG..GAAACCCCCCCCCccccc                                                                               >  2:2315635/1‑101 (MQ=255)
                        catcaCCAGCCGCCGCCCGGTCAGCGCGACCTCGCGCCCGGG.CTCGACGCCCCGGGCCCGGCGCAGGGTGAAATGGC.GGGCG..AAATCCGCCCCCTccgccc                                                                              <  2:3374998/101‑1 (MQ=255)
                        catcaCCAGCCGCCGCCCGGTCAGCGCGACCTCGCGCCCGGG.CTCGACGCCCCGGGCCCGGCGCAGGGTGAAATGGC.GGGCG..AAATCCGCCCCCTccgccc                                                                              <  2:3836046/101‑1 (MQ=255)
                         atcaCCAGCCGCCGCCCGGTCAGCGCGACCTCGCGCCCGGG.CTCGACGCCCCGGGCCCGGCGCAGGGTGAAATGGC.GGGCG..AAATCCGCCCCCTccgcccg                                                                             <  1:270087/101‑1 (MQ=255)
                          tcaCCAGCCGCCGCCCGGTCAGCGCGACCTCGCGCCCGGG.CTCGACGCCCCGGGCCCGGCGCAGGGTGAAATGGC.GGGCG..GAAAACGCCCCCCccgccccc                                                                            >  2:2487757/1‑99 (MQ=255)
                           caCCAGCCGCCGCCCGGTCAGCGCGACCTCGCGCCCGGG.CTCGACGCCCCGGGCCCGGCGCAGGGTGAAATGGC.GGGCG..AAATCCGCCCCCTCCGCCCGCt                                                                           <  1:2246314/101‑1 (MQ=255)
                           caCCAGCCGCCGCCCGGTCAGCGCGACCTCGCGCCCGGG.CTCGACGCCCCGGGCCCGGCGCAGGGTGAAATGGC.GGGCG..AAAACCGCCCCCCccccccccc                                                                           >  2:3623202/1‑98 (MQ=255)
                            cccAGCCGCCGCCCGGTCAGCGCGACCTCGCGCCCGGG.CTCGACGCCCCGGGCCCGGCGCAGGGTGAAATGGC.GGGCG..AAATCCGCCCCCTCCGCCCGCTc                                                                          <  1:2813791/100‑1 (MQ=255)
                                  cgccgccCGGTCAGCGCGACCTCGCGCCCGGG.CTCGACGCCCCGGGCCCGGCGCAGGGTGAAATGGC.GGGCG..AAATCCGCCCCCTCCGCCCGCTCCGGccc                                                                    <  1:1836275/101‑1 (MQ=255)
                                   gccgccCGGTCAGCGCGACCTCGCGCCCGGG.CTCGACGCCCCGGGCCCGGCGCAGGGTGAAATGGC.GGGCG..AAATCCGCCCCCTCCGCCCGCTCCGGcccc                                                                   <  1:3546917/101‑1 (MQ=255)
                                    ccgccCGGTCAGCGCGACCTCGCGCCCGGG.CTCGACGCCCCGGTCCCGGCGCAGGGTGAAATGGC.GGGCG..AAATCCGCCCCCTCCGCCCGCTCCGGCCCCt                                                                  <  1:323928/101‑1 (MQ=255)
                                    ccgccCGGTCAGCGCGACCTCGCGCCCGGG.CTCGACGCCCCGGGCCCGGCGCAGGGTGAAATGGC.GGGCG..AAATCCGCCCCCTCCGCCCGCTCCGGCCCCt                                                                  <  1:1332910/101‑1 (MQ=255)
                                       cccGGTCAGCGCGACCTCGCGCCCGGG.CTCGACGCCCCGGGCCCGGCGCAGGGTGAAATGGC.GGGCG..AAATCCGCCCCCTCCGCCCGCTCCGGCCCCTGAc                                                               <  1:4511689/101‑1 (MQ=255)
                                        ccGGTCAGCGCGACCTCGCGCCCGGG.CTCGACGCCCCGGGCCCGGCGCAGGGTGAAATGGC.GGGGG..GAAAGCCCCCCCCCCCCCCGCCCCCGccccccgcc                                                              >  1:583949/1‑96 (MQ=255)
                                        ccGGCCAGCGCGCCCTCGCGCCCGGG.CTCGACGCCCCGGGCCCGGCGCAGGGTGAAATGGC.GGGCA..AAATCCGCCCCCTCCGCCCGCTCCGGCCCCTGAcc                                                              <  1:1608491/101‑1 (MQ=255)
                                         cGGTCAGCGCGACCTCGCGCCCGGG.CTCGACGCCCCGGGCCCGGCGCAGGGTGAAATGGC.GGGCG..AAATCCGCCCCCTCCGCCCGCTCCGGCCCCTGACcg                                                             <  1:1546304/101‑1 (MQ=255)
                                         cGGTCAGCGCGACCTCGCGCCCGGG.CTCGACGCCCCGGGCCCGGCGCAGGGTGAAATGGC.GGGCG..AAATCCGCCCCCTCCGCCCGCTCCGGCCCCTGACcg                                                             <  2:4351889/101‑1 (MQ=255)
                                         cGGTCAGCGCGACCTCGCGCCCGGG.CTCGACGCCCCGGGCCCGGCGCAGGGTGAAATGGC.GGGCG..AAATCCGCCCCCTCCGCCCGCTCCGGCCCCTGACcg                                                             <  2:329165/101‑1 (MQ=255)
                                            tCAGCGCGACCTCGCGCCCGGG.CTCGACGCCCCGGGCCCGGCGCAGGGTGAAATGGC.GGGCG..AAATCCGCCCCCTCCGCCCGCTCCGGCCCCTGACcgacg                                                          <  2:1897411/101‑1 (MQ=255)
                                             cAGCGCGACCTCGCGCCCGGG.CTCGACGCCCCGGGCCCGGCGCAGGGTGAAATGGC.GGGGG..GAAAACGCCCCCCCCCCCCCCCCCCGcccccgccccccca                                                         >  2:2467853/1‑91 (MQ=255)
                                             cAGCGCGACCTCGCGCCCGGG.CTCGACGCCCCGGGCCCGGCGCAGGGTGAAATGGA.GGGCG..GAAAACGCCCCCTCCCCCCCCCCCGGccccccccccaccc                                                         >  1:1278629/1‑91 (MQ=255)
                                              cgcgcGACCTCGCGCCCGGG.CTCGACGCCCCGGGCCCGGCGCAGGGTGAAATGGC.GGGCG..AAATCCGCCCCCTCCGCCCGCTCCGGCCCCTGACCGACGAt                                                        <  2:1561498/100‑1 (MQ=255)
                                              aGCGCGACCTCGCGCCCGGG.CTCGACGCCCCGGGCCCGGCGCAGGGTGAAATGGC.GGGCG..AAATCCGCCCCCTCCGCCCGCTCCGGCCCCTGACCGACGAt                                                        <  2:3538558/101‑1 (MQ=255)
                                               gcgcGACCTCGCGCCCGGG.CTCGACGCCCCGGGCCCGGCGCAGGGTGAAATGGC.GGGCG..AAATCCGCCCCCTCCGCCCGCTCCGGCCCCTGACCGACGATc                                                       <  1:4725252/101‑1 (MQ=255)
                                                cgcgACCTCGCGCCCGGG.CTCGACGCCCCGGGCCCGGCGCAGGGTGAAATGGC.GGGCG..AAATCCGCCCCCTCCGCCCGCTCCGGCCCCTGACCGACGATcc                                                      <  2:736257/101‑1 (MQ=255)
                                                    aCCTCGCGCCCGGG.CTCGACGCCCCGGGCCCGGCGCAGGGTGAAATGGC.GGGCG..AAAACCGCCCCCCCCCCCCCCCCCCGCCCCCCACCCACCACCCgcgc                                                  >  2:1826309/1‑101 (MQ=255)
                                                      cTCGCGCCCGGG.CTCGACGCCCCGGGCCCGGCGCAGGGTGAAATGGC.GGGCG..AAATCCGCCCCCTCCGCCCGCTCCGGCCCCTGACCGACGATCCGCGCCa                                                <  2:4305491/101‑1 (MQ=255)
                                                       tCGCGCCCGGG.CTCGACGCCCCGGGCCCGGCGCAGGGTGAAATGGC.GGGGG..GAAACCGCCCCCCCCGCCCCCTCCCGCCCCCCACCCACCAGCCGCGCCAg                                               >  2:2568851/1‑101 (MQ=255)
                                                         gcgcCCGGG.CTCGACGCCCCGGGCCCGGCGCAGGGTGAAATGGC.GGGCG..AAATCCGCCCCCTCCGCCCGCTCCGGCCCCTGACCGACGATCCGCGCCAGCg                                             <  1:4199981/101‑1 (MQ=255)
                                                           gcCCGGG.CTCGACGCCCCGGGCCCGGCGCAGGGTGAAATGGC.GGGGG..AAAACCGCCCCCCCCCCCCCCCCCCGCCCCCGCCCCACCATCCGCGCCCGCgcc                                           >  1:1512671/1‑99 (MQ=255)
                                                           gcCCGGG.CTCGACGCCCCGGGCCCGGCGCAGGGTGAAATGGC.GGGCG..AAATCCGCCCCCTCCGCCCGCTCCGGCCCCTGACCGACGATCCGCGCCAGCGTc                                           <  2:2622448/101‑1 (MQ=255)
                                                              cGGG.CTCGACGCCCCGGGCCCGGCGCAGGGTGAAATGGC.GGGGG..GAAACCGCCCCCCCCCCCCCCCCCCGCCCCCCACCCACGATCCCCGCCCCCGTCGGc                                        >  2:3917067/1‑101 (MQ=255)
                                                               ggg.CTCGACGCCCCGGGCCCGGCGCAGGGTGAAATGGC.GGGGGGAAAAACCGCCCCC.CC.CCCCCTTCGGCCCCCCACCCCCAAACCCCGCCAGCAGCGGcc                                       >  1:3900018/1‑101 (MQ=255)
                                                               ggg.CTCGACGCCCCGGGCCCGGCGCAGGGTGAAATGGC.GGGCG..AAATCCGCCCCCTCCGCCCGCTCCGGCCCCTGACCGACGATCCGCGCCAGCGTCGGcc                                       <  1:2487901/101‑1 (MQ=255)
                                                                 g.CTCGACGCCCCGGGCCCGGCGCAGGGTGAAATGGC.GGGCG..AAATCCGCCCCCTCCGCCCGCTCCGGCCCCTGACCGACGATCCGCGCCAGCGTCGGcccc                                     <  1:29422/101‑1 (MQ=255)
                                                                   cTCGACGCCCCGGGCCCGGCGCAGGGTGAAATGGC.GGGGG..AAAAACCCCCCCCCCCCCCCCCCCGGCCCCCGACCCCCGATCCGCGCCCCCCGCGGccccc                                    >  1:103252/1‑100 (MQ=255)
                                                                    tCGACGCCCCGGGCCCGGCGCAGGGTGAAATGGC.GGGCG..AAATCCGCCCCCTCCGCCCGCTCCGGCCCCTGACCGACGATCCGCGCCAGCGTCGGCCCCAt                                   <  2:3709293/101‑1 (MQ=255)
                                                                     cGACGCCCCGGGCCCGGCGCAGGGTGAAATGGC.GGGCG..AAATCCGCCCCCTCCGCCCGCTCCGGCCCCTGACCGACGATCCGCGCCAGCGTCGGCCCCATg                                  <  1:3380863/101‑1 (MQ=255)
                                                                       aCGCCCCGGGCCCGGCGCAGGGTGAAATGGC.GGGCG..AAATCCGCCCCCTCCGCCCGCTCCGGCCCCTGACCGACGATCCGCGCCAGCGTCGGCCCCATGgc                                <  2:1000473/101‑1 (MQ=255)
                                                                        cGCCCCGGGCCCGGCGCAGGGTGAAATGGC.GGGCG..AAATCCGCCCCCTCCGCCCGCTCCGGCCCCTGACCGACGATCCGCGCCAGCGTCGGCCCCATGgcg                               <  2:1596689/101‑1 (MQ=255)
                                                                         gCCCCGGGCCCGGCGCAGGGTGAAATGGC.GGGGG..GAAACCGCCCCCCCCCCCCCCCCCCGCCCCCCACCCCCGACACCCGCCACCGCCGGCCCCCCGgcgc                              >  2:1640089/1‑101 (MQ=255)
                                                                             cGGGCCCGGCGCAGGGTGAAATGGC.GGGGG..AAAACCCCCCCCCCCCCCCCCTCCCGCCCCCGACCGACGATTCCCCCCCGCCTCCGCCCCCTCGCCCgggc                          >  2:1430199/1‑99 (MQ=255)
                                                                             cGGGCCCGGCGCAGGGTGAAATGGC.GGGCG..GAAACCCCCCCCCCCCCCCCCCCCGGCCCCCTCCCCCCGACCCCCGCCCGCGGCGGCCCCCTCGGGCggct                          >  1:3428453/1‑99 (MQ=255)
                                                                             cGGGCCCGGCGCAGGGTGAAAGGGC.GGGCG..AAATCCGCCCCCCCCGCCCGCTCCGGCCCCTGACCGACGATCCGCGCCAGCGTCGGCCCCATGGCGCGGAt                          <  1:891147/101‑1 (MQ=255)
                                                                              gggCCCGGCTCACGGTGAAATGGC.GGGCG..AAATCCGCCCCCTCCGCCCGCTCCGGCCCCTGACCGACGATCCGCGCCAGCGTCGGCCCCATGGCGCGGATc                         <  2:1468686/101‑1 (MQ=255)
                                                                              gggCCCGGCGCAGGGTGAAATGGC.GGGCG..GAAACCGCCCCCCCCCCCCCCCCCCGCCCCCGACCGACGATCCCCGCCCGCGGCGGCCCCCCGGGGCggggc                         >  1:4144838/1‑98 (MQ=255)
                                                                              gggCCCGGCGCAGGGTGAAATGGC.GGGCG..AAATCCGCCCCCTCCGCCCGCTCCGGCCCCTGACCGACGATCCGCGCCAGCGTCGGCCCCATGGCGCGGATc                         <  1:4309523/101‑1 (MQ=255)
                                                                              gggCCCGGCGCAGGGTGAAATGGC.GGGCG..AAATCCCCCCCCCCCCCCCCCCCCCGCCCCCCACCCACGATCCGCGCCCCCGCCCGCCCCCCGGCGCgggcc                         >  1:2225451/1‑98 (MQ=255)
                                                                              gggCCCGGCGCAGGGTGAAATGGC.GGGCG..AAAAACGCCCCCCCCCCCCCCCCCCGCCCCCGACCGACGAACCCCGCCCGCGCCGGCCCCCCGGGGCGGagc                         >  1:3584115/1‑99 (MQ=255)
                                                                                  ccGGCGCAGGGTGAAATGGC.GGGCG..AAAACCCCCCCCCCCCCCCCCCCCCGCCCCCGACCCACGATCCGCGCCAGCGTCGGCCCCCCCGCGCggggccccc                     >  1:31671/1‑94 (MQ=255)
                                                                                  ccGGCGCAGGGTGAAATGGC.GGGCG..AAAAACGCCCCCCCCGCCCGCCCCGGCCCCCGACCGACGATCCGCGCCCGCGTCGGCCCCCAGGCGCGGACccccc                     >  1:1220371/1‑101 (MQ=255)
                                                                                   cGGCGCAGGGTGAAATGGC.GGGCG..AAATCCGCCCCCTCCGCCCGCTCCGGCCCCTGACCGACGATCCGCGCCAGCGTCGGCCCCATGGCGCGGATccgccc                    <  1:2710822/101‑1 (MQ=255)
                                                                                   cGGCGCAGGGTGAAATGGC.GGGCG..AAATCCGCCCCCTCCGCCCGCTCCGGCCCCTGACCGACGATCCGCGCCAGCGTCGGCCCCATGGCGCGGATccgccc                    <  1:2134524/101‑1 (MQ=255)
                                                                                     gcgcAGGGTGAAATGGC.GGGCG..AAAACCGCCCCCCCCCCCCCCCCCCGCCCCCGACCGACGATCCGCGCCCGCGGCGGCCCCCCGGCGCGGGGccgccccc                  >  1:4278009/1‑99 (MQ=255)
                                                                                      cgcAGGGTGAAATGGC.GGGCG..AAATCCGCCCCCTCCGCCCGCTCCGGCCCCTGACCGACGATCCGCGCCAGCGTCGGCCCCATGGCGCGGATCCGCCcgcg                 <  2:3101965/101‑1 (MQ=255)
                                                                                       gcAGGGTGAAATGGC.GGGCG..AAATCCGCCCCCTCCGCCCGCTCCGGCCCCTGACCGACGATCCGCGCCAGCGTCGGCCCCATGGCGCGGATCCGCCcgcgc                <  1:778366/101‑1 (MQ=255)
                                                                                        cAGGGTGAAATGGC.GGGCG..AAATCCGCCCCCTCCCCCCGCTCCGGCCCCCGACCCACGATCCGCGCCAGCGTCGGCCCCCTGGCGCGGGTCCGCCcgcgcg               >  2:1923454/1‑101 (MQ=255)
                                                                                         aGGGTGAAATGGC.GGGCG..AAATTCCCCCCCTCCGCCCCCTCCCGGCCCCGACCGGCGATCCGCGCCAGCGCCGGCCCCATGGCGCGGGACCGCCcgcgccc              >  2:2373472/1‑99 (MQ=255)
                                                                                          gggTGAAATGGC.GGGCG..AAATCCGCCCCCTCCGCCCGCTCCGGCCCCTGACCGACGATCCGCGCCAGCGTCGGCCCCATGGCGCGGATCCGCCCGCGCGCa             <  1:1731112/101‑1 (MQ=255)
                                                                                           gttGAAATGGC.GGCCT..AAATCCCCCCCCTCCGCCCGCCCCGGCCCCTGACCGACGACCCGCGCCAGCGTCGGCCCCATGGCGCGGATCCGCCCGCGCGCAt            <  2:412943/99‑1 (MQ=255)
                                                                                             tGAAATGGC.GGGCG..AAATCCGCCCCCTCCGCCCGCTCCGGCCCCTGACCGACGATCCGCGCCAGCGTCGGCCCCATGGCGCGGATCCGCCCGCGCGCATcc          <  1:4344775/101‑1 (MQ=255)
                                                                                             tGAAATGGC.GGGCG..AAATCCGCCCCCCCCCCCCGCTCCGGCCCCCCACCCACCATCCCCGCCCGCGTCGGCCCCCTTGCGCGGGTccgccccccccccccc          >  2:3393876/1‑91 (MQ=255)
                                                                                              gAAATGGC.GGGCG..AAATCCGCCCCCTCCGCCCGCTCCGGCCCCTGACCGACGATCCGCGCCAGCGTCGGCCCCATGGCGCGGATCCGCCCGCGCGCATCCa         <  2:3337805/101‑1 (MQ=255)
                                                                                              gAAATGGC.GGGCG..AAATCCGCCCCCTCCGCCCGCTCCGGCCCCTGACCGACGATCCGCGCCAGCGTCGGCCCCATGGCGCGGATCCGCCCGCGCGCATCCa         <  1:1740581/101‑1 (MQ=255)
                                                                                              gAAATGGC.GGGCG..AAATCCGCCCCCTCCGCCCGCTCCGGCCCCTGACCGACGATCCGCGCCAGCGTCGGCCCCATGGCGCGGATCCGCCCGCGCGCATCCa         <  2:4406141/101‑1 (MQ=255)
                                                                                              gAAATGGC.GGGCG..AAAACCGCCCCCCCCCCCCCCTCCGGCCCCCGAGCGACGATCCGCGCCAGCGTCGGGCCCCTGGCGCGGGTCCCCCCgcgcgccgcca         >  2:2845658/1‑96 (MQ=255)
                                                                                               aaaTGGC.GGGCG..AAATCCGCCCCCTCCGCCCGCTCCGGCCCCTGACCGACGATCCGCGCCAGCGTCGGCCCCATGGCGCGGATCCGCCCGCGCGCATCCAg        <  2:907640/101‑1 (MQ=255)
                                                                                                aaTGGC.GGGCG..AAATCCGCCCCCTCCGCCCGCTCCGGCCCCTGACCGACGATCCGCGCCAGCGTCGGCCCCATGGCGCGGATCCGCCCGCGCGCATCCAGc       <  2:2798634/101‑1 (MQ=255)
                                                                                                 aTGGC.GGGCG..AAATCCGCCCCCTCCGCCCGCTCCGGCCCCTGACCGACGATCCGCGCCAGCGTCGGCCCCATGGCGCGGATCCGCCCGCGCGCATCCAGcc      <  1:3869973/101‑1 (MQ=255)
                                                                                                  tGGC.GGGCG..AAATCCGCCCCCTCCGCCCGCTCCGGCCCCTGACCGACGATCCGCGCCAGCGTCGGCCCCATGGCGCGGATCCGCCCGCGCGCATCCAGCca     <  2:12619/101‑1 (MQ=255)
                                                                                                  tGGC.GGGCG..AAATCCGCCCCCTCCGCCCGCTCCGGCCCCTGACCGACGATCCGCGCCAGCGTCGGCCCCATGGCGCGGATCCGCCCGCGCGCATCAAGCca     >  2:1238607/1‑101 (MQ=255)
                                                                                                   ggc.gggcg..AAACCCGCCCCCTCCGCCCGCTCCGGCCCCTGACCGACGATCCGCGCCAGCGTCGGCCCCATGGCGCGGATCCGCCCGCGCGCATCCAGCcac    <  2:4378887/101‑1 (MQ=255)
                                                                                                     c.gggcg..AAATCCGCCCCCTCCGCCCGCTCCGGCCCCTGACCGACGATACGCGCCAGCGTCAGCCCCATGGCGCGGATCCGCCCGCGCGCATCCAGCCACAg  >  2:3925896/1‑101 (MQ=255)
                                                                                                      |                                                                                                      
GAACATCGGATAGCGCACCAGCGTCATCACCAGCCGCCGCCCGGTCAGCGCGACCTCGCGCCCGGG.CTCGACGCCCCGGGCCCGGCGCAGGGTGAAATGGC.GGGCG..AAATCCGCCCCCTCCGCCCGCTCCGGCCCCTGACCGACGATCCGCGCCAGCGTCGGCCCCATGGCGCGGATCCGCCCGCGCGCATCCAGCCACAG  >  NC_014034/2821348‑2821548

Base quality scores:  ATCG  < 3 ≤  ATCG  < 15 ≤  ATCG  < 24 ≤  ATCG  < 33 ≤  ATCG  < 39 ≤  ATCG