Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_014034 1,253,264 +G 69.3% coding (240/2421 nt) copA2 ← copper‑transporting P‑type ATPase

Read alignment evidence...
  seq id position change freq score reads annotation genes product
*NC_0140341,253,2641.→G69.3% 80.6 ‑1.6 ‑0.175coding (240/2421 nt)copA2copper‑transporting P‑type ATPase
Reads supporting (aligned to +/- strand):  new base (18/34):  ref base (12/11):  total (30/45)
Fisher's exact test strand distribution p-value = 2.03e-01
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 7.86e-01

Only 100 of 147 total aligned reads displayed.

CCGTGCCCTGCGCATAGGTGACGCGGGCGCTGCCCTCGGGCAAAAGCAGGACCGAGAGCACCCCGGGCGTCGCCTGCAGCGCGGACACGACCCTGT..CCGC.C..GGGGGGGGGGCGTCG.AAGCCCAGACAGGTCTCGGCCTCGGCAATCGGATAGCCCGCCCGCG.CCAGCGCCTCGGCCAGGGTTTCCGCCCGCGCCGG  >  NC_014034/1253164‑1253359
                                                                                                        |                                                                                                  
ccGTGCCCTGCGCATAGGTGACGCGGGCGCTGCCCTCGGGCAAAAGCAGGACCGAGAGCACCCCGGGCGTCGCCTGCAGCGCGGACACGACCCTGT..ccgc.c                                                                                                     >  1:1048442/1‑101 (MQ=255)
  gTGCCCTGCGCATAGGTGACGCGGGCGCTGCCCTCGGGCAAAAGCAGGACCGAGAGCACCCCGGGCGTCGCCTGCAGCGCGGACACGACCCTGT..CCGC.C..gg                                                                                                 >  2:945643/1‑101 (MQ=255)
    gCCCTGCGCATAGGTGACGCGGGCGCTGCCCTCGGGCAAAAGCAGGACCGAGAGCACCCCGGGCGTCGCCTGCAGCGCGGACACGACCCTGT..CCGC.C..gggg                                                                                               >  1:2699068/1‑101 (MQ=255)
       cTGCGCATAGGTGACGCCGGCGCTGCCCTCGGGCAAAAGCAGGACCGAGACCACCCCGGGCGTCGCCTGCTGCGCGGCCACGCCCCTGT..CCGC.C..ggggggg                                                                                            >  2:761984/1‑101 (MQ=255)
        tGCGCATAGGTGACGCGGGCGCTGCCCTCGGGCAAAAGCAGGACCGAGAGCACCCCGGGCGTCGCCTGCAGCGCGGACACGACCCTGT..CCGC.C..gggggggg                                                                                           <  2:3066324/101‑1 (MQ=255)
         gcgcATAGGTGACGCGGGCGCTGCCCTCGGGCAAAAGCAGGACCGAGAGCACCCCGGGCGTCGCCTGCAGCGCGGACACGACCCTGT..CCGC.C..ggggggggg                                                                                          <  1:1089248/101‑1 (MQ=255)
         gcgcATAGGTGACGCGGGCGCTGCCCTCGGGCAAAAGCAGGACCGAGAGCACCCCGGGCGTCGCCTGCAGCGCGGACACGACCCTGT..CCGC.C..ggggggggg                                                                                          <  2:1866837/101‑1 (MQ=255)
          cgcATAGGTGACGCGGGCGCTGCCCTCGGGCAAAAGCAGGACCGAGAGCACCCCGGGCGTCGCCTGCAGCGCGGACACGACCCTGT..CCGC.C..gggggggggg                                                                                         >  2:3918262/1‑101 (MQ=255)
          cgcATAGGTGACGCGGGCGCTGCCCTCGGGCAAAAGCAGGACCGAGAGCACCCCGGGCGTCGCCTGCAGCGCGGACACGACCCTGT..CCGC.C..gggggggggg                                                                                         >  2:2951678/1‑101 (MQ=255)
             atgggTGACGCGGGCGCCGCCCTCGGGCAAAAGCAGGACCGAGAGCACCCCGGGCGTCCCCTGCAGCGCGGACACGACCCTGT..CCGC.C..gggggggggggcg                                                                                      <  2:1083957/98‑4 (MQ=255)
             aTAGGTGACGCGGGCGCTGCCCTCGGGCAAAAGCAGGACCGAGAGCACCCCGGGCGTCGCCTGCAGCGCGGACACGACCCTGT..CCGC.C..gggggggggggcg                                                                                      >  1:3473378/1‑98 (MQ=255)
              tAGGTGACGCGGGCGCTGCCCTCGGGCAAAAGCAGGACCGAGAGCACCCCGGGCGTCGCCTGCAGCGCGGACACGACCCTGT..CCGC.C..GGGGGGGGGGCGTc                                                                                     >  1:113755/1‑101 (MQ=255)
                 ggggcgcgGGCGCTGCCCTCGGGCAAAAGCAGGACCGAGGGCCCCCCGGGCGTCGCCTGCAGCGCGGACACGCCCCTGT..CCGC.CG.GGGGGGGGGGCGTCG.a                                                                                  <  2:662038/97‑1 (MQ=255)
                  tGACGCGGGCGCTGCCCTCGGGCAAAAGCAGGACCGAGAGCACCCCGGGCGTCGCCTGCAGCGCGGACACGACCCTGT..CCGC.CG.GGGGGGGGGGCGTCgga                                                                                  >  2:1994265/1‑99 (MQ=255)
                  tGACGCGGGCGCTGCCCTCGGGCAAAAGCAGGACCGAGAGCACCCCGGGCGTCGCCTGCAGCGCGGACACGACCCTGT..CCGC.CG.GGGGGGGGGGCGTCG.aa                                                                                 >  1:1970522/1‑101 (MQ=255)
                  ggACGGGGGCGCTGCCCCCGGGCAAAAGCAGGACCGAGAGCCCCCCGGGCGTCGCCTGCAGCCCGGACACGACCCTGG..CCGC.CG.GGGGGGGGGGCGTCG.aa                                                                                 <  2:429459/100‑1 (MQ=255)
                   gACGCGGGCGCTGCCCTCGGGCAAAAGCAGGACCGAGAGCACCCCGGGCGTCGCCTGCAGCGCGGACACGACCCTGT..CCGC.CG.GGGGGGGGGGCGTCG.AAg                                                                                >  2:2838162/1‑101 (MQ=255)
                   gACGCGGGCGCTGCCCTCGGGCAAAAGCAGGACCGAGAGCACCCCGGGCGTCGCCTGCAGCGCGGACACGACCCTGT..CCGC.C..gggggggggggcggagaag                                                                                >  1:3086177/1‑92 (MQ=255)
                   gACGCGGGCGCTGCCCTCGGGCAAAAGCAGGACCGAGAGCACCCCGGGCGTCGCCTGCAGCGCGGACACGACCCTGT..CCGC.C..GGGGGGGGGGCGTCG.AAGc                                                                               >  1:2849270/1‑101 (MQ=255)
                    gcgcgGGCGCTCCCCTGGGGCAAAAGGGGGGCCGAGAGCCCCCGGGGCGGCGCCGCCGGCGCGGACACCCCCCTGTCCCCGG.G..GGGGGGGGGGCGTCG.AAg                                                                                <  2:1468364/100‑1 (MQ=255)
                    aCGCGGGCGCTGCCCTCGGGCAAAAGCAGGACCGAGAGCACCCCGGGCGTCGCCTGCAGCGCGGACACGACCCTGT..CCGC.C..ggggggggggggggcgaggc                                                                               >  2:1124626/1‑91 (MQ=255)
                    aCGCGGGCGCTGCCCTCGGGCAAAAGCAGGACCGAGAGCACCCCGGGCGTCGCCTGCAGCGCGGACACGACCCTGT..CCGC.C..ggggggggggggcgtcgagg                                                                               >  1:3392185/1‑91 (MQ=255)
                     cgcgGGCTCGGCCCGCGGGCAAAAGCAGGGCCGAGAGTCCCCCGTGCGTCGCCTGCAGGGCGGGCCCCCCCCTGT..CCTC.CGGGGGGGGGGGGCGTCG.AAGc                                                                               <  1:1393414/101‑1 (MQ=255)
                      ggggggcgtGCCCTGGGGAAAAAGGGGGACCGGGAGCCCCCCGGGCGTCGCCTGCAGGGGGGACCCGACCCTGT..CCCC.CG.GGGGGGGGGGCGTCG.AAGccc                                                                             <  1:3090723/93‑1 (MQ=255)
                      gcgGGCGCTGCCCTCGGGCAAAAGCAGGACCGAGAGCACCCCGGGCGTCGCCTGCAGCGCGGACACGACCCTGT..CCGC.CG.GGGGGGGGGGCGCCG.AAGccc                                                                             >  2:390068/1‑101 (MQ=255)
                      gcgGGCGCTGCCCTCGGGCAAAAGCAGGACCGAGAGCACCCCGGGCGTCGCCTGCAGCGCGGACACGACCCTGT..CCGC.C..GGGGGGGGGGCGTCG.AAGCCCa                                                                            >  2:2305355/1‑101 (MQ=255)
                       cgGGCGCTGCCCTCGGGCAAAAGCAGGACCGAGAGCACCCCGGGCGTCGCCTGCAGCGCGGACACGACCCTGT..CCGC.CG.GGGGGGGGGGCGTCG.AAGCCCa                                                                            >  1:2461709/1‑101 (MQ=255)
                        gggCGCGCCCCCCGGGCAAAAGCAGGGCCGGGAGCCCCCCGGGGGTCCCCTGCCGCGCGGACCCGACCCTGT..CCCC.CG.GGGGGGGGGGCGTCG.AAGCCCAg                                                                           <  1:1568586/101‑1 (MQ=255)
                         ggccccgCCGTCGGGAAAAAGCAGGACCGAGAGCCCCCCGGGCGTCGCCTGCAGCGCGGGCAGGACCCTGC..CCGC.C..GGGGGGGGGGCGTCG.AAGCCCAGAc                                                                         <  2:1038690/95‑1 (MQ=255)
                          gggcTGCCCTGGGGCAAAAGGAGGCCCGAGAGCCCCCCGGGCGTCGCCTGCAGCGCGGACACGACCCTGT..CCGC.CG.GGGGGGGGGGCGTCT.AAGCCCAGAc                                                                         <  2:2567513/99‑1 (MQ=255)
                            gcTGCCCTCGGGCAAAAGCAGGACCGAGAGCACCCCGGGCGTCGCCTGCAGCGCGGACACGACCCTGT..CCGC.C..GGGGGGGGGGCGTCG.AAGCCCAGACAgg                                                                      <  1:615010/101‑1 (MQ=255)
                             cTGCCCTCGGGCAAAAGCAGGACCGAGAGCACCCCGGGCGTCGCCTGCAGCGCGGACACGACCCTGT..CCGC.CG.GGGGGGGGGGGGTCG.AAGCCCAGACCgg                                                                      >  2:1104944/1‑101 (MQ=255)
                              gggccTCGGGGAAAAGCAGGACCGAGAGCACCCCGGGCGTCGCCTGCAGCGCGGACACGACCCTGT..CCGC.CG.GGGGGGGGGGCGTCG.AAGCCCAGACAgct                                                                     <  2:2258305/98‑3 (MQ=255)
                               gCCCTCGGGCAAAAGCAGGACCGAGAGCACCCCGGGCGTCGCCTGCAGCGCGGACACGACCCTGT..CCGC.CG.GGGGGGGGGGCGGCG.AAGCCCAGCCAGGtc                                                                    >  2:350501/1‑101 (MQ=255)
                               gCCCTCGGGCAAAAGCAGGACCGAGAGCACCCCGGGCGTCGCCTGCAGCGCGGACACGACCCTGT..CCGC.C..GGGGGGGGGGCGTCG.AAGCCCAGACAGGcct                                                                   >  1:2697819/1‑101 (MQ=255)
                                ccccggggCAAAAGCAGGCCCGAGAGCACCCCGGGGCTCCCCTGCAGCGGGGACACCACCCTGT..CCCC.CG.GGGGGGGGGGCGTCG.AAGCCCAGACAGGtct                                                                   <  2:1060645/96‑1 (MQ=255)
                                 ccTCGGGCAAAAGCAGGACCGAGGGCCCCCCGGGCGTCGCCTGCAGCGCGGACACGCCCCTGT..CCGC.CG.GGGGGGGGGGCGTCG.AAGCCCAGACAGGtctc                                                                  <  1:1051637/101‑1 (MQ=255)
                                 ccTCGGGCAAAAGCAGGACCGAGAGCCCCCCGGGCGTCGCCTGCGGCGCGGACACGACCCTGT..CCGC.CG.GGGGGGGGGGCGTCG.AAGCCCAGACAGGtctc                                                                  <  1:3890520/101‑1 (MQ=255)
                                  ccggggCAAAAGCAGGACCGAGAGCCCCCCGGGGGTCGCCTGCAGCGCGGACACCACCCTGT..CCGC.CG.GGGGGGGGGGCGTCG.AAGCCCAGACAGGTctcg                                                                 <  2:4115714/98‑1 (MQ=255)
                                      ggCAAAAGCAGGACCGAGAGCACCCCGGGCGTCGCCTGCAGCGCGGACACGACCCTGT..CCGC.CG.GGGGGGGGGGCGTCG.AAGCCCAGACAGGTctggggcc                                                             >  1:1001641/1‑98 (MQ=255)
                                      ggCAAAAGCAGGACCGAGAGCACCCCGGGCGTCGCCTGCAGCGCGGACACGACCCTGT..CCGC.CG.GGGGGGGGGGCGTCG.AAGCCCAGACAGGGcccggccc                                                             >  2:4151364/1‑100 (MQ=255)
                                        cAAAAGCAGGCCCGAGGGCCCCCCGGGCGTCCCCTGCAGCGCGGACCCGCCCCTGT..CCGC.G..GGGGGGGGGGCGTCG.AAGCCCAGACAGGTctcggcctcgg                                                          <  2:3581508/101‑1 (MQ=255)
                                        aaaagGCAGGCCCGAGAGCCCCCCGGGCGTCCCCTGCAGCGCGGCCCCGACCCTGT..CCGC.CG.GGGGGGGGGGCGTCG.AAGCCCAGACAGGTctcggcctcg                                                           <  1:828990/100‑1 (MQ=255)
                                         aaaaGCAGGACCGAGAGCACCCCGGGCGTCGCCTGCAGCGCGGACACGACCCTGT..CCGC.CG.GGGGGGGGGGGGTGG.AAGCCCAGCCAGGCcccggccccgg                                                          >  1:3575129/1‑101 (MQ=255)
                                         aaaaGCAGGACCGAGAGCACCCCGGGCGTCGCCTGCAGCGCGGACACGACCCTGT..CCGC.C..GGGGGGGGGGCGTCG.AAGCCCAGACAGGTctcggcctcggc                                                         >  2:509635/1‑101 (MQ=255)
                                           aaGCAGGACCGGGACCCCCCCGGGCGGCGCCCGCAGCGCGGACCCGGCCCTGT..CCCC.CG.GGGGGGGGGGCGTCG.AAGCCCAGACAGGTCTCGGCCTCGGCa                                                        <  2:3472094/101‑1 (MQ=255)
                                           aaGCAGGACCGAGAGCACCCCGGGCGTCGCCTGCAGCGCGGACACGACCCTGT..CCGC.CG.GGGGGGGGGGCGGCG.AAGCCCAGCCAGCCCTGGGCCTCCGCa                                                        >  2:2785680/1‑101 (MQ=255)
                                               gggACGGAGGGCACCCCGGGGTGCGCCTGCGGCCGGGCCACGCCCCTTT..CCCC.C..GGGGGGGGGGCGTCG.AAGCCCAGACAGGTCTCGGCCTCGGCAATCgg                                                   <  1:537369/100‑1 (MQ=255)
                                               aGGACCGAGAGCACCCCGGGCGTCGCCTGCAGCGCGGACACGACCCTGT..CCGC.CG.GGGGGGGGGGGGTCG.AAGCCCAGACAGGTCTCGGCCTCGGCAATCg                                                    >  1:855370/1‑101 (MQ=255)
                                               aGGACCGAGAGCACCCCGGGCGTCGCCTGCAGCGCGGACACGACCCTGT..CCGC.CG.GGGGGGGGGGCGTCG.AGCCCCAGACAGGCCTCTGCCTCGCCAATCg                                                    >  1:2075135/1‑101 (MQ=255)
                                                 gACCGGGAGCCCCCCGGGGGTCCCCCTCAGCGCGGGCACCGCCCCGT..CCCC.CG.GGGGGGGGGGCGTCG.AAGCCCAGACAGGTCTCGGCCTCGGCAATCGGa                                                  <  2:2709618/101‑1 (MQ=255)
                                                 gACCGAGAGCCCCCCGGGCGTCGCCTGCAGCGCGGACACGACCCTGT..CCGC.C..GGGGGGGGGGCGTCG.AAGCCCAGACAGGTCTCGGCCTCGGCAATCGGAt                                                 <  2:3144938/101‑1 (MQ=255)
                                                    cgggagCCCCCCGGGCGTCGCCTGCGGCGCGGACCCGCCCCTGT..CCCC.C..GGGGGGGGGGCGTCG.AAGCCCAGACAGGTCTCGGCCTCGGCAATCGGATAgc                                              <  1:2305523/98‑1 (MQ=255)
                                                     gaggGCCCCCCGGGCGTCCCCGGCAGCGCGGACCCGACCCTGT..CCGC.CG.GGGGGGGGGGCGTCG.AAGCCCAGACAGGTCTCGGCCTCGGCAATCGGATAgc                                              <  2:2287107/101‑1 (MQ=255)
                                                     gagaGCACCCCGGGCGTCGCCTGCAGCGCGGACACGACCCTGT..CCGC.CG.GGGGGGGGGGGGCCG.AAGTCCAGCCAGGCCCCGGCCCCGGCAACCGGaaacc                                              >  1:3669881/1‑97 (MQ=255)
                                                     gagaGCACCCCGGGCGTCGCCTGCAGCGCGGACACGACCCTGT..CCGC.CG.GGGGGGGGGGCGTCG.AAGCCCAGACAGGTCTCGGCCTCGGCAATCGGATAgc                                              >  1:4107508/1‑101 (MQ=255)
                                                       gggcccccgGGGCGGCCCCGGCAGCGCGGACCCGCCCCGGC..CCGCGG..GGGGGGGGGGCGTCG.AAGCCCAGACAGGTCTCGGCCTCGGCAATCGGATAgccc                                            <  1:3722225/96‑1 (MQ=255)
                                                       gggccccccGGGCGTCGCCTGCAGCGCGGACACGGCCCTGT..CCGC.CG.GGGGGGGGGGCGTCG.AAGCCCAGACAGGTCTCGGCCTCGGCAATCGGATAgccc                                            <  1:390119/96‑1 (MQ=255)
                                                       gagCACCCCGGGCGTCGCCTGCAGCGCGGACACGACCCTGT..CCGC.CG.GGGGGGGGGGCCTCC.AAGCCCGAACCGGACTCGCCCCCGGCAATCGGAAAgccc                                            >  2:3247098/1‑101 (MQ=255)
                                                            ccccGGGCGGCCCCTGCAGCGCGGACCCGACCCTGT..CCCC.CG.GGGGGGGGGGCGTCG.AAGCCCAGACAGGTCTCGGCCTCGGCAATCGGATAGcccgcccg                                       <  1:97576/101‑1 (MQ=255)
                                                             cccGGGCGTCGCCTGCAGCGCGGACACGACCCTGT..CCGC.CG.GGGGGGGGGGGGTGG.AAGCCCAAACAGGTCTCGGCCTCGGCAATCGGAAAacccgccccc                                      >  2:2575789/1‑99 (MQ=255)
                                                             cccGGGCGTCGCCTGCAGCGCGGACACGACCCTGT..CCGC.C..GGGGGGGGGGCGTCG.AAGCCCAGACAGGTCTCGGCCTCGGCAATCGGATAGCCCGCCcgcg                                     >  1:2249029/1‑101 (MQ=255)
                                                             cccGGGCGGCCCCTGCCGCGGGGACCCGCCCCGGC..CCGCGG..GGGGGGGGGGCGTCG.AAGCCCAGACAGGTCTCGGCCTCGGCAATCGGATAGCccgcccgc                                      <  1:773742/101‑1 (MQ=255)
                                                              ccGGTCGTGGCCTGCAGCGCGGACACGACCCTGT..CCGC.CG.GGGGGGGGGGCGTCA.AAGCCCGGCCAGTCCCCGCCCTCGCCACT.GGGTAgcccgcgcacgg                                    >  2:1592303/1‑95 (MQ=255)
                                                              ccGGGCGTCGCCTGCAGCGCGGACACGACCCTGT..CCGC.CG.GGGGGGGGGGCGTCG.AAGCCCAGACAGGTCTCGGCCTCGGCAATCGGATAGCccccccgcc                                     >  1:1639808/1‑100 (MQ=255)
                                                               cGGGCGTCGCCTGCAGCGCGGACACGACCCTGG..CCGC.CG.GGGGGGGGGGCGTCG.AAGCCCAGACAGGTCTCGGCCTCGGCAATCGGATAGCCCGCCCgcg.c                                   <  2:2361095/101‑1 (MQ=255)
                                                                 ggCGTCGCCTGCAGGGCGGACACGCCCCTGT..CCGCGG..GGGGGGGGGGCGTCG.AAGCCCAGACAGGTCTCGGCCTCGGCAATCGGATAGCCCGCCCGCG.CCa                                 <  2:1535915/101‑1 (MQ=255)
                                                                 ggCGCCGCCTCCGGCGGGGCCCCGCCCCTGC..CCGC.CG.GGGGGGGGGGCGTCG.AAGCCCAGACAGGTCTCGGCCTCGGCAATCGGATAGCCCGCCCGCG.CCa                                 <  2:4031964/101‑1 (MQ=255)
                                                                   cggcGCCTGCAGCCCGGACACGACCCTGT..CCGC.C..GGGGGGGGGGCGTCG.AAGCCCAGACAGGTCTCGGCCTCGGCAATCGGATAGCCCGCCCGCG.CCAgcg                              <  1:2564642/98‑1 (MQ=255)
                                                                    gTCGCCGGCAGGCCGGACCCGACCCTGT..CCGC.CG.GGGGGGGGGGCGTCG.AAGCCCAGACAGGTCTCGGCCTCGGCAATCGGATAGCCCGCCCGCG.CCAgcg                              <  2:1186929/101‑1 (MQ=255)
                                                                     tCGCCTGCAGCGCGGACACGACCCTGT..CCGC.CG.GGGGGGGGGGGGCCG.AAGCCCAGACAGGTCTCGGCCCCGGCAACCGGAAAGCCCCCCCCCC.CCAccgc                             >  2:1972470/1‑101 (MQ=255)
                                                                       gCCTGCAGCGCGGACACGACCCTGT..CCGC.CG.GGGGGGGGGGCGGCG.AAGCCCAGGCAGGTCTCGGCCTCGGCAAACGGATAGCCCGCCCGCG.CCAGCCccc                           >  1:4128688/1‑100 (MQ=255)
                                                                        ccTGCAGCGCGGACACGACCCTGT..CCGC.CG.GGGGGGGGGGCGTCG.AAGCCCAGACAGGTCTCGGCCTCGGCAATCGGATAGCCCGCCCGCG.CCAGCGCCTc                          <  2:2263929/101‑1 (MQ=255)
                                                                         cggCAGCGCGGACAGGACCCAGT..CCGCGG..GGGGGGGGGGCGTCG.AAGCCCAGACAGGTCTCGGCCTCGGCAATCGGATAGCCCGCCCGCG.CCAGCGCCTCg                         <  1:3247339/99‑1 (MQ=255)
                                                                         ccgCAGCGGGGACCCGCCCCTGT..CCCC.CG.GGGGGGGGGGCGTCG.AAGCCCAGACAGGTCTCGGCCTCGGCAATCGGATAGCCCGCCCGCG.CCAGCGCCTCg                         <  1:1322415/99‑1 (MQ=255)
                                                                          tGCAGCGCGGACACGACCCTGT..CCGC.CG.GGGGGGGGGGGGTCC.AAGCCCAAAAAGGTCTTGGCCTCGGCAATCGGATACGCCCCCCGCG.CCCGGGCCCCgg                        >  2:4202381/1‑101 (MQ=255)
                                                                           gCAGCGCGGACACGACCCTGT..CCGC.C..GGGGGGGGGGCGTCG.AAGCCCAGACAGGTCTCGGCCTCGGCAATCGGATAGCCCGCCCGCG.CCAGCGCCTCGGcc                      >  2:2360769/1‑101 (MQ=255)
                                                                            cAGCGCGGACACGCCCCTGT..CCGC.C...GGGGGGGGGCGTCG.AAGCCCAGACAGGTCTCGGCCTCGGCAATCGGATAGCCCGCCCGCG.CCAGCGCCTCGGCCAg                    <  1:3593164/101‑1 (MQ=255)
                                                                             aGCGCGGACACGACCCTGT..CCGC.CG.GGGGGGGGGGCGTCG.AAGCCCAGACAGGTCTCGGCCTCGGCAATCGGATAGCCCGCCCGCG.CCAGCGCCTCGGCCa                     >  2:3231845/1‑101 (MQ=255)
                                                                              gcgcGGACACGACCCTGC..CCGC.CG.GGGGGGGGGGCGTCG.AAGCCCAGACAGGTCTCGGCCTCGGCAATCGGATAGCCCGCCCGCG.CCAGCGCCTCGGCCAg                    <  1:1680887/101‑1 (MQ=255)
                                                                               cgcgGACACGGCCCTGT..CCGC.CG.GGGGGGGGGGCGTCG.AAGCCCAGACAGGTCTCGGCCTCGGCAATCGGATAGCCCGCCCGCG.CCAGCGCCTCGGCCAgg                   <  2:3011283/101‑1 (MQ=255)
                                                                               cgcgGACACGACCCTGT..CCGC.CG.GGGGGGGGGGCGTCG.AAGCCCAGACAGGTCTCGGCCTCGGCAATCGGATAGCCCGCCCCCG.CCACCGCCCCGGCCCgg                   >  1:3350246/1‑101 (MQ=255)
                                                                                 cgGACACGACCCTGT..CCGC.C..GGGGGGGGGGCGTCG.AAGCCCAGACAGGTCTCGGCCTCGGCAATCGGCTAGCCCGCCCGCG.CCCGCGCCTCGGCCCggggt                >  1:3756152/1‑99 (MQ=255)
                                                                                  ggACCCGACCCTGT..CCGCGG..GGGGGGGGGGCGTCG.AAGCCCAGACAGGTCTCGGCCTCGGCAATCGGATAGCCCGCCCGCG.CCAGCGCCTCGGCCAGGGtt                <  2:2461525/101‑1 (MQ=255)
                                                                                  ggACACGACCCTGT..CCGC.CG.GGGGGGGGGGGGTCG.AAGCCCAGACAGGTCTCGGCCTCGGCAATCGGATAGCCCGCCCGCG.CCCGCGCCTCGGCCAGGGtt                >  1:431416/1‑101 (MQ=255)
                                                                                   gACACGACCCTGTC.CCGC.G..GGGGGGGGGGCGTCG.AAGCCCAGACAGGTCTCGGCCTCGGCAATCGGATAGCCCGCCCGCG.CCAGCGCCTCGGCCAGGGttt               <  2:1787845/101‑1 (MQ=255)
                                                                                   gACACGACCCTGT..CCGC.CG.GGGGGGGGGGGGTGG.AACCCCAAACAGGTCCCGGCCTGGGCAATCGGATAGCCCCCCCGCC.CCCGCCCCCCGCCCAGGGttc               >  1:1860116/1‑100 (MQ=255)
                                                                                    gcccgccccTGG..CCCC.CG.GGGGGGGGGGCGTCG.AAGCCCAGACAGGTCTCGGCCTCGGCAATCGGATAGCCCGCCCGCG.CCAGCGCCTCGGCCAGGGTTTc              <  2:855289/95‑1 (MQ=255)
                                                                                    acacGACCCTGT..CCGC.CG.GGGGGGGGGGGGTCG.AAGCCCAGACAGGTCTCGGCCTCGGCAATCGGATAGCCCGCCCCCG.CCCGGGCCCCCGCCAGGGTTTc              >  2:4068490/1‑101 (MQ=255)
                                                                                      acGACCCTGT..CCGC.CG.GGGGGGGGGGGGCCG.AAGCCCAGAAAGGTTCCGGCCCCGGAAATCGGAAAGCCCCCCCGCC.CCAGCCCCCCGCCCGGGGTTTccc            >  2:2638742/1‑100 (MQ=255)
                                                                                      acGACCCTGT..CCGC.C..GGGGGGGGGGCGTCG.AAGCCCAGACAGGTCTCGGCCTCGGCAATCGGATAGCCCGCCCGCG.CCAGCGCCTCGGCCAGGGTTTccgc           >  1:1925150/1‑101 (MQ=255)
                                                                                       cGACCCCGT..CCCC.CG.GGGGGGGGGGCGTCG.AAGCCCAGACAGGTCTCGGCCTCGGCAATCGGATAGCCCGCCCGCG.CCAGCGCCTCGGCCAGGGTTTccgc           <  1:4236238/101‑1 (MQ=255)
                                                                                        gACCCTGT..CCGC.C..GGGGGGGGGGCGTCG.AAGCCCAGACAGGTCTCGGCCTCGGCAATCGGATAGCCCGCCCGCG.CCAGCGCCTCGGCCAGGGTTTccgccc         >  1:2364155/1‑101 (MQ=255)
                                                                                             tGT..CCGC.CG.GGGGGGGGGGGGTCG.AAGCCCAGACAGGTCTCGGCCCCGGCAATCGGATAGCCCCCCCGCC.CCCGCCCCTGGGCCAGGGTTTCCCCCCccgc     >  2:1920335/1‑101 (MQ=255)
                                                                                             tGT..CCGC.CG.GGGGGGGGGGCGTCG.AAGCCCAGACAGGTCTCGGCCTCGGCAATCGGATAGCCCGCCCGCG.CCAGCGCCCCGGCCAGGGTTTccccccgccc     >  2:4017957/1‑99 (MQ=255)
                                                                                             tGT..CCGC.C..GGGGGGGGGGGGGCGAAAACCCAGAAAGTTTTGGGCCCCGCCAACGGGGAAACCCCCCCGGG.CCCGCCCCCCGCCCCGGGTTTccgcccgccc     >  2:2967454/1‑99 (MQ=255)
                                                                                              gT..CCGC.CG.GGGGGGGGGGGGTCG.AAGCCAAGACAGGCCTCGGCCCGGCCAATCGGA.AACCCGCCCCCGCCCACCGCCGCGGCCGGGGTTT.CGCCCCCGccc   >  2:1731272/1‑100 (MQ=255)
                                                                                                 gccgc.cg.GGGGGGGGGGCGTCG.AAGCCCAGACAGGTCTCGGCCTCGGCAATCGGATAGCCCGCCCGCG.CCAGCGCCTCGGCCAGGGTTTCCGCCCGCGCCg   <  1:2838381/100‑1 (MQ=255)
                                                                                                 ccccc.c..gGGGGGGGGGCGTCG.AAGCCCAGACAGGTCTCGGCCTCGGCAATCGGATAGCCCGCCCGCG.CCAGCGCCTCGGCCAGGGTTTCCGCCCGCGCCgg  <  2:3755895/97‑1 (MQ=255)
                                                                                                  ccgc.cg.GGGGGGGGGGCGTCG.AAGCCCAGACAGGTCTCGGCCTCGGCAATCGGATAGCCCGCCCGCG.CCAGCGCCTCGGCCAGGGTTTCCGCCCGCGCCgg  <  1:2423286/101‑1 (MQ=255)
                                                                                                        |                                                                                                  
CCGTGCCCTGCGCATAGGTGACGCGGGCGCTGCCCTCGGGCAAAAGCAGGACCGAGAGCACCCCGGGCGTCGCCTGCAGCGCGGACACGACCCTGT..CCGC.C..GGGGGGGGGGCGTCG.AAGCCCAGACAGGTCTCGGCCTCGGCAATCGGATAGCCCGCCCGCG.CCAGCGCCTCGGCCAGGGTTTCCGCCCGCGCCGG  >  NC_014034/1253164‑1253359

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