Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_014034 1,396,956 +C 11.9% coding (838/1578 nt) RCAP_rcc01307 → hypothetical protein

Read alignment evidence...
  seq id position change freq score reads annotation genes product
*NC_0140341,396,9561.→C11.9% 159.5 ‑0.1 ‑1.8402coding (838/1578 nt)RCAP_rcc01307hypothetical protein
Reads supporting (aligned to +/- strand):  new base (24/23):  ref base (173/181):  total (198/204)
Fisher's exact test strand distribution p-value = 8.77e-01
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 1.48e-02

Only 100 of 446 total aligned reads displayed.

GTGGATCGAAACCGGTCCAGAGATAGACCGCGCGGACCTCCTCAGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACA.GC.CT.GCGCG.GC.CGCA..TCGGCT.TCCGC.CTG..CGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCAGCTCGTCGTTGGTCAGCCGGACGATGTCCGTCGCTCCCCCCGCGGGGGCGTG  >  NC_014034/1396856‑1397056
                                                                                                            |                                                                                                        
gTGGATCGAAACCGGTCCAGAGATAGACCGCGCGGACCTCCTCAGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACA.GC.CT.GCGCG.GC.CGCA..TCGGCt                                                                                                           >  2:2364252/1‑101 (MQ=255)
 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)
        aaaCCGGTCCAGAGATAGACCGCGCGGACCTCCTCAGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACA.GC.CT.GCGCG.GC.CGCA..TCGGCT.TCCGC.CTg                                                                                                 <  2:3108831/101‑1 (MQ=255)
         aaCCGGTCCAGAGATAGACCGCGCGGACCTCCTCAGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACA.GC.CT.GCGCG.GC.CGCA..TCGGCT.TCCGC.CTG..c                                                                                              >  1:4145564/1‑101 (MQ=255)
         aaCCGGTCCAGAGATAGACCGCGCGGACCTCCTCAGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACA.GC.CT.GCGCG.GC.CGCA..TCGGCT.TCCCC.CCG..c                                                                                              >  1:3920367/1‑101 (MQ=255)
            cGGTCCAGAGATAGACCGCGCGGACCTCCTCAGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACA.GC.CT.GCGCG.GC.CGCA..TCGGCT.TCCGC.CTG..CGTg                                                                                           >  1:3662621/1‑101 (MQ=255)
            cGGTCCAGAGATAGACCGCGCGGACCTCCTCAGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACA.GC.CT.GCGCG.GC.CGCA..TCGGCT.TCCGC.CTG..CGTg                                                                                           <  2:3053645/101‑1 (MQ=255)
                  agagATAGACCGCGCGGACCTCCTCAGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACA.GC.CT.GCGCG.GC.CGCA..TCGGCT.TCCGC.CTG..CGTG.GTCGCt                                                                                    <  2:3512569/101‑1 (MQ=255)
                   gagaTAGACCGCGCGGACCTCCTCAGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACA.GC.CT.GCGCG.GC.CGCA..TCGGCT.TCCGC.CTG..CGTG.GTCGCTc                                                                                   >  2:1116184/1‑101 (MQ=255)
                        aGACCGCGCGGACCTCCTCAGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACA.GC.CT.GCGCG.GC.CGCA..TCGGCT.TCCGC.CTG..CGTG.GTCGCTcccccc                                                                              <  2:3374768/101‑1 (MQ=255)
                        aGACCGCGCGGACCTCCTCAGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACA.GC.CT.GCGCG.GC.CGCA..TCGGCT.TCCGC.CTG..CGTG.GTCGCTcccccc                                                                              <  2:228429/101‑1 (MQ=255)
                         gACCGCGCGGACCTCCTCAGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACA.GC.CT.GCGCG.GC.CGCA..TCGGCT.TCCGC.CTG..CGTG.GTCGCTCCCCCcg                                                                             <  1:2897212/101‑1 (MQ=255)
                          cccGCGCGGACCTCCTCAGGCGTCGCTCCCCCCGCGGGGGCGTGGAGCGAAACA.GC.CT.GCGCG.GC.CGCA..TCGGCT.TCCGC.CTG..CGTG.GTCGCGCCCGCcgc                                                                            <  2:1878481/100‑1 (MQ=38)
                              cgcgGACCTCCTCAGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACA.GC.CT.GCGCG.GC.CGCA..TCGGCT.TCCGC.CTG..CGTG.GTCGCTCCCCCCGCgggg                                                                        <  1:1988471/101‑1 (MQ=255)
                                cgGACCTCCTCAGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACA.GC.CT.GCGCG.GC.CGCA..TCGGCT.TCCGC.CTC..CGTG.GTCGCTCCCCCCGCGGGGGc                                                                      >  1:2130045/1‑101 (MQ=255)
                                cgGACCTCCTCAGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACA.GC.CT.GCGCG.GC.CGCA..TCGGCT.TCCGC.CTC..CGTG.GTCGCTCCCCCCGCGGGGGc                                                                      >  1:4994181/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                                                                  <  1:4696864/92‑1 (MQ=0)
                                     tcgcggttgcGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATGG.CTATGGCG.GC.AGTG..TCAGGT.T.TTC.CAG..CGGC.GTCGCTCCCCCCGCGGGGGCGTgg                                                                  <  2:409452/92‑1 (MQ=0)
                                     attggccccggTCGCTCCCCCCGCGGGGGCGTGGATCGAAACA..A.CT.GACC...C.GGCA..ACGGCC.TATGCGCTGGACGTG.GTCGCTCCCCCCGCGGGGGCGTGGa                                                                 <  2:1157589/91‑1 (MQ=0)
                                      ttggccccggTCGCTCCCCCCGCGGGGGCGTGGATCGAAACA..A.CT.GACC...C.GGCA..ACGGCC.TATGCGCTGGACGTG.GTCGCTCCCCCCGCGGGGGCGTGGAt                                                                >  1:4088813/10‑101 (MQ=2)
                                      ccctcAGACGTCGCCCCCCCCGGGGGGGCGTGGATCGAAACA.GC.CT.GCGCG.GC.CGCA..TCGGCT.TCCGC.CTG..CGTG.GTCGCTCCCCCCGC.GGGGCGTGGATc                                                               <  2:41905/100‑1 (MQ=38)
                                       tggccccggTCGCTCCCCCCGCGGGGGCGTGGATCGAAACA..A.CT.GACC...C.GGCA..ACGGCC.TATGCGCTGGACGTG.GTCGCTCCCCCCGCGGGGGCGTGGATc                                                               <  1:3217404/93‑1 (MQ=2)
                                       tggccccggTCGCTCCCCCCGCGGGGGCGTGGATCGAAACA..A.CT.GACC...C.GGCA..ACGGCC.TATGCGCTGGACGTG.GTCGCTCCCCCCGCGGGGGCGTGGATc                                                               <  1:629030/93‑1 (MQ=2)
                                       cctcAGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACA.GC.CT.GCGCG.GC.CGCA..TCGGCT.TCCGC.CTG..CGTG.GTCGCTCCCCCCGCGGGGGCGTGGATc                                                               >  2:2606539/1‑101 (MQ=255)
                                        ctcAGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACA.GC.CT.GCGCG.GC.CGCA..TCGGTT.TCCGC.CTG..CGTG.GTCGCTCCCCCCGCGGGGGCGTGGAGCg                                                              >  1:3701711/1‑101 (MQ=255)
                                        ctcAGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACA.GC.CT.GCGCG.GC.CGCA..TCGGCT.TCCGC.CTG..CGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCg                                                              <  2:1054391/101‑1 (MQ=255)
                                         tcAGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACA.GC.CT.GCGCG.GC.CGCA..TCGGCT.TCCGC.CTG..CGTG.GTCGCTCCCCCCGCGGGGGCGGGGATCGa                                                             >  1:3802797/1‑101 (MQ=255)
                                         aacactgTCGCTCCCCCCGCGGGGGCGTGGATCGAAAC..GC.TG.GCGCGCGC.CGCG..CCGGATCTTTGT.CTG..CGCG.GTCGCTCCCCCCGCGGGGGCGTGGATCg                                                              >  1:4290267/7‑101 (MQ=9)
                                           aGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACA.GC.CT.GCGCG.GC.CGCA..TCGGCT.TCCGC.CTG..CGTG.GTCGCTCCCCCCGCGGGGGGGGGGGTCGaaa                                                           >  1:20549/1‑101 (MQ=255)
                                            gACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACA.GC.CT.GCGCG.GC.CGCA..TCGGCT.TCCGC.CTG..CGTG.GTTGCTCCCCCCGCGGGGGCGGGGATCGAAAc                                                          >  1:747289/1‑101 (MQ=255)
                                              tgTCGCTCCCCCCGCGGGGGCGTGGATCGAAAC..GC.TG.GCGCGCGC.CGCG..CCGGATCTTTGT.CTG..CGCG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAAct                                                         >  1:1780662/2‑100 (MQ=9)
                                              cGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GC.GT.GCGCG.CCATGAAT.TCATCG.AGCGC.CAT..CGAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGaaa                                                           >  2:477790/1‑101 (MQ=2)
                                              cGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GC.GT.GCGCG.CCATGAAT.TCATCG.AGCGC.CAT..CGAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGaaa                                                           <  1:2838087/101‑1 (MQ=2)
                                              cGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GC.GT.GCGCG.CCATGAAT.TCATCG.AGCGC.CAT..CGAGCGGCGCTCCCCCCGCGGGGCCGTGGATCGaaa                                                           >  2:1677022/1‑101 (MQ=2)
                                               gTCGCTCCCCCCGCGGGGGCGTGGATCGAAAC..GC.TG.GCGCGCGC.CGCG..CCGGATCTTTGT.CTG..CGCG.GTCGCTCCCCCCGCGGGGGCGTGGGTCGAAActc                                                        >  1:1057255/1‑99 (MQ=9)
                                               gTCGCTCCCCCCGCGGGGGCGTGGATCGAAAC..GC.TG.GCGCGCGC.CGCG..CCGGATCTTTGT.CTG..CGCG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAActc                                                        <  1:989258/101‑3 (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)
                                                tCGCTCCCCCCGCGGGGGCGTGGATCGAAACA.GC.CT.GCGCG.GC.CGCA..TCGGCT.TCCGC.CTG..CGTG.GTCGCTCCCCCCGCGGGGGCGTGGGTCGAAACATCa                                                      >  1:3577883/1‑101 (MQ=255)
                                                tCGCTCCCCCCGCGGGGGCGTGGATCGAAACA.GC.CT.GCGCG.GC.CGCA..TCGGCT.TCCGC.CTG..CGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCa                                                      <  1:2667094/101‑1 (MQ=255)
                                                tCGCTCCCCCCGCGGGGGCGTGGATCGAAACA.GC.CT.GCGCG.GC.CGCA..TCGGCT.TCCGC.CTG..CGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCa                                                      >  1:4860951/1‑101 (MQ=255)
                                                tCGCTCCCCCCGCGGGGGCGTGGATCGAAACA..A.CT.GACC...C.GGCA..ACGGCC.TATGCGCTGGACGTG.GTCGCTCCCACCGCGGGGGCGTGGATCGAAActaca                                                      >  1:4238981/1‑97 (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.GC.CT.GCGCG.GC.CGCA..TCGGCT.TCCGC.CTG..CGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCAg                                                     >  2:1351571/1‑101 (MQ=255)
                                                 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)
                                                  gCTCCCCCCGCGGGGGCGTGGATCGAAACT.GC.GT.GCGCG.CCATGAAT.TCATCG.AGCGC.CAT..CGAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACacc                                                       >  1:4938516/1‑99 (MQ=2)
                                                  gCTCCCCCCGCGGGGGCGTGGATCGAAACA..A.CT.GACC...C.GGCA..ACGGCC.TATGCGCTGGACGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAActacacc                                                    >  1:3393917/1‑95 (MQ=2)
                                                  gCTCCCCCCGCGGGGGCGTGGATCGAAACA..A.CT.GACC...C.GGCA..ACGGCC.TATGCGCTGGACGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAActacacc                                                    <  1:44290/101‑7 (MQ=2)
                                                  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)
                                                   ccccccccGCGGGGGCGTGGATCGAAACT.GC.GT.GCGCG.CCATGAAT.TCATCG.AGCGC.CAT..CGAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACACCa                                                      <  1:4471075/99‑1 (MQ=2)
                                                   cTCCCCCCGCGGGGGCGTGGATCGAAACT.GC.GT.GCGCG.CCATGAAT.TCATCG.AGCGC.CAT..CGAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACACCa                                                      <  1:4542812/101‑1 (MQ=2)
                                                   cTCCCCCCGCGGGGGCGTGGATCGAAACA.GC.CT.GCGCG.GC.CGCA..TCGGCT.TCCGC.CTG..CGTG.GTCGCTCCCTCCGCGGGGGCGTGGATCGAAACATCAGCt                                                   <  2:468931/101‑1 (MQ=255)
                                                   cTCCCCCCGCGGGGGCGTGGATCGAAAC..GC.TG.GCGCGCGC.CGCG..CCGGATCTTTGT.CTG..CGCG.GTCGCTCCCCCCGCGGGGGCGTGGATGGAAActcggcc                                                    >  2:319463/1‑95 (MQ=9)
                                                    tCCCCCCGCGGGGGCGTGGATCGAAACT.GC.GT.GCGCG.CCATGAAT.TCATCG.AGCGC.CAT..CGAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACACCaa                                                     <  1:2944103/101‑2 (MQ=2)
                                                    cccccccGCGGGGGCGTGGATCGAAACA.GC.CT.GCGCG.GC.CGCA..TCGGCT.TCCGC.CTG..CGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCAGCtc                                                  <  2:2759261/100‑1 (MQ=255)
                                                     ccccccGCGGGGGCGTGGATCGAAACT.GC.GT.GCGCG.CCATGAAT.TCATCG.AGCGC.CAT..CGAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACACCaat                                                    <  1:1519451/101‑3 (MQ=2)
                                                     ccccccGCGGGGGCGTGGATCGAAACATGG.CTATGGCG.GC.AGTG..TCAGGT.T.TTC.CAG..CGGC.GTCGCTCCCCCCGCGGGGGCGTGGATCGaaaaccgccccg                                                  >  2:1446/1‑92 (MQ=2)
                                                     ccccccGCGGGGGCGTGGATCGAAACA..A.CT.GACC...C.GGCA..ACGGCC.TATGCGCTGGACGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAActacaccggg                                                 <  2:4469526/101‑10 (MQ=2)
                                                     ccccccGCGGGGGCGTGGATCGAAAC..GC.TG.GCGCGCGC.CGCG..CCGGATCTTTGT.CTG..CGCG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAActcggcctc                                                  >  1:903897/1‑93 (MQ=9)
                                                      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)
                                                         ccGGGGGGGCGTGGATCGAAACA.GC.CT.GCGCG.GC.CGCA..TCGGCT.TCCGC.CTG..CGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCAGCTcgtcgt                                             <  1:682592/101‑1 (MQ=255)
                                                          cgcgGGGGCGTGGATCGAAACA.GC.CT.GCGCG.GC.CGCA..TCGGCT.TCCGC.CTG..CGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCAGCTCGTCGtt                                            <  1:3333269/101‑1 (MQ=255)
                                                          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)
                                                            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)
                                                             gggggCGTGGATCGAAACA.GC.CT.GCGCG.GC.CGCA..TCGGCT.TCCGC.CTG..CGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCAGCTCGTCGTTGGt                                         >  1:4004196/1‑101 (MQ=255)
                                                                  cGTGGATCGAAACA.GC.CT.GCGCG.GC.CGCA..TCGGCT.TCCGC.CTG..CGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCAGCTCGTCGTTGGTCAGcc                                    <  1:3496113/101‑1 (MQ=255)
                                                                   gTGGATCGAAACA.GC.CT.GCGCG.GC.CGCA..TCGGCT.TCCGC.CTG..CGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCAGCTCCTCGGTGGGCAGCCg                                   >  2:4335797/1‑101 (MQ=255)
                                                                       aTCGAAACA.GC.CT.GCGCG.GC.CGCA..TCGGCT.TCCGC.CTG..CGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCAGCTCGTCGTTGGTCAGCCGGACg                               <  1:2350078/101‑1 (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:2536770/1‑101 (MQ=255)
                                                                          gAAACA.GC.CT.GCGCG.GC.CGCA..TCGGCT.TCCGC.CTG..CGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCAGCTCGTCGTTGGTCAGCCGGACGATg                            >  1:309201/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                         >  1:377430/1‑101 (MQ=255)
                                                                             aCA.GC.CT.GCGCG.GC.CGCA..TCGGCT.TCCGC.CTG..CGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCAGCTCGTCGTTGGTCAGCCGGACGATGTcc                         <  2:196544/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)
                                                                                 gC.CT.GCGCG.GC.CGCA..TCGGCT.TCCGC.CTG..CGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCAGCTCGTCGTTGGTCAGCCGGACGATGTCCGTc                      <  1:2753155/101‑1 (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)
                                                                                             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)
                                                                                                  cA..TCGGCT.TCCGC.CTG..CGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCAGCTCGTCGTTGGTCAGCCGGACGATGTCCGTCGCTCCCCCCGCgg         >  1:2300340/1‑101 (MQ=255)
                                                                                                        ggCT.TCCGC.CTG..CGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCAGCTCGTCGTTGGTCAGCCGGACGATGTCCGTCGCTCCCCCCGCGGGGGc     >  2:4591499/1‑101 (MQ=255)
                                                                                                        ggCT.TCCGC.CTG..CGTG.GTCGCTCCCCCCGCGGGAGCGTGGATCGAAACATCAGCTCGTCGTTGGTCAGCCGGACGATGTCCGTCGCTCCCCCCGCGGGGGc     <  1:4368236/101‑1 (MQ=255)
                                                                                                            |                                                                                                        
GTGGATCGAAACCGGTCCAGAGATAGACCGCGCGGACCTCCTCAGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACA.GC.CT.GCGCG.GC.CGCA..TCGGCT.TCCGC.CTG..CGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCAGCTCGTCGTTGGTCAGCCGGACGATGTCCGTCGCTCCCCCCGCGGGGGCGTG  >  NC_014034/1396856‑1397056

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