Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_014034 221,958 A→G 4.3% D233D (GAT→GAC RCAP_rcc00185 ← hypothetical protein

Read alignment evidence...
  seq id position change freq score reads annotation genes product
*NC_014034221,9580A→G 4.3% 40.3 ‑4.4 ‑1.4559D233D (GAT→GACRCAP_rcc00185hypothetical protein
Reads supporting (aligned to +/- strand):  new base (19/6):  ref base (264/270):  total (283/276)
Fisher's exact test strand distribution p-value = 1.27e-02
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 3.74e-02

Only 100 of 1180 total aligned reads displayed.

CGCCGGTGCCCGACACCGTCGCCGAGGCTTCGGTGTTGGTGGTCAGCGTCGCCGGATCGCGCTCATCGACGGTGGGCGTCGGATGCGTGGGCTCGGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGT  >  NC_014034/221858‑222058
                                                                                                    |                                                                                                    
  tccgtgggatcggtcggatccgtgggatcggtcgggtccgtgggatcggtcggatccgtgggatcggtcgggtccggggggtcggtggggtccgtgggatc                                                                                                    <  2:886430‑M1/5‑1 (MQ=255)
      tGCCCGACACCGTCGCCGAGGCTTCGGTGTTGGTGGTCAGCGTCGCCGGATCGCGCTCATCGACGGTGGGCGTCGGATGCGTGGGCTCGGTCggatcggtc                                                                                                >  2:1217782/1‑101 (MQ=255)
      tGCCCGACACCGTCGCCGAGGCTTCGGTGTTGGTGGTCAGCGTCGCCGGATCGCGCTCATCGACGGTGGGCGTCGGATGCGTGGGCTCGGTCggatcggtc                                                                                                <  1:3958424/101‑1 (MQ=255)
        ggatcggtcggatccgtgggatcggtcgggtccgtgggatcggtcggatccgtgggatcggtcggatccgtgggatcggtcggatccgtgggatcggtcgg                                                                                              >  2:2264653‑M1/91‑101 (MQ=255)
          gtccgtcggatccgtgggatcggtcggatccgtgggatcggtcggatccgtgggatcggtcgggtccgtgggatcggtcggatccgtgggatcggtcgggt                                                                                            <  1:3475728‑M1/13‑1 (MQ=255)
          cGACACCGGCGCCGAGGCTTCGGTGTTGGTGGCCAGCGTCGGCGGATCGCGCTCATCGACGGTGGGCGTCGGATGCGTGGGCTCGGTCggatcggtcgggt                                                                                            <  1:3749904/101‑1 (MQ=255)
            cggtcggatccgtgggatcggtcgggtccgtgggatcggtcgggtccgtgggatcggtggggtccgtgggatcggtcgggtccgtgggatcggtcggctcg                                                                                          <  1:571998‑M1/15‑5 (MQ=255)
             cgtcggatccgtgggatcggtcggatccgtgggatcggtcggatccgtgggatcggtcgggtccgtgggatcggtcggatccgtgggatcggtcgggtccg                                                                                         >  2:3641978‑M1/86‑101 (MQ=255)
              gtcggatccgtgggatcggtcgggtccgggggatcggtcgggtccgtgggatcggtcgggtccgtggcatcggtcgggtccgtgggatcggtcgggtccgt                                                                                        <  1:2796018‑M1/17‑1 (MQ=255)
              gtcggatccgtgggatcggtcggatccgtgggatcggtcggatccgtgggatcggtcgggtccgtgggatcggtcggatccgtgggatcggtcgggtccgt                                                                                        >  1:726129‑M1/85‑101 (MQ=255)
              gtcggatccgtgggatcggtcggatccgtgggatcggtcggatccgtgggatcggtcgggtccgtgggatcggtcggatccgtgggatcggtcgggtccgt                                                                                        <  1:124666‑M1/17‑1 (MQ=255)
                 gTCGCCGAGGCTTCGGTGTTGGTGGTCAGCGTCGCCGGATCGCGCTCATCGACGGTGGGCGTCGGATGCGTGGGCTCGGTCggatcggtcgggtccgtggg                                                                                     <  1:4138812/101‑1 (MQ=255)
                   gtccgtgggatcggtccgatccgtgggatcggtcgggtccgtgggatcggtcgggtccgtgggatcggtcgggtccgtgggatcggtcgggtccgtgggat                                                                                   >  2:1356802‑M1/80‑101 (MQ=255)
                   ggccgtggggtggggcgggtccgtgggatcggtcggatccgtgggatcggtcgggtccgtgggatcggtcgggtccgtgggatcggtcgggtccgtgggat                                                                                   >  1:633025‑M1/80‑101 (MQ=255)
                    tctgtgggagcggtcgggtccgtgggatcggtcgggtccgtgggatcggtcgggtccgtgggatcggtcgggtccgtgggatcggtcgggtccgtgggatc                                                                                  >  2:2082558‑M1/79‑101 (MQ=255)
                    gCCGAGGCTTCGGTGTTGGTGGTGAGCGTCGCCGGATCGCGCTCATCGACGGTGGGCGTCGGATGCGTGGGCTCGGTCggatcggtcgggtccgtgggatc                                                                                  <  1:1389327/101‑1 (MQ=255)
                         gggatcggtcgggtccgtgggatcggtcgggtccgtgggatcggtcgggtccgtgggatcggtcgggtccgtgggatcggtcgggtccgtggggtccgtcg                                                                             <  2:3544811‑M1/28‑1 (MQ=255)
                         gggatcggtcgggtccgtcggatcggtgggatccgtgggatcggtcgggtccgtgggatcggtcgggtccgtcggatcggtcgggtccgtgggatcggtcg                                                                             >  1:4248120‑M1/74‑101 (MQ=255)
                           gatcggtcgggtccgtgggatcggtcgggtccgtgggatcggtcgggtccgtgggatcggtcgggtccgtgggatcggtcggttccgtgggatcgggcggt                                                                           <  2:4138366‑M1/30‑2 (MQ=255)
                             tcggtcggatccgtgggatcggtcgggtccgtgggatcggtcgggtccgtgggatcggtcgggtccgtgggatcggtcgggtccgtgggctcggtcgggtc                                                                         <  1:4067896‑M1/32‑1 (MQ=255)
                              cggtcgggtccgtgggatcggtcggatccgtgggatcggtcgggtccgtgggatcggtcggatccgtgggatcggtcggatccgtgggatcggtcggatcc                                                                        >  1:2468478‑M1/69‑101 (MQ=255)
                              cggtcgggtccgtgggatcggtcggatccgtgggatcggtcgggtccgtgggatcggtcggatccgtgggatcggtcggatccgtgggatcggtcggatcc                                                                        >  2:2572851‑M1/69‑101 (MQ=255)
                               ggtcgggtccgtggggtccgtcggatccgtgggatcggtcggatccgtgggatcggtcgggtccgtgggatcggtcgggtccgtgggatcggtcgggtccg                                                                       <  2:3351831‑M1/34‑1 (MQ=255)
                                  gtTGGTGGTCAGCGTCGCCGGATCGCGCTCATCGACGGTGGGCGTCGGATGCGTGGGCTCGGTCGgatcggtcgggtccgtgggatcggtcgggtccgtgg                                                                    >  1:2030326/1‑101 (MQ=255)
                                  ggggtccgtgggatcggtcgggtccgtgggatcggtcggatccgtgggatcggtcgggtccgtgggatcggtcgggtccgtgggatcggtcgggtccgtgg                                                                    >  1:904153‑M1/65‑101 (MQ=255)
                                  gggatccgtgggatcggtcgggtccgtgggatcggtcgggtccgtcggatcggtcgggtccgtgggatcggtcgggtccgtgggatcggtcgggtccgtgg                                                                    <  1:4309991‑M1/37‑1 (MQ=255)
                                     ggtggtCAGCGTCGCCGGATCGCGCTCATCGACGGTGGGCGTCGGATGCGTGGGCTCGGTCGGATcggtcgggtccgtgggatcggtcgggtccgtgggat                                                                 <  1:4515745/101‑1 (MQ=255)
                                      tccgtaggatcggtgggatccgtgggatcggtcgggtccgtgggatcggtcgggtccgtcggatcggtcgggtccgtgggatcggtcgggtccgtgggatc                                                                >  1:795188‑M1/61‑101 (MQ=255)
                                        cgtgggatcggtcggatccgtgggatcggtcgggtccgtgggatcgggcgggtccgtgggatcgggcgggtccgtgggatcggtcgggtccgtgggattgg                                                              <  2:3202744‑M1/43‑4 (MQ=255)
                                                  gtcggatccgtgggatcggtcgggtccgtgggatcggtcggatccgtgGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGt                                                    >  2:4067218‑M1/49‑101 (MQ=255)
                                                    cgggtccgtgggatcggtcgggtccgtgggatcggtcgggtccgtgGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTgg                                                  <  2:1669997‑M1/55‑1 (MQ=255)
                                                    cGGATCGCGCTCATCGACGGTGGGCGTCGGATGCGTGGGCTCGGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCTGTGGGATCCGTgg                                                  >  2:201296/1‑101 (MQ=255)
                                                    cGGATCGCGCTCATCGACGGTGGGCGTCGGATGCGTGGGCTCGGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTgg                                                  <  1:3819727/101‑1 (MQ=255)
                                                      ggtccgtgggatcggtcgggtccgtgggatcggtcgggtccgtgGGATCGGTCGGGTCCGTGGGGTCCGTCGGATCCGTGGGATCGGTCGGATCCGTGGGa                                                <  1:1196788‑M1/57‑1 (MQ=255)
                                                       gtccgtgggatcggtcgggtccgtgggatcggtcgggtccgtgGGATCGGTCGGGTCCGTGGGGTCCGTCGGATCCGTGGGATCGGTCGGATCCGTGGgaa                                               <  2:1102170‑M1/58‑2 (MQ=255)
                                                          cgtgggatcggtcgggtccgtgggatcggtcgggtccgtgGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGGTCCGTCGGATCCGTGGGCtcgg                                            <  1:2584018‑M1/61‑1 (MQ=255)
                                                          cgtgggatcggtcgggtccgtgggatcggtcgggtccgtgGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTcggatccgtagatcgga                                            <  2:1590721‑M1/61‑9 (MQ=255)
                                                          cgtgggatcggtcgggtccgtgggatcggtcggatccgtgGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGCGGGAtcgg                                            <  2:1017996‑M1/61‑1 (MQ=255)
                                                           gtgggatcggtcgggtccgtgggatcggtcgggtccgtgGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGatccgtggggtccgt                                           >  1:4435381‑M1/40‑95 (MQ=255)
                                                             gggatcggtcgggtccgtgggatcggtcgggtccgtgGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGGTCGGTCGGGTCCGTGGGAtcggtcg                                         <  2:1871703‑M1/64‑1 (MQ=255)
                                                             cTCATCGACGGTGGGCGTCGGATGCGTGGGCTCGGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGAtcggtcg                                         <  2:2462851/101‑1 (MQ=255)
                                                               cATCGACGGTGGGCGTCGGATGCGTGGGCTCGGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCggg                                       >  2:3711007/1‑101 (MQ=255)
                                                                 tcggtcgggtccgtggggtcggtcgggtccgtgGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTc                                     >  2:1775029‑M1/34‑101 (MQ=255)
                                                                 tCGACGGTGGGCGTCGGATGCGTGGGCTCGGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTc                                     >  1:2264815/1‑101 (MQ=255)
                                                                  cggtcggttccgtgggatcggtcggatccgtgGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTcc                                    >  2:3797830‑M1/33‑101 (MQ=255)
                                                                  cggtcgggtccgtgggatcggtcgggtccgtgGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGGTCCGTCGGATCCGTGGGATCGGTCGGATcc                                    >  2:4427076‑M1/33‑101 (MQ=255)
                                                                     tcgggtccgtgggatcggtcgggtccgtgGGGTCCGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTg                                 >  1:3394024‑M1/30‑101 (MQ=255)
                                                                     tcgggtccgtgggatcggtcggatccgtgGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTg                                 <  2:4502122‑M1/72‑1 (MQ=255)
                                                                       gggtccgtgggatcggtcggatccgtgGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGCTCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTggg                               <  1:4657873‑M1/74‑1 (MQ=255)
                                                                       ggatccgtgggatcggtcggatccgtgGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTggg                               <  2:2725255‑M1/74‑1 (MQ=255)
                                                                         gttcgtgggatcggtcgggtccgtgGGGTCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGAt                             >  2:3289672‑M1/26‑101 (MQ=255)
                                                                          tccgtgggatcggtcgggtccgtgGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGGTCCGTCGGATCCGTGGGATCGGTCGGATCCGTGGGAtc                            >  1:3545069‑M1/25‑101 (MQ=255)
                                                                          tccgtgggatcggtcggatccgtgGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGAtc                            >  1:1657491‑M1/25‑101 (MQ=255)
                                                                           gCGTCGGATGCGTGGGCTCGGTCGGATCGGGCGGGTCCGTGGGATCGGTGGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGAtcg                           <  2:2157761/101‑1 (MQ=255)
                                                                           ccgtgggatcggtcgggtccgtgGGATCGGTCGGGTCCGTGGGGTCCGTCGGATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGAtcg                           >  1:1968992‑M1/24‑101 (MQ=255)
                                                                           ccgtgggatcggtcggatccgtgGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGGCGGGTCCGTGGGAtcg                           <  1:1111259‑M1/78‑1 (MQ=255)
                                                                            cgtgggatcggtcgggtccgtgGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGGtcgg                          <  2:3969873‑M1/79‑1 (MQ=255)
                                                                             gtgggatcggtcgggtccgtgGGATCGGTCGGGTCCGTGGGGTCCGTCGGATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGatcggg                         <  2:3025502‑M1/80‑2 (MQ=255)
                                                                             gtgggatcggtcgggtccgtcGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTCGGATCGGTGGGATCCGTGGGATCGGTCGGGTCCGTGGGAtcggt                         >  1:1467774‑M1/22‑101 (MQ=255)
                                                                             gTCGGATGCGTGGGCTCGGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGAtcggt                         >  1:3323969/1‑101 (MQ=255)
                                                                              tgggatcggtcgggtccgtgGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGGTCCGTCGGATCCGTGGGATCGGTCGGATCCGTGGGAtcggtc                        <  2:799164‑M1/81‑1 (MQ=255)
                                                                               gggatcggtcgggtccgtgGGATCGGTCGGGTCCGTGGGGTCCGTCGGATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGAtcggtcg                       <  1:1009700‑M1/82‑1 (MQ=255)
                                                                               gggatcggtcgggtccgtgGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGGTCCGTCGGATCCGTGGGATCGGTCGGATCCGTGGGAtcggtcg                       <  1:1583882‑M1/82‑1 (MQ=255)
                                                                                 gatcggtcggggccgtgGGATCGGTCGGGTCCGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGatcggtcggg                     >  1:3107983‑M1/18‑100 (MQ=255)
                                                                                 gATGCGTGGGCTCGGTCGGATCGGTCGGGTCCGTGGGAGCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGa                     <  2:713164/101‑1 (MQ=255)
                                                                                  aTGCGTGGGCTCGGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGAt                    >  1:3022066/1‑101 (MQ=255)
                                                                                   tcggtcggtgccgtgGGATCGGTCGGGTCCGTGGGATCCGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATc                   >  2:4351057‑M1/16‑101 (MQ=255)
                                                                                   tcggtcgggtccgtgGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGGTCCGTCGGATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGATc                   >  1:4586144‑M1/16‑101 (MQ=255)
                                                                                    gCGTGGGCTCGGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATcc                  <  1:4394106/101‑1 (MQ=255)
                                                                                    cggtcgggtccgtgGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGGTCCGTCGGATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGATcc                  <  1:3895357‑M1/87‑1 (MQ=255)
                                                                                    ccGTGGGTTCGGTCTGATCCGTAGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTAGGATCGGTCGGGTCCGTCGGGTCCGTGGGATcg                  <  2:571950/100‑2 (MQ=14)
                                                                                     ggtcgggtccgtgGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGGTCCGTCGGATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGATCCg                 <  2:173136‑M1/88‑1 (MQ=255)
                                                                                      gtcggatccgtgGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGATCcgg                <  1:323552‑M1/89‑2 (MQ=255)
                                                                                       tcgggtccgtgGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTCGGATCGGTGGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCgtc               >  2:346192‑M1/12‑100 (MQ=255)
                                                                                        gggGTCGGTCGGGTCCGTGGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTAGGATCGGTCGGGTCCGTCGGGTCCGTGGGATCGGtcg              <  1:3662204/101‑3 (MQ=17)
                                                                                        gggCTCGGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTgg              >  2:1491998/1‑101 (MQ=255)
                                                                                        cgggtccgtgGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGGTCCGTCGGATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGATCCGTgg              >  2:2952777‑M1/11‑101 (MQ=255)
                                                                                          tccgatctggatcggtgggatccgtgggatcggtcgggtccgtgGGATCGGTCGGGTCCGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGa            <  1:2654211/93‑1 (MQ=11)
                                                                                          ggtccgtcggATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTGGGATCCGTGGGATCGGTCGGGGCCGTGGGATCGGTCGGGGCCGTGGGa            >  2:459549/6‑101 (MQ=11)
                                                                                          gCTCGGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGa            <  2:1262480/101‑1 (MQ=255)
                                                                                           gtgcgtccgATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTCGGATCAGTGGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTCGGAt           <  1:1823649/97‑1 (MQ=11)
                                                                                            gctcggtgattggtggggtccgtgggagcggtcgggtccgtggCGTCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGAtc          <  1:289508/94‑1 (MQ=11)
                                                                                            gccgtcggGTCGGCGGGATCCGGGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGAtc          <  1:2150083/98‑1 (MQ=11)
                                                                                             tcgtgggatcggtcgggcccgtggcatcggtcgggcccgtgGCATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGGtcc         <  2:2816303/96‑2 (MQ=11)
                                                                                             cggtcggATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTCGGATCGGTGGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTCGGAtcg         <  1:1715433/101‑1 (MQ=14)
                                                                                             ccgggggatcggtcggatccgtgggatcggtcggatccgtaGGATCGGTCGGGTCCGTAGGATCGGTCGGGTCCGTCGGATCGGTCGGATCCGTGGGAtcg         <  2:3103215/96‑1 (MQ=11)
                                                                                               gtcggATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGAtcggt       <  1:511233/101‑1 (MQ=255)
                                                                                                tgggatcggtcgggtccgtgggatcggtcgggtccgtgGGATCGGTCGGATCCGTGGGATCGGTGGGGTCCGTGGGATAGGTGGGGTCTGTGGGAtaggag      >  2:4266455/3‑99 (MQ=11)
                                                                                                tgGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGGTCCGTCGGATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGATCCGTGGGAtcggtc      >  2:389867‑M1/3‑101 (MQ=255)
                                                                                                tgGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGgatccgtggggtccgtg      >  2:2213516‑M1/3‑94 (MQ=255)
                                                                                                tcggAGCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGGTCCGTCGGATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGAtcggtc      <  2:2396674/101‑1 (MQ=11)
                                                                                                 gggatcggtcgggtcggtggggtctgtcgggtccgtgGGTTCGGTCGGATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGAtcggtcg     <  2:508319/100‑1 (MQ=11)
                                                                                                 gggatcggtcggatccgtgggatcggtcggatccgtaGGATCGGTCGGGTCCGTAGGATCGGTCGGGTCCGTCGGATCGGTCGGATCCGTGGGAtcggtcg     >  1:2953315/2‑101 (MQ=11)
                                                                                                 gGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGAtcggtcg     <  1:2213678‑M1/100‑1 (MQ=255)
                                                                                                 cggATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTCGGATCGGTGGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTCGGAtcggtcg     <  1:1327296/101‑1 (MQ=11)
                                                                                                 cggATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTCGGAGCGGTGGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTCGGAtcggtcg     <  1:2732976/101‑1 (MQ=11)
                                                                                                  ggatcggtgggatccgtgggatcggtcgggtccgtgGGATCGGTCGGGTCCGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGaagatcgga    >  2:2654026/1‑93 (MQ=11)
                                                                                                  ggaacgggcgggtgcgtggggacggtaggatccgggGGGTCGGCCGGGTCCGTGGGATCGGTCGGGTCCGTGGGGTCCGTCGGATCCGTGGGAtcggtcgg    <  1:2089415/101‑1 (MQ=2)
                                                                                                   gatcggtcgggtccgtgggatcggtcgggtccgtggGGTCCGTCGGATCCGTGGGAAAGGTCGGATCCGTGGGATCGGTAGGTTCCGTGGGATCGGTCggg   >  1:2434410/1‑101 (MQ=11)
                                                                                                    atcggtcggatcggggggatcggtcgggtgcgtgggATCGGTCGGGTCCGTGGGGTCCGTCGGATCCGTGGGATCGGTCGGATCCGTGGGAtcggtcggat  <  1:2340583/101‑3 (MQ=11)
                                                                                                    |                                                                                                    
CGCCGGTGCCCGACACCGTCGCCGAGGCTTCGGTGTTGGTGGTCAGCGTCGCCGGATCGCGCTCATCGACGGTGGGCGTCGGATGCGTGGGCTCGGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGT  >  NC_014034/221858‑222058

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