Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_014034 1,396,955 C→A 11.8% G279G (GGC→ATA RCAP_rcc01307 → hypothetical protein

Read alignment evidence...
  seq id position change freq score reads annotation genes product
*NC_0140341,396,9550C→A11.8% 109.5 ‑0.3 0.0451G279G (GGC→GGARCAP_rcc01307hypothetical protein
Reads supporting (aligned to +/- strand):  new base (27/22):  ref base (191/173):  total (237/214)
Fisher's exact test strand distribution p-value = 7.62e-01
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 1.00e+00

Only 100 of 476 total aligned reads displayed.

CGTGGATCGAAACCGGTCCAGAGATAGACCGCGCGGACCTCCTCAGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACA.GCCT.GCGCG.GC.CGCA.TCGG.CT.TCCGC.CTG..CGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCAGCTCGTCGTTGGTCAGCCGGACGATGTCCGTCGCTCCCCCCGCGGGGGCG  >  NC_014034/1396855‑1397054
                                                                                                          |                                                                                                        
 gTGGATCGAAACCGGTCCAGAGATAGACCGCGCGGACCTCCTCAGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACA.GCCT.GCGCG.GC.CGCA.TCGG.Ct                                                                                                         >  2:1579379/1‑101 (MQ=255)
   ggATCGAAACCGGTCCAGAGATAGACCGCGCGGACCTCCTCAGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACA.GCCT.GCGCG.GC.CGCA.TCGG.CT.Tc                                                                                                      >  1:1639557/1‑101 (MQ=255)
    gATCGAAACCGGTCCAGAGATAGACCGCGCGGACCTCCTCAGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACA.GCCT.GCGCG.GC.CGCA.TCGG.CT.Tcc                                                                                                     >  2:4818366/1‑101 (MQ=255)
       cGAAACCGGTCCAGAGATAGACCGCGCGGACCTCCTCAGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACA.GCCT.GCGCG.GC.CGCA.TCGG.CT.TCCGc.c                                                                                                 >  1:5010490/1‑101 (MQ=255)
       cGAAACCGGTCCAGAGATAGACCGCGCGGACCTCCTCAGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACA.GCCT.GCGCG.GC.CGCA.TCGG.CT.TCCGc.c                                                                                                 >  2:3269746/1‑101 (MQ=255)
        gAAACCGGTCCAGAGATAGACCGCGCGGACCTCCTCAGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACA.GCCT.GCGCG.GC.CGCA.TCGG.CT.TCCGC.Ct                                                                                                >  1:3632926/1‑101 (MQ=255)
         aaaCCGGTCCAGAGATAGACCGCGCGGACCTCCTCAGATGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACA.GCCT.GCGCG.GC.CGCA.TCGG.CT.TCCGC.CTg                                                                                               >  2:3572558/1‑101 (MQ=255)
         aaaCCGGTCCAGAGATAGACCGCGCGGACCTCCTCAGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACA.GCCT.GCGCG.GC.CGCA.TCGG.CT.TCCGC.CTg                                                                                               <  1:3935478/101‑1 (MQ=255)
         aaaCCGGTCCAGAGATAGACCGCGCGGACCTCCTCAGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACA.GCCT.GCGCG.GC.CGCA.TCGG.CT.TCCGC.CTg                                                                                               <  2:2011086/101‑1 (MQ=255)
          aaCCGGTCCAGAGATAGACCGCGCGGACCTCCTCAGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACA.GCCT.GCGCG.GC.CGCA.TCGG.CT.TCCGC.CTG..c                                                                                            <  2:4452353/101‑1 (MQ=255)
           aCCGGTCCAGAGATAGACCGCGCGGACCTCCTCAGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACA.GCCT.GCGCG.GC.CGCA.TCGG.CT.TCCGC.CTG..Cg                                                                                           <  1:4402696/101‑1 (MQ=255)
             cGGTCCAGAGATAGACCGCGCGGACCTCCTCAGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACA.GCCT.GCGCG.GC.CGCA.TCGG.CT.TCCGC.CTG..CGTg                                                                                         >  2:2975354/1‑101 (MQ=255)
              ggTCCAGAGATAGACCGCGCGGACCTCCTCAGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACA.GCCT.GCGCG.GC.CGCA.TCGG.CT.TCCGC.CTG..CGTg.g                                                                                       <  1:1164154/101‑1 (MQ=255)
                  cAGAGATAGACCGCGCGGACCTCCTCAGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACA.GCCT.GCGCG.GC.CGCA.TCGG.CT.TCCGC.CTG..CGTG.GTCGc                                                                                   >  1:5232273/1‑101 (MQ=255)
                    gagaTAGACCGCGCGGCCCTCCTCAGACGTCGCGCCCCCCGCGGGGGCGTGGATCGAAACA.GCCT.GCGCG.GC.CGCA.TCGG.CT.TCCGC.CTG..CGTG.GTCGCTc                                                                                 <  1:3080646/101‑1 (MQ=255)
                     agaTAGACCGCGCGGACCTCCTCAGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACA.GCCT.GCGCG.GC.CGCA.TCGG.CT.TCCGC.CTG..CGTG.GTCGCTcc                                                                                >  1:4571478/1‑101 (MQ=255)
                       aTAGACCGCGCGGACCTCCTCAGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACA.GCCT.GCGCG.GC.CGCA.TCGG.CT.TCCGC.CTG..CGTG.GTCGCTcccc                                                                              <  1:2961852/101‑1 (MQ=255)
                        tAGACCGCGCGGACCTCCTCAGACGTCGCTCCCACCGCGGGGGCGTGGATCGAAACA.GCCT.GCGCG.GC.CGCA.TCGG.CT.TCCGC.CTG..CGTG.GTCGCTccccc                                                                             >  2:5213149/1‑101 (MQ=255)
                          gACCGCGCGGACCTCCTCAGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACA.GCCT.GCGCG.GC.CGCA.TCGG.CT.TCCGC.CTG..CGTG.GTCGCTCCCCCcg                                                                           <  2:1793242/101‑1 (MQ=255)
                           aCCGCGCGGACCTCCTCAGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACA.GCCT.GCGCG.GC.CGCA.TCGG.CT.TCCGC.CTG..CGTG.GTCGCTCCCCCcgc                                                                          >  2:3474830/1‑101 (MQ=255)
                                   gACCTCCTCAGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACA.GCCT.GCGCG.GC.CGCA.TCGG.CT.TCCGC.CTG..CGTG.GTCGCTCCCCCCGCGGGGGCgg                                                                  >  2:3082602/1‑100 (MQ=255)
                                     cctcctCAGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACA.GCCT.GCGCG.GC.CGCA.TCGG.CT.TCCGC.CTG..CGTG.GTCGCTCCCCCCGCGGGGGCGTgg                                                                <  2:1324988/101‑1 (MQ=255)
                                      ctcctcAGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACA.GCCT.GCGCG.GC.CGCA.TCGG.CT.TCCGC.CTG..CGTG.GTCGCTCCCCCCGCGGGGGCGTGGa                                                               <  1:852556/101‑1 (MQ=255)
                                       tcctcAGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACA.GCCT.GCGCG.GC.CGCA.TCGG.CT.TCCGC.CTG..CGTG.GTCGCTCCCCCCGCGGGGGCGTGGAt                                                              <  1:2499455/101‑1 (MQ=255)
                                       tcctcAGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACA.GCCT.GCGCG.GC.CGCA.TCGG.CT.TCCGC.CTG..CGTG.GTCGCTCCCCCCGCGGGGGCGTGGAt                                                              <  2:3358717/101‑1 (MQ=255)
                                       ccctcAGATGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACA.GCCT.GCGCG.GC.CGCA.TCGG.CT.TCCGC.CTG..CGTG.GTCGCTCCCCCCGCGGGGGCGTGGAt                                                              <  2:213819/100‑1 (MQ=255)
                                        tggccccggTCGCTCCCCCCGCGGGGGCGTGGATCGAAACA..ACT.GACC...C.GGCA.ACGG.CC.TATGCGCTGGACGTG.GTCGCTCCCCCCGCGGGGGCGTGGATc                                                             >  2:5744647/9‑101 (MQ=2)
                                        tcaacactgTCGCTCCCCCCGCGGGGGCGTGGATCGAAAC..GCTG.GCGCGCGC.CGCG.CCGG.ATCTTTGT.CTG..CGCG.GTCGCTCCCCCCGCGGGGGCGTGGAt                                                              <  2:3737644/93‑1 (MQ=9)
                                        cctcAGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACA.GCCT.GCGCG.GC.CGCA.TCGG.CT.TCCGC.CTG..CGTG.GTCGCTCCCCCCGCGGGGGCGTGGATc                                                             <  2:854770/101‑1 (MQ=255)
                                        cctcAGACGTCGCCCCCCCCGCGGGGGCGTGGATCGAAACA.GCCT.GCGCG.GC.CGCA.TCGG.CT.TCCGC.CTG..CGTG.GTCGCTCCCCCCGCGGGGGCGTGGATc                                                             <  1:2744148/101‑1 (MQ=255)
                                         ggccccggTCGCTCCCCCCGCGGGGGCGTGGATCGAAACA..ACT.GACC...C.GGCA.ACGG.CC.TATGCGCTGGACGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCg                                                            >  1:3719194/8‑101 (MQ=2)
                                         ggccccagTCGCTCCCCCCGCGGGGGCGTGGATCGAAACA..ACT.GACC...C.GGCA.ACGGCCTATGCGC.TTGA.CGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCg                                                            <  1:3268649/94‑1 (MQ=2)
                                          tcAGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACA.GCCT.GCGCG.GC.CGCA.TCGG.CT.TCCGC.CTG..CGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGa                                                           >  2:1721406/1‑101 (MQ=255)
                                          tcAGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACA.GCCT.GCGCG.GC.CGCA.TCGG.CT.TCCGC.CTG..CGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGa                                                           <  2:3938533/101‑1 (MQ=255)
                                          ggttgcGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATGGCTATGGCG.GC.AGTG.TCAG.GT.T.TTC.CAG..CGGC.GTCGCTCCCCCCGCGGGGGCGTGGATCg                                                            >  1:4526889/6‑101 (MQ=2)
                                           cgcccGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GCGT.GCGCG.CCATGAATTCAT.CG.AGCGC.CAT..CGAGCGTCGCTCCCCCCGCGGGGGCGGGGATc                                                             >  2:767479/5‑101 (MQ=0)
                                           acactgTCGCTCCCCCCGCGGGGGCGTGGATCGAAAC..GCTG.GCGCGCGC.CGCG.CCGG.ATCTTTGT.CTG..CGCG.GTCGCTCCCCCCGCGGGGGCGTGGATCGa                                                           >  1:493656/6‑101 (MQ=9)
                                           acactgTCGCTCCCCCCGCGGGGGCGTGGATCGAAAC..GCTG.GCGCGCGC.CGCG.CCGG.ATCTTTGT.CTG..CGCG.GTCGCTCCCCCCGCGGGGGCGTGGATCGa                                                           >  2:946426/6‑101 (MQ=9)
                                            cccggTCGCTCCCCCCGCGGGGGCGTGGATCGAAACA..ACT.GACC...C.GGCA.ACGG.CC.TATGCGCTGGACGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGaaa                                                         <  2:4361005/97‑1 (MQ=2)
                                            aGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACA.GCCT.GCGCG.GC.CGCA.TCGG.CT.TCCGC.CTG..CGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGaaa                                                         >  2:683994/1‑101 (MQ=255)
                                             gACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACA.GCCT.GCGCG.GC.CGCA.TCGG.CT.TCCGC.CTG..CGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAAc                                                        <  1:411050/101‑1 (MQ=255)
                                             ccggTCGCTCCCCCCGCGGGGGCGTGGATCGAAACA..ACT.GACC...C.GGCA.ACGG.CC.TATGCGCTGGACGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAAc                                                        >  1:2873016/4‑101 (MQ=2)
                                              ctgTCGCTCCCCCCGCGGGGGCGTGGATCGAAAC..GCTG.GCGCGCGC.CGCG.CCGG.ATCTTTGT.CTG..CGCG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAAc                                                        >  1:1803395/3‑101 (MQ=9)
                                               cGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GCGT.GCGCG.CCATGAATTCAT.CG.AGCGC.CAT..CGAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGaaa                                                         <  1:5316910/101‑1 (MQ=2)
                                                gTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GCGT.GCGCG.CCATGAATTCAT.CG.AGCGC.CAT..CGAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAAc                                                        >  1:3565486/1‑101 (MQ=2)
                                                gTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GCGT.GCGCG.CCATGAATTCAT.CG.AGCGC.CAT..CGAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAAc                                                        >  2:3618041/1‑101 (MQ=2)
                                                gTCGCTCCCCCCGCGGGGGCGTGGATCGAAAC..GCTG.GCGCGCGC.CGCG.CCGG.ATCTTTGT.CTG..CGCG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAActc                                                      <  2:2254809/101‑3 (MQ=9)
                                                 tCGCTCCCCCCGCGGGGGCGTGGATCGAAACATGGCTATGGCG.GC.AGTG.TCAG.GT.T.TTC.CAG..CGGC.GTCGCTCCCCCCGCGGGGGCGGGGATCGAAAcccg                                                     >  2:4208788/1‑98 (MQ=2)
                                                 tCGCTCCCCCCGCGGGGGCGTGGATCGAAACA.GCCT.GCGCG.GC.CGCA.TCGG.CT.TCCGC.CTG..CGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCa                                                    <  2:2763548/101‑1 (MQ=255)
                                                 tCGCTCCCCCCGCGGGGGCGTGGATCGAAACA.GCCT.GCGCG.GC.CGCA.TCGG.CT.TCCGC.CTG..CGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCa                                                    >  1:4289973/1‑101 (MQ=255)
                                                 tCGCCCCCCCCGCGGGGGCGTGGATCGAAACA.GCCT.GCGCG.GC.CGCA.TCGG.CT.TCCGC.CTG..CGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCa                                                    <  2:4279339/101‑1 (MQ=255)
                                                  cGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GCGT.GCGCG.CCATGAATTCAT.CG.AGCGC.CAT..CGAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACac                                                      <  2:632506/101‑2 (MQ=2)
                                                  cGCTCCCCCCGCGGGGGCGTGGATCGAAACA.GCCT.GCGCG.GC.CGCA.TCGG.CT.TCCGC.CTG..CGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCAg                                                   <  2:4639410/101‑1 (MQ=255)
                                                  cGCTCCCCCCGCGGGGGCGTGGATCGAAAC..GCTG.GCGCGCGC.CGCG.CCGG.ATCTTTGT.CTG..CGCG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAActcgg                                                    >  2:724523/1‑97 (MQ=9)
                                                  cGCTCCCCCCGCGGGGGCGTGGATCGAAAC..GCTG.GCGCGCGC.CGCG.CCGG.ATCTTTGT.CTG..CGCG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAActcgg                                                    >  1:4029964/1‑97 (MQ=9)
                                                  cGCTCCCCCCGCGGGGGCGTGGATCGAAAC..GCTG.GCGCGCGC.CGCG.CCGG.ATCTTTGT.CTG..CGCG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAActcgg                                                    >  2:431342/1‑97 (MQ=9)
                                                   gCTCCCCCCGCGGGGGCGTGGATCGAAACA.GCCT.GCGCG.GC.CGCA.TCGG.CT.TCCGC.CTG..CGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCAGc                                                  <  2:239843/101‑1 (MQ=255)
                                                   gCTCCCCCCGCGGGGGCGTGGATCGAAACA..ACT.GACC...C.GGCA.ACGG.CC.TATGCGCTGGACGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAActacacc                                                  >  2:5157848/1‑95 (MQ=2)
                                                    cTCCCCCCGCGGGGGCGTGGATCGAAACT.GCGT.GCGCG.CCATGAATTCAT.CG.AGCGC.CAT..CGAGCGTCGCTCCCCCCGCGGGGGCGGGGATCGAAACACCa                                                    >  1:5382939/1‑101 (MQ=2)
                                                    cTCCCCCCGCGGGGGCGTGGATCGAAACATGGCTATGGCG.GC.AGTG.TCAG.GT.T.TTC.CAG..CGGC.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAAcccgcca                                                  >  1:2210409/1‑95 (MQ=2)
                                                    cTCCCCCCGCGGGGGCGTGGATCGAAACA.GCCT.GCGCG.GC.CGCA.TCGG.CT.TCCGC.CTG..CGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCAGCt                                                 <  1:5760898/101‑1 (MQ=255)
                                                    cTCCCCCCGCGGGGGCGTGGATCGAAACA..ACT.GACC...C.GGCA.ACGG.CC.TATGCGCTGGACGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAActacaccg                                                 <  2:2829983/101‑8 (MQ=2)
                                                     tCCCCCCGCGGGGGCGTGGATCAAAACATGGCTATGGCG.GC.AGTG.TCAG.GT.T.TTC.CAG..CGGC.GTCCCCCCCCCCGCGGGGGCGTGGATCAAAAcccgccac                                                 <  2:5586089/101‑8 (MQ=2)
                                                      ccccccGCGGGGGCGTGGATCGAAACT.GCGT.GCGCG.CCATGAATTCAT.CG.AGCGC.CAT..CGAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACACCaat                                                  <  1:22563/101‑3 (MQ=2)
                                                      ccccccGCGGGGGCGTGGATCGAAACATGGCTATGGCG.GC.AGTG.TCAG.GT.T.TTC.CAG..CGGG.GTCGCCCCCCCCGCGGGGGCGTGGATCGAAAcccgccacg                                                <  2:490523/101‑9 (MQ=2)
                                                      ccccccGCGGGGGCGTGGATCGAAACATGGCTATGGCG.GC.AGTG.TCAG.GT.T.TTC.CAG..CGGC.GTCGCTCCCCCCGCGGGGGCGGGGATCGAAAcccgccacg                                                >  1:1232911/1‑93 (MQ=0)
                                                      ccccccGCGGGGGCGTGGATCGAAAC..GCTG.GCGCGCGC.CGCG.CCGG.ATCTTTGT.CTG..CGCG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAActcggcctc                                                <  1:4154103/101‑9 (MQ=9)
                                                      ccccccGCGGGGGCGTGGATCGAAAC..GCTG.GCGCGCGC.CGCG.CCGG.ATCTTTGT.CTG..CGCG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAActcggcctc                                                <  1:946643/101‑9 (MQ=9)
                                                       cccccGCGGGGGCGTGGATCGAAACA.GCCT.GCGCG.GC.CGCA.TCGG.CT.TCCGC.CTG..CGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCAGCtcgt                                              <  1:1579473/101‑1 (MQ=255)
                                                       cccccGCGGGGGCGTGGATCGAAACA.GCCT.GCGCG.GC.CGCA.TCGG.CT.TCCGC.CTG..CGTG.GTCACTCCCCCCGCGGGGGCGTGGATCGAAACATCAGCtcgt                                              <  1:1988384/101‑1 (MQ=255)
                                                       cccccGCGGGGGCGTGGATCGAAACA..ACT.GACC...C.GGCA.ACGG.CC.TATGCGCTGGACGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAActacaccgggg                                              <  2:1110004/101‑11 (MQ=2)
                                                       cccccGCGGGGGCGTGGATCGAAAC..GCTG.GCGCGCGC.CGCG.CCGG.ATCTTTGT.CTG..CGCG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAActcggcctcg                                               <  1:212554/101‑10 (MQ=9)
                                                             cgGGGGCGTGGATCGAAACA.GCCT.GCGCG.GC.CGCA.TCGG.CT.TCCGC.CTG..CGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCAGCTCGTCGTTgg                                        <  2:1442549/101‑1 (MQ=255)
                                                             cgGGGGCGTGGATCGAAACA.GCCT.GCGCG.GC.CGCA.TCGG.CT.TCCGC.CTG..CGTG.CTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCAGCTCGTCGTTgg                                        <  2:5010209/101‑1 (MQ=255)
                                                              gggggCGTGGATCGAAACA.GCCT.GCGCG.GC.CGCA.TCGG.CT.TCCGC.CTG..CGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCAGCTCGTCGTTGGt                                       >  2:5065095/1‑101 (MQ=255)
                                                               ggggCGTGGATCGAAACA.GCCT.GCGCG.GC.CGCA.TCGG.CT.TCCGC.CTG..CGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCAGCTCGTCGTTGGTc                                      <  1:1753461/101‑1 (MQ=255)
                                                               ggggCGTGGATCGAAACA.GCCT.GCGCG.GC.CGCA.TCGG.CT.TCCGC.CTG..CGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCAGCTCGTCGTTGGTc                                      <  1:470214/101‑1 (MQ=255)
                                                                   cGTGGATCGAAACA.GCCT.GCGCG.GC.CGCA.TCGG.CT.TCCGC.CTG..CGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCAGCTCGTCGTTGGTCAGcc                                  >  2:1250221/1‑101 (MQ=255)
                                                                    gTGGATCGAAACA.GCCT.GCGCG.GC.CGCA.TCGG.CT.TCCGC.CTG..CGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCAGCTCGTCGTTGGTCAGCCg                                 >  1:1678737/1‑101 (MQ=255)
                                                                     tGGATCGAAACA.GCCT.GCGCG.GC.CGCA.TCGG.CT.TCCGC.CTG..CGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCAGCTCGTCGTTGGTCAGCCgg                                >  2:3202754/1‑101 (MQ=255)
                                                                     tGGATCGAAACA.GCCT.GCGCG.GC.CGCA.TCGG.CT.TCCGC.CTG..CGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCAGCTCGTCGTTGGTCAGCCgg                                >  1:184949/1‑101 (MQ=255)
                                                                      ggATCGAAACA.GCCT.GCGCG.GC.CGCA.TCGG.CT.TCCGC.CTG..CGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCAGCTCGTCGTTGGTCAGCCGGa                               <  2:1862340/101‑1 (MQ=255)
                                                                       gATCGAAGCA.GCCT.GCGCG.GC.CGCA.TCGG.CT.TCCGC.CTG..CGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCAGCTCGTCGTTGGTCAGCCGGAc                              <  1:4818647/101‑1 (MQ=255)
                                                                       gATCGAAACA.GCCT.GCGCG.GC.CGCA.TCGG.CT.TCCGC.CTG..CGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCAGCTCGTCGTTGGTCAGCCGGAc                              <  1:2316727/101‑1 (MQ=255)
                                                                       gATCGAAACA.GCCT.GCGCG.GC.CGCA.TCGG.CT.TCCGC.CTG..CGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCAGCTCGTCGTTGGTCAGCCGGAc                              <  2:5357170/101‑1 (MQ=255)
                                                                           gAAACA.GCCT.GCGCG.GC.CGCA.TCGG.CT.TCCGC.CTG..CGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCAGCTCGTCGTTGGTCAGCCGGACGATg                          >  2:456703/1‑101 (MQ=255)
                                                                           gAAACA.GCCT.GCGCG.GC.CGCA.TCGG.CT.TCCGA.CTG..CGTG.GTTGCCTCCCCCGCGGGGGCGTGAATCGAAACATCAGCTCGTCGTTTGGCAGCCGGACGATg                          >  2:3817537/1‑101 (MQ=255)
                                                                            aaaCA.GCCT.GCGCG.GC.CGCA.TCGG.CT.TCCGC.CTG..CGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCAGCTCGTCGTTGGTCAGCCGGACGATGt                         >  1:1460664/1‑101 (MQ=255)
                                                                             aaCA.GCCT.GCGCG.GC.CGCA.TCGG.CT.TCCGC.CTG..CGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCAGCTCGTCGTTGGTCAGCCGGACGATGTc                        <  2:4751067/101‑1 (MQ=255)
                                                                                  gCCT.GCGCG.GC.CGCA.TCGG.CT.TCCGC.CTG..CGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCAGCTCGTCGTTGGTCAGCCGGACGATGTCCGTc                    <  1:1370919/101‑1 (MQ=255)
                                                                                      cgcgcG.GC.CGCA.TCGG.CT.TCCGC.CTG..CGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCAGCTCGTCGTTGGTCAGCCGGACGATGTCCGTCGCt                 <  1:4244751/100‑1 (MQ=255)
                                                                                         gcg.GC.CGCA.TCGG.CT.TCCGC.CTG..CGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCAGCTCGTCGTTGGTCAGCCGGACGATGTCCGTCGCTccc              >  2:2329306/1‑101 (MQ=255)
                                                                                           g.gC.CGCA.TCGG.CT.TCCGC.CTG..CGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCAGCTCGTCGTTGGTCAGCCGGACGATGTCCGGCGCTccccc            >  2:2560670/1‑101 (MQ=255)
                                                                                             gC.CGCA.TCGG.CT.TCCGC.CTG..CGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCAGCTCGTCGTTGGTCAGCCGGACGATGTCCGTCGCTcccccc           >  1:1048022/1‑101 (MQ=255)
                                                                                                cGCA.TCGG.CT.TCCGC.CTG..CGTG.GTCGCCCCCCCCGCGGGGGCGTGGATCGAAACATCAGCTCGTCGTTGGTCAGCCGGACGATGTCCGTCGCTCCCCCcgc         <  1:674880/101‑1 (MQ=255)
                                                                                                   a.TCGG.CT.TCCGC.CTG..CGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCAGCTCGTCGTTGGTCAGCCGGACGATGTCCGTCGCTCCCCCCGCggg      >  1:5568686/1‑101 (MQ=255)
                                                                                                     tCGG.CT.TCCGC.CTG..CGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCAGCTCGTCGTTGGTCAGCCGGACGATGTCCGTCGCTCCCCCCGCgggg     >  1:1014489/1‑101 (MQ=255)
                                                                                                       gg.CT.TCCGC.CTG..CGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCAGCTCGTCGTTGGTCAGCCGGACGATGTCCGTCGCTCCCCCCGCGGGGGc   >  1:1308751/1‑101 (MQ=255)
                                                                                                       gg.CT.TCCGC.CTG..CGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCAAAACATCAGCTCGTCGTTGGTCAGCCGGACGATGTCCGTCGCTCCCCCCGCGGGGGc   <  2:4088154/101‑1 (MQ=255)
                                                                                                        g.CT.TCCGC.CTG..CGTG.GTCGCACCCCCCGCGGGGGCGTGGATCGAAACATCAGCTCGTCGTTGGTCAGCCGGACGATGTCCGTCGCTCCCCCCGCGGGGGCg  <  2:5687584/101‑1 (MQ=255)
                                                                                                          |                                                                                                        
CGTGGATCGAAACCGGTCCAGAGATAGACCGCGCGGACCTCCTCAGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACA.GCCT.GCGCG.GC.CGCA.TCGG.CT.TCCGC.CTG..CGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCAGCTCGTCGTTGGTCAGCCGGACGATGTCCGTCGCTCCCCCCGCGGGGGCG  >  NC_014034/1396855‑1397054

Base quality scores:  ATCG  < 3 ≤  ATCG  < 18 ≤  ATCG  < 25 ≤  ATCG  < 33 ≤  ATCG  < 40 ≤  ATCG