Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_014034 1,396,953 G→A 17.3% G279S (GGC→AGC)  RCAP_rcc01307 → hypothetical protein

Read alignment evidence...
  seq id position change freq score reads annotation genes product
*NC_0140341,396,9530G→A17.3% 0.0 ‑2.6 ‑1.5698G279S (GGC→AGC) RCAP_rcc01307hypothetical protein
Reads supporting (aligned to +/- strand):  new base (50/70):  ref base (293/282):  total (344/354)
Fisher's exact test strand distribution p-value = 7.09e-02
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 3.36e-02

Only 100 of 727 total aligned reads displayed.

GGCGTGGATCGAAACCGGTCCAGAGATAGACCGCGCGGACCTCCTCAGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACA.GCCT.G.CGCG.GC.CGCA.TCGGCT.TCCGC.CTG..CGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCAGCTCGTCGTTGGTCAGCCGGACGATGTCCGTCGCTCCCCCCGCGGGGGC  >  NC_014034/1396853‑1397053
                                                                                                          |                                                                                                         
     ggATCGAAACCGGTCCAGAGATAGACCGCGCGGTCCTCCTCAGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACA.GCCT.G.CGCG.GC.CGCA.TCGGCT.Tc                                                                                                     >  1:6438986/1‑101 (MQ=255)
     ggATCGAAACCGGTCCAGAGATAGACCGCGCGGACCTCCTCAGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACA.GCCT.G.CGCG.GC.CGCA.TCGGCT.Tc                                                                                                     >  1:5822777/1‑101 (MQ=255)
     ggATCGAAACCGGTCCAGAGATAGACCGCGCGGACCTCCTCAGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACA.GCCT.G.CGCG.GC.CGCA.TCGGCT.Tc                                                                                                     <  2:6575066/101‑1 (MQ=255)
      gATCGAAACCGGTCCAGAGATAGACCGCGCGGACCTCCTCAGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACA.GCCT.G.CGCG.GC.CGCA.TCGGCT.Tcc                                                                                                    <  2:5871810/101‑1 (MQ=255)
       aTCGAAACCGGTCCAGAGATAGACCGCGCGGACCTCCTCAGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACA.GCCT.G.CGCG.GC.CGCA.TCGGCT.Tccc                                                                                                   >  2:8133332/1‑100 (MQ=255)
       aTCGAAACCGGTCCAGAGATAGACCGCGCGGACCTCCTCAGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACA.GCCT.G.CGCG.GC.CGCA.TCGGCT.TCCg                                                                                                   >  1:7570068/1‑101 (MQ=255)
         cGAAACCGGTCCAGAGATAGACCGCGCGGACCTCCTCAGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACA.GCCT.G.CGCG.GC.CGCA.TCGGCT.TCCGc.c                                                                                                >  2:67661/1‑101 (MQ=255)
          gAAACCGGTCCAGAGATAGACCGCGCGGACCTCCTCAGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACA.GCCT.G.CGCG.GC.CGCA.TCGGCT.TCCGC.Ct                                                                                               <  2:101209/101‑1 (MQ=255)
           aaaCCGGTCCAGAGATAGACCGCGCGGACCTCCTCAGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACA.GCCT.G.CGCG.GC.CGCA.TCGGCT.TCCGC.CTg                                                                                              >  2:6863829/1‑101 (MQ=255)
            aaCCGGTCCAGAGATAGACCGCGCGGACCTCCTCAGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACA.GCCT.G.CGCG.GC.CGCA.TCGGCT.TCCGC.CTG..c                                                                                           <  2:7496479/101‑1 (MQ=255)
             aCCGGTCCAGAGATAGACCGCGCGGACCTCCTCAGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACA.GCCT.G.CGCG.GC.CGCA.TCGGCT.TCCGC.CTG..Cg                                                                                          >  1:4839554/1‑101 (MQ=255)
                   ccAGATATAGGCCGCGCGGGCCTCCTCAGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACA.GCCT.G.CGCG.GC.CGCA.TCGGCT.TCCGC.CTG..CGTG.GTCg                                                                                   <  2:1416312/101‑1 (MQ=255)
                   ccAGAGATAGACCGCGCGGACCTCCTCAGGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACA.GCCT.G.CGCG.GC.CGCA.TCGGCT.TCCGC.CTG..CGTG.GTCg                                                                                   <  2:4875917/101‑1 (MQ=255)
                   ccAGAGATAGACCGCGCGGACCTCCTCAGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACA.GCCT.G.CGCG.GC.CGCA.TCGGCT.TCCGC.CTG..CGTG.GTCg                                                                                   >  2:8731556/1‑101 (MQ=255)
                     agagATAGACCGCGCGGACCTCCTCAGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACA.GCCT.G.CGCG.GC.CGCA.TCGGCT.TCCGC.CTG..CGTG.GTCGcc                                                                                 >  2:3052073/1‑100 (MQ=255)
                     agagATAGACCGCGCGGACCTCCTCAGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACA.GCCT.G.CGCG.GC.CGCA.TCGGCT.TCCGC.CTG..CGTG.GTCGCt                                                                                 <  1:8189984/101‑1 (MQ=255)
                           aGACCGCGCGGACCTCCTCAGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACA.GCCT.G.CGCG.GC.CGCA.TCGGCT.TCCGC.CTG..CGTG.GTCGCTcccccc                                                                           >  1:2164476/1‑101 (MQ=255)
                             aCCGCGCGGACCTCCTCAGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACA.GCCT.G.CGCG.GC.CGCA.TCGGCT.TCCGC.CTG..CGTG.GTCGCTCCCCCcgc                                                                         >  1:7414099/1‑101 (MQ=255)
                              ccGCGCGGACCTCCTCAGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACA.GCCT.G.CGCG.GC.CGCA.TCGGCT.TCCGC.CTG..CGTG.GTCGCTCCCCCcgcg                                                                        <  1:1974533/101‑1 (MQ=255)
                                  gcgGACCTCCTCAGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACA.GCCT.G.CGCG.GC.CGCA.TCGGCT.TCCGC.CTG..CGTG.GTCGCTCCCCCCGCggggg                                                                    <  2:8299118/101‑1 (MQ=255)
                                      aCCTCCTCAGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACA.GCCT.G.CGCG.GC.CGCA.TCGGCT.TCCGC.CTG..CGTG.GTCGCTCCCCCCGCGGGGGCGTg                                                                <  1:3977567/101‑1 (MQ=255)
                                       cctcctCAGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACA.GCCT.G.CGCG.GC.CGCA.TCGGCT.TCCGC.CTG..CGTG.GTCGCTCCCCCCGCGGGGGCGTgg                                                               <  2:4054113/101‑1 (MQ=255)
                                       cctcctCAGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACA.GCCT.G.CGCG.GC.CGCA.TCGGCT.TCCGC.CTG..CGTG.GTCGCTCCCCCCGCGGGGGCGTgg                                                               <  2:3940352/101‑1 (MQ=255)
                                        ggtcaacactgTCGCTCCCCCCGCGGGGGCGTGGATCGAAAC..GCTG.G.CGCGCGC.CGCG.CCGGATCTTTGT.CTG..CGCG.GTCGCTCCCCCCGCGGGGGCGTgg                                                               <  2:4642391/91‑1 (MQ=9)
                                        ctcctcAGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACA.GCCT.G.CGCG.GC.CGCA.TCGGCT.TCCGC.ATG..CGTG.GTCGCTCCCCCCGCGGGGGCGTGGa                                                              >  1:4990518/1‑101 (MQ=255)
                                         ttggccccggTCGCTCCCCCCGCGGGGGCGTGGATCGAAACA..ACT.G.ACC...C.GGCA.ACGGCC.TATGCGCTGGACGTG.GTCGCTCCCCCCGCGGGGGCGTGGAt                                                             <  2:7752115/92‑1 (MQ=2)
                                         ttggccccggTCGCTCCCCCCGCGGGGGCGTGGATCGAAACA..ACT.G.ACC...C.GGCA.ACGGCC.TATGCGCTGGACGTG.GTCGCTCCCCCCGCGGGGGCGTGGAt                                                             >  2:2841783/10‑101 (MQ=2)
                                         tcctcAGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACA.GCCT.G.CGCG.GC.CGCA.TCGGCT.TCCGC.CTG..CGTG.GTCGCTCCCCCCGCGGGGGCGTGGAt                                                             >  2:7050003/1‑101 (MQ=255)
                                          tggccccggTCGCTCCCCCCGCGGGGGCGTGGATCGAAACA..ACT.G.ACC...C.GGCA.ACGGCC.TATGCGCTGGACGTG.GTCGCTCCCCCCGCGGGGGCGTGGATc                                                            >  1:1236750/9‑101 (MQ=2)
                                          tggccccggTCGCTCCCCCCGCGGGGGCGTGGATCGAAACA..ACT.G.ACC...C.GGCA.ACGGCC.TATGCGCTGGACGTG.GTCGCTCCCCCCGCGGGGGCGTGGATc                                                            <  1:9159490/93‑1 (MQ=2)
                                          cctcAGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACA.GCCT.G.CGCG.GC.CGCA.TCGGCT.TCCGC.CTG..CGTG.GTCGCTCCCCCCGCGGGGGCGTGGATc                                                            >  2:128897/1‑101 (MQ=255)
                                           cggttgcGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATGGCTAT.GGCG.GC.AGTG.TCAGGT.T.TTC.CAG..CGGC.GTCGCTCCCCCCGCGGGGGCGTGGATc                                                            >  1:840943/7‑101 (MQ=2)
                                            tcAGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACA.GCCT.G.CGCG.GC.CGCA.TCGGCT.TCCGC.CTG..CGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGa                                                          >  1:4425440/1‑101 (MQ=255)
                                            ggttgcGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATGGCTAT.GGCG.GC.AGTG.TCAGGT.T.TTC.CAG..CGGC.GTCGCTCCCCCCGCGGGGGCGTGGATCg                                                           >  1:3065514/6‑101 (MQ=2)
                                            ccgcccGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GCGT.G.CGCG.CCATGAATTCATCG.AGCGC.CAT..CGAGCGTCGCTCCCCCCGCGGGGGCGTGGAt                                                             <  1:3376433/96‑1 (MQ=2)
                                            ccgcccGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GCGT.G.CGCG.CCATGAATTCATCG.AGCGC.CAT..CGAGCGTCGCTCCCCCCGCGGGGGCGTGGAt                                                             <  2:6397868/96‑1 (MQ=2)
                                            aacactgTCGCACCCCCCGCGGGGGCGTGGATCGAAAC..GCTG.G.CGCGCGC.CGCG.CCGGATCTTTGT.CTG..CGCG.GTCGCTCCCCCCGCGGGGGCGTGGATCg                                                           >  2:416270/7‑101 (MQ=9)
                                             cgcccGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GCGT.G.CGCG.CCATGAATTCATCG.AGCGC.CAT..CGAGCGTCGCTCCCCCCGCGGGGGCGTGGATc                                                            >  2:348861/5‑101 (MQ=2)
                                             cAGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACA.GCCT.G.CGCG.GC.CGCA.TCGGCT.TCCGC.CTG..CGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGaa                                                         >  2:8538471/1‑101 (MQ=255)
                                              ttgcGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATGGCTAT.GGCG.GC.AGTG.TCAGGT.T.TTC.CAG..CGGC.GTCGCTCCCCCCGCGGGGGCGTGGATCGaa                                                         <  2:3626005/98‑1 (MQ=2)
                                                cggTCGCTCCCCCCGCGGGGGCGTGGATCGAAACA..ACT.G.ACC...C.GGCA.ACGGCC.TATGCGCTGGACGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAAct                                                      <  2:5659945/99‑2 (MQ=2)
                                                 tgTCGCTCCCCCCGCGGGGGCGTGGATCGAAAC..GCTG.G.CGCGCGC.CGCG.CCGGATCTTTGT.CTG..CGCG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAAct                                                      >  1:7955316/2‑100 (MQ=9)
                                                  gTCGCTCCCCCCGCGGGGGCGTGGATCGTAACA..ACT.G.ACC...C.GGCA.ACGGCC.TATGCGCTGGACGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAActac                                                    <  1:1752065/101‑4 (MQ=2)
                                                  gTCGCTCCCCCCGCGGGGGCGTGGATCGAAACA..ACT.G.ACC...C.GGCA.ACGGCC.TATGCGCTGGACGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAActac                                                    <  2:2233003/101‑4 (MQ=2)
                                                  gTCGCTCCCCCCGCGGGGGCGTGGATCGAAACA..ACT.G.ACC...C.GGCA.ACGGCC.TATGCGCTGGACGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAActac                                                    <  2:2564852/101‑4 (MQ=2)
                                                   tctccccccccGCGGTGGCGTGGATCGAAACA.GCCT.G.CGCG.GC.CGCA.TCGGCT.TCCGC.CTG..CGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCa                                                   <  2:5160511/96‑1 (MQ=255)
                                                   tCGCTCCCCCCGCGGGGGCGTGGATCGAAACA.GCCT.G.CGCG.GC.CGCA.TCGGCT.TCCGC.CTG..CGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCa                                                   >  1:5566310/1‑101 (MQ=255)
                                                   tCGCCCCCCCCGCGGGGGCGTGGATCGAAACA.GCCT.G.CGCG.GC.CGCA.TCGGCT.TCCGC.CTG..CGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCa                                                   <  2:5389740/101‑1 (MQ=255)
                                                    cGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GCGT.G.CGCG.CCATGAATTCATCG.AGCGC.CAT..CGAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACac                                                     <  2:891981/101‑2 (MQ=2)
                                                    cGCTCCCCCCGCGGGGGCGTGGATCGAAACA..ACT.G.ACC...C.GGCA.ACGGCC.TATGCGCTGGACGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAActacac                                                  >  1:3523817/1‑96 (MQ=2)
                                                    cGCTCCCCCCGCGGGGGCGTGGATCGAAACA..ACT.G.ACC...C.GGCA.ACGGCC.TATGCGCTGGACGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAActacac                                                  >  1:8508117/1‑96 (MQ=2)
                                                    cGCTCCCCCCGCGGGGGCGTGGATCGAAAC..GCTG.G.CGCGCGC.CGCG.CCGGATCTTTGT.CTG..CGCG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAActcgg                                                   >  2:7185895/1‑97 (MQ=9)
                                                    cGCTCCCCCCGCGGGGGCGTGGATCGAAAC..GCTG.G.CGCGCGC.CGCG.CCGGATCTTTGT.CTG..CGCG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAActcgg                                                   >  1:4689965/1‑97 (MQ=9)
                                                     gCTCCCCCCGCGGGGGCGTGGATCGAAACA..ACT.G.ACC...C.GGCA.ACGGCC.TATGCGCTGGACGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAActacacc                                                 >  1:2377156/1‑95 (MQ=2)
                                                     ccTCCCCCCGCGGGGGCGTGGATCGAAACT.GCGT.G.CGCG.CCATGAATTCATCG.AGCGC.CAT..CGAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACacc                                                    <  2:9108524/100‑3 (MQ=2)
                                                      cTCCCCCCGCGGGGGCGTGGGTCGAAACT.GCGT.G.CGCG.CCATGAATTCATCG.AGCGC.CAT..CGAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAAcccca                                                   <  1:8808696/101‑5 (MQ=2)
                                                      cTCCCCCCGCGGGGGCGTGGATCGAAACT.GCGT.G.CGCG.CCATGAATTCATCG.AGCGC.CAT..CGAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACACCa                                                   <  2:7265493/101‑1 (MQ=2)
                                                      cTCCCCCCGCGGGGGCGTGGATCGAAACATGGCTAT.GGCG.GC.AGTG.TCAGGT.T.TTC.CAG..CGGC.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAAcccgccc                                                 >  1:8977033/1‑95 (MQ=2)
                                                      cTCCCCCCGCGGGGGCGTGGATCGAAACA.GCCT.G.CGCG.GC.CGCA.TCGGCT.TCCGC.CTG..CGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCAGCt                                                <  2:5789628/101‑1 (MQ=255)
                                                      cTCCCCCCGCGGGGGCGTGGATCGAAAC..GCTG.G.CGCGCGC.CGCG.CCGGATCTTTGT.CTG..CGCG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAActcggcc                                                 >  1:1643562/1‑95 (MQ=9)
                                                      cTCCCCCCGCGGGGGCGTGGATCGAAAC..GCTG.G.CGCGCGC.CGCG.CCGGATCTTTGT.CTG..CGCG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAActcggcc                                                 >  1:2771859/1‑95 (MQ=9)
                                                        ccccccGCGGGGGCGTGGATCGAAAC..GCTG.G.CGCGCGC.CGCG.CCGGATCTTTGT.CTG..CGCG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAActcggcctc                                               >  1:2659623/1‑93 (MQ=9)
                                                         cccccGCGGGGGCGTGGATCGAAACT.GCGT.G.CGCG.CCATGAATTCATCG.AGCGC.CAT..CGAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACACCaatg                                                <  1:5695275/101‑4 (MQ=2)
                                                         cccccGCGGGGGCGTGGATCGAAACATGGCTAT.GGCG.GC.AGTG.TCAGGT.T.TTC.CAG..CGGC.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAAcccgccacgt                                              <  2:3255287/101‑10 (MQ=2)
                                                         cccccGCGGGGGCGTGGATCGAAACA.GCCT.G.CGCG.GC.CGCA.TCGGCT.TCCGC.CTG..CGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCAGCtcgt                                             <  1:5340711/101‑1 (MQ=255)
                                                          ccccGCGGGGGCGTGGATCGAAACATGGCT.GTGGCG.GC.AGTG.TCAGGT.T.TTC.CAG..CGGC.GTCGCCCCCCCCGCGGGGGCGTGGATCGAAAcccgccacgtc                                             <  1:3793565/101‑11 (MQ=2)
                                                          ccccGCGGGGGCGTGGATCGAAACA.GCCT.G.CGCG.GC.CGCA.TCGGCT.TCCGC.CTG..CGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCAGCtcgtc                                            <  1:6785191/101‑1 (MQ=255)
                                                              gcgGGGGGGTGGATCGAAACA.GCCT.G.CGCG.GC.CGCA.TCGGCT.TCCGC.CTG..CGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCAGCTCGTCGTTg                                        <  1:4248635/101‑1 (MQ=255)
                                                               cgtgggCGTGGATCGAAACA.GCCT.G.CGCG.GC.CGCA.TCGGCT.TCCGC.CTG..CGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCAGCTCGTCGTTgg                                       <  2:7577439/98‑1 (MQ=255)
                                                                gggggCGTGGATCGAAACA.GCCT.G.CGCG.GC.CGCA.TCGGCT.TCCGC.CTG..CGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCAGCTCGTCGTTGGt                                      <  1:158091/101‑1 (MQ=255)
                                                                gggggCGTGGATCGAAACA.GCCT.G.CGCG.GC.CGAA.TCGGCT.TCCGC.CTG..CGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCAGCTCGTCGTTGGt                                      >  1:6989282/1‑101 (MQ=255)
                                                                 ggggCGTGGATCGAAACA.GCCT.G.CGCG.GC.CGCA.TCGGCT.TCCGC.CTG..CGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCAGCTCGTCGTTGGTc                                     <  2:2394063/101‑1 (MQ=255)
                                                                  gggCGTGGATCGAAACA.GCCT.G.CGCG.GC.CGCA.TCGGCT.TCCGC.CTG..CGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCAGCTCGTCGTTGGTCa                                    <  1:1013270/101‑1 (MQ=255)
                                                                   ggCGTGGATCGAAACA.GCCT.G.CGCG.GC.CGCA.TCGGCT.TCCGC.CTG..CGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCAGCTCGTCGTTGGTCAg                                   <  1:6210885/101‑1 (MQ=255)
                                                                      gTGGATCGAAACA.GCCT.G.CGCG.GC.CGCA.TCGGCT.TCCGC.CTG..CGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCAGCTCGTCGTTGGTCAGCCg                                <  1:8060737/101‑1 (MQ=255)
                                                                        ggATCGAAACA.GCCT.G.CGCG.GC.CGCA.TCGGCT.TCCGC.CTG..CGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCAGCTCGTCGTTGGTCCGCCGGa                              >  2:6146303/1‑101 (MQ=255)
                                                                        ggATCGAAACA.GCCT.G.CGCG.GC.CGCA.TCGGCT.TCCGC.CTG..CGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCAGCTCGTCGTTGGTCAGCCGGa                              >  1:6358391/1‑101 (MQ=255)
                                                                        ggATCGAAACA.GCCT.G.CGCG.GC.CGCA.TCGGCT.TCCGC.CTG..CGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCAGCTCGTCGTTGGGCAGCCGGa                              >  2:6696335/1‑101 (MQ=255)
                                                                         ggtCGAAACA.GCCT.G.CGCC.GC.CGCA.TCGGCC.TCCGC.CTG..CGTG.GTCGCGCCCCCCGCGGGGGCGTGGATCGAAACATCAGCTCGTCGTTGGTCAGCCGGAc                             <  2:4966164/99‑1 (MQ=255)
                                                                            cGAACCA.GCCT.G.CGCG.GC.CGCA.TCGGCT.TCCGC.CTG..CGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCAGCTCGTCGTTGGTCAGCCGGACGAt                          <  2:2055873/101‑1 (MQ=255)
                                                                            cGAAACA.GCCT.G.CGCG.GC.CGCA.TCGGCT.TCCGC.CTG..CGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCAGCTCGTCGTTGGTCAGCCGGACGAt                          >  1:1848443/1‑101 (MQ=255)
                                                                            cGAAACA.GCCT.G.CGCG.GC.CGCA.TCGGCT.TCCGC.CTG..CGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCAGCTCGTCGTTGGTCAGCCGGACGAt                          <  1:7777412/101‑1 (MQ=255)
                                                                             gAAACA.GCCT.G.CGCG.GC.CGCA.TCGGCT.TCCGC.CTG..CGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCAGCTCGTCGTTGGTCAGCCGGACGATg                         >  1:4092171/1‑101 (MQ=255)
                                                                             gAAACA.GCCT.G.CGCG.GC.CGCA.TCGGCT.TCCGC.CTG..CGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCAGCTCGTCGTTGGTCAGCCGGACGATg                         >  2:6016472/1‑101 (MQ=255)
                                                                               aaCA.GCCT.G.CGCG.GC.CGCA.TCGGCT.TCCGC.CTG..CGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCAGCTCGTCGTTGGTCAGCCGGACGATGTc                       >  1:8688806/1‑101 (MQ=255)
                                                                                  a.GCCT.G.CGCG.GC.CGCA.TCGGCT.TCCGC.CTG..CGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCAGCTCGTCGTTGGTCAGCCGGACGATGTCCGt                    <  2:7876987/101‑1 (MQ=255)
                                                                                      cT.G.CGCG.GC.CGCC.TCGGCT.TCCGC.CTG..CGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAAGATCAGCTCGTCGTTGGTCAGCCGGACGATGTCCGTCGc                 >  1:6550690/1‑101 (MQ=255)
                                                                                      cT.G.CGCG.GC.CGCA.TCGGCT.TCCGC.CTG..CGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCAGCTCGTCGTTGGTCAGCCGGACGATGTCCGTCGc                 >  1:6634942/1‑101 (MQ=255)
                                                                                      cT.G.CGCG.GC.CGCA.TCGGCT.TCCGC.CTG..CGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCAGCTCGTCGTTGGTCAGCCGGACGATGTCCGTCGc                 >  1:2467130/1‑101 (MQ=255)
                                                                                       ggg.cgcG.GC.CGCA.TCGGCT.TCCGC.CTG..CGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCAGCTCGTCGTTGGTCAGCCGGACGATGTCCGTCGc                 >  2:3367182/3‑101 (MQ=255)
                                                                                         g.cgcG.GC.CGCA.TCGGCT.TCCTC.CTG..CGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCAGCTCGTCGTTGGTCAGCCGGACGATGTCCGTCGCTc               >  1:7259857/1‑101 (MQ=255)
                                                                                         g.cgcG.GC.CGCA.TCGGCT.TCCGC.CTG..CGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCAGCTCGTCGTTGGTCAGCCGGACGATGTCCGTCGCTc               <  2:1108040/101‑1 (MQ=255)
                                                                                                gC.CGCA.TCGGCT.TCCGC.CTG..CGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCAGCTCGTCGTTGGTCAGCCGGACGATGTCCGTCGCTcccccc          <  2:1586497/101‑1 (MQ=255)
                                                                                                 c.cGCA.TCGGCT.TCCGC.CTG..CGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCAGCTCGTCGTTGGTCAGCCGGACGATGTCCGTCGCTCCCCCcg         <  2:1846291/101‑1 (MQ=255)
                                                                                                   cGCA.TCGGCT.TCCGC.CTG..CGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCAGCTCGTCGTTGGTCAGCCGGACGATGTCCGTCGCTCCCCCcgc        <  1:6016707/101‑1 (MQ=255)
                                                                                                    gCA.TCGGCT.TCCGC.CTG..CGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCAGCTCGTCGTTGGTCAGCCGGACGATGTCCGTCGCTCCCCCcgcg       >  1:3236605/1‑101 (MQ=255)
                                                                                                    gCA.TCGGCT.TCCGC.CTG..CGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCAGCTCGTCGTTGGTCAGCCGGACGATGTCCGTCGCTCCCCCcgcg       >  1:5215369/1‑101 (MQ=255)
                                                                                                    gCA.TCGGCT.TCCGC.CTG..CGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCAGCTCGTCGTTGGTCAGCCGGACGATGTCCGTCGCTCCCCCcgcg       <  1:7139032/101‑1 (MQ=255)
                                                                                                         cGGCT.TCCGC.CTG..CGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCAGCTCGTCGTTGGTCAGCCGGACGATGTCCGTCGCTCCCCCCGCggggg   <  2:7375655/101‑1 (MQ=255)
                                                                                                         cGGCT.TCCGC.CTG..CGTG.GTCGCCCCCCCCGCGGGGGCGTGGATCGAAACATCAGCTCGTCGTTGGTCAGCCGGACGATGTCCGTCGCTCCCCCCGCggggg   <  2:4425313/101‑1 (MQ=255)
                                                                                                          |                                                                                                         
GGCGTGGATCGAAACCGGTCCAGAGATAGACCGCGCGGACCTCCTCAGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACA.GCCT.G.CGCG.GC.CGCA.TCGGCT.TCCGC.CTG..CGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCAGCTCGTCGTTGGTCAGCCGGACGATGTCCGTCGCTCCCCCCGCGGGGGC  >  NC_014034/1396853‑1397053

Base quality scores:  ATCG  < 3 ≤  ATCG  < 20 ≤  ATCG  < 27 ≤  ATCG  < 33 ≤  ATCG  < 40 ≤  ATCG