Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_014034 1,396,968 G→C 10.4% G284R (GGT→CGT)  RCAP_rcc01307 → hypothetical protein

Read alignment evidence...
  seq id position change freq score reads annotation genes product
*NC_0140341,396,9680G→C10.4% 108.0 ‑1.4 ‑0.0527G284R (GGT→CGT) RCAP_rcc01307hypothetical protein
Reads supporting (aligned to +/- strand):  new base (21/34):  ref base (219/250):  total (240/287)
Fisher's exact test strand distribution p-value = 2.54e-01
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 9.72e-01

Only 100 of 596 total aligned reads displayed.

CGGTCCAGAGATAGACCGCGCGGACCTCCTCAGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACA.GCCT.GCGCG.GC.CGCA.TCGGCT.TCCGC.CTG..CGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCAGCTCG.TCGTTGGTCAGCCGGAC.GATGTCCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGC  >  NC_014034/1396868‑1397068
                                                                                                             |                                                                                                       
cGGTACAGAGATAGACCGCGCGGACCTCCTCAGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACA.GCCT.GCGCG.GC.CGCA.TCGGCT.TCCGC.CTG..CGTg                                                                                                         >  2:1644415/1‑101 (MQ=255)
   tCCAGAGATAGACCGCGCGGACCTCCTCAGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACA.GCCT.GCGCG.GC.CGCA.TCGGCT.TCCGC.CTG..CGTG.GTc                                                                                                     <  2:2187321/101‑1 (MQ=255)
      cgagaTAGCCCGCGCGGCCCTCCTCAGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACA.GCCT.GCGCG.GC.CGCA.TCGGCT.TCCGC.CTG..CGTG.GTCGCt                                                                                                  <  2:3160027/100‑1 (MQ=255)
          aTAGACCGCGCGGACCTCCTCAGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACA.GCCT.GCGCG.GC.CGCA.TCGGCT.TCCGC.CTG..CGTG.GTCGCTcccc                                                                                              >  1:3089302/1‑101 (MQ=255)
           tagacccgccgGACCTCCTCAGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACA.GCCT.GCGCG.GC.CGCA.TCGGCT.TCCGC.CTG..CGTG.GTCGCTccccc                                                                                             <  2:3794856/92‑1 (MQ=255)
           tAGACCGCGCGGACCTCCTCAGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACA.GCCT.GCGCG.GC.CGCA.TCGGCT.TCCGC.CTG..CGTG.GTCGCTccccc                                                                                             <  2:2515112/101‑1 (MQ=255)
             gACCGCGCGGCCCTCCTCAGACGTCGCCCCCCCCGCGGGGGCGTGGATCGAAACA.GCCT.GCGCG.GC.CGCA.TCGGCT.TCCGC.CTG..CGTG.GTCGCTCCCCCcg                                                                                           <  2:5670868/101‑1 (MQ=255)
              aCCGCGCGGCCCGCCGCAGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACA.GCCT.GCGCG.GC.CGCA.TCGGCT.TCCGC.CTG..CGTG.GTCGCTCCCCCcgc                                                                                          <  2:6443827/101‑1 (MQ=255)
                    cgGACCTCCTCAGACGTCGCCCCCCCCGCGGGGGCGTGGATCGAAACA.GCCT.GCGCG.GC.CGCA.TCGGCT.TCCGC.CTG..CGTG.GTCGCTCCCCCCGCCGGGGc                                                                                    <  2:3514505/101‑1 (MQ=255)
                         ggtcaacactgTCGCTCCCCCCGCGGGGGCGTGGATCGAAAC..GCTG.GCGCGCGC.CGCG.CCGGATCTTTGT.CTG..CGCG.GTCGCTCCCCCCGCGGGGGCGTgg                                                                                >  1:5711524/11‑101 (MQ=9)
                         attggccccggTCGCTCCCCCCGCGGGGGCGTGGATCGAAACA..ACT.GACC...C.GGCA.ACGGCC.TATGCGCTGGACGTG.GTCGCTCCCCCCGCGGGGGCGTGGa                                                                               >  1:4989459/11‑101 (MQ=0)
                         attggccccggTCGCTCCCCCCGCGGGGGCGTGGATCGAAACA..ACT.GACC...C.GGCA.ACGGCC.TATGCGCTGGACGTG.GTCGCTCCCCCCGCGGGGGCGTGGa                                                                               <  2:2468786/91‑1 (MQ=0)
                          ttggccccggTCGCTCCCCCCGCGGGGGCGTGGATCGAAACA..ACT.GACC...C.GGCA.ACGGCC.TATGCGCTGGACGTG.GTCGCTCCCCCCGCGGGGGCGTGGAt                                                                              >  2:3867606/10‑101 (MQ=2)
                          cgcggttgcGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATGGCTATGGCG.GC.AGTG.TCAGGG.T.TTC.CAG..CGGC.GTCGCTCCCCCCGCGGGGGCGTGGa                                                                               >  1:1994423/9‑101 (MQ=0)
                           tcaacactgTCGCTCCCCCCGCGGGGGCGTGGATCGAAAC..GCTG.GCGCGCGC.CGCG.CCGGATCTTTGT.CTG..CGCG.GTCGCTCCCCCCGCGGGGGCGTGGAt                                                                              >  1:1279405/9‑101 (MQ=9)
                           cctcGGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACA.GCCT.GCGCG.GC.CGCA.TCGGCT.TCCGC.CTG..CGTG.GTCGCTCCCCCCGCGGGGGCGTGGATc                                                                             <  1:2742438/101‑1 (MQ=255)
                            ggccccggTCGCTCCCCCCGCGGGGGCGTGGATCGAAACA..ACT.GACC...C.GGCA.ACGGCC.TATGCGCTGGACGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCg                                                                            <  2:3216650/94‑1 (MQ=2)
                            ctcAGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACA.GCCT.GCGCG.GC.CGCA.TCGGCT.TCCGC.CTG..CGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCg                                                                            <  1:1904131/101‑1 (MQ=255)
                            ctcAGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACA.GCCT.GCGCG.GC.CGCA.TCGGCT.TCCGC.CTG..CGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCg                                                                            <  1:4822804/101‑1 (MQ=255)
                            caacactgTCGCTCCCCCCGCGGGGGCGTGGATCGAAAC..GCTG.GCGCGCGC.CGCG.CCGGATCTTTGT.CTG..CGCG.GTCGCTCCCCCCGCGGGGGCGTGGATc                                                                             <  2:1101817/94‑1 (MQ=9)
                            caacactgTCGCTCCCCCCGCGGGGGCGTGGATCGAAAC..GCTG.GCGCGCGC.CGCG.CCGGATCTTTGT.CTG..CGCG.GTCGCTCCCCCCGCGGGGGCGTGGATc                                                                             <  1:5108832/94‑1 (MQ=9)
                             gccccggTCGCTCCCCCCGCGGGGGCGTGGATCGAAACA..ACT.GACC...C.GGCA.ACGGCC.TATGCGCTGGACGTG.GTCGCTCCCCCCCGGGGGGCGCGGATCGa                                                                           >  2:5975680/7‑101 (MQ=2)
                             ccAGACGTCGCGCCCCCCGCGGGGGCGTGGATCGAAACA.GCCT.GCGCG.GC.CGCA.TCGGCT.TCCGC.CTG..CGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGa                                                                           <  2:231277/100‑1 (MQ=255)
                              gttgcGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATGGCTATGGCG.GC.AGTG.TCAGGT.T.TTC.CAG..CGGC.GTCGCTCCCCCCGCGGGGGCGTGGATCGa                                                                           <  1:6306006/97‑1 (MQ=2)
                              cAGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACA.GCCT.GCGCG.GC.CGCA.TCGGCT.TCCGC.CTG..CGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGaa                                                                          >  1:3175547/1‑101 (MQ=255)
                              cAGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACA.GCCT.GCGCG.GC.CGCA.TCGGCT.TCCGC.CTG..CGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGaa                                                                          >  2:6848057/1‑101 (MQ=255)
                               ttgcGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATGGCTATGGCG.GC.AGTG.TCAGGT.T.TTC.CAG..CGGC.GTCGCTCCCCCCGCGGGGGCGTGGATCGaa                                                                          <  1:161251/98‑1 (MQ=2)
                               aGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACA.GCCT.GCGCG.GC.CGCA.TCGGCT.TCCGC.CTG..CGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGaaa                                                                         <  1:702851/101‑1 (MQ=255)
                                actgTCGCTCCCCCCGCGGGGGCGTGGATCGAAAC..GCTG.GCGCGCGC.CGCG.CCGGATCTTTGT.CTG..CGCG.GTCGCTCCCCCCGCGGGGGCGTGGATCGaaa                                                                         >  1:5521182/4‑101 (MQ=9)
                                 tcGTCGTTCCCCCCGCGGGTTCCTGATTCTAAACA.TCAG.GCGCG.GC.CGCA.TCGGCT.TCCGC.CTG..TGTG.GTCGCACCCCCCGCGGGGGCGTGGATCGAAACa                                                                       <  2:405880/100‑1 (MQ=21)
                                 gcGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATGGCTATGGCG.GC.AGTG.TCAGGT.T.TTC.CAG..CGGC.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAAc                                                                        <  1:5854134/100‑1 (MQ=2)
                                 gcGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATGGCTATGGCG.GC.AGTG.TCAGGT.T.TTC.CAG..CGGC.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAAc                                                                        <  2:5003180/100‑1 (MQ=2)
                                 ctgTCGCTCCCCCCGCGGGGGCGTGGATCGAAAC..GCTG.GCGCGCGC.CGCG.CCGGATCTTTGT.CTG..CGCG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAAc                                                                        >  2:3368841/3‑101 (MQ=9)
                                  cGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GCGT.GCGCG.CCATGAATTCATCG.AGCGC.CAT..CGAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGaaa                                                                         >  1:6413968/1‑101 (MQ=2)
                                  cGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACA.GCCT.GCGCG.GC.CGCA.TCGGCT.TCCGC.CTG..CGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACac                                                                      >  2:2687303/1‑100 (MQ=255)
                                   gTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GCGT.GCGCG.CCATGAATTCATCG.AGCGC.CAT..CGAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAAc                                                                        <  2:213600/101‑1 (MQ=2)
                                   gTCGCTCCCCCCGCGGGGGCGTGGATCGAAACA.GCCT.GCGCG.GC.CGCA.TCGGCT.TCCGC.CTG..CGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATc                                                                     <  2:1278016/101‑1 (MQ=255)
                                   gTCGCTCCCCCCGCGGGGGCGTGGATCGAAACA..ACT.GACC...C.GGCA.ACGGCC.TATGCGCTGGACGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAActac                                                                     <  2:1457818/101‑4 (MQ=2)
                                   gTCGCTCCCCCCGCGGGGGCGTGGATCGAAAC..GCTG.GCGCGCGC.CGCG.CCGGATCTTTGT.CTG..CGCG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAActc                                                                      <  2:6528541/101‑3 (MQ=9)
                                   gTCGCCCCCCCCGCGGGGGCGTGGATCGAAACT.GCGT.GCGCG.CCATGAATTCATCG.AGCGC.CAT..CGAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAAc                                                                        <  1:1169098/101‑1 (MQ=2)
                                   gTCGCCCCCCCCGCGGGGGCGTGGATCGAAACA.GCCT.GCGCG.GC.CGCA.TCGGCT.TCCGC.CTG..CGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATc                                                                     <  1:2892583/101‑1 (MQ=255)
                                    tCGCTCCCCCCGCGGGGGCGTGGATCGAAACA.GCCT.GCGCG.GC.CGCA.TCGGCT.TCCGC.CTG..CGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCa                                                                    <  2:4592094/101‑1 (MQ=255)
                                    tCGCTCCCCCCGCGGGGGCGTGGATCGAAACA.GCCT.GCGCG.GC.CGCA.TCGGCT.TCCGC.CTG..CGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCa                                                                    <  1:222824/101‑1 (MQ=255)
                                    tCGCTCCCCCCGCGGGGGCGTGGATCGAAACA.GCCT.GCGCG.GC.CGCA.TCGGCT.TCCGC.CTG..CGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCa                                                                    <  1:23042/101‑1 (MQ=255)
                                    tCGCTCCCCCCGCGGGGGCGTGGATCGAAAC..GCTG.GCGCGCGC.CGCG.CCGGATCTTTGT.CTG..CGCG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAActcg                                                                     >  2:1037758/1‑98 (MQ=9)
                                    tCGCCCCCCCCGCGGGGGCGTGGATCGAAACA.GCCT.GCGCG.GC.CGCA.TCGGCT.TCCGC.CTG..CGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCa                                                                    <  1:6047529/101‑1 (MQ=255)
                                    cccgcccccccGCGGGGGCGTGGATCGAAACA.GCCT.GCGCG.GC.CGCA.TCGGCT.TCCGC.CTG..CGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCa                                                                    <  2:533246/96‑1 (MQ=255)
                                     cGCTCCCCCCGCGGGGGCGTGGATCGAAAC..GCTG.GCGCGCGC.CGCG.CCGGATCTTTGT.CTG..CGCG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAActcgg                                                                    >  2:5863637/1‑97 (MQ=9)
                                      gCTCCCCCCGCGGGGGCGTGGATCGAAACT.GCGT.GCGCG.CCATGAATTCATCG.AGCGC.CAT..CGAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACacc                                                                     <  2:6973197/101‑3 (MQ=2)
                                      gCTCCCCCCGCGGGGGCGTGGATCGAAACA.GCCT.GCGCG.GC.CGCA.TCGGCT.TCCGC.CTG..CGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCAGc                                                                  <  1:2640560/101‑1 (MQ=255)
                                      gCTCCCCCCGCGGGGGCGTGGATCGAAACA.GCCT.GCGCG.GC.CGCA.TCGGCT.TCCGC.CTG..CGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCAGc                                                                  >  1:3389383/1‑101 (MQ=255)
                                      gCTCCCCCCGCGGGGGCGTGGATCGAAAC..GCTG.GCGCGCGC.CGCG.CCGGATCTTTGT.CTG..CGCG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAActcggc                                                                   >  2:6216243/1‑96 (MQ=9)
                                       cTCCCCCCGCGGGGGCGTGGATCGAAACATGGCTATGGCG.GC.AGTG.TCAGGT.T.TTC.CAG..CGGC.GTCGCTCCCCCCGCGGGGGCGTGGATCGaaaaccgccc                                                                  >  2:3071441/1‑94 (MQ=2)
                                       cTCCCCCCGCGGGGGCGTGGATCGAAACA.GCCT.GCGCG.GC.CGCA.TCGGCT.TCCGC.CTG..CGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCAGcc                                                                 >  2:1852078/1‑100 (MQ=255)
                                       cTCCCCCCGCGGGGGCGTGGATCGAAACA..ACT.GACC...C.GGCA.ACGGCC.TATGCGCTGGACGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAActacaccg                                                                 >  1:1683084/1‑94 (MQ=2)
                                         ccccccGCGGGGGCGTGGATCGAAACT.GCGT.GCGCG.CCATGAATTCATCG.AGCGC.CAT..CGAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACACCaat                                                                  <  1:3819063/101‑3 (MQ=2)
                                         ccccccGCGGGGGCGTGGATCGAAACT.GCGT.GCGCG.CCATGAATTCATCG.AGCGC.CAT..CGAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACACCaat                                                                  <  2:4688640/101‑3 (MQ=2)
                                         ccccccGCGGGGGCGTGGATCGAAACATGGCTATGGCG.GC.AGTG.TCAGGT.T.TTC.CAG..CGGC.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAAcccgccacg                                                                <  2:3606681/101‑9 (MQ=0)
                                         ccccccGCGGGGGCGTGGATCGAAACA.GCCT.GCGCG.GC.CGCA.TCGGCT.TCCGC.CTG..CGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCAGCtcg                                                               >  1:4502997/1‑101 (MQ=255)
                                         ccccccGCGGGGGCGTGGATCGAAACA..ACT.GACC...C.GGCA.ACGGCC.TATGCGCTGGACGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAActacaccggg                                                               >  2:893597/1‑92 (MQ=2)
                                          cccccGCGGGGGCGTGGATCGAAACA..ACT.GACC...C.GGCA.ACGGCC.TATGCGCTGGACGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAActacaccgggg                                                              <  1:2828917/101‑11 (MQ=2)
                                            cccGCGGGGGCGTGGATCGAAACT.GCGT.GCGCG.CCATGAATTCATCG.AGCGC.CAT..CGAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACACCaatgcc                                                               <  1:2098123/101‑6 (MQ=0)
                                              cgcgGGGGCGTGGATCGAAACT.GCGT.GCGCG.CCATGAATTCATCG.AGCGC.CAT..CGAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACACCaatgccga                                                             <  1:3556455/101‑8 (MQ=0)
                                               gcgGGGGCGTGGATCGAAACA.GCCT.GCGCG.GC.CGCA.TCGGCT.TCCGC.CTG..CGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCAGCTCG.TCGTTg                                                        <  1:4008226/101‑1 (MQ=255)
                                                  ggggCGTGGATCGAAACA.GCCT.GCGCG.GC.CGCA.TCGGCT.TCCGC.CTG..CGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCAGCTCG.TCGTTGGTc                                                     <  1:4194239/101‑1 (MQ=255)
                                                   gggCGTGGATCGAAACA.GCCT.GCGCG.GC.CGCA.TCGGCT.TCCGC.CTG..CGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCAGCTCG.TCGTTGGTCa                                                    <  2:6670335/101‑1 (MQ=255)
                                                      cGTGGATCGAAACA.GCCT.GCGCG.GC.CGCA.TCGGCT.TCCGC.CTG..CGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCAGCTCG.TCGTTGGTCAGcc                                                 <  2:1381927/101‑1 (MQ=255)
                                                       gTGGATCGAAACA.GCCT.GCGCG.GC.CGCA.TCGGCT.TCCGC.CTG..CGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCCGCTCG.TCGTTGGTCAGCCg                                                >  2:5755535/1‑101 (MQ=255)
                                                        tGGATCGAAACA.GCCT.GCGCG.GC.CGCA.TCGGCT.TCCGC.CTG..CGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCAGCTCG.TCGTTGGTCAGCCgg                                               >  2:5227036/1‑101 (MQ=255)
                                                         ggATCGAAACA.GCCT.GCGCG.GC.CGCA.TCGGCT.TCCGC.CTG..CGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCAGCTCG.TCGTTGGTCAGCCGGa                                              >  2:2983762/1‑101 (MQ=255)
                                                           aTCGAAACA.GCCT.GCGCG.GC.CGCA.TCGGCT.TCCGC.CTG..CGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCAGCTCG.TCGTTGGTCAGCCGGAC.g                                           <  1:2167717/101‑1 (MQ=255)
                                                             cGAAACA.GCCT.GCGCG.GC.CGCA.TCGGCT.TCCGC.CTG..CGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCAGCTCG.TCGTTGGTCAGCCGGAC.GAt                                         >  1:3745500/1‑101 (MQ=255)
                                                                  cA.GCCT.GCGCG.GC.CGCA.TCGGCT.TCCGC.CTG..CGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCAGCTCG.TCGTTGGTCAGCCGGAC.GATGTCCg                                    >  1:2465466/1‑101 (MQ=255)
                                                                   a.GCCT.GCGCG.GC.CGCA.TCGGCT.TCCGC.CTG..CGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCAGCTCG.TCGTTAGTCAGCCGGAC.GATGTCCGt                                   >  1:5631395/1‑101 (MQ=255)
                                                                      ccT.GCGCG.GC.CGCA.TCGGCT.TCCGC.CTG..CGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCAGCTCG.TCGTTGGTCAGCCGGAC.GATGTCCGTCg                                 >  1:4640683/1‑101 (MQ=255)
                                                                        t.GCGCG.GC.CGCA.TCGGCT.TCCGC.CTG..CGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCAGCTCG.TCGTTGGTCAGCCGGAC.GATGTCCGTCGCt                               <  2:4905796/101‑1 (MQ=255)
                                                                            gcg.GC.CGCA.TCGGCT.TCCGC.CTG..CGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCAGCTCG.TCGTTGGTCAGCCGGAC.GATGTCCGTCGCTccc                            <  1:3234189/101‑1 (MQ=255)
                                                                                gC.CGCA.TCGGCT.TCCGC.CTG..CGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCAGCTCG.TCGTTGGTCAGCCGGAC.GATGTCCGTCGCTcccccc                         <  1:6197700/101‑1 (MQ=255)
                                                                                 c.cGCA.TCGGCT.TCCGC.CTG..CGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCAGCTCG.TCGTTGGTCAGCCGGAC.GATGTCCGTCGCTCCCCCcg                        <  2:1983307/101‑1 (MQ=255)
                                                                                    gCA.TCGGCT.TCCGC.CTG..CGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCAGCTCG.TCGTTGGTCAGCCGGAC.GATGTCCGTCGCTCCCCCcgcg                      >  2:3976264/1‑101 (MQ=255)
                                                                                           gCT.TCCGC.CTG..CGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCAGCTCG.TCGTTGGTCAGCCGGAC.GATGTCCGTCGCTCCCCCCGCGGGGGCg                >  1:695906/1‑101 (MQ=255)
                                                                                            cT.TCCGC.CTG..CGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCAGCTCG.TCGTTGGTCAGCCGGAC.GATGTCCGTCGCTCCCCCCGCGGGGGCgg               >  2:4409490/1‑100 (MQ=255)
                                                                                            cT.TCCGC.CTG..CGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCAGCTCG.TCGTTGGTCAGCCGGAC.GATGTCCGTCGCTCCCCCCGCGGGGGCGt               >  2:1521011/1‑101 (MQ=255)
                                                                                                 cGC.CTG..CGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCAGCTCG.TCGTTGGTCAGCCGGAC.GATGTCCGTCGCTCCCCCCGCGGGGGCGTGGAt           <  1:1490653/101‑1 (MQ=255)
                                                                                                  gC.CTG..CGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCAGCTCG.TCGTTGGTCAGCCGGAC.GATGTCCGTCGCTCCCCCCGCGGGGGCGTGGATc          <  2:4219874/101‑1 (MQ=255)
                                                                                                  gC.CTG..CGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCAGCTCG.TCGTTGGTCAGCCGGAC.GATGTCCGTCGCTCCCCCCGCGGGGGCGTGGATc          <  1:6518074/101‑1 (MQ=255)
                                                                                                    ggtcccgccg.gTCGCGCCCCCCGCGGGGGCGTGGATCGAAACATCAGCTTGACCGCATCGAAGCTGGTCACATCGCCGTCGCTCCCCCCGCGGGGGCGTgg             <  2:2024939/92‑1 (MQ=2)
                                                                                                     atcccgccg.gTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCAGCTTGACCGCATCGAAGCTGGTCACATCGCCGTCGCTCCCCCCGCGGGGGCGTGGa            <  1:944340/93‑1 (MQ=2)
                                                                                                      tcccgccg.gTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCAGCTTGACCGCATCGAAGCTGGTCACATCGCCGTCGCTCCCCCCGCGGGGGCGTGGAt           >  1:6728917/8‑101 (MQ=2)
                                                                                                      tG..CGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCAGCTCG.TCGTTGGTCAGCCGGAC.GATGTCCGTCGCTCCCCCCGCGGGGGCGTGGATCGaa       <  1:5259699/101‑1 (MQ=255)
                                                                                                      tG..CGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCAGCTCG.TCGTTGGTCAGCCGGAC.GATGTCCGTCGCTCCCCCCGCGGGGGCGTGGATCGaa       <  2:6151447/101‑1 (MQ=255)
                                                                                                        ccgccg.gTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCAGCTTGACCGCATCGAAGCTGGTCACATCGCCGTCGCTCCCCCCGCGGGGGCGTGGATCg         <  1:3918233/96‑1 (MQ=2)
                                                                                                        ccgccg.gTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCAGCTTGACCGCATCGAAGCTGGTCACATCGCCGTCGCTCCCCCCGCGGGGGCGTGGATCg         <  1:5414395/96‑1 (MQ=2)
                                                                                                         cgccg.gTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCAGCTTGACCGCATCGAAGCTGGTCACATCGCCGTCGCTCCCCCCGCGGGGGCGTGGATCGa        <  2:2869386/97‑1 (MQ=2)
                                                                                                          gccg.gTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCAGCTTGACCGCATCGAAGCTGGTCACATCGCCGTCGCTCCCCCCGCGGGGGCGTGGATCGaa       <  2:5548641/98‑1 (MQ=2)
                                                                                                          gccg.gTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCAGCTTGACCGCATCGAAGCTGGTCACATCGCCGTCGCTCCCCCCGCGGGGGCGTGGATCGaa       <  2:7090206/98‑1 (MQ=2)
                                                                                                           gTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCAGCTCG.TCGTTGGTCAGCCGGAC.GATGTCCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACt    >  1:1003564/1‑101 (MQ=255)
                                                                                                            tG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCAGCTCG.TCGTTGGTCAGCCGGAC.GATGTCCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTg   >  1:4754998/1‑101 (MQ=255)
                                                                                                            cg.gTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCAGCTTGACCGCATCGAAGCTGGTCACATCGCCGTCGCTCCCCCTGCGGGGGCGTGGATCGAAAc     >  1:4372443/2‑101 (MQ=2)
                                                                                                             g.gTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCAGCTCG.TCGTTGGTCAGCCGGAC.GATGTCCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGc  <  1:741681/101‑1 (MQ=255)
                                                                                                             |                                                                                                       
CGGTCCAGAGATAGACCGCGCGGACCTCCTCAGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACA.GCCT.GCGCG.GC.CGCA.TCGGCT.TCCGC.CTG..CGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCAGCTCG.TCGTTGGTCAGCCGGAC.GATGTCCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGC  >  NC_014034/1396868‑1397068

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