Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_014034 222,021 G→A 24.8% D212D (GAC→GAT RCAP_rcc00185 ← hypothetical protein

Read alignment evidence...
  seq id position change freq score reads annotation genes product
*NC_014034222,0210G→A24.8% 0.0 ‑3.4 0.01454D212D (GAC→GATRCAP_rcc00185hypothetical protein
Reads supporting (aligned to +/- strand):  new base (149/212):  ref base (522/569):  total (673/781)
Fisher's exact test strand distribution p-value = 3.30e-02
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 1.00e+00

Only 100 of 1790 total aligned reads displayed.

CATCGACGGTGGGCGTCGGATGCGTGGGCTCGGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTAGGATCGGTCGGGTCCGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGT  >  NC_014034/221921‑222121
                                                                                                    |                                                                                                    
  tcggtcgggtccgtgggatcggtcgggtccgtgGGATCGGTCGGGTCCGTGGGGTCCGTCGGATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTc                                                                                                    >  2:4116470‑M1/34‑101 (MQ=255)
  tcggtcgggtccgtgggatcggtcggatccgtgGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGtcggggc                                                                                                    <  2:231927‑M1/68‑3 (MQ=255)
      tggggtctgtgggatcggtcgggtcggtgGGATCGGTCGGGTCCGGGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTg                                                                                                >  1:584891‑M1/30‑101 (MQ=255)
       cgggtccgtgggatcggtcgggtccgtgGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGGTCCGTCGGATCCGTGGGATCGGTCGGATCCGTgg                                                                                               <  2:1013565‑M1/73‑1 (MQ=255)
           tccgtggggtcggtcgggtccgtgGGATCGGTCGGGTCCGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGAtc                                                                                           >  2:1529534‑M1/25‑101 (MQ=255)
             cgtgggatcggtcgggtccgtgGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGGTCCGTCGGATCCGTGGGATCGGTCGGATCCGTGGGAtcgg                                                                                         >  1:2925269‑M1/23‑101 (MQ=255)
             cgtgggatcggtcgggtccgtgGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGAtcgg                                                                                         <  1:4434320‑M1/79‑1 (MQ=255)
              gtgggatcggtcggatccgtgGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGAtcggt                                                                                        <  1:12000‑M1/80‑1 (MQ=255)
              gTCGGATGCGTGGGCTCGGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGAtcggt                                                                                        <  2:2047280/101‑1 (MQ=255)
               tCGGATGCGTGGGCTCGGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGAtcggtc                                                                                       <  1:1116024/101‑1 (MQ=255)
                   atccGTGGGATCGGTCGGGTCCGTCGGATCCGTGGGGTCCGTCGGGTCCGTGGGATCGGTCGGGTCCGTCGGATCGGTCGGGGCCGTGGGAtcggtcgggt                                                                                   >  2:4072354/4‑99 (MQ=11)
                     gcggtggCTCGGTGGGGGCGGTGGGGGGGGTGGGGCCGGGAGGGGCGGTCGGGTCCGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTcc                                                                                 <  1:3050787/96‑1 (MQ=9)
                       gTGGGATCGGTCGGGTCCGTCGGATCCGTGGGGTCCGTCGGGTCCGTGGGATCGGTCGGGTCCGTCGGATCGGTCGGGTCCGGGGGATCGGTCGGGCCCGt                                                                               >  2:2409906/1‑101 (MQ=11)
                         cgggtccgtgGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGGTCCGTCGGATCCGTgg                                                                             <  1:3779145‑M1/91‑1 (MQ=255)
                          ggCTCGGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTggg                                                                            <  2:3822963/101‑1 (MQ=255)
                           ggtccgtgGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGGTCCGTCGGATCCGTGGGa                                                                           <  1:282336‑M1/93‑1 (MQ=255)
                           ggtccgtgGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGGTCCGTCGGATCCGTGGGa                                                                           <  1:2352430‑M1/93‑1 (MQ=255)
                            gtccgtgGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGGTCCGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGAt                                                                          <  2:4325670‑M1/94‑1 (MQ=255)
                             tccttgGGATCGGTCGGCTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGAtc                                                                         >  2:2558805‑M1/7‑101 (MQ=255)
                             tccgtgGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGGTCCGTCGGATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGATCCGTGGGAtc                                                                         <  1:3612132‑M1/95‑1 (MQ=255)
                             tccgtcggATCGGTCGGCTCTGTGGGATCGGTCGGGTCCGTCGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGAtc                                                                         <  2:2139993/98‑1 (MQ=11)
                              cggtcggGTCCGTCGGATCCGTGGGGTCCGTCGGGTCCGTGGGATCGGTCGGGTCCGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGGCGGAtcg                                                                        >  2:2637041/1‑101 (MQ=11)
                               cgtgggatcggtcgggtccgtggggtccgtcggatccgtgGGACCGGTCGGATCCGTGGGGTCCGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGaacag                                                                       >  1:915925/5‑97 (MQ=11)
                               cgtgGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGtccgtggggtctt                                                                       <  2:682130‑M1/97‑6 (MQ=255)
                               cgtcggATCGGTGGGATCCGTGGGGTCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGAtcgg                                                                       <  1:1999200/100‑1 (MQ=11)
                               cgtcggATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTCGGATCGGTGGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTCGGAtcgg                                                                       <  1:3141996/100‑1 (MQ=11)
                                gtggggtccgtcaggtccgtgggatcggtcgggtccgtcggatcgGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGCGGGAGCCGTGGGAtcggc                                                                      >  1:2418537/10‑100 (MQ=11)
                                gtgggatcggtcgggtccgtgggatcggtcgggtccgtgGGATCGGTCGGGTCCGGGGGCGCGGTGGGGCCCGGGGGA..GGTCGGGGCCGTGGGGtcggccg                                                                    >  2:3148337/4‑101 (MQ=2)
                                gtgggatcggtcgggtccgtgggatcggtcgggtccgtgGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATGGGTAGGGTCCGTGGGagggga                                                                      >  2:3601999/4‑96 (MQ=11)
                                 tgggatcggtcgggtccgtgggatcggtcgggtccgtgGGATCGGTGGGATACCTGGGATCGGGGGGGGCGGTGGGATCTGACGGATAGGAGGGAtgtgaa                                                                     >  2:850437/3‑96 (MQ=2)
                                 tgggatcggtcgggtccgtcggatcggtcgggtccgtgGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGTAtcagtc                                                                     >  2:602321/3‑101 (MQ=11)
                                  gggatggggggggtccgtgggatcggtcgggtccgtgggatcggtcGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGGtccgtcg                                                                    <  1:4072614/91‑1 (MQ=11)
                                  gggatcggtcgggtccgtgggatcggtcgggtccgtgGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGGGGGATCGGTCGGGTCGCGGGGAtcggtcg                                                                    >  2:2723461/2‑101 (MQ=11)
                                  cggATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCTTGGGATCGGTCGGATCCGTGGGAtcggtcg                                                                    >  1:318561/1‑101 (MQ=255)
                                    gatcggtcgggtccgtgggatcggtcgggtccgtggGGTCCGTCGGATCCGTGGGATCGGTCGGCTCCGTGGGGTCGGTAGGATCGGTGGGATCGGGCggg                                                                  >  2:3876194/1‑101 (MQ=14)
                                     gtcggtgcggtgcgtgggatcggtggggtccgtgggaTCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGt                                                                 <  1:4181055/100‑1 (MQ=11)
                                     gtccgtcgggtccgtgggatcggtcgggtccgtcggatcgGTCGGGTCCGTGGGATCGGTCGGGTCCGTTGGATCGGGGGGGTGCGTGGGGTCGGGCgggg                                                                 >  1:2081666/5‑100 (MQ=11)
                                     ctcgggcggggcggtggagtcggggggggcggtgggaTCGGTCGGGTCGGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGt                                                                 <  2:2781852/100‑1 (MQ=11)
                                     atcggtcgggtccgtgggatcggtcgggtccgtgggGTCCGTCGGATCCGTGGGATCGGTCGGGGCCGTGGGATCGGTCGGGGCCGTGGGCTCGGCCGGGt                                                                 >  2:3264076/1‑101 (MQ=14)
                                      tcggtcgggtccgtgggatcggtcgggtccggggggGCGGTGGGGTCCGTGGGGTCGGTGGGGTCCGTTGGGTGGGTGGGGTCCGCGTGGtcggaggcgga                                                                >  2:3050566/1‑94 (MQ=0)
                                      ccggtcggggccgtgggggcgatcgggtccctgggccCGGCCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTc                                                                <  2:925465/100‑1 (MQ=11)
                                       cggtcgggtccgtgggatcggtcgggtccgtggggtCCGTCGGATCCGTGGGATCGGTCGGATCCGTGGGCTCGGTCGGGTCCGCGGGATCGGTCGGGGcc                                                               >  2:780194/1‑101 (MQ=11)
                                          tcggatccgtgggatcggtcggatccgtaggatcggTCGGGTCCGTAGGATCGGTCGGGTCCGTCGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGtc                                                            <  1:134543/101‑2 (MQ=11)
                                          tcggatccgtgggatcggtcggatccgtaggatcggTCGGGTCCGTAGGATCGGTCGGGTCCGTCGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGtc                                                            >  2:3319301/1‑100 (MQ=11)
                                           cggatccgtgggatcggtcgggtccgtgggatcggtCGGGTCCGTGGGATCGGTCGGGTCCGTGGGGTCGGTCGGGGCCGTGGGGTCGGTCGGGTCCgcgg                                                           >  2:2249609/1‑98 (MQ=11)
                                            ggatccgtgggatcggtcggatccggaggatcggtcggGTCCGTAGGATCGGTCGGGTCCGTCGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTCgg                                                          <  2:2955295/98‑1 (MQ=11)
                                             ggtccgtgggatcggtcgggtccgtggggtccgtcgGATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGCCCGTGGGa                                                         >  2:3397751/1‑101 (MQ=14)
                                             gatccgtgggatcggtcggatccgtgggatcggtcggATCCGTAGGATCGGTCGGGTCCGTGGGATCGGGCGGGGCCGGGGGATCGGTCGGGTCGGTGggg                                                         >  2:3719193/3‑100 (MQ=11)
                                              gtccgtggcatcggtcgggtccgtgggctcggtcggGTCCCTTGCATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGAt                                                        <  2:2484548/101‑1 (MQ=11)
                                               tccgtgggatcggtcgggtccgtgggatcggtcggGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTgggatc                                                       >  1:3990264‑M2/1‑95 (MQ=255)
                                                ccgtgggatcggtcgggtccgtggggtccgtcggATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTgggatcg                                                      >  1:2496926‑M2/1‑94 (MQ=255)
                                                ccgtgggatcggtcgggtccgtgggatcggtcggGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTggggtcc                                                      >  1:3325490‑M2/1‑94 (MQ=255)
                                                ccgtgggatcggtcgggtccgtcggatcggtgggATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTgggatcg                                                      <  1:1417161‑M2/101‑8 (MQ=255)
                                                ccgtgggatcggtcgggtccgtcggatcggtcggGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTGGGATCCGTGGGATCGGTCGGGTCCGTgggatcg                                                      <  2:4416156‑M2/101‑8 (MQ=255)
                                                ccgtgggatcggtcggatccgtgggatcggtcggGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTgggatcg                                                      <  1:2569474‑M2/101‑8 (MQ=255)
                                                 cgtggggtccgtcgggtccgtgggatcggtcggGTCCGTGGGGTCGGTCGGGTCCGTGGGATCGGTCGGGTCAGTGGGATCGGTGGGGACCggggggcggg                                                     >  2:4236056/1‑92 (MQ=11)
                                                 cgtgggatcggtcgggtccgtgggatcggtcggATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGGGGGATCGGTCGGGTCGGTgggatcgg                                                     <  2:1182068‑M2/101‑9 (MQ=255)
                                                 cgtgggatcggtcggatccgtgggatcggtcggATCCGTAGGATCGGTCGGGTCCGTAGGATCGGTCGGGTCCGTCGGATCGGTCGGATCCGTGGGAtcgg                                                     >  2:1989559/1‑101 (MQ=11)
                                                    gggatcggtcgggtccgtgggatcggtcggATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGATCCGTgggatcggtcg                                                  <  1:4108174‑M2/101‑12 (MQ=255)
                                                    gggagcggtcgggtccgtggggtccgtcggATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGAtcggtcg                                                  <  1:50344/101‑1 (MQ=11)
                                                    cggctcggtcgggcccgttttgtccgtcggATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGAtcggtcg                                                  <  2:3202219/97‑1 (MQ=11)
                                                     ggatcggtcggatccgtgggatcggtcggGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTgggatcggtcgg                                                 >  1:674494‑M2/1‑89 (MQ=255)
                                                      gatcggtcggatccgtgggatcggtcggGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTgggatcggtcgga                                                <  1:1482123‑M2/101‑14 (MQ=255)
                                                       gtaggtcgggtccgtggggttggtcggGGCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGgtcgggt                                               <  1:3413746/98‑1 (MQ=14)
                                                       ggcgggcgggtccgtcggatcggtgggATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTCGGATCGGTCGGGTCCGTGGGATCGgtcgggt                                               <  1:1551528/99‑1 (MQ=11)
                                                       atcggtcgggtccgtggggtccgtcggATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTgggatcggtcggat                                               <  1:435797‑M2/101‑15 (MQ=255)
                                                       atcggtcgggtccgtcggatcggtgggATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTCGGATCGGTCGGGTCCGGGGGATCGgtcgggc                                               >  1:1050274/1‑100 (MQ=11)
                                                        tcggtcgcgtcgctgggctcgggcggGGCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGgtcgggtc                                              <  2:1260469/101‑1 (MQ=11)
                                                         cggtcggatccgtgggatcggtcggATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGGTCCGTCGGCTCCGTGGGATCGgtcggatcc                                             >  2:3948639/1‑101 (MQ=11)
                                                            tcgggtccgtcggatcggtcggGTCCGTGGGATCGGTCGGGTCCGTCGGATCGGTGGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGgtcgggtccgtc                                          <  1:1284586/101‑1 (MQ=11)
                                                              ggatccgtgggatcggtcggATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGATCCGTAGGATTGTTCGGGTCCGTAGGATCGgtcggggccgtcgg                                        >  1:3483795/4‑101 (MQ=14)
                                                               ggtccgtggggtccgtcggATCCGTGGGATCGGTCGGATCCGTGGGATCTGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGgtcgggtccgtgggg                                       >  1:3703656/1‑101 (MQ=11)
                                                                atccgtgggatcggtcggATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGATCCGTAGGATCGGTCGGGTCCGTAGGATCGgtcgggtccgtcggat                                      <  1:1001929/100‑3 (MQ=255)
                                                                atccgtgggatcggtcggATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGATCCGTAGGATCGGTCGGGTCCGTAGGATCGgtcgggtccgtcggat                                      >  2:132870/2‑99 (MQ=255)
                                                                 tccgtggggtccgccggATCAGTGGGATCGGTCGTGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGgtcgggtccgtgggatc                                     <  2:51803/101‑1 (MQ=11)
                                                                   cgtgggatcggtcggATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGATCCGTAGGATCGGTCGGGTCCGTAGGATCGgtcgggtccgtcggatcgg                                   >  2:1726511/1‑99 (MQ=11)
                                                                   cgtgggatcggtcggATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGATCCGTAGGATCGGTCGGGTCCGTAGGATCGgtcgggtccgtcggatcgg                                   >  2:1633837/1‑99 (MQ=11)
                                                                     tgggatcggtcggGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTgggatcggtcgggtccgtgggatcggtc                                 <  2:1466770‑M2/101‑29 (MQ=255)
                                                                     tgggatcggtcggGTCCGTGGGATCGGTCGGGTCCGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTgggatcggtcgggtccgtgggatcggtc                                 <  1:2166776‑M2/101‑29 (MQ=255)
                                                                     taggatcggtcggGTCCGTAGGATCGGTCGGGTCCGTCGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTCGGATCCGTGGGGTCCGTCGGGTCCGTg                                 >  1:970207/3‑101 (MQ=14)
                                                                          tccgtgggATCGGTCGGATCCGTGGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTAGGATCGGTCGGGTCCGTCGGGTCCGTGGGAtc                            <  1:2666295/95‑1 (MQ=14)
                                                                            cgtcggATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGGTCCGTCGGAtccg                          <  2:1805137/100‑3 (MQ=11)
                                                                             gtcgggggcGTGGGATGGGTCGGGGCCGTGGGGTCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGGTCCGTGGGAtcggt                         <  1:3876448/93‑1 (MQ=11)
                                                                             gtcggATCCGTGGGGTCCGTCGGGTCCGTGGGATCGGTCGGGTCCGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTCGGATCGGTGGGAtccgt                         <  1:1794868/101‑1 (MQ=11)
                                                                               cggATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGATCCGTgggatcggtggggtccgtgggatcggtggggtccgtgg                       <  2:3995249‑M2/101‑39 (MQ=255)
                                                                                ggATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGGGCCGTCGGAtccgtggg                      >  1:756352/1‑101 (MQ=11)
                                                                                 ggtCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGCCGGGTCCGTGGGGAGCGTCGGAtcagtagaa                     >  1:2743936/3‑97 (MQ=11)
                                                                                  gtcCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTggggtccgtcgggtccgtgggatcggtcgggtccgtgggat                    <  1:4276627‑M2/100‑42 (MQ=255)
                                                                                   tcggtcggATCCGTGGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTAGGATCGGTCGGGTCCGTCGGGTCCGTGGGATCGGTCGGATc                   <  2:1089327/95‑1 (MQ=17)
                                                                                      gtaggATCGGTCGGGTCCGTCGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTCGGATCCGTGGGGTCCGTCGGGTCCGTGGGATCGGTCGGGTCCGt                <  1:1726622/98‑1 (MQ=14)
                                                                                      gTGGGATCGGTCGGGTCCGTGGGGTCCGTCGGATCCGTGGGATCGGTCGGATCCGTgggatcggtcgggtccgtgggatcggtcgggtccgtgggatcggc                <  1:3261733‑M2/101‑46 (MQ=255)
                                                                                      gTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGGTCCGTGGGATCGGTCGGGTCCGt                >  1:190448/1‑101 (MQ=14)
                                                                                        gggATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATTGGTCGGGTCCGTgggatcggtcggatccgtgggatcggtcggatccgtgggatcggtcg              >  1:2275891‑M2/1‑54 (MQ=255)
                                                                                        gggATCGGTCGGATCCGTGGGATCGGTCGGATCCGTAGGATCGGTCGGGTCCGTAGGATCGGTCGGGTCCGTCGGATCGGTGGGATCCGTGGGATCGGtcg              >  2:2875497/1‑99 (MQ=11)
                                                                                           aTCGGTCGGATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTgggatcggtcggatccgtgggatcggtcgggtccgtgggatcggtcgggt           >  1:3943289‑M2/1‑51 (MQ=255)
                                                                                               gtcggGTCCGTGGGGTCGGTCGGATCCGTGGGATCGGTCGGATCCGTgggatcggtcggatccgtgggatcggtcgggtccgtgggatcggtcgggtccgt       >  2:3483409‑M2/1‑47 (MQ=255)
                                                                                               gtcggGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTgggatcggtcggatccgtgggatcggtcgggtccgtgggatcggtcggatccgt       >  2:1395055‑M2/1‑47 (MQ=255)
                                                                                               gtcggGTCCGTCGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTCGGATCCGTGGGGTCCGTCGGGTCCGTGGGATCGGTCGGGTCCGGGGGAtcggt       >  2:1879785/1‑101 (MQ=14)
                                                                                                  gggTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTgggatcggtcgggtccgtgggatcggtcgggtccgtggggtcggtcggatccgtggg    <  1:839472‑M2/101‑58 (MQ=255)
                                                                                                  gggTCCGTCGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTCGGATCCGTGGGGTCCGTCGGGTCCGTGGGATCGGTCGGGTCCGTCGGAtcggtcgg    >  1:1657094/1‑101 (MQ=17)
                                                                                                    |                                                                                                    
CATCGACGGTGGGCGTCGGATGCGTGGGCTCGGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTAGGATCGGTCGGGTCCGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGT  >  NC_014034/221921‑222121

Base quality scores:  ATCG  < 3 ≤  ATCG  < 15 ≤  ATCG  < 23 ≤  ATCG  < 31 ≤  ATCG  < 40 ≤  ATCG