Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_014034 1,397,016 G→C 19.8% V300L (GTC→CTA)  RCAP_rcc01307 → hypothetical protein

Read alignment evidence...
  seq id position change freq score reads annotation genes product
*NC_0140341,397,0160G→C19.8% 97.1 ‑1.0 ‑0.0219V300L (GTC→CTC) RCAP_rcc01307hypothetical protein
Reads supporting (aligned to +/- strand):  new base (15/23):  ref base (79/84):  total (112/107)
Fisher's exact test strand distribution p-value = 3.69e-01
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 9.59e-01

Only 100 of 233 total aligned reads displayed.

CGGGGGCGTGGATCGAAACAGCCTGCGCGGCCGCATCGGCTTCCGCCTGCGTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCAGCTCG.TCGTTGGTCAGCCGGAC.GATGTCC...GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCCGTCTTGAGTGATGGGGCGCGAGCACATGTCGTCGCTCCCCCCGCGGG  >  NC_014034/1396916‑1397116
                                                                                                     |                                                                                                        
 gggggCGTGGATCGAAACAGCCTGCGCGGCCGCATCGGCTTCCGCCTGCGTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCAGCTCG.TCGTTGGt                                                                                                         >  1:391132/1‑101 (MQ=255)
 gggggCGTGGATCGAAACAGCCTGCGCGGCCGCATCGGCTTCCGCCTGCGTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCAGCTCG.TCGTTGGt                                                                                                         >  2:1269658/1‑101 (MQ=255)
   gggCGTGGATCGAAACAGCCTGCGCGGCCGCATCGGCTTCCGCCTGCGTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCAGCTCG.TCGTTGGTCa                                                                                                       >  2:3068961/1‑101 (MQ=255)
    ggCGTGGATCGAAACAGCCTGCGCGGCCGCATCGGCTTCCGCCTGCGTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCAGCTCG.TCGTTGGTCAg                                                                                                      <  1:1579571/101‑1 (MQ=255)
     gCGTGGATCGAACCAGCCTGCGCGGCCGCATCGGCTTCCGCCTGCGTGGTCGCTGGCCCCGCGGGGGCGTGGATCGAAACATCAGCTCG.TCGTTGGTCAGc                                                                                                     <  1:4178471/101‑1 (MQ=255)
     gCGTGGATCGAAACAGCCTGCGCGGCCGCATCGGCTTCCGCCTGCGTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCAGCTCG.TCGTTGGTCAGc                                                                                                     <  2:1797123/101‑1 (MQ=255)
        tGGATCGAAACAGCCTGCGCGGCCGCATCGGCTTCCGCCTGCGTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCAGCTCG.TCGTTGGTCAGCCgg                                                                                                  <  1:3361384/101‑1 (MQ=255)
         ggATCGAAACAGCCTGCGCGGCCGCATCGGCTTCCGCCTGCGTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCAGCTCG.TCGTTGGTCAGCCGGa                                                                                                 >  1:641595/1‑101 (MQ=255)
          gATCGAAACCGCCTGCGCGGCCGCATCGGCTTCCGCCTGCGTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCAGCTCG.TCGTTGGTCAGCCGGAc                                                                                                <  1:4132035/101‑1 (MQ=255)
          gATCGAAACAGCCTGCGCGGCCGCATCGGCTTCCGCCTGCGTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCAGCTCG.TCGTTGGTCAGCCGGAc                                                                                                <  2:2067419/101‑1 (MQ=255)
           aTCGAAACAGCCTGCGCGGCCGCATCGGCTTCCGCCTGCGTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCAGCTCG.TCGTTGGTCAGCCGGAC.g                                                                                              <  2:3891940/101‑1 (MQ=255)
               aaaCAGCCTGCGCGGCCGCATCGGCTTCCGCCTGCGTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCAGCTCG.TCGTTGGTCAGCCGGAC.GATGt                                                                                          >  1:769828/1‑101 (MQ=255)
               aaaCAGCCTGCGCGGCCGCATCGGCTTCCGCCTGCGTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCAGCTCG.TCGTTGGTCAGCCGGAC.GATGt                                                                                          >  1:2899939/1‑101 (MQ=255)
                  cAGCCTGCGCGGCCGCATCGGCTTCCGCCTGCGTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCAGCTCG.TCGTTGGTCAGCCGGAC.GATGTAC...g                                                                                    <  2:3518471/101‑1 (MQ=255)
                   aGCCTGCGCGGCCGCATCGGCTTCCGCCTGCGTGGTCGCTCCCCCCGCGGGGGCGGGGATCGAAACATCAGCTCG.TCGTTGGTCAGCCGGAC.GATGTCC...Gt                                                                                   >  2:1068396/1‑101 (MQ=255)
                    gCCTGCGCGGCCGCATCGCCTTCCGCCTGCGTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCAGCTCG.TCGTTGGTCAGCCGGAC.GATGTCC...GTc                                                                                  <  2:1934989/101‑1 (MQ=255)
                     ccTGCGCGGCCGCATCGGCTTCCGCCTGCGTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCAGCTCG.TCGTTGGTCAGCCGGAC.GATGTCC...GTCg                                                                                 <  2:641518/101‑1 (MQ=255)
                      cTGCGCGGCCGCATCGGCTTCCGCCTGCGTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCAGCTCG.TCGTTGGTCAGCCGGAC.GGTGTCC...GTCGc                                                                                >  1:3720226/1‑101 (MQ=255)
                       tGCGCGGCCGCATCGGCTTCCGCCTGCGTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCAGCTCG.TCGTTGGTCAGCCGGAC.GATGTCC...GTCGCt                                                                               <  2:3257704/101‑1 (MQ=255)
                       tGCGCGGCCGCATCGGCTTCCGCCTGCGTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCAGCTCG.TCGTTGGTCAGCCGGAC.GATGTCC...GTCGCt                                                                               <  1:3125621/101‑1 (MQ=255)
                       tGCGCGGCCGCATCGGCTTCCGCCTGCGTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCAGCTCG.TCGTTGGTCAGCCGGAC.GATGTCC...GTCGCt                                                                               <  2:1108690/101‑1 (MQ=255)
                            ggCCGCATCGGCTTCCGCCTGCGTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCAGCTCG.TCGTTGGTCAGCCGGAC.GATGTCC...GTCGCTccccc                                                                          >  1:1038889/1‑101 (MQ=255)
                               cGCATCGGCTTCCGCCTGCGTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCAGCTCG.TCGTTGGTCAGCCGGAC.GATGTCC...GTCGCTCCCCCcgc                                                                       <  2:559397/101‑1 (MQ=255)
                                  aTCGGCTTCCGCCTGCGTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCAGCTCG.TCGTTGGTCAGCCGGAC.GATGTCC...GTCGCTCCCCCCGCggg                                                                    >  2:1591196/1‑101 (MQ=255)
                                  aTCGGCTTCCGCCTGCGTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCAGCTCG.TCGTTGGCCAGCCGGAC.GATGTCC...GTCCCTCCCCCCGCggg                                                                    >  1:3644290/1‑101 (MQ=255)
                                      gCTTCCGCCTGCGTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCAGCTCG.TCGTTGGTCAGCCGGAC.GATGTCC...GTCGCTCCCCCCGCGGGGGCg                                                                >  1:1353852/1‑101 (MQ=255)
                                       cTTCCGCCTGCGTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCAGCTCG.TCGTTGGTCAGCCGGAC.GATGTCC...GTCGCTCCCCCCGCGGGGGCgg                                                               >  2:544704/1‑100 (MQ=255)
                                        ttCCGCCTGCGTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCAGCTCG.TCGTTGGTCAGCCGGAC.GATGTCC...GTCGCTCCCCCCGCGGGGGCGTg                                                              >  2:317413/1‑101 (MQ=255)
                                          ccGCCTGCGTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCAGCTCG.TCGTTGGTCAGCCGGAC.GATGTCC...GTCGCTCCCCCCGCGGGGGCGTGGa                                                            >  2:316410/1‑101 (MQ=255)
                                           gatcccgccggTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCAGCTTGACCGCATCGAAGCTGGTCACATCGCC...GTCGCTCCCCCCGCGGGGGCGTgg                                                             <  2:2729424/92‑1 (MQ=2)
                                           cGCCTGCGTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCAGCTCG.TCGTTGGTCAGCCGGAC.GATGTCC...GTCGCTCCCCCCGCGGGGGCGTGGAt                                                           <  1:1222848/101‑1 (MQ=255)
                                            gCCTGCGTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCAGCTCG.TCGTTGGTCAGCCGGAC.GATGTCC...GTCGCTCCCCCCGCGGGGGCGTGGATc                                                          <  1:3566638/101‑1 (MQ=255)
                                             gcgtgccggTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCAGCTTGACCGCATCGAAGCTGGTCACATCGCC...GTCGCTCCCCCCGCGGGGGCGTGGAt                                                           >  2:1424000/8‑101 (MQ=2)
                                             ccTGCGTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCAGCTCG.TCGTTGGTCAGCCGGAC.GATGTCC...GTCGCTCCCCCCGCGGGGGCGTGGATCg                                                         <  2:281081/101‑1 (MQ=255)
                                              cccgccggTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCAGCTTGACCGCATCGAAGCTGGTCACATCGCC...GTCGCTCCCCCCGCGGGGGCGTGGATc                                                          <  1:2422780/95‑1 (MQ=2)
                                              cccgccggTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCAGCTTGACCGCATCGAAGCTGGTCACATCGCC...GTCGCTCCCCCCGCGGGGGCGTGGATc                                                          >  1:247355/7‑101 (MQ=2)
                                              cagcggcgTCGCTCCCCCCGCGGGGGCGTGGATCGAAAC.CCGCCACG.TCTTGGAACAGCGTGCC.ATTGGCCCCGGTCGCTCCCCCCGCGGGGGCGTGGATc                                                          >  2:41432/8‑101 (MQ=2)
                                              cTGCGTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCAGCTCG.TCGTTGGTCAGCCGGAC.GATGTCC...GTCGCTCCCCCCGCGGGGGCGTGGATCGa                                                        <  2:1038798/101‑1 (MQ=255)
                                              cTGCGTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCAGCTCG.TCGTTGGTCAGCCGGAC.GATGTCC...GTCGCTCCCCCCGCGGGGGCGTGGATCGa                                                        >  2:2504107/1‑101 (MQ=255)
                                               tGCGTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCAGCTCG.TCGTTGGTCAGCCGGAC.GATGTCC...GTCGCTCCCCCCGCGGGGGCGTGGATCGaa                                                       <  2:1340577/101‑1 (MQ=255)
                                               tGCGTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCAGCTCG.TCGTTGGTCAGCCGGAC.GATGTCC...GTCGCTCCCCCCGCGGGGGCGTGGATCGaa                                                       <  2:801482/101‑1 (MQ=255)
                                               tGCGTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCAGCTCG.TCGTTGGTCAGCCGGAC.GATGTCC...GTCGCTCCCCCCGCGGGGGCGTGGATCGaa                                                       <  2:3450813/101‑1 (MQ=255)
                                               ccgccggTCGCTCCCCCCGCGGGGGCGTGGATCCAAACATCAGCTTGACCGCATCGAAGCTGGTCACATCGCC...GTCGCTCCCCCCGCGGGGGCGTGGATCg                                                         <  1:4208987/96‑1 (MQ=2)
                                               agcggcgTCGCTCCCCCCGCGGGGGCGTGGATCGAAAC.CCGCCACG.TCTTGGAACAGCGTGCC.ATTGGCCCCGGTCGCTCCCCCCGCGGGGGCGTGGATCg                                                         >  1:594992/7‑101 (MQ=2)
                                                gCGTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCAGCTCG.TCGTTGGTCAGCCGGAC.GATGTCC...GTCGCTCCCCCCGCGGGGGCGTGGATCGaaa                                                      >  1:3544371/1‑101 (MQ=255)
                                                 gccggTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCAGCTTGACCGCATCGAAGCTGGTCACATCGCC...GTCGCTCCCCCCGCGGGGGCGTGGATCGaa                                                       >  2:900211/4‑101 (MQ=2)
                                                 cggcgTCGCTCCCCCCGCGGGGGCGTGGATCGAAAC.CCGCCACG.TCTTGGAACAGCGTGCC.ATTGGCCCCGGTCGCTCCCCCCGCGGGGGCGTGGATCGaa                                                       >  1:104632/5‑101 (MQ=2)
                                                  ggcgTCGCTCCCCCCGCGGGGGCGTGGATCGAAAC.CCGCCACG.TCTTGGAACAGCGTGCC.ATTGGCCCCGGTCGCTCCCCCCGCGGGGGCGTGGATCGaaa                                                      >  2:1618386/4‑101 (MQ=0)
                                                  gTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCAGCTCG.TCGTTGGTCAGCCGGAC.GATGTCC...GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACt                                                    <  1:1088998/101‑1 (MQ=255)
                                                  gTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCAGCTCG.TCGTTGGTCAGCCGGAC.GATGTCC...GTCCCTCCCCCCGCGGGGGCGCGGGCCGaaagt                                                    >  1:3958755/1‑99 (MQ=37)
                                                  ccggTCCCTCCCCCCGCGGGGGCGTGGATCGAAACATCAGCTTGACCGCATCGAAGCTGGTCACATCGCC...GTCGCTCCCCCCGCGGGGGCGTGGATCGaaa                                                      <  1:457378/99‑1 (MQ=2)
                                                   tGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCAGCTCG.TCGTTGGTCAGCCGGAC.GATGTCC...GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTg                                                   <  2:287412/101‑1 (MQ=255)
                                                   tGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCAGCTCG.TCGTTGGTCAGCCGGAC.GATGTCC...GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTg                                                   >  1:1909/1‑101 (MQ=255)
                                                   tGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCAGCTCG.TCGTTGGTCAGCCGGAC.GATGTCC...GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTg                                                   >  2:1560910/1‑101 (MQ=255)
                                                   cggTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCAGCTTGACCGCATCGAAGCTGGTCACATCGCC...GTCGCTCCCCCCGCGGGGGCGTGGATCGAAAc                                                     <  1:1394292/100‑1 (MQ=2)
                                                   cggTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCAGCTTGACCGCATCGAAGCTGGTCACATCGCC...GTCGCTCCCCCCGCGGGGGCGTGGATCGAAAc                                                     <  2:3104666/100‑1 (MQ=2)
                                                    ggTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCAGCTTGACCGCATCGAAGCTGGTCACATCGCC...GTCGCTCCCCCCGCGGGGGCGTGGATCGAAAcg                                                    >  2:2272359/1‑100 (MQ=2)
                                                    ggTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCAGCTCG.TCGTTGGTCAGCCGGAC.GATGTCC...GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGc                                                  <  1:3419684/101‑1 (MQ=255)
                                                    ggTCGCGCCCCCCGCGGGGGCGTGGATCGAAACATCAGCTCG.TCGTTGGTCAGCCGGAC.GATGTCC...GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGc                                                  <  1:1717563/101‑1 (MQ=255)
                                                     gTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCAGCTCG.TCGTTGGTCAGCCGGAC.GATGTCC...GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGcc                                                 <  1:166502/101‑1 (MQ=255)
                                                     gTCGCTCCCCCCGCGGGGGCGTGGATCGAAAC.CCGCCACG.TCTTGGAACAGCGTGCC.ATTGGCCCCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAAcaa                                                   >  2:2588455/1‑99 (MQ=2)
                                                      tCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCAGCTCG.TCGTTGGTCAGCCGGAC.GATGTCC...GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCCg                                                <  1:3845500/101‑1 (MQ=255)
                                                      tCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCAGCTCG.TCGTTGGTCAGCCGGAC.GATGTCC...GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCCg                                                <  1:2248836/101‑1 (MQ=255)
                                                       cGCTCCCCCCGCGGGGGCGTGGATCTAAACATCAGCTCG.TCGTTTGTCAGCCGGAC.GATGTCC...GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCCGt                                               <  2:2423515/101‑1 (MQ=255)
                                                       cGCTCCCCCCGCGGGGGCGTGGATCGAAAC.CCGCCACG.TCTTGGAACAGCGTGCC.ATTGGCCCCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAAcaact                                                 >  1:2800531/1‑97 (MQ=2)
                                                        gCTCCCCCCGCGGGGGCGTGGATCGAAACATCAGCTTGACCGCATCGAAGCTGGTCACATCTCC...GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGGcgt                                                <  2:2527389/101‑3 (MQ=2)
                                                        gCTCCCCCCGCGGGGGCGTGGATCGAAAC.CCGCCACG.TCTTGGAACAGCGTGCC.ATTGGCCCCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAAcaactg                                                >  2:2378021/1‑96 (MQ=2)
                                                        gCTCCCCCCGCGGGGGCGTGGATCGAAAC.CCGCCACG.TCTTGGAACAGCGTGCC.ATTGGCCCCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAAcaactg                                                >  1:12822/1‑96 (MQ=2)
                                                         cTCCCCCCGCGGGGGCGTGGATCGAAACATCAGCTTGACCGCATCGAAGCTGGTCACATCGCC...GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGGcgta                                               >  2:4025181/1‑98 (MQ=2)
                                                         cTCCCCCCGCGGGGGCGTGGATCGAAACATCAGCTCG.TCGTTGGTCAGCCGGAC.GATGTCC...GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCCGTCt                                             >  1:2498002/1‑101 (MQ=255)
                                                         cTCCCCCCGCGGGGGCGTGGATCGAAAC.CCGCCACG.TCTTGGAACAGCGTGCC.ATTGGCCCCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAAcaactga                                               >  2:2933242/1‑95 (MQ=2)
                                                          cccccccGCGGGGGCGTGGATCGAAACATCAGCTCG.TCGTTGGTCAGCCGGAC.GATGTCC...GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCCGTCtt                                            <  2:691260/100‑1 (MQ=255)
                                                           ccccccGCGGGGGCGTGGATCGAAACATCAGCTTGACCGCATCGAAGCTGGTCACATCGCC...GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGGcgtaga                                             >  1:2193943/1‑96 (MQ=2)
                                                           ccccccGCGGGGGCGTGGATCGAAACATCAGCTCG.TCGTTGGTCAGCCGGAC.GATGTCC...GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCCGTCTTg                                           >  1:1419940/1‑101 (MQ=255)
                                                           ccccccGCGGGGGCGTGGATCGAAACATCAGCTCG.TCGTTGGTCAGCCGGAC.GATGTCC...GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCCGTCTTg                                           >  1:3084131/1‑101 (MQ=255)
                                                             ccccGCGGGGGCGTGGATCGAAACATCAGCTCG.TCGTTGGTCAGCCGGAC.GATGTCC...GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCCGTCTTGAg                                         <  2:3319454/101‑1 (MQ=255)
                                                              cccGGGGGGGCGTGGATCGAAACATCAGCTCG.TCGTTGGTCAGCCGGAC.GATGTCC...GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCCGTCTTGAGt                                        <  2:1162674/101‑1 (MQ=255)
                                                               ccGTGGGGGCGTGGATCGAAACATCAGCTTGACCGCATCGAAGCTGGTCACATCGCC...GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGGcgtagatctc                                         <  1:3797443/101‑10 (MQ=2)
                                                                  cgGGGGCGTGGATCGAAACATCAGCTCG.TCGTTGGTCAGCCGGAC.GATGTCC...GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCCGTCTTGAGTGATg                                    >  1:260904/1‑101 (MQ=255)
                                                                   gggggCGTGGATAGAAACATCAGCTCG.TCGTTGGTCAGCCGGAC.GATGTCC...GTCGCTCCCCCCGCGGGGGCGTGGATCGCAACTGGCGGCTTGAGTGATgg                                   >  1:3350459/1‑101 (MQ=38)
                                                                    ggggCGTGGATCGAAACATCAGCTCG.TCGTTGGTCAGCCGGAC.GATGTCC...GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCCGTCTTGAGTGATggg                                  >  1:1971143/1‑101 (MQ=255)
                                                                    ggggCGTGGATCGAAACATCAGCTCG.TCGTTGGTCAGCCGGAC.GATGTCC...GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCCGTCTTGAGGGATggg                                  >  2:3578131/1‑101 (MQ=255)
                                                                         gTGGATCGAAACATCAGCTCG.TCGTTGGTCAGCCGGAC.GATGTCC...GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCCGTCTTGAGTGATGGGGcgcg                             >  2:2147484/1‑101 (MQ=255)
                                                                          tGGATCGAAACATCAGCTCG.TCGTTGGTCAGCCGGAC.GATGTCC...GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCCGTCTTGAGTGATGGGGCGCGa                            >  1:3694879/1‑101 (MQ=255)
                                                                           ggATCGAAACATCAGCTCG.TCGTTGGTCAGCCGGAC.GATGTCC...GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCCGTCTTGAGTGATGGGGCGCGAg                           >  1:2053963/1‑101 (MQ=255)
                                                                               cGAAACATCAGCTCG.TCGTTGGTCAGCCGGAC.GATGTCC...GTCGCTCCCCCCGCGGGGGCGTGGATCGGAACTGCCGTCTTGAGTGATGGGGCGCGAGcaca                       >  1:3604314/1‑101 (MQ=255)
                                                                                 aaaCATCAGCTCG.TCGTTGGTCAGCCGGAC.GATGTCC...GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCCGTCTTGAGTGATGGGGCGCGAGCACATg                     >  2:2645615/1‑101 (MQ=255)
                                                                                  aaCATCAGCTCG.TCGTTGGTCAGCCGGAC.GATGTCC...GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCCGTCTTGAGTGATGGGGCGCGAGCACATgt                    >  1:2972715/1‑101 (MQ=255)
                                                                                   aCATCAGCTCG.TCGTTGGTCAGCCGGAC.GATGTCC...GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCCGTCTTGAGTGATGGGGCGCGAGCCCATgtc                   >  2:3356885/1‑101 (MQ=255)
                                                                                   aCATCAGCTCG.TCGTTGGTCAGCCGGAC.GATGTCC...GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCCGTCTTGAGTGATGGGGCGCGAGCACATgtc                   <  1:1525287/101‑1 (MQ=255)
                                                                                    cATCAGCTCG.TCGTTGGTCAGCCGGAC.GATGTCC...GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCCGTCTTGAGTGATGGGGCGCGAGCACATgtcg                  >  1:3308422/1‑101 (MQ=255)
                                                                                     aTCAGCTCG.TCGTTGGTCAGCCGGAC.GATGTCC...GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCCGCCTTGAGTGATGGGGCGCGCGCCCATgtcgc                 >  1:1860660/1‑100 (MQ=255)
                                                                                         gCTCG.TCGTTGGTCAGCCGGAC.GATGTCC...GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCCGTCTTGAGTGATGGGGCGCGAGCACATGTCGTCGCt             <  2:3705442/101‑1 (MQ=255)
                                                                                           tcg.tcgTTGGTCAGCCGGAC.GATGTCC...GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCCGTCTTGAGTGATGGGGCGCGAGCACATGTCGTCGCTcc           >  2:4231140/1‑101 (MQ=255)
                                                                                             g.tcgtTGGTCAGCCGGAC.GATGTCC...GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCCGTCTTGAGTGATGGGGCGCGAGCACATGTCGTCGCTcccc         <  2:2928050/101‑1 (MQ=255)
                                                                                               tcgtTGGTCAGCCGGAC.GATGTCC...GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCCGTCTTGAGTGATGGGGCGCGAGCACATGTCGTCGCTccccc        >  1:883300/1‑101 (MQ=255)
                                                                                               tcgtTGGTCAGCCGGAC.GATGTCC...GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCCGTCTTGAGTGATGGGGCGCGAGCACATGTCGTCGCTccccc        <  2:391081/101‑1 (MQ=255)
                                                                                               tcgtTGGTCAGCCGGAC.GATGTCC...GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCCGTCTTGAGTGATGGGGCGCGAGCACATGTCGTCGCTccccc        >  1:1645164/1‑101 (MQ=255)
                                                                                                cgtTGGTCAGCCGGAC.GATGTCC...GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCCGTCTTGAGTGATGGGGCGCGAGCACATGTCGTCGCTcccccc       <  2:3749567/101‑1 (MQ=255)
                                                                                                     gTCAGCGGGGC.GATGTCC...GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGGCATCTCCAGTAATGGCGCGCGCACACGTGTCGTCGCTCCCTCcgcgcg  >  2:4208026/1‑99 (MQ=16)
                                                                                                     |                                                                                                        
CGGGGGCGTGGATCGAAACAGCCTGCGCGGCCGCATCGGCTTCCGCCTGCGTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCAGCTCG.TCGTTGGTCAGCCGGAC.GATGTCC...GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCCGTCTTGAGTGATGGGGCGCGAGCACATGTCGTCGCTCCCCCCGCGGG  >  NC_014034/1396916‑1397116

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