Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_014034 1,396,957 T→A 9.7% F280Y (TTC→CAT)  RCAP_rcc01307 → hypothetical protein

Read alignment evidence...
  seq id position change freq score reads annotation genes product
*NC_0140341,396,9570T→A 9.7% 104.6 ‑3.2 ‑0.2402F280Y (TTC→TAC) RCAP_rcc01307hypothetical protein
Reads supporting (aligned to +/- strand):  new base (13/26):  ref base (185/177):  total (198/204)
Fisher's exact test strand distribution p-value = 4.27e-02
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 6.75e-01

Only 100 of 445 total aligned reads displayed.

TGGATCGAAACCGGTCCAGAGATAGACCGCGCGGACCTCCTCAGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACA.GC.CT.GCGCG.GC.CGCA..TCGGCT.TCCGC.CTG..CGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCAGCTCGTCGTTGGTCAGCCGGACGATGTCCGTCGCTCCCCCCGCGGGGGCGTGG  >  NC_014034/1396857‑1397057
                                                                                                            |                                                                                                        
tGGATCGAAACCGGTCCAGAGATAGACCGCGCGGACCTCCTCAGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACA.GC.CT.GCGCG.GC.CGCA..TCGGCt.t                                                                                                          >  2:958840/1‑101 (MQ=255)
tGGATCGAAACCGGTCCAGAGATAGACCGCGCGGACCTCCTCAGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACA.GC.CT.GCGCG.GC.CGCA..TCGGCt.t                                                                                                          <  1:1484668/101‑1 (MQ=255)
 ggATCGAAACCGGTCCAGAGATAGACCGCGCGGACCTCCTCAGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACA.GC.CT.GCGCG.GC.CGCA..TCGGCT.Tc                                                                                                         >  2:1468206/1‑101 (MQ=255)
  gATCGAAACCGGTCCAGAGATAGACCGCGCGGACCTCCTCAGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACA.GC.CT.GCGCG.GC.CGCA..TCGGCT.Tca                                                                                                        >  1:3671096/1‑100 (MQ=255)
     cGAAACCGGTCCAGAGATAGACCGCGCGGACCTCCTCAGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACA.GC.CT.GCGCG.GC.CGCA..TCGGCT.TCCGc.c                                                                                                    <  1:3593500/101‑1 (MQ=255)
      gAAACCGGTCCAGAGATAGACCGCGCGGACCTCCTCAGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACA.GC.CT.GCGCG.GC.CGCA..TCGGCT.TCCGC.Ct                                                                                                   >  2:1694137/1‑101 (MQ=255)
       aaaCCGGTCCAGAGATAGACCGCGCGGACCTCCTCAGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACA.GC.CT.GCGCG.GC.CGCA..TCGGCT.TCCGC.CTg                                                                                                  <  2:3108831/101‑1 (MQ=255)
        aaCCGGTCCAGAGATAGCCCGCGCGGACCTCCTCAGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACA.GC.CT.GCGCG.GC.CGCA..TCGGCT.TCCGC.CTG..c                                                                                               <  2:3208583/101‑1 (MQ=255)
                 agagATAGACCGCGCGGACCTCCTCAGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACA.GC.CT.GCGCG.GC.CGCA..TCGGCT.TCCGC.CTG..CGTG.GTCGCt                                                                                     <  2:3512569/101‑1 (MQ=255)
                 agagATAGACCGCGCGGACCTCCTCAGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACA.GC.CT.GCGCG.GC.CGCA..TCGGCT.TCCGC.CTG..CGTG.GTCGCt                                                                                     >  2:4417943/1‑101 (MQ=255)
                  gagaTAGACCGCGCGGACCTCCTCAGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACA.GC.CT.GCGCG.GC.CGCA..TCGGCT.TCCGC.CTG..CGTG.GTCGCTc                                                                                    >  2:1116184/1‑101 (MQ=255)
                   agaTAGACCGCGCGGACCTCCTCAGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACA.GC.CT.GCGCG.GC.CGCA..TCGGCT.TCCGC.CTG..CGTG.GTCGCTcc                                                                                   <  1:2909133/101‑1 (MQ=255)
                      tAGACCGCGCGGACCTCCTCAGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACA.GC.CT.GCGCG.GC.CGCA..TCGGCT.TCCCC.CTG..CGTG.GTCGCTccccc                                                                                >  1:3955416/1‑101 (MQ=255)
                       agcccGCGCGGACCTCCTCAGACGTCGCCCCCCCCGCGTGGGCGTGGATCGAAACA.GC.CG.GCGCG.GC.CGCA..TCGGCT.TCCGC.CTG..CGTG.GTCGCTgccccc                                                                               <  2:225060/98‑1 (MQ=38)
                       aGACCGCGCGGACCTCCTCAGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACA.GC.CT.GCGCG.GC.CGCA..TCGGCT.TCCGC.CTG..CGTG.GTCGCTcccccc                                                                               >  2:3407247/1‑101 (MQ=255)
                               cgGACCTCCTCAGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACA.GC.CT.GCGCG.GC.CGCA..TCGGCT.TCCGC.CTC..CGTG.GTCGCTCCCCCCGCGGGGGc                                                                       >  1:2130045/1‑101 (MQ=255)
                                 gccctcctCAGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACA.GC.CT.GCGCG.GC.CGCA..TCGGCT.TCCGC.CTG..CGTG.GTCGTTCCCCCCGCGGGGGCGt                                                                     <  1:1854325/99‑1 (MQ=255)
                                 gccctcctCAGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACA.GC.CT.GCGCG.GC.CGCA..TCGGCT.TCCGC.CTG..CGTG.GTCGCTCCCCCCGCGGGGGCGt                                                                     <  2:3633955/99‑1 (MQ=255)
                                 gACCTCCTCAGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACA.GC.CT.GCGCG.GC.CGCA..TCGGCT.TCCGC.CTG..CGTG.GTCGCTCCCCCCGCGGGGGCGt                                                                     <  2:2958992/101‑1 (MQ=255)
                                 gACCTCCTCAGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACA.GC.CT.GCGCG.GC.CGCA..TCGGCT.TCCGC.CTG..CGTG.GTCGCTCCCCCCGCGGGGGCGt                                                                     <  2:1809747/101‑1 (MQ=255)
                                   cctcctCAGACGTCGCCCCCCCCGCGGGGGCGTGGATCGAAACA.GC.CT.GCGCG.GC.CCCA..TCGGCT.TCCGC.CTG..CGTG.GTCGCTCCCCCCGCGGGGGCGTgg                                                                   <  2:4993851/101‑1 (MQ=255)
                                    tcgcggttgcGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATGG.CTATGGCG.GC.AGTG..TCAGGT.T.TTC.CAG..CGGC.GTCGCTCCCCCCGCGGGGGCGTgg                                                                   <  2:409452/92‑1 (MQ=0)
                                    ctcctcAGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACA.GC.CT.GCGCG.GC.CGCA..TCGGCT.TCCGC.CTG..CGTG.GTCGCTCCCCCCGCGGGGGCGTGGa                                                                  <  1:2808515/101‑1 (MQ=255)
                                     ccctcAGACGTCGCCCCCCCCGGGGGGGCGTGGATCGAAACA.GC.CT.GCGCG.GC.CGCA..TCGGCT.TCCGC.CTG..CGTG.GTCGCTCCCCCCGC.GGGGCGTGGATc                                                                <  2:41905/100‑1 (MQ=38)
                                      cctcAGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACA.GC.CT.GCGCG.GC.CGCA..TCGGCT.TCCGC.CTG..CGTG.GTCGCTCCCCCCGCGGGGGCGTGGATc                                                                <  1:936471/101‑1 (MQ=255)
                                      cctcAGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACA.GC.CT.GCGCG.GC.CGCA..TCGGCT.TCCGC.CTG..CGTG.GTCGCTCCCCCCGCGGGGGCGTGGATc                                                                >  2:2606539/1‑101 (MQ=255)
                                        ggttgcGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATGG.CTATGGCG.GC.AGTG..TCAGGT.T.TTC.CAG..CGGC.GTCGCTCCCCCCGCGGGGGCGGGGACCg                                                               >  2:2397628/6‑101 (MQ=0)
                                         cAGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACA.GC.CT.GCGCG.GC.CGCA..TCGGCT.TCCGC.CTG..CGTG.GTCGCTCCCCCCGCGGGGGCGGGGATCGaa                                                             >  1:635910/1‑101 (MQ=255)
                                          aGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACA.GC.CT.GCGCG.GC.CGCA..TCGGCT.TCCGC.CTG..CGTG.GTCGCTCCCCCCGCGGGGGGGGGGGTCGaaa                                                            >  1:20549/1‑101 (MQ=255)
                                           cccGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GC.GT.GCGCG.CCATGAAT.TCATCG.AGCGC.CAT..CGAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGa                                                              <  1:4788053/99‑1 (MQ=2)
                                            gcGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATGG.CTATGGCG.GC.AGTG..TCAGGT.T.TTC.CAG..CGGC.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAAc                                                           <  1:716399/100‑1 (MQ=2)
                                            cggTCGCTCCCCCCGCGGGGGCGTGGATCGAAACA..A.CT.GACC...C.GGCA..ACGGCC.TATGCGCTGGACGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAAct                                                          >  2:3058242/3‑100 (MQ=2)
                                             tgTCGCTCCCCCCGCGGGGGCGTGGATCGAAAC..GC.TG.GCGCGCGC.CGCG..CCGGATCTTTGT.CTG..CGCG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAAct                                                          >  1:1780662/2‑100 (MQ=9)
                                             cGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GC.GT.GCGCG.CCATGAAT.TCATCG.AGCGC.CAT..CGAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGaaa                                                            <  1:2838087/101‑1 (MQ=2)
                                             cGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GC.GT.GCGCG.CCATGAAT.TCATCG.AGCGC.CAT..CGAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGaaa                                                            >  2:477790/1‑101 (MQ=2)
                                              gTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GC.GT.GCGCG.CCATGAAT.TCATCG.AGCGC.CAT..CGAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAAc                                                           >  1:440269/1‑101 (MQ=2)
                                              gTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATGG.CTATGGCG.GC.AGTG..TCAGGT.T.TTC.CAG..CGGC.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAAccc                                                         <  2:1090364/101‑3 (MQ=2)
                                              gTCGCTCCCCCCGCGGGGGCGTGGATCGAAACA.GC.CT.GCGCG.GC.CGCA..TCGGCT.TCCGC.CTG..CGTG.GTCGCTCCCCCCGCGGGGGCGTGGGTCGAAACATc                                                        >  2:1004889/1‑101 (MQ=255)
                                              gTCGCTCCCCCCGCGGGGGCGTGGATCGAAAC..GC.TG.GCGCGCGC.CGCG..CCGGATCTTTGT.CTG..CGCG.GTCGCTCCCCCCGCGGGGGCGTGGGTCGAAActc                                                         >  1:1057255/1‑99 (MQ=9)
                                              gTCGCTCCCCCCGCGGGGGCGTGGATCGAAAC..GC.TG.GCGCGCGC.CGCG..CCGGATCTTTGT.CTG..CGCG.GTCGCTCCCCCCGCGGGGGCGCGGATCGAAActc                                                         >  1:1664245/1‑99 (MQ=9)
                                               tCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GC.GT.GCGCG.CCATGAAT.TCATCG.AGCGC.CAT..CGAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACa                                                          <  2:1261210/101‑1 (MQ=2)
                                               tCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GC.GT.GCGCG.CCATGAAT.TCATCG.AGCGC.CAT..CGAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACa                                                          >  1:4031554/1‑101 (MQ=2)
                                               tCGCTCCACCCGCGGGGTCGTGGGCCGAAACA.ACGCG.GCGCG.GG.CCTA..T.GGCTGTCGCC.CCC..CCCC.GTCGCTCCCCCCCCGGGGGCGTGGATCGAAActgc                                                        <  2:710065/101‑4 (MQ=2)
                                               gcGCGCCCCCCGCGGGGGCGTGGATCGAAACA.GC.CT.GCGCG.GC.CGCA..TCGGCT.TCCGC.CTG..CGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCa                                                       <  2:624089/100‑1 (MQ=255)
                                                cGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GC.GT.GCGCG.CCATGAAT.TCATCG.AGCGC.CAT..CGAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACac                                                         >  1:1221739/1‑100 (MQ=2)
                                                cGCTCCCCCCGCGGGGGCGTGGATCGAAACA..A.CT.GACC...C.GGCA..ACGGCC.TATGCGCTGGACGTG.GTCGCTCCCCCCGCGGGTGCGTGGATCGAAActaaac                                                      >  2:4009691/1‑96 (MQ=2)
                                                cGCTCCCCCCGCGGGGGCGTGGATCGAAAC..GC.TG.GCGCGCGC.CGCC..CCGGATCTTTGT.CTG..CGCG.GGCGCTCCCCCCGCGGGGGCGTGGATCGAAActcgg                                                       >  1:2235468/1‑97 (MQ=9)
                                                 gCTCCCCCCGCGGGGGCGTGGATCGAAACATGG.CTATGGCG.GC.AGTG..TCAGGT.T.TTC.CAG..CGGC.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAAcccgcc                                                      <  1:3386612/101‑6 (MQ=2)
                                                 gCTCCCCCCGCGGGGGCGTGGATCGAAACA..A.CT.GACC...C.GGCA..ACGGCC.TATGCGCTGGACGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAActacacc                                                     <  2:1447142/101‑7 (MQ=2)
                                                 gCTCCCCCCGCGGGGGCGTGGATCGAAAC..GC.TG.GCGCGCGC.CGCG..TCGGATCTTTGT.CTG..CGCG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAActcggc                                                      <  1:4247461/101‑6 (MQ=9)
                                                 gCTCCCCCCGCGGGGGCGTGGATCGAAAC..GC.TG.GCGCGCGC.CGCG..CCGGATCTTTGT.CTG..CGCG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAActcggc                                                      >  1:1927665/1‑96 (MQ=9)
                                                 ccTCCCCCCGCGGGGGCGGGGATCGAAACA.GC.CT.GCGCG.GC.CGCA..TCGGCT.TCCGC.CTG..CGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCacc                                                     <  2:4909765/100‑3 (MQ=255)
                                                  cacccccgGCGGGGGCGTGGATCGAAACA.GC.CT.GCGCG.GC.CGCA..TCGGCT.TCCGC.CTG..CGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCAGCt                                                    <  1:4995760/99‑1 (MQ=255)
                                                  cTCCCCCCGCGGGGGCGTGGATCGAAACT.GC.GT.GCGCG.CCATGAAT.TCATCG.AGCGC.CAT..CGAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACACCa                                                       <  1:3411648/101‑1 (MQ=2)
                                                  cTCCCCCCGCGGGGGCGTGGATCGAAACA.GC.CT.GCGCG.GC.CGCA..TCGGCT.TCCGC.CTG..CGTG.GTCGCTCCCTCCGCGGGGGCGTGGATCGAAACATCAGCt                                                    <  2:468931/101‑1 (MQ=255)
                                                   tCCCCCCGCGGGGGCGTGGATCGAAACT.GC.GT.GCGCG.CCATGAAT.TCATCG.AGCGC.CAT..CGAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACACCaa                                                      <  1:2944103/101‑2 (MQ=2)
                                                   tCCCCCCGCGGGGGCGTGGATCGAAACA..A.CT.GACC...C.GGCA..ACGGCC.TATGCGCTGGACGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAActacaccgg                                                   >  1:322135/1‑93 (MQ=2)
                                                   cccccccGCGGGGGCGTGGATCGAAACA.GC.CT.GCGCG.GC.CGCA..TCGGCT.TCCGC.CTG..CGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCAGCtc                                                   <  2:2759261/100‑1 (MQ=255)
                                                   cccccccGCGGGGGCGTGGATCGAAAC..GC.TG.GCGCGCGC.CGCG..CCGGATCTTTGT.CTG..CGCG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAActcggcct                                                    <  1:4517929/100‑8 (MQ=9)
                                                    ccccccGCGGGGGCGTGGATCGAAACT.GC.GT.GCGCG.CCATGAAT.TCATCG.AGCGC.CAT..CGAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACACCaat                                                     <  1:313218/101‑3 (MQ=2)
                                                    ccccccGCGGGGGCGTGGATCGAAACT.GC.GT.GCGCG.CCATGAAT.TCATCG.AGCGC.CAT..CGAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACACCaat                                                     <  1:1519451/101‑3 (MQ=2)
                                                    ccccccGCGGGGGCGTGGATCGAAACA..A.CT.GACC...C.GGCA..ACGGCC.TATGCGCTGGACGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAActacaccggg                                                  <  2:4469526/101‑10 (MQ=2)
                                                     cccccGCGGGGGCGTGGATCGAAAC..GC.TG.GCGCGCGC.CGCG..CCGGATCTTTGT.CTG..CGCG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAActcggcctcg                                                  <  1:1301550/101‑10 (MQ=9)
                                                     cccccGCGGGGGCGTGGATCGAAAC..GC.TG.GCGCGCGC.CGCG..CCGGATCTTTGT.CTG..CGCG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAActcggcctcg                                                  <  2:1579696/101‑10 (MQ=9)
                                                      ccccGCGGGGGCGTGGATCGAAACT.GC.GT.GCGCG.CC.ATGAGTTCATCG.AGCGC.CAT..CGAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACACCaatgc                                                   <  2:4299223/101‑5 (MQ=0)
                                                       cccGCGGGGGCGTGGATCGAAACA.GC.CT.GCGCG.GC.CGCA..TCGGCT.TCCGC.CTG..CGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCAGCtcgtcg                                               <  2:2885350/101‑1 (MQ=255)
                                                         cgcgGGGGCGTGGATCGAAACT.GC.GT.GCGCG.CCATGAAT.TCATCG.AGCGC.CAT..CGAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACACCaatgccga                                                <  1:4192660/101‑8 (MQ=0)
                                                         cgcgGGGGCGTGGAGCGAACCA.GC.CT.GCGCG.GC.CGCA..TCGGCT.TCCGC.CTG..CGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCAGCTCGTCGtt                                             <  2:799502/101‑1 (MQ=255)
                                                          gggggggCGTGGATCGAAACA.GC.CT.GCGCG.GC.CCCA..TCGGCT.TCCGC.CTG..CGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCAGCTCGTCGTTg                                            <  1:2195171/99‑1 (MQ=255)
                                                          gcgGGGGCGTGGATCGAAACA.GC.CT.GCGCG.GC.CGCA..TCGGCT.TCCGC.CTG..CGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCAGCTCGTCGTTg                                            >  1:3748730/1‑101 (MQ=255)
                                                           cgGGGGCGTGGATCGAAACA.GC.CT.GCGCG.GC.CGCA..TCGGCT.TCCGC.CTG..CGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCAGCTCGTCGTTgg                                           <  1:2689497/101‑1 (MQ=255)
                                                           cgGGGGCGTGGATCGAAACA.GC.CT.GCGCG.GC.CGCA..TCGGCT.TCCGC.CTG..CGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCAGCTCGTCGTTgg                                           <  1:1394581/101‑1 (MQ=255)
                                                            gggggCGTGGATCGAAACA.GC.CT.GCGCG.GC.CGCA..TCGGCT.TCCGC.CTG..CGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCAGCTCGTCGTTGGt                                          <  2:2651635/101‑1 (MQ=255)
                                                             ggggGGTGGATCCAAACA.GC.CT.GCGCG.GC.CGCA..TCGGTT.TCCGC.CTG..CGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCAGCTCGTCGTTGGTc                                         <  1:4875229/101‑1 (MQ=255)
                                                             ggggCGTGGATCGAAACA.GC.CT.GCGCG.GC.CGCA..TCGGCT.TCCGC.CTG..CGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAAAATCAGCTCGTCGTTGGTc                                         >  1:3835259/1‑101 (MQ=255)
                                                                  gTGGATCGAAACA.GC.CT.GCGCG.GC.CGCA..TCGGCT.TCCGC.CTG..CGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCAGCTCGTCGTTGGTCAGCCg                                    >  1:1970120/1‑101 (MQ=255)
                                                                    ggATCGAAACA.GC.CT.GCGCG.GC.CGCA..TCGGCT.TCCGC.CTG..CGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCAGCTCGTCGTTGGTCAGCCGGa                                  >  1:3777720/1‑101 (MQ=255)
                                                                      aTCGAAACA.GC.CT.GCGCG.GC.CGCA..TCGGCT.TCCGC.CTG..CGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCAGCTCGTCGTTGGTCAGCCGGACg                                <  1:2350078/101‑1 (MQ=255)
                                                                       tCGAAACA.GC.CT.GCGCG.GC.CGCA..TCGGCT.TCCGC.CTG..CGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCAGCTCGTCGTTGGTCAGCCGGACGa                               >  2:1353300/1‑101 (MQ=255)
                                                                        cGAAACA.GC.CG.GCGCG.GC.CGCA..TCGGCT.TCCGC.CTG..CGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCAGCTCGTCGTTGGTCAGCCGGACGAt                              <  2:2069792/101‑1 (MQ=255)
                                                                         gAAACA.GC.CT.GCGCG.GC.CGCA..TCGGCT.TCCGC.CTG..CGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCAGCTCGTCGTTGGTCAGCCGGACGATg                             <  1:771636/101‑1 (MQ=255)
                                                                          aaaCA.GC.CT.GCGCG.GC.CGCA..TCGGCT.TCCGC.CTG..CGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCAGCTCGTCGTTGGTCAGCCGGACGATGt                            >  1:4963719/1‑101 (MQ=255)
                                                                            aCA.GC.CT.GCGCG.GC.CGCA..TCGGCT.TCCGC.CTG..CGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCAGCTCGTCGTTGGTCAGCCGGACGATGTcc                          <  1:1439366/101‑1 (MQ=255)
                                                                            aCA.GC.CT.GCGCG.GC.CGCA..TCGGCT.TCCGC.CTG..CGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCAGCTCGTCGTTGGTCAGCCGGACGATGTcc                          <  2:196544/101‑1 (MQ=255)
                                                                             cA.GC.CT.GCGCG.GC.CGCA..TCGGCT.TCCGC.CTG..CGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCAGCTCGTCGTTGGTCAGCCGGACGATGTCCg                         <  1:2595785/101‑1 (MQ=255)
                                                                              cccc.cT.GCGCG.GC.CGCA..TCGGCT.TCCGC.CTG..CGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCAGCTCGTCGTTGGTCAGCCGGACGATGTCCg                         <  2:3292606/98‑1 (MQ=255)
                                                                                gC.CT.GCGCG.GC.CGCA..TCGGCT.TCCGC.CTG..CGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCAGCTCGTCGTTGGTCAGCCGGACGATGTCCGTc                       <  2:2825696/101‑1 (MQ=255)
                                                                                   cT.GCGCG.GC.CGCA..TCGGCT.TCCGC.CTG..CGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCAGCTCGTCGTTGGTCAGCCGGACGATGTCCGTCGc                     >  1:4091246/1‑101 (MQ=255)
                                                                                   cT.GCGCG.GC.CGCA..TCGGCT.TCCGC.CTG..CGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCAGCTCGTCGTTGGTCAGCCGGACGATGTCCGTCGc                     >  2:3896688/1‑101 (MQ=255)
                                                                                      gggcg.GC.CGCA..TCGGCT.TCCGC.CTG..CGGG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCAGCTCGTCGTTGGTCAGCCGGACGATGTCCGTCGCTc                   <  2:3748452/99‑1 (MQ=255)
                                                                                       cgcg.GC.CGCA..TCGGCT.TCCGC.CTG..CGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCAGCTCGTCGTTGGTCAGCCGGACGATGTCCGTCGCTcc                  >  1:3527274/1‑101 (MQ=255)
                                                                                          g.gC.CGCA..TCGGCT.TCCGC.CTG..CGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCAGCTCGTCGTTGGTCAGCCGGACGATGTCCGTCGCTccccc               >  2:121791/1‑101 (MQ=255)
                                                                                            gC.CGCA..TCGGCT.TCCGC.CTG..CGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCAGCTCGTCGTTGGTCAGCCGGACGATGTCCGTCGCTcccccc              <  1:2427407/101‑1 (MQ=255)
                                                                                                gCA..TCGGCT.TCCGC.CTG..CGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCAGCTCGTCGTTGGTCAGCCGGACGATGTCCGTCGCTCCCCCcgcg           >  2:1396510/1‑101 (MQ=255)
                                                                                                gCA..TCGGCT.TCCGC.CTG..CGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCAGCTCGTCGTTGGTCAGCCGGACGATGTCCGTCGCTCCCCCcgcg           <  1:4434090/101‑1 (MQ=255)
                                                                                                 cA..TCGGCT.TCCGC.CTG..CGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCAGCTCGTCGTTGGTCAGCCGGACGATGTCCGTCGCTCCCCCCGCgg          >  1:2300340/1‑101 (MQ=255)
                                                                                                     tCGGCT.TCCGC.CTG..CGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCAGCTCGTCGTTGGTCAGCCGGACGATGTCCGTCGCTCCCCCCGCgggg        <  2:1894191/101‑1 (MQ=255)
                                                                                                      cGGCT.TCCGC.CTG..CGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCAGCTCGTCGTTGGTCAGCCGGACGATGTCCGTCGCTCCCCCCGCggggg       >  1:4604600/1‑101 (MQ=255)
                                                                                                       ggCT.TCCGC.CTG..CGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCAGCTCGTCGTTGGTCAGCCGGACGATGTCCGTCGCTCCCCCCGCGGGGGc      >  2:4591499/1‑101 (MQ=255)
                                                                                                          t.tCCGC.CTG..CGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCAGCTCGTCGTTGGTCAGCCTGACGATGTCCGTCGCTCCCGCCGCGGGGGCGTg   >  1:4344578/1‑101 (MQ=255)
                                                                                                            |                                                                                                        
TGGATCGAAACCGGTCCAGAGATAGACCGCGCGGACCTCCTCAGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACA.GC.CT.GCGCG.GC.CGCA..TCGGCT.TCCGC.CTG..CGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCAGCTCGTCGTTGGTCAGCCGGACGATGTCCGTCGCTCCCCCCGCGGGGGCGTGG  >  NC_014034/1396857‑1397057

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