Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_014034 1,397,030 G→C 17.6% M304I (ATG→ATC RCAP_rcc01307 → hypothetical protein

Read alignment evidence...
  seq id position change freq score reads annotation genes product
*NC_0140341,397,0300G→C17.6% 97.1 ‑2.6 ‑0.6216M304I (ATG→ATCRCAP_rcc01307hypothetical protein
Reads supporting (aligned to +/- strand):  new base (15/23):  ref base (100/78):  total (115/101)
Fisher's exact test strand distribution p-value = 7.37e-02
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 2.38e-01

Only 100 of 235 total aligned reads displayed.

GAAACAGCCTGCGCGGCCGCATCGGCTTCCGCCTGCGTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCAGCTCG.TCGTTGGTCAGCCGGAC.GATGTCC...GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCCGTCTTGAGTGATGGGGCGCGAGCACATGTCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAA  >  NC_014034/1396930‑1397130
                                                                                                      |                                                                                                       
gAAACAGCCTGCGCGGCCGCATCGGCTTCCGCCTGCGTGGTCTCTCCCCCCGCGGGGGCGTGGATCGAAACATCAGCTCG.TCGTTCGTCAGCCGGAC.GATg                                                                                                         >  1:476185/1‑101 (MQ=255)
   aCAGCCTGCGCGGCCGCATCGGCTTCCGCCTGCGTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCAGCTCG.TCGTTGGTCAGCCGGAC.GATGTcc                                                                                                      <  1:3699062/101‑1 (MQ=255)
    cAGCCTGCGCGGCCGCATCGGCTTCCGCCTGCGTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCAGCTCG.TCGTTGGTCAGCCGGGC.GATGTCC...g                                                                                                  >  2:381142/1‑101 (MQ=255)
    cAGCCTGCGCGGCCGCATCGGCTTCCGCCTGCGTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCAGCTCG.TCGTTGGTCAGCCGGAC.GATGTAC...g                                                                                                  <  2:3518471/101‑1 (MQ=255)
      gCCTGCGCGGCCGCATCGGCTTCCGCCTGCGTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCAGCTCG.TCGTTGGTCAGCCGGAC.GATGTCC...GTc                                                                                                <  2:3251375/101‑1 (MQ=255)
       ccTGCGCGGCCGCATCGGCTTCCGCCTGCGTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCAGCTCG.TCGTTGGTCAGCCGGAC.GATGTCC...GTCg                                                                                               <  2:641518/101‑1 (MQ=255)
        cTGCGCGGCCGCATCGGCTTCCGCCTGCGTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCAGCTCG.TCGTTGGTCAGCCGGAC.GATGTCC...GTCGc                                                                                              >  2:2782294/1‑101 (MQ=255)
              ggCCGCATCGGCTTCCGCCTGCGTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCAGCTCG.TCGTTGGTCAGCCGGAC.GATGTCC...GTCGCTccccc                                                                                        >  1:1038889/1‑101 (MQ=255)
               gCCGCATCGGCTTCCGCCTGCGTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCAGCTCG.TCGTTGGTCAGCCGGAC.GATGTCC...GTCGCTcccccc                                                                                       <  1:1410445/101‑1 (MQ=255)
                ccGCATCGGCTTCCGCCTGCGTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCAGCTCG.TCGTTGGTCAGCCGGAC.GATGTCC...GTCGCTCCCCCcg                                                                                      >  1:2729697/1‑101 (MQ=255)
                 cGCATCGGCTTCCGCCTGCGTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCAGCTCG.TCGTTGGTCAGCCGGAC.GATGTCC...GTCGCTCCCCCcgc                                                                                     <  2:559397/101‑1 (MQ=255)
                  gCATCGGCTTCCGCCTGCGTGGTCGCTCCCCCCGCGGGGGCGTGGATCGGAACATAAGCTCG.TCGTTGGTCAGCAGGAC.GATGTCC...GTCGCTCCCCCcgcg                                                                                    >  1:3555133/1‑101 (MQ=255)
                      cGGCTTCCGCCTGCGTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCAGCTCG.TCGTTGGTCAGCCGGAC.GATGTCC...GTCGCCCCCCCCGCggggg                                                                                >  2:334141/1‑101 (MQ=255)
                       ggCTTCCGCCTGCGTGGTCGCCCCCCCCGCGGGGGCGTGGATCGAAACATCAGCTCG.TCGTTGGTCAGCCGGAC.GATGTCC...GTCGCTCCCCCCGCGGGGGc                                                                               <  2:3714353/101‑1 (MQ=255)
                          ttCCGCCTGCGTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCAGCTCG.TCGTTGGTCAGCCGGAC.GATGTCC...GTCGCTCCCCCCGCGGGGGCGTg                                                                            >  1:4141334/1‑101 (MQ=255)
                           tCCGCCTGCGTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCAGCTCG.TCGTTGGTCAGCCGGAC.GATGTCC...GTCGCTCCCCCCGCGGGGGCGTgg                                                                           >  2:1991795/1‑101 (MQ=255)
                           tCCGCCTGCGTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCAGCTCG.TCGTTGGTCAGCCGGAC.GATGTCC...GTCGCTCCCCCCGCGGGGGCGTgg                                                                           <  1:1959975/101‑1 (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:1784257/1‑101 (MQ=255)
                              gCCTGCGTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCAGCTCG.TCGTTGGTCAGCCGGAC.GATGTCC...GTCGCTCCCCCCGCGGGGGCGTGGATc                                                                        <  1:1126987/101‑1 (MQ=255)
                              gCCTGCGTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCAGCTCG.TCGTTGGTCAGCCGGAC.GATGTCC...GTCGCTCCCCCCGCGGGGGCGCGGATc                                                                        >  1:2423691/1‑101 (MQ=255)
                              cccTGGGTGGTCGCCCCCCCCGCGGGGGGGTGGATCGAAACATCAGCTCG.TCGTTGGTCAGCCGGAC.GATGTCC...GTCGCTCCCCCCGCGGGGGCGTGGATc                                                                        <  2:538218/100‑1 (MQ=255)
                              atcccgccggTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCAGCTTGACCGCATCGAAGCTGGTCACATCGCC...GTCGCTCCCCCCGCGGGGGCGTGGa                                                                          <  1:2586171/93‑1 (MQ=2)
                               gcgtgccggTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCAGCTTGACCGCATCGAAGCTGGTCACATCGCC...GTCGCTCCCCCCGCGGGGGCGTGGAt                                                                         >  2:1424000/8‑101 (MQ=2)
                                cccgccggTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCAGCTTGACCGCATCGAAGCTGGTCACATCGCC...GTCGCTCCCCCCGCGGGGGCGTGGATc                                                                        >  1:247355/7‑101 (MQ=2)
                                cTGCGTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCAGCTCG.TCGTTGGTCAGCCGGAC.GATGTCC...GTCGCTCCCCCCGCGGGGGCGTGGATCGa                                                                      >  2:2504107/1‑101 (MQ=255)
                                 tGCGTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCAGCTCG.TCGTTGGTCAGCCGGAC.GATGTCC...GTCGCTCCCCCCGCGGGGGCGTGGATCGaa                                                                     <  1:2830518/101‑1 (MQ=255)
                                 tGCGTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCAGCTCG.TCGTTGGTCAGCCGGAC.GATGTCC...GTCGCTCCCCCCGCGGGGGCGTGGATCGaa                                                                     <  2:1340577/101‑1 (MQ=255)
                                 tGCGTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCAGCTCG.TCGTTGGTCAGCCGGAC.GATGTCC...GTCGCTCCCCCCGCGGGGGCGTGGATCGaa                                                                     >  2:1271464/1‑101 (MQ=255)
                                 tGCGTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCAGCTCG.TCGTTGGTCAGCCGGAC.GATGTCC...GTCGCTCCCCCCGCGGGGGCGTGGATCGaa                                                                     <  2:3450813/101‑1 (MQ=255)
                                 tGCGTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCAGCTCG.TCGTTGGTCAGCCGGAC.GATGTCC...GTCGCTCCCCCCGCGGGGGCGTGGATCGaa                                                                     <  1:3472301/101‑1 (MQ=255)
                                  cgccggTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCAGCTTGACCGCATCGAAGCTGGTCACATCGCC...GTCGCTCCCCCCGCGGGGGCGTGGATCGa                                                                      <  1:1280086/97‑1 (MQ=2)
                                   gccggTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCAGCTTGACCGCATCGAAGCTGGTCACATCGCC...GTCGCTCCCCCCGCGGGGGCGTGGATCGaa                                                                     >  2:900211/4‑101 (MQ=2)
                                   cggcgTCGCTCCCCCCGCGGGGGCGTGGATCGAAAC.CCGCCACG.TCTTGGAACAGCGTGCC.ATTGGCCCCGGTCGCTCCCCCCGCGGGGGCGTGGATCGaa                                                                     >  1:104632/5‑101 (MQ=2)
                                   cGTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCAGCTCG.TCGTTGGTCAGCCGGAC.GATGTCC...GTCGCTCCCCCCGCGGGGGCGTGGATCGAAAc                                                                   >  2:3573389/1‑101 (MQ=255)
                                   cGTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCAGCTCG.TCGTTGGTCAGCCGGAC.GATGTCC...GTCGCTCCCCCCGCGGGGGCGTGGATCGAAAc                                                                   >  1:85659/1‑101 (MQ=255)
                                    ggcgTCGCTCCCCCCGCGGGGGCGTGGATCGAAAC.CCGCCACG.TCTTGGAACAGCGTGCC.ATTGGCCCCGGTCGCTCCCCCCGCGGGGGCGTGGATCGaca                                                                    >  1:4196142/4‑99 (MQ=2)
                                    ggcgTCGCTCCCCCCGCGGGGGCGTGGATCGAAAC.CCGCCACG.TCTTGGAACAGCGTGCC.ATTGGCCCCGGTCGCTCCCCCCGCGGGGGCGTGGATCGaaa                                                                    >  1:1645675/4‑101 (MQ=0)
                                    ggcgTCGCTCCCCCCGCGGGGGCGTGGATCGAAAC.CCGCCACG.TCTTGGAACAGCGTGCC.ATTGGCCCCGGTCGCTCCCCCCGCGGGGGCGTGGATCGaaa                                                                    >  2:1618386/4‑101 (MQ=0)
                                    gTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCAGCTCG.TCGTTGGTCAGCCGGAC.GATGTCC...GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACt                                                                  <  1:3069186/101‑1 (MQ=255)
                                    gTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCAGCTCG.TCGTTGGTCAGCCGGAC.GATGTCC...GTCCCTCCCCCCGCGGGGGCGCGGGCCGaaagt                                                                  >  1:3958755/1‑99 (MQ=37)
                                    ccggTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCAGCTTGACCGCATCGAAGCTGGTCACATCGCC...GTCGCTCCCCCCGCGGGGGCGTGGATCGaaa                                                                    <  2:1112195/99‑1 (MQ=2)
                                     cggTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCAGCTTGACCGCATCGAAGCTGGTCACATCGCC...GTCGCTCCCCCCGCGGGGGCGTGGATCGAAAc                                                                   <  2:3104666/100‑1 (MQ=2)
                                      ggTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCAGCTTGACCGCATCGAAGCTGGTCACATCGCC...GTCGCTCCCCCCGCGGGGGCGTGGATCGAAAcg                                                                  <  1:1337728/101‑2 (MQ=2)
                                      ggTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCAGCTTGACCGCATCGAAGCTGGTCACATCGCC...GTCGCTCCCCCCGCGGGGGCGTGGATCGAAAcg                                                                  >  2:2272359/1‑100 (MQ=2)
                                          gCTCCCCCCGCGGGGGCGTGGATCGAAACATCAGCTTGACCGCATCGAAGCTGGTCACATCTCC...GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGGcgt                                                              <  2:2527389/101‑3 (MQ=2)
                                          gCTCCCCCCGCGGGGGCGTGGATCGAAACATCAGCTTGACCGCATCGAAGCTGGTCACATCGCC...GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGGcgt                                                              >  1:1424694/1‑99 (MQ=2)
                                          gCTCCCCCCGCGGGGGCGTGGATCGAAAC.CCGCCACG.TCTTGGAACAGCGTGCC.ATTGGCCCCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAAcaactg                                                              >  2:2378021/1‑96 (MQ=2)
                                           ccccccccGCGGGGGCGTGGATCGAAACATCAGCTTGCCCGCATCGAAGCTGGTCACATCGCC...GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGGcgta                                                             <  1:3622218/99‑4 (MQ=2)
                                           cTCCCCCCGCGGGGGCGTGGATCGAAACATCAGCTTGACCGCATCGAAGCTGGTCACATCGCC...GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGGcgta                                                             >  2:4025181/1‑98 (MQ=2)
                                           cTCCCCCCGCGGGGGCGTGGATCGAAACATCAGCTTGACCGCATCGAAGCTGGTCACATCGCC...GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGGcgta                                                             >  2:1805460/1‑98 (MQ=2)
                                           cTCCCCCCGCGGGGGCGTGGATCGAAAC.CCGCCACG.TCTTGGAACAGCGTGCC.ATTGGCCCCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAAcaactga                                                             >  2:2933242/1‑95 (MQ=2)
                                            cccccccGCGGGGGCGTGGATCGAAACATCAGCTCG.TCGTTGGTCAGCCGGAC.GATGTCC...GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCCGTCtt                                                          <  1:3774031/100‑1 (MQ=255)
                                             ccccccGCGGGGGCGTGGATCGAAACATCAGCTTGACCGCATCGAAGCTGGTCACATCGCC...GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGGcgtaga                                                           <  1:3804552/101‑6 (MQ=2)
                                             ccccccGCGGGGGCGTGGATCGAAACATCAGCTCG.TCGTTGGTCAGCCGGAC.GATGTCC...GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCCGTCTTg                                                         >  2:2846190/1‑101 (MQ=255)
                                              cccccGCGGGGGCGTGGATCGAAACATCAGCTTGACCGCATCGAAGCTGGTCACATCGCC...GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGGcgtagat                                                          <  1:3683906/101‑7 (MQ=2)
                                              cccccGCGGGGGCGTGGATCGAAACATCAGCTCG.TCGTTGGTCAGCCGGAC.GATGTCC...GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCCGTCTTGa                                                        <  1:933071/101‑1 (MQ=255)
                                               ccccGCGGGGGCGTGGATCGAAACATCAGCTTGACCGCATCGAAGCTGGTCACATCGCC...GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGGcgtagatc                                                         >  1:2367063/1‑94 (MQ=2)
                                               ccccGCGGGGGCGTGGATCGAAACATCAGCTCG.TCGTTGGTCAGCCGGAC.GATGTCC...GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCCGTCTTGAg                                                       <  2:3319454/101‑1 (MQ=255)
                                                cccGCGGGGGCGTGGATCGAAACATCAGCTTGACCGCATCGAAGCTGGTCACATCGCC...GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGGcgtagatct                                                        <  1:2307054/101‑9 (MQ=2)
                                                cccGCGGGGGCGTGGATCGAAACATCAGCTTGACCGCATCGAAGCTGGTCACATCGCC...GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGGcgtagatcg                                                        >  1:3109828/1‑93 (MQ=2)
                                                    cgGGGGCGTGGATCGAAACATCAGCTCG.TCGTTGGTCAGCCGGAC.GATGTCC...GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCCGTCTTGAGTGATg                                                  >  2:65895/1‑101 (MQ=255)
                                                     gggggCGTGGATCGAAACATCAGCTCG.TCGTTGGTCAGCCGGAC.GATGTGC...GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCCGTCTTGAGTGATgg                                                 >  1:3578817/1‑101 (MQ=255)
                                                       gggCGTGGATCGAAACATCAGCTCG.TCGTTGGTCAGCCGGAC.GATGTCC...GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCCGTCTTGAGTGATgggg                                               <  2:1586523/101‑1 (MQ=255)
                                                         gCGTGGATCGAAACATCAGCTCG.TCGTTGGTCAGCCGGAC.GATGTCC...GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCCGTCTTGAGTGATGGGgcg                                             <  2:3644037/101‑1 (MQ=255)
                                                           gTGGATCGAAACATCAGCTCG.TCGTTGGTCAGCCGGAC.GATGTCC...GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCCGTCTTGAGTGATGGGGcgcg                                           >  2:2147484/1‑101 (MQ=255)
                                                           gTGGATCGAAACATCAGCTCG.TCGTTGGTCAGCCGGAC.GATGTCC...GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCCGTCTTGAGTGATGGGGcgcg                                           <  1:639678/101‑1 (MQ=255)
                                                             ggATCGAAACATCAGCTCG.TCGTTGGTCAGCCGGAC.GATGTCC...GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCCGTCTTGAGTGATGGGGCGCGAg                                         >  1:2053963/1‑101 (MQ=255)
                                                                 cGAAACATCAGCTCG.TCGTTGGTCAGCCGGAC.GATGTCC...GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCCGTCTTGAGTGATGGGGCGCGAGcaca                                     <  1:1681646/101‑1 (MQ=255)
                                                                   aaaCATCAGCTCG.TCGTTGGTCAGCCGGAC.GATGTCC...GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCCGTCTTGAGTGATGGGGCGCGAGCACATg                                   >  1:2018896/1‑101 (MQ=255)
                                                                   aaaCATCAGCTCG.TCGTTGGTCAGCCGGAC.GATGTCC...GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCCGTCTTGAGTGATGGGGCGCGAGCACATg                                   >  1:2928272/1‑101 (MQ=255)
                                                                   aaaCATCAGCTCG.TCGTTGGTCAGCCGGAC.GATGTCC...GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCCGTCTTGAGTGATGGGGCGCGAGCACATg                                   >  2:833732/1‑101 (MQ=255)
                                                                   aaaCATCAGCTCG.TCGTTGGTCAGCCGGAC.GATGTCC...GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCCGTCTTGAGTGATGGGGCGCGAGCACATg                                   >  2:2645615/1‑101 (MQ=255)
                                                                      cATCAGCTCG.TCGTTGGTCAGCCGGAC.GATGTCC...GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCCGTCTTGAGTGATGGGGCGCGAGCACATgtcg                                >  1:3308422/1‑101 (MQ=255)
                                                                      aaTCAGCTCG.TCCTTGATCAGCCGGAC.GATGACC...GTCAGAGCCCCCGCGGGGGCGTAGATCGAAACTGCCGTCTTGAGTGATGGGGCGCGAGCACATgtcg                                <  1:1391324/100‑1 (MQ=255)
                                                                          aGCTCG.TCGTTGGTCAGCCGGAC.GATGTCC...GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCCGTCTTGAGTGATGGGGCGCGAGCACATGTCGTCGc                            <  2:1045144/101‑1 (MQ=255)
                                                                            cTCG.TCGTTGGTCAGCCGGAC.GATGTCC...GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCCGTCTTGAGTGATGGGGCGCGAGCACATGTCGTCGCTc                          <  1:1591305/101‑1 (MQ=255)
                                                                             tcg.tcgTTGGTCAGCCGGAC.GATGTCC...GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCCGTCTTGAGTGATGGGGCGCGAGCACATGTCGTCGCTcc                         >  2:4231140/1‑101 (MQ=255)
                                                                              cg.tcgtTGGTCAGCCGGAC.GATGTCC...GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCCGTCTTGAGTGATGGGGCGCGAGCACATGTCGTCGCTccc                        >  1:3356190/1‑101 (MQ=255)
                                                                               g.tcgtTGGTCAGCCGGAC.GATGTCC...GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCCGTCTTGAGTGATGGGGCGCGAGCACATGTCGTCGCTcccc                       <  2:2928050/101‑1 (MQ=255)
                                                                                 tcgtTGGTCAGCCGGAC.GATGTCC...GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCCGTCTTGAGTGATGGGGCGCGAGCACATGTCGTCGCTccccc                      <  2:391081/101‑1 (MQ=255)
                                                                                   gtTTGTCAGCCGGAC.GATGTCC...GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCCGTCTTGAGTGATGGGGCGCGAGCACATGTCGTCGCTCCCCCcg                    <  1:2217244/101‑1 (MQ=255)
                                                                                   gtTGGTCAGCCGGAC.GATGTCC...GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCCGTCTTGAGTGATGGGGCGCGAGCACATGTCGTCGCTCCCCCcg                    >  2:4114158/1‑101 (MQ=255)
                                                                                    ttGGTCAGCCGGAC.GATGTCC...GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCCGTCTTGAGTGATGGGGCGCGAGCACATGTCGTCGCTCCCCCcgc                   >  2:2472911/1‑101 (MQ=255)
                                                                                       gTCAGCGGGGC.GATGTCC...GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGGCATCTCCAGTAATGGCGCGCGCACACGTGTCGTCGCTCCCTCcgcgcg                >  2:4208026/1‑99 (MQ=16)
                                                                                        tCAGCCGGAC.GATGTCC...GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCCGTCTTGAGTGATGGGGCGCGAGCACATGTCGTCGCTCCCCCCGCgggg               >  1:3873286/1‑101 (MQ=255)
                                                                                         cAGCCGGAC.GATGTCC...GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCCGTCTTGAGTGATGGGGCGCGAGCACATGTCGTCGCTCCCCCCGCggggg              <  2:1840585/101‑1 (MQ=255)
                                                                                          aGCCGGAC.GATGTCC...GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCCGTCTTGAGTGATGGGGCGCGAGCACATGTCGTCGCTCCCCCCGCGGGGGc             >  1:1908003/1‑101 (MQ=255)
                                                                                          aGCCGGAC.GATGTCC...GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCCGTCTTGAGTGATGGGGCGCGAGCACATGTCGTCGCTCCCCCCGCCGGGGc             >  1:822985/1‑101 (MQ=255)
                                                                                             cGGAC.GATGTCC...GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCCGTCTTGAGTGATGGGGCGCGAGCACATGTCGTCGCTCCCCCCGCGGGGGCGTg          <  1:4231402/101‑1 (MQ=255)
                                                                                              ggAC.GATGTCC...GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCCGTCTTGAGTGATGGGGCGCGAGCACATGTCGTCGCTCCCCCCGCGGGGGCGTgg         >  2:1573271/1‑101 (MQ=255)
                                                                                               gAC.GATGTCC...GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCCGTCTTGAGTGATGGGGCGCGAGCACATGTCGTCGCTCCCCCCGCGGGGGCGTGGa        <  2:1072026/101‑1 (MQ=255)
                                                                                                aC.GATGTCC...GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCCGTCTTGAGTGATGGGGCGCGAGCACATGTCGTCGCTCCCCCCGCGGGGGCGTGGAt       <  1:3900099/101‑1 (MQ=255)
                                                                                                aC.GATGTCC...GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCCGTCTTGAGTGATGGGGCGCGAGCACATGTCGTCGCTCCCCCCGCGGGGGCGTGGAt       <  1:571925/101‑1 (MQ=255)
                                                                                                   gATGTCC...GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCCGTCTTGAGTGATGGGGCGCGAGCACATGTCGTCGCTCCCCCCGCGGGGGCGTGGATCg     >  1:1525658/1‑101 (MQ=255)
                                                                                                   gATGTCC...GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCCGTCTTGAGTGATGGGGCGCGAGCACATGTCGTCGCTCCCCCCGCGGGGGCGTGGATCg     >  2:2213322/1‑101 (MQ=255)
                                                                                                   gATGTCC...GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCCGTCTTGAGTGATGGGGCGCGAGCACATGTCGTCGCTCCCCCCGCGGGGGCGTGGATCg     <  1:1865235/101‑1 (MQ=255)
                                                                                                   gATGTCC...GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCCGTCTTGAGTGATGGGGCGCGAGCACATGTCGTCGCTCCCCCCGCGGGGGCGTGGATCg     >  2:2175096/1‑101 (MQ=255)
                                                                                                      gTCC...GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCCGTCTTGAGTGATGGGGCGCGAGCACATGTCGTCGCTCCCCCCGCGGGGGGGTGGATCGaaa  >  1:3343433/1‑101 (MQ=255)
                                                                                                      |                                                                                                       
GAAACAGCCTGCGCGGCCGCATCGGCTTCCGCCTGCGTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCAGCTCG.TCGTTGGTCAGCCGGAC.GATGTCC...GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCCGTCTTGAGTGATGGGGCGCGAGCACATGTCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAA  >  NC_014034/1396930‑1397130

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