Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_014034 1,396,961 C→T 12.1% R281R (CGC→TGT RCAP_rcc01307 → hypothetical protein

Read alignment evidence...
  seq id position change freq score reads annotation genes product
*NC_0140341,396,9610C→T12.1% 115.1 ‑0.8 ‑0.1380R281R (CGC→CGTRCAP_rcc01307hypothetical protein
Reads supporting (aligned to +/- strand):  new base (26/20):  ref base (166/168):  total (192/188)
Fisher's exact test strand distribution p-value = 4.33e-01
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 8.02e-01

Only 100 of 405 total aligned reads displayed.

TCGAAACCGGTCCAGAGATAGACCGCGCGGACCTCCTCAGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACA.GCCT.GCGCG.GC.CGCA.TCGGCT.TCCGC.CTG..CGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCAGCTCGTCGTTGGTCAGCCGGACGATGTCCGTCGCTCCCCCCGCGGGGGCGTGGATC  >  NC_014034/1396861‑1397060
                                                                                                          |                                                                                                       
tCGAAACCGGTCCAGAGATAGACCGCGCGGACCTCCTCAGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACA.GCCT.GCGCG.GC.CGCA.TCGGCT.TCCGc                                                                                                         <  1:134819/101‑1 (MQ=255)
 cGAAACCGGTCCAGAGATAGACCGCGCGGACCTCCTCAGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACA.GCCT.GCGCG.GC.CGCA.TCGGCT.TCCGc.c                                                                                                       <  2:2115465/101‑1 (MQ=255)
  gAAACCGGTCCAGAGATAGACCGCGCGGACCTCCTCAGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACA.GCCT.GCGCG.GC.CGCA.TCGGCT.TCCGC.Ct                                                                                                      >  2:2447377/1‑101 (MQ=255)
   aaaCCGGTCCAGAGATAGACCGCGCGGACCTCCTCAGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACA.GCCT.GCGCG.GC.CGCA.TCGGCT.TCCGC.CTg                                                                                                     >  2:1263224/1‑101 (MQ=255)
      ccGGTCCAGAGATAGACCGCGCGGACCTCCTCAGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACA.GCCT.GCGCG.GC.CGCA.TCGGCT.TCCGC.CTG..CGt                                                                                                >  1:712688/1‑101 (MQ=255)
      ccGGTCCAGAGATAGACCGCGCGGACCTCCTCAGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACA.GCCT.GCGCG.GC.CGCA.TCGGCT.TCCGC.CTG..CGt                                                                                                >  1:3947848/1‑101 (MQ=255)
      ccGGTCCAGAGATAGACCGCGCGGACCTCCTCAGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACA.GCCT.GCGCG.GC.CGCA.TCGGCT.TCCGC.CTG..CGt                                                                                                >  1:281182/1‑101 (MQ=255)
       cGGTCCAGAGATAGACCGCGCGGACCTCCTCAGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACA.GCCT.GCGCG.GC.CGCA.TCGGCT.TCCGC.CTG..CGTg                                                                                               >  2:2828122/1‑101 (MQ=255)
           ccAGAGATAGCCCGCGCGGACCTCCTCAGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACA.GCCT.GCGCG.GC.CGCA.TCGGCT.TCCGC.CTG..CGTG.GTCg                                                                                          <  2:4586546/101‑1 (MQ=255)
            cAGAGATAGACCGCGCGGACCTCCTCAGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACA.GCCT.GCGCG.GC.CGCA.TCGGCT.TCCGC.CTG..CGTG.GTCGc                                                                                         >  2:3206088/1‑101 (MQ=255)
             agagATGGGCCGCGCGGCCCTCCTCAGACGTCTGGCCCCCCGCGGGGGCGTGGATCGAAACA.GCCT.GCGCG.GC.CGCA.TCGGCT.TCCGC.CTG..CGTG.GTCGCt                                                                                        <  2:1701737/101‑1 (MQ=38)
             agagATAGACCGCGCGGACCTCCTCAGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACA.GCCT.GCGCG.GC.CGCA.TCGGCT.TCCGC.CTG..CGTG.GTCGCt                                                                                        <  1:857447/101‑1 (MQ=255)
             agagATAGACCGCGCGGACCTCCTCAGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACA.GCCT.GCGCG.GC.CGCA.TCGGCT.TCCGC.CTG..CGTG.GTCGCt                                                                                        <  1:117957/101‑1 (MQ=255)
              gagaTAGACCGCTCTGACCTCCTCAGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACA.GCCT.GCGCT.GC.CGCA.TCGGCT.TCCGC.CTG..CGTG.GTCGCTc                                                                                       <  2:2006466/101‑1 (MQ=255)
                gaTAGACCGCGCGGACCTCCTCAGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACA.GCCT.GCGCG.GC.CGCA.TCGGCT.TCCGC.CTG..CGTG.GTCGCTccc                                                                                     >  1:1625200/1‑101 (MQ=255)
                   aGACCGCGCGGACCTCCTCAGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACA.GCCT.GCGCG.GC.CGCA.TCGGCT.TCCGC.CTG..CGTG.GTCGCTcccccc                                                                                  <  1:2737324/101‑1 (MQ=255)
                           cgGACCTCCTCAGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACA.GCCT.GCGCG.GC.CGCA.TCGGCT.TCCGC.CTG..CGTG.GTCGCTCCCCCCGCGGGGGc                                                                          >  1:3891049/1‑101 (MQ=255)
                             gACCTCCTCAGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACA.GCCT.GCGCG.GC.CGCA.TCGGCT.TCCGC.CTG..CGTG.GTCGCTCCCCCCGCGGGGGCGt                                                                        <  2:2270000/101‑1 (MQ=255)
                              aCCTCCTCAGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACA.GCCT.GCGCG.GC.CGCA.TCGGCT.TCCGC.CTG..CGTG.GTCGCTCCCCCCGCGGGGGCGTg                                                                       <  1:1263315/101‑1 (MQ=255)
                               cctcctCAGACGTCGCTCCCCCCGCGGGGTCGTTGATCGATACA.GCCT.GCGCT.GC.CGCA.TCGGCT.TCCGC.CTG..CGTT.GTCGCCCCCCCCGCGGGGGCGTgg                                                                      <  1:4712738/101‑1 (MQ=38)
                               cctcctCAGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACA.GCCT.GCGCG.GC.CGCA.TCGGCT.TCCGC.CTG..CGTG.GTCGCTCCCCCCGCGGGGGCgcat                                                                      <  2:4388010/101‑4 (MQ=255)
                                cgtggatccggTCGCTCCCCCCGCGGGGGCGTGGATCGAAACA..ACT.GACC...C.GGCA.ACGGCC.TATGCGCTGGACGTG.GTCGCTCCCCCCGCGGGGGCGTGGa                                                                     >  2:3892539/11‑101 (MQ=0)
                                attggccccggTCGCTCCCCCCGCGGGGGCGTGGATCGAAACA..ACT.GACC...C.GGCA.ACGGCC.TATGCGCTGGACGTG.GTCGCTCCCCCCGCGGGGGCGTGGa                                                                     >  2:1594426/11‑101 (MQ=0)
                                 ttggccccggTCGCTCCCCCCGCGGGGGCGTGGATCGAAACA..ACT.GACC...C.GGCA.ACGGCC.TATGCGCTGGACGTG.GTCGCTCCCCCCGCGGGGGCGTGGAt                                                                    <  1:4183820/92‑1 (MQ=2)
                                 tcctcAGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACA.GCCT.GCGCG.GC.CGCA.TCGGCT.TCCGC.CTG..CGTG.GTCGCTCCCCCCGCGGGGGCGTgggt                                                                    >  2:4038780/1‑99 (MQ=255)
                                 cgcggttgcGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATGGCTATGGCG.GC.AGTG.TCAGGT.T.TTC.CAG..CGGC.GTCGCTCCCCCCGCGGGGGCGTGGa                                                                     <  2:624531/93‑1 (MQ=2)
                                  gcggttgcGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATGGCTATGGCG.GC.AGTG.TCAGGT.T.TTC.CAG..CGGC.GTCGCTCCCCCCGCGGGGGCGTGGAt                                                                    <  2:1403685/94‑1 (MQ=2)
                                  ccccgcccGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GCGT.GCGCG.CCATGAATTCATCG.AGCGC.CAT..CGAGCGTCGCTCCCCCCGCGGGGGCGTgg                                                                      <  1:1730043/94‑1 (MQ=0)
                                  ccccAGACCTCGCTCCCCCGGGGGGGGCGTGGATGGAAACG.GCCT.GCCCG.GC.CGCA.GCGGCT.TCCGC.CTG..CGTG.GTCGCTCCCCCCGCGGGGGCGTGGATc                                                                   <  1:2945353/98‑1 (MQ=25)
                                   ggccccggTCGCTCCCCCCGCGGGGGCGTGGATCGAAACA..ACT.GACC...C.GGCA.ACGGCC.TATGCGCTGGACGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCg                                                                  <  2:479052/94‑1 (MQ=2)
                                   ctcaagcggcGCGCCCCCCGCGGGGGCGTGGATCGAAACA.GCCT.GCGCG.GC.CGCA.TCGGCT.TCCGC.CTG..CGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCg                                                                  <  2:266986/92‑1 (MQ=37)
                                   caacactgTCGCTCCCCCCGCGGGGGCGTGGATCGAAAC..GCTG.GCGCGCGC.CGCG.CCGGATCTTTGT.CTG..CGCG.GTCGCTCCCCCCGCGGGGGCGTGGATc                                                                   <  2:1669496/94‑1 (MQ=9)
                                    ggttgcGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATGGCTATGGCG.GC.AGTG.TCAGGT.T.TTC.CAG..CGGC.GTCGCTCCCCCCGCGGGGGCGTGGATCg                                                                  <  2:3985152/96‑1 (MQ=2)
                                    gccccggTCGCTCCCCCCGCGGGGGCGTGGATCGAAACA..ACT.GACC...C.GGCA.ACGGCC.TATGCGCTGGACGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGa                                                                 >  2:3398284/7‑101 (MQ=2)
                                    gccccggTCGCTCCCCCCGCGGGGGCGTGGATCGAAACA..ACT.GACC...C.GGCA.ACGGCC.TATGCGCTGGACGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGa                                                                 >  2:2799351/7‑101 (MQ=2)
                                    gccccggTCGCTCCCCCCGCGGGGGCGTGGATCGAAACA..ACT.GACC...C.GGCA.ACGGCC.TATGCGCTGGACGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGa                                                                 >  2:3681033/7‑101 (MQ=2)
                                    gccccggTCGCCCCCCCCGCGGGGGCGTGGATCGAAACA..ACT.GACC...C.GGCA.ACGGCC.TATGCGCTGGACGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGa                                                                 <  1:3354320/95‑1 (MQ=2)
                                     cgcccGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GCGT.GCGCG.CCATGAATTCATCG.AGCGC.CAT..CGAGCGTCGCTCCCCCCGCGGGGGCGTGGATc                                                                   <  2:2391373/97‑1 (MQ=2)
                                     cAGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACA.GCCT.GCGCG.GC.CGCA.TCGGCT.TCCGC.CTG..CGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGaa                                                                <  2:102252/101‑1 (MQ=255)
                                     cAGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACA.GCCT.GCGCG.GC.CGCA.TCGGCT.TCCGC.CTG..CGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGaa                                                                >  2:3001107/1‑101 (MQ=255)
                                     cAGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACA.GCCT.GCGCG.GC.CGCA.TCGGCT.TCCGC.CTG..CGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGaa                                                                >  2:1137759/1‑101 (MQ=255)
                                     cAGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACA.GCCT.GCGCG.GC.CGCA.TCGGCT.TCCGC.CTG..CGTG.GTCGCTCCCCCCGCGGGGGCGGGGATCGaa                                                                >  1:4031864/1‑101 (MQ=255)
                                     cAGACGTCGCTCCCCCCGCGGGGGCGGGGATCGAAACA.GCCT.GCGCG.GC.CGCA.TCGGCT.TCCGC.CTG..CGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGaa                                                                <  2:1334934/101‑1 (MQ=255)
                                      gcccGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GCGT.GCGCG.CCATGAATTCATCG.AGCGC.CAT..CGAGCGTCGCTCCCCCCGCGGGGGCGTGGATCg                                                                  >  1:158061/4‑101 (MQ=2)
                                       tgcGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATGGCTATGGCG.GC.AGTG.TCAGGT.T.TTC.CAG..CGGC.GTCGCTCCCCCCGCGGGGGCGTGGATCGaaa                                                               <  2:148623/99‑1 (MQ=2)
                                       tgcGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATGGCTATGGCG.GC.AGTG.TCAGGT.T.TTC.CAG..CGGC.GTCGCTCCCCCCGCGGGGGCGTGGATCGaaa                                                               <  1:3168204/99‑1 (MQ=2)
                                       gACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACA.GCCT.GCGCG.GC.CGCA.TCGGCT.TCCGC.CTG..CGTG.GTCGCTCCCCCCGCGGGGGCGTGGGTCGAAAc                                                              >  1:4618125/1‑101 (MQ=255)
                                       ccggTCGCTCCCCCCGCGGGGGCGTGGATCGAAACA..ACT.GACC...C.GGCA.ACGGCC.TATGCGCTGGACGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAAc                                                              <  1:4547005/98‑1 (MQ=2)
                                       ccggTCGCTCCCCCCGCGGGGGCGTGGATCGAAACA..ACT.GACC...C.GGCA.ACGGCC.TATGCGCTGGACGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAAc                                                              <  2:539112/98‑1 (MQ=2)
                                        ctgTCGCTCCCCCCGCGGGGGCGTGGATCGAAAC..GCTG.GCGCGCGC.CGCG.CCGGATCTTTGT.CTG..CGCG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAAc                                                              >  1:3794378/3‑101 (MQ=9)
                                        cggTCGCTCCCCCCGCGGGGGCGTGGATCGAAACA..ACT.GACC...C.GGCA.ACGGCC.TATGCGCTGGACGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAAct                                                             <  1:783361/99‑2 (MQ=2)
                                         tgTCGCTCCCCCCGCGGGGGCGTGGATCGAAAC..GCTG.GCGCGCGC.CGCG.CCGGATCTTTGT.CTG..CGCG.GTCGCTCCCCCCGCGGGGGCGTGGATCGaaaga                                                             >  1:2091579/2‑99 (MQ=9)
                                         cGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GCGT.GCGCG.CCATGAATTCATCG.AGCGC.CAT..CGAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGaaa                                                               <  2:374910/101‑1 (MQ=2)
                                         cGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACA.GCCT.GCGCG.GC.CGCA.TCGGCT.TCCGC.CTG..CGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAt                                                            >  2:1791219/1‑101 (MQ=255)
                                          gTCGCTCCCCCCGCGGGGGCGTGGATCGAAAC..GCTG.GCGCGCGC.CGCG.CCGGATCTTTGT.CTG..CGCG.GTCGCTCCCCCCGCGGGGGGGGGGGGCGAAActc                                                            >  1:4097007/1‑99 (MQ=9)
                                          gTCGCGCCCCCCGCGGGGGCGTGGATCGAAACA..ACT.GACC...C.GGCA.ACGGCC.TATGCGCTGGACGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAActac                                                           <  1:444856/101‑4 (MQ=2)
                                          gTCGCCCCCCCCGCGGGGGCGTGGATCGAAACATGGCTATGGCG.GC.AGTG.TCAGGT.T.TTC.CAG..CGGC.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAAccc                                                            <  2:4067860/101‑3 (MQ=2)
                                           tCGCTCCCCCCGCGGGGGCGTGGATCGAAACA.GCCT.GCGCG.GC.CGCA.TCGGCT.TCCGC.CTG..CGTG.GTCGCTCCCCCCGCGGGGGCGGGGATCGAAACATCa                                                          >  1:614202/1‑101 (MQ=255)
                                           tCGCTCCCCCCGCGGGGGCGTGGATCGAAACA..ACT.GACC...C.GGCA.ACGGCC.TATGCGCTGGACGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAActaca                                                          <  2:2572451/101‑5 (MQ=2)
                                            cGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GCGT.GCGCG.CCATGAATTCATCG.AGCGC.CAT..CGAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACac                                                            >  1:3064434/1‑100 (MQ=2)
                                            cGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GCGT.GCGCG.CCATGAATTCATCG.AGCGC.CAT..CGAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACac                                                            >  1:1931556/1‑100 (MQ=2)
                                             gCTCCCCCCGCGGGGGCGTGGATCGAAACA..ACT.GACC...C.GGCA.ACGGCC.TATGCGCTGGACGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAActacacc                                                        >  2:360535/1‑95 (MQ=2)
                                              cTCCCCCCGCGGGGGCGTGGATCGAAACA..ACT.GACC...C.GGCA.ACGGCC.TATGCGCTGGACGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAActacaccg                                                       <  1:3614395/101‑8 (MQ=2)
                                              cTCCCCCCGCGGGGGCGTAGATCGAAACA..ACT.GACC...C.GGCA.ACGGCC.TATGCGCTGGACGTG.GTCGCTCCCACCGCGGGGGCGTGGATCGAAActacaccg                                                       >  2:802017/1‑94 (MQ=2)
                                               tCCCCCCGCGGGGGCGTGGATCGAAACT.GCGT.GCGCG.CCATGAATTCATCG.AGCGC.CAT..CGAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACACCaa                                                         >  1:79411/1‑100 (MQ=2)
                                               tCCCCCCGCGGGGGCGTGGATCGAAACT.GCGT.GCGCG.CCATGAATTCATCG.AGCGC.CAT..CGAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACACCaa                                                         >  2:3496730/1‑100 (MQ=2)
                                               tCCCCCCGCGGGGGCGTGGATCGAAACA.GCCT.GCGCG.GC.CGCA.TCGGCT.TCCGC.CTG..CGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCAGCtc                                                      >  2:1741829/1‑101 (MQ=255)
                                               tCCCCCCGCGGGGGCGTGGATCGAAAC..GCTG.GCGCGCGC.CGCC.CCGGATCTTTGT.CTG..CGCG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAActcggcct                                                       >  1:1468121/1‑94 (MQ=9)
                                               cccccccGCGGGGGCGTGGATCGAAAC..GCTG.GCGCGCGC.CGCG.CCGGATCTTTGT.CTG..CGCG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAActcggcct                                                       <  1:1000110/100‑8 (MQ=9)
                                                ccccccGCGGGGGCGTGGATCGAAACT.GCGT.GCGCG.CCATGAATTCATCG.AGCGC.CAT..CGAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACACCaat                                                        <  1:2406552/101‑3 (MQ=2)
                                                ccccccGCGGGGGCGTGGATCGAAACA.GCCT.GCGCG.GC.CGCA.TCGGCT.TCCGC.CTG..CGTG.GTCGCTCCCCCCGCGGGGGCGGGGATCGAAACATCAGCtcg                                                     >  2:2165642/1‑101 (MQ=255)
                                                ccccccGCGGGGGCGTGGATCGAAACA..ACT.GACC...C.GGCA.ACGGCC.TATGCGCTGGACGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAActacaccggg                                                     >  1:2689073/1‑92 (MQ=2)
                                                 cccccGCGGGGGCGTGGATCGAAACT.GCGT.GCGCG.CCATGAATTCATCG.AGCGC.CAT..CGAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACACCaatg                                                       <  2:3274638/101‑4 (MQ=2)
                                                 cccccGCGGGGGCGTGGATCGAAACA.GCCT.GCGCG.GC.CGCA.TCGGCT.TCCGC.CTG..CGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCAGCtcgt                                                    >  1:2283097/1‑101 (MQ=255)
                                                 cccccGCGGGGGCGTGGATCGAAACA..ACT.GACC...C.GGCA.ACGGCC.TATGCGCTGGACGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAActacaccgggg                                                    <  1:2188209/101‑11 (MQ=2)
                                                 cccccGCGGGGGCGTGGATCGAAAC..GCTG.GCGCGCGC.CGCG.CCGGATCTTTGT.CTG..CGCG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAActcggcctcg                                                     <  1:2505007/101‑10 (MQ=9)
                                                  ccccGCGGGGGCGTGGATCGAAAC..GCTG.GCGCGCGC.CGCG.CCGGATCTTTGT.CTG..CGCG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAActcggcctcga                                                    >  2:684213/1‑91 (MQ=2)
                                                  ccccGCGGGGGCGTGGATCGAAAC..GCTG.GCGCGCGC.CGCG.CCGGATCTTTGT.CTG..CGCG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAActcggcctcga                                                    >  1:2357086/1‑91 (MQ=2)
                                                    ccGCGGGGGCGTGGATCGAAACT.GCGT.GCGCG.CCATGAATTCATCG.AGCGC.CAT..CGAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACACCaatgccg                                                    >  1:1094457/1‑95 (MQ=0)
                                                      gcgGGGGCGTGGATCGAAACA.GCCT.GCGCG.GC.CGCA.TCGGCT.TCCGC.CTG..CGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCAGCTCGGCGTTg                                               >  2:4195614/1‑101 (MQ=255)
                                                             cGTGGATCGAAACA.GCCT.GCGCG.GC.CGCA.TCGGCT.TCCGC.CTG..CGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCAGCTCGTCGTTGGTCAGcc                                        <  2:1893078/101‑1 (MQ=255)
                                                              gTGGATCGAAACA.GCCT.GCGCG.GC.CGCA.TCGGCT.TCCGC.CTG..CGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCAGCTCGTCGTTGGTCAGCCg                                       <  1:3967017/101‑1 (MQ=255)
                                                              gTGGATCGAAACA.GCCT.GCGCG.GC.CGCA.TCGGCT.TCCGC.CTG..CGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCAGCTCGTCGTTGGTCAGCCg                                       >  1:262244/1‑101 (MQ=255)
                                                                ggATCGAAACA.GCCT.GCGCG.GC.CGCA.TCGGCT.TCCGC.CTG..CGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCAGCTCGTCGTTGGTCAGCCGGa                                     >  2:332009/1‑101 (MQ=255)
                                                                     gAAACA.GCCT.GCGCG.GC.CGCA.TCGGCT.TCCGC.CTG..CGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCAGCTCGTCGTTGGTCAGCCGGACGATg                                >  1:1396218/1‑101 (MQ=255)
                                                                      aaaCA.GCCT.GCGCG.GC.CGCA.TCGGCT.TCCGC.CTG..CGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCAGCTCGTCGTTGGTCAGCCGGACGATGt                               <  2:2717255/101‑1 (MQ=255)
                                                                       aaCA.GCCT.GCGCG.GC.CGCA.TCGGCT.TCCGC.CTG..CGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCAGCTCGTCGTTGGTCAGCCGGACGATGTc                              <  1:2594532/101‑1 (MQ=255)
                                                                          cacccT.GCGCG.GC.CGCA.TCGGCT.TCCGC.CTG..CGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCAGCTCGTCGTTGGTCAGCCGGACGATGTCCg                            <  2:773028/98‑1 (MQ=255)
                                                                           ggCCC.GCGCG.GC.AGCA.CCGTCT.CCCGC.TTG..CGTG.GTCACTCCCCCCGCTGTGTCGTGGATCGAAACATCAGCTCGTCGTTGGTCAGCCGGACGATGTCCGt                           <  1:898098/100‑1 (MQ=255)
                                                                              cT.GCGCG.GC.CGCA.TCGGCT.TCCGC.CTG..CGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCAGCTCGTCGTTGGTCAGCCGGACGATGTCCGTCGc                        <  2:3851075/101‑1 (MQ=255)
                                                                                 gcgcG.GC.CGCA.TCGGCT.TCCGC.CTG..CGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCAGCTCGTCGTTGGTCAGCCGGACGATGTCCGTCGCTc                      <  2:4337311/101‑1 (MQ=255)
                                                                                    cg.GC.CGCA.TCGGCT.TCCGC.CTG..CGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCAGCTCGTCGTTGGTCAGCCGGACGATGTCCGTCGCTcccc                   <  1:3238880/101‑1 (MQ=255)
                                                                                    cg.GC.CGCA.TCGGCT.TCCGC.CTG..CGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCAGCTCGTCGTTGGTCAGCCGGACGATGTCCGTCGCTcccc                   <  2:466608/101‑1 (MQ=255)
                                                                                     g.gC.CGCA.TCGGCT.TCCGC.CTG..CGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATAAGCTCGTCGTTGGTCAGCCGGACGATGTCCGTCGCTccccc                  >  1:3489041/1‑101 (MQ=255)
                                                                                       gC.CGCA.TCGGCT.TCCGC.CTG..CGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCAGCTCGTCGTTGGTCAGCCGGACGATGTCCGTCGCTcccccc                 <  2:866586/101‑1 (MQ=255)
                                                                                        c.cGCA.TCGGCT.TCCGC.CTG..CGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCAGCTCGTCGTTGGTCAGCCGGACGATGTCCGTCGCTCCCCCcg                <  1:2039043/101‑1 (MQ=255)
                                                                                        c.cGCA.TCGGCT.TCCGC.CTG..CGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCAGCTCGTCGTTGGTCAGCCGGACGATGTCCGTCGCTCCCCCcg                <  1:4363033/101‑1 (MQ=255)
                                                                                               tCGGCT.TCCGC.CTG..CGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCAGCTCGTCGTTGGTCAGCCGGACGATGTCCGTCGCTCCCCCCGCgggg           >  1:3422955/1‑101 (MQ=255)
                                                                                                        cGC.CTG..CGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCAGCTCGTCGTTGGTCAGCCGGACGATGTCCGTCGCTCCCCCCGCGGGGGCGTGGAt   <  1:3001295/101‑1 (MQ=255)
                                                                                                         gC.CTG..CGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCAGCTCGTCGTTGGTCAGCCGGACGATGTCCGTCGCTCCCCCCGCGGGGGCGTGGATc  <  2:3061606/101‑1 (MQ=255)
                                                                                                          |                                                                                                       
TCGAAACCGGTCCAGAGATAGACCGCGCGGACCTCCTCAGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACA.GCCT.GCGCG.GC.CGCA.TCGGCT.TCCGC.CTG..CGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCAGCTCGTCGTTGGTCAGCCGGACGATGTCCGTCGCTCCCCCCGCGGGGGCGTGGATC  >  NC_014034/1396861‑1397060

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