Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_014034 221,979 G→C 8.9% T226T (ACC→ACG)  RCAP_rcc00185 ← hypothetical protein

Read alignment evidence...
  seq id position change freq score reads annotation genes product
*NC_014034221,9790G→C 8.9% 92.5 ‑0.3 ‑0.6540T226T (ACC→ACG) RCAP_rcc00185hypothetical protein
Reads supporting (aligned to +/- strand):  new base (25/24):  ref base (235/254):  total (261/279)
Fisher's exact test strand distribution p-value = 7.65e-01
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 2.78e-01

Only 100 of 1525 total aligned reads displayed.

CCGAGGCTTCGGTGTTGGTGGTCAGCGTCGCCGGATCGCGCTCATCGACGGTGGGCGTCGGATGCGTGGGCTCGGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTAGGATCGGTCGGGTCCGTCGGGTCC  >  NC_014034/221879‑222087
                                                                                                    |                                                                                                            
   tgggatcggtcgggtccgtggggtccgtcggatccgtgggatcggtcggatccgtgggatcggtcgggtccgtgggatcggtcgggtccgtgggatcggtc                                                                                                           <  2:349286‑M1/27‑1 (MQ=255)
     ggatcggtcgggtccgtgggatcggtcgggtccgtgggatcggtcgggtccgtgggatcggtcggatccgtgggatcggtcgggtccgtgggatcggtcgg                                                                                                         >  1:1219842‑M1/73‑101 (MQ=255)
     gCTTCGGTGTTGGTGGTCAGCGTCGCCGGATCGCGCTCATCGACGGTGGGCGTCGGATGCGTGGGCTCGGTCggatcggtcgggtccgtgggatcggtcgg                                                                                                         >  2:4380582/1‑101 (MQ=255)
     gCTTCGGAGTTGGGGGGCAGCGTCGCCGGATCGCGCTGATCGCCGGGGGGCGTCGGATGCGTGGGCTCGGTCggatcggtcggggccgtgggatcggtcgg                                                                                                         >  2:1624758/1‑101 (MQ=255)
      gatcggtcgggtccgtgggatcggtcgggtccgtgggatcggtcgggtccgtgggatcggtcggatccgtgggatcggtcgggtccgtgggatcggtcggg                                                                                                        >  2:1106692‑M1/72‑101 (MQ=255)
      gatcggtcggatccgtgggatcggtcggatccgtgggatcggtcgggtccgtgggatcggtcggatccgtgggatcggtcgggtccgtgggatcggtcggg                                                                                                        <  1:3444209‑M1/30‑1 (MQ=255)
        tcggtcgggtccgtcggatcggtgggatccgtgggatcggtcgggtccgtgggatcggtcgggtccgtgggatcggtcgggtccgtgggatcgggcgggcc                                                                                                      <  1:1615692‑M1/32‑3 (MQ=255)
         gggttgggtccgtggggtcggtcggatccgtgggatcggtcggatccgtgggatcggtcgggtccgtgggatcggtcgggtccgtgggatcggtcgggtcc                                                                                                     >  1:3995505‑M1/69‑101 (MQ=255)
           gtcgggtccgtgggatcggtcgggtccgtgggatcggtcgggtccgtgggatcggtcggatccgtgggatcggtcgggtccgtgggatcggtcggatccgt                                                                                                   >  2:1395055‑M1/67‑101 (MQ=255)
            tgtTGGTGGTCAGCGTCGCCTGATCGCGCTCATCGACGGTGGGCGTCGGATGCGTGGGCTCGGGCggatcggtcgggtccgtgggatcggtcggatccgtg                                                                                                  <  2:146380/101‑1 (MQ=255)
             gtTGGTGGTCAGCGTCGCCGGATCGCGCTCATCGACGGTGGGCGTCGGATGCGTGGGCTCGGTCGgatcggtcgggtccgtgggatcggtcgggtccgtgg                                                                                                 >  1:2299375/1‑101 (MQ=255)
              gggtccgtgggagcggtcgggtccgtcggatcggtcgggtccgtgggatcggtcgggtccgtgggatcggtcgggtccgtgggatcggtcgggtccgtggg                                                                                                <  1:2427296‑M1/38‑1 (MQ=255)
              ggatccgtgggatcggtcgggtccgtgggatcggtcgggtccgtgggatcggtcgggtccgtgggatcggtcgggtccgtgggatcggtcgggtccggggg                                                                                                <  2:2494217‑M1/38‑1 (MQ=255)
                atccgtgggatcggtcgggtccgtgggatcggtcgggtccgtgggatcggtcgggtccgtgggatcggtcgggtccgtgggatcggtcgggcccgggggat                                                                                              <  2:4114831‑M1/40‑1 (MQ=255)
                 ggggtCAGCGTCGCCGGATCGCGCTCATCGACGGTGGGCGTCGGATGCGTGGGCTCGGTCGGATCggtcgggtccgtgggatcggtcgggtccgtgggatc                                                                                             <  2:3416810/99‑1 (MQ=255)
                  tggtCAGCGTCGCCGGATCGCGCTCATCGACGGTGGGCGTCGGATGCGTGGGCTCGGTCGGATCGgtcgggtccgtgggatcggtcgggtccgtgggatcg                                                                                            <  2:2328442/101‑1 (MQ=255)
                    gggggatcggtcggatccgtgggatcggtcgggtccgtgggatcggtcggatccgggggatcggtcgggtccgtgggatcggtcgggtccgtgggatcggt                                                                                          >  1:3285926‑M1/58‑101 (MQ=255)
                     cttccgatctcggatccgtgggatcggtcgggtccgtgggatcggtcgggtccgtgggatcggtcgggtccgtgggatcggtcgggtccgtgggatcggtc                                                                                         >  2:291534‑M1/57‑101 (MQ=255)
                       aGCGTCGCCGGATCGCGCTCATCGACGGTGGGCGTCGGATGCGTGGGCTCGGTCGGATCGGccgggtccgtgggctcggtcgggtccgcgggatcggcggg                                                                                       >  2:2652984/1‑97 (MQ=255)
                         atcggtcgggtccgtgggatcggtcgggtccgtcggatcggtcgggtccgtgggatcggtcgggtccgtgggatcggtcgggtccgtgggatcggccgggt                                                                                     <  1:4438951‑M1/49‑3 (MQ=255)
                             gtcggatccgtgggatcggtcgggtccgtgggaacggtcggatccgtgGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGt                                                                                 >  1:3815086‑M1/49‑101 (MQ=255)
                                      gtgggatcggtcggatccgtgggatcggtcgggtccgtgGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGAtcggt                                                                        >  1:511576‑M1/40‑101 (MQ=255)
                                       tgggatcggtcgggtccgtgggatcggtcgggtccgtgGGATCGGTCGGGTCCGTGGGGTCCGTCGGATCCGTGGGATCGGTCGGATCCGTGGGAtcggtc                                                                       >  1:2494380‑M1/39‑101 (MQ=255)
                                        gggatcggtcggatccgtgggagcggtcgggtccgtgGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGAtcggtcg                                                                      >  2:731549‑M1/38‑101 (MQ=255)
                                           accGACGGTTGCCCTCGGATGCGTGGCCTCGTCCGGTCCGGTCGGGCCCTTGGAATCGGCCGGG.CCGGTGGATCGGTCGGATCCGTGGGATCGGTCGGGTc                                                                  <  1:2653156/99‑1 (MQ=9)
                                            tCGACGGTGGGCGTCGGATGCGTGGGCTCGGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTc                                                                  <  1:2576961/101‑1 (MQ=255)
                                             cggtcgggtcggtgggatcggtcggatccgtgGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTcc                                                                 >  1:1958859‑M1/33‑101 (MQ=255)
                                              ggtcgggtccgtgggatcggtcggatccgtgGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCg                                                                >  1:1877272‑M1/32‑101 (MQ=255)
                                              gACGGTGGGCGTCGGATGCGTGGGCTCGGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCg                                                                <  2:4059339/101‑1 (MQ=255)
                                                tggggtctgtgggatcggtcgggtcggtgGGATCGGTCGGGTCCGGGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTg                                                              >  1:584891‑M1/30‑101 (MQ=255)
                                                  gTGGGCGTCGGATGCGTGGGCTCGGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTggg                                                            >  1:3352891/1‑101 (MQ=255)
                                                    gtccgtgggatcggtcgggtccgtgGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGAt                                                          <  1:2864132‑M1/76‑1 (MQ=255)
                                                    atccgtgggatcggtcggatccgtgGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCgtggggt                                                          <  2:1919591‑M1/76‑3 (MQ=255)
                                                     tccgtcggagcggtcgggtccgtgGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGAtc                                                         >  2:3123421‑M1/25‑101 (MQ=255)
                                                     ggCGCCGGATGCGTGGGCTCGGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGAtc                                                         <  2:1323347/101‑1 (MQ=255)
                                                       cgtgggatcggtcgggtccgtcGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTCGGATCGGTGGGATCCGTGGGATCGGTCGGGTCCGTGGGAtcgg                                                       >  2:1263656‑M1/23‑101 (MQ=255)
                                                        gtgggatcggtcgggtccgtgGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGAtcggt                                                      <  1:1441598‑M1/80‑1 (MQ=255)
                                                         tCGGATGCGTGGGCTCGGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGAtcggtc                                                     <  1:1116024/101‑1 (MQ=255)
                                                          gggatcggtcgggtccgtgGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCGGTGGGAtcggtcg                                                    <  1:2950617‑M1/82‑1 (MQ=255)
                                                          gggatcggtcggatccgtgGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGAtcggtcg                                                    <  2:3324359‑M1/82‑1 (MQ=255)
                                                           ggATGCGTGGGCTCGGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGCGCGAACGCACGGATCCGAGGGATCGGTCGGGTCCGTGGGAtcggtcgg                                                   >  2:3681979/1‑101 (MQ=18)
                                                            gatcggtcgggtccgtcGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTCGGATCGGTGGGATCCGTGGGATCGGTCGGGTCCGTGGGatcggtcggg                                                  <  1:1355440‑M1/84‑2 (MQ=255)
                                                             atcggtgggatccgtgGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGgatcggtcgggt                                                 <  2:3938319‑M1/85‑3 (MQ=255)
                                                                 gtcgggtccgtcGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCcgg                                             <  2:1841816‑M1/89‑2 (MQ=255)
                                                                 gtcggatccgtgGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGt                                             <  1:3153483‑M1/89‑1 (MQ=255)
                                                                 gTGGGATCGGTCGGGTCCGTCGGATCCGTGGGGTCCGTCGGGTCCGTGGGATCGGTCGGGTCCGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGt                                             >  1:2108375/1‑101 (MQ=1)
                                                                  tcgggtccgtgGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGGTCGGGCGGATccgcg                                            <  1:4326483‑M1/90‑3 (MQ=255)
                                                                   gggATCGGTCGGGTCCGTCGGATCCGTGGGGTCCGTCGGGTCCGTGGGATCGGTCGGGTCCGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGtcg                                           >  1:3651038/1‑99 (MQ=11)
                                                                    gtgtccttgggatcggtcgggtccgtggggtccgtcgggtccgtgGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTggg                                          <  2:86152/92‑1 (MQ=11)
                                                                    gggtccgtgggatcggtcgggcccgtgggatcggtcggggccgtgGGATCGGGCGGGTCAGACGGATCGGACGGGTACGAGGGGTCCATCGAATACGGggg                                          >  1:2231372/10‑101 (MQ=2)
                                                                     gttccgtgggatcgggcgggtccgtgggatcggtcgggtccgtgGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGGTCCGTCGGATCCGTGGGa                                         <  2:557347/93‑1 (MQ=11)
                                                                     gatccgtgGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCTGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGa                                         >  2:2278853‑M1/9‑101 (MQ=255)
                                                                      gtccgtgggatcggtcgggtccgtgggatcggtcgggtccgtgGGGTCCGTCGGATCCGTGGGATCGGTCGGGTAAGTGGGATCGGCGGGGTCCGTGGGAt                                        >  2:2210026/8‑101 (MQ=11)
                                                                        cggtcggGTCCGTCGGATCCGTGGGGTCCGTCGGGTCCGTGGGATCGGTCGGGTCCGTCGGATCGGTCGGGTCCGTGGGATCGGTGGGGGCCGTGGGAtcg                                      >  1:847279/1‑101 (MQ=11)
                                                                        ccgtgGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGGTCCGTCGGATCCGTGGGAtcg                                      >  2:2579650‑M1/6‑101 (MQ=255)
                                                                        ccgtcggATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTCGGATCGGCGGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTCGGAtag                                      >  2:2742604/3‑99 (MQ=11)
                                                                         cgtgGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGGTCCGTCGGATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGATCCGTGGGAtcgg                                     >  1:131580‑M1/5‑101 (MQ=255)
                                                                         cgtcggATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGGGGGATCGGTCGGGTCCGGGGGAtcgg                                     >  2:1739849/2‑101 (MQ=11)
                                                                          gtgggatcggtcgggtccgtgggatcggtcgggtccgtgGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATGGGTAGGGTCCGTGGGagggga                                    >  2:3601999/4‑96 (MQ=11)
                                                                          gtgggatcgggcgggtccgtgggatcggtcgggtccgtgGGGTCCGTCGGATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGATCCGTGGGAtcggt                                    <  1:3959890/98‑1 (MQ=11)
                                                                          gtgGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGGTCCGTCGGATCCGTGGGAtcggt                                    <  1:4123204‑M1/98‑1 (MQ=255)
                                                                          gtcggATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCgaccgatccca                                    <  2:4322508/101‑11 (MQ=255)
                                                                          ggcggATCGGTCGGGTCCGGGGGATCGGGCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGAtcggt                                    <  1:3217596/99‑1 (MQ=255)
                                                                           tggggtccgtcgggtccgtgggatcggtcgggtccgtcggatcgGTCGGGTCCGTGGGATCGGTCGGGGCCGTCGGATCGGTGGGGGCCGTGGGAtcgggg                                   >  2:1999071/9‑99 (MQ=11)
                                                                           tgggatcggtcgggtccgtgggatcggtcgggtccgtgGGGTCCGTCGGATCCGTGGGATGGGTAGGATCGGCGGGCTCTGTCGGCCCCGCGGGGccggtc                                   >  1:1710313/3‑101 (MQ=2)
                                                                           tgggatcggtcgggtccgtcggatcggtcgggtccgtgGGATCGGTCGGGTCCGTCGGATCGGTGGGGTCCGTGGGATCGGTCGGGGCCGTGGGGtcggtc                                   >  2:3871707/3‑101 (MQ=11)
                                                                           tgGGGTCGGTGGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGAtcggtc                                   >  1:3724135‑M1/3‑101 (MQ=255)
                                                                           ggGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGgatccgtggggtccgtg                                   >  1:4373936‑M1/3‑94 (MQ=255)
                                                                            gggatcggtcgggtccgtgggatcggtcgggtccgtgGGATCGGTCGGGTCCGTGGGGTCGGGCGGATCGGTGGGGTCGGGCGGATGCGGGGGGtcgggcg                                  >  2:803499/2‑101 (MQ=11)
                                                                            gGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGGTCCGTCGGATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGATCCGTGGGAtcggtcg                                  <  2:1039697‑M1/100‑1 (MQ=255)
                                                                             gggtccgtcgggtccgtgggatcggtcgggtccgtcggatcgGTCGGGTCCGTGGGATCGGTCGGGTCCGTCGGATCGGTGGGATCCGTGGGAtcggtcgg                                 <  1:917936/95‑1 (MQ=14)
                                                                              gatcggtggtgtgcgtgggatgggtgggatccgtggGATCGGTCGGATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGATCCGTGGGAtcggtcgga                                <  2:2087427/101‑2 (MQ=11)
                                                                               gtcggtcggatccgtggggtcggtcggatccgtgggaTCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGGTCCGTCGGATCCGTGGGAtcggtcggat                               <  1:151645/100‑3 (MQ=11)
                                                                                 cggtcgggtccgggggatcggtcgggtccgtgggatCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGGTCCGTCGGATCCGTGGGATCGGTCGGATcc                             <  1:3081528/101‑1 (MQ=14)
                                                                                 cggccgggtcggtgtggtcggtggggtccggcggctCCGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTcc                             <  1:2356222/101‑1 (MQ=14)
                                                                                  cgtcgggtccgtgggatcggtcgggtccgtcggatcgGTCGGGTCCGTGGGATCGGTCGGGTCCGTCGGATCGGGGGGGTGCGTGGGATCGGGAGGGTCCg                            >  2:94349/2‑101 (MQ=11)
                                                                                   gtcgggtccgtgggatcggtcgggtccgtggggtccGTCGGATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTGGGGTcgga                           >  1:3515457/1‑98 (MQ=11)
                                                                                    tcgggtccgtgggatcggtcgggtccgtgggatcggTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGATCCGTg                          >  2:423797‑M2/1‑100 (MQ=255)
                                                                                    tcggatccgtgggatcggtctgatccgtgggatcggTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGGTCCGTCGGATCCGTGGGATCGGTCGGATCCGtg                          <  2:2823256/101‑2 (MQ=11)
                                                                                    tcggatccgtgggatcggtcggatccgtgggatcggTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGGTCCGTCGGATCCGTGGGATCGGTCGGATCCGtg                          <  1:3878177/101‑2 (MQ=11)
                                                                                     tggggccgtgggatcggtcgggtccgtcgggtcggtgGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTCGGATCGGTCGGGTCCGtgg                         <  1:2027845/100‑3 (MQ=11)
                                                                                     cgggtcctggtgatccgtcggggccgcgggaccggtCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGtgg                         <  2:2418380/101‑3 (MQ=1)
                                                                                     cgggtccgtggggtccgtcggatccgtgggatcggtCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGGCCGtgg                         >  2:3264247/1‑99 (MQ=11)
                                                                                     cgggtccgtgggatcggtcgggtccgtcgggtcggtGGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTCGGATCGGTCGGGTCTGtgg                         <  2:3766932/101‑3 (MQ=11)
                                                                                     cggatccgtgggatcggtcggatccgtgggatcggtCGGGTCCGTGGGGTCGGTCGGGTCCGTGGGGTCCGTCGGATCCGTGGGGTCGGGGGGATCCGtgg                         >  1:1004321/1‑99 (MQ=11)
                                                                                       ggtccggggggtgcgtcggatccgtggggtcggtcgGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGa                       <  1:979490/101‑1 (MQ=11)
                                                                                        gtccgtgggatcggtcgggtccgtgggatcggtcggGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTgggat                      >  1:77440‑M2/1‑96 (MQ=255)
                                                                                        gtccatgggatcggtcggggccgtgggatcggtcggGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGCGTCCGTGGGATCGGTCGGGTCCGTGGGac                      >  2:3988721/1‑100 (MQ=11)
                                                                                         tcggtgggttcggtcggttccgtcggatccgtcggGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGAtc                     <  2:1192136/98‑1 (MQ=11)
                                                                                         tccgtgggatcggtcggatccgtgggatcggtcggGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTgggatc                     >  1:4219789‑M2/1‑95 (MQ=255)
                                                                                         tccgtgggatcggtcggatccgtgggatcggtcggATCCGTAGGATCGGTCGGGTCCGTAGGATCGGTCGGGTCCGTCGGATCGGTCGGATCCGTGGGAtc                     <  1:1788139/101‑1 (MQ=11)
                                                                                             tggggtccgtcggatccgtgggatcggtcggATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGACCGTGGGAtcggtc                 >  2:635182/1‑101 (MQ=11)
                                                                                              tggaacggtcgggtccttgcgatcgttcggGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGAtcggtcg                <  2:3017325/100‑1 (MQ=11)
                                                                                              gggatcggtcgggtccgtggggtccgtcggATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTACGTGGGAtcggtcg                >  2:487622/1‑101 (MQ=11)
                                                                                              gggatcggtcgggtccgtgggatcggtcggGTCCGTGGGATCGGTCGGGTCCGTGGGGTCCGTCGGATCCGTGGGATCGGTCGGATCCGTGGGCtcggtcg                >  1:2764250/1‑101 (MQ=14)
                                                                                              cggctcggtcgggcccgttttgtccgtcggATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGAtcggtcg                <  2:3202219/97‑1 (MQ=11)
                                                                                               ggatcggtcggatccgtgggatcggtcggGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTgggatcggtcgg               >  1:674494‑M2/1‑89 (MQ=255)
                                                                                                 gtcggtcgggtccgtggggtccgtcggATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGAtcggtcggat             <  1:4079928/100‑3 (MQ=11)
                                                                                                 atcggtcgggtccgtcggatcggtgggATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTCGGATCGGTCGGGTCCGGGGGATCGgtcgggc             >  1:1050274/1‑100 (MQ=11)
                                                                                                   cggtcgggtgcggggggtcgggcggATCCGTGGGACCGGGCGGGTCCCGTGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGgtcgggtcc           <  2:1527888/101‑1 (MQ=11)
                                                                                                    |                                                                                                            
CCGAGGCTTCGGTGTTGGTGGTCAGCGTCGCCGGATCGCGCTCATCGACGGTGGGCGTCGGATGCGTGGGCTCGGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTAGGATCGGTCGGGTCCGTCGGGTCC  >  NC_014034/221879‑222087

Base quality scores:  ATCG  < 3 ≤  ATCG  < 12 ≤  ATCG  < 20 ≤  ATCG  < 31 ≤  ATCG  < 40 ≤  ATCG