Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_014034 221,949 C→A 10.9% E236D (GAG→GAT RCAP_rcc00185 ← hypothetical protein

Read alignment evidence...
  seq id position change freq score reads annotation genes product
*NC_014034221,9490C→A10.9% 46.2 ‑2.2 ‑0.5201E236D (GAG→GATRCAP_rcc00185hypothetical protein
Reads supporting (aligned to +/- strand):  new base (15/7):  ref base (87/92):  total (102/99)
Fisher's exact test strand distribution p-value = 1.13e-01
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 3.21e-01

Only 100 of 245 total aligned reads displayed.

AGCCGGTTTCGCCGGTGCCCGACACCGTCGCCGAGGCTTCGGTGTTGGTGGTCAGCGTCGCCGGATCGCGCTCATCGACGGTGGGCGTCGGATGCGTGGGCTCGGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGAT  >  NC_014034/221849‑222049
                                                                                                    |                                                                                                    
aGCCGGTTTCGCCGGTGCCCGACACCGTCGCCGAGGCTTCGGTGTTGGTGGTCAGCGTCGCCGGATCGCGCTCATCGACGGTGGGCGTCGGATGCGTGGGc                                                                                                      >  1:521915/1‑101 (MQ=255)
    ggTTTCGCCGGTGCCCGACACCGTCGCCGAGGCTTCGGTGTTGGTGGTCAGCGTCGCCGGATCGCGCCCATCGACGGTGGGCGTCGGGTGGGTGGGCtcgg                                                                                                  >  1:4343487/1‑101 (MQ=255)
     gTTTCGCCGGTGCCCGACACCGTCGCCGAGGCTTCGGTGTTGGTGGTCAGCGTCGCCGGATCGCGCTCATCGACGGTGGGGGTCGGATGCGTGGGCtcggt                                                                                                 >  1:2277324/1‑101 (MQ=255)
      tttCGCCGGTGCCCGACACCGTCGCCGAGGCTTCGGTGTTGGTGGTCAGCGTCGCCGGATCGCGCTCATCGACGGTGGGCGTCGGATGCGTGGGCtcggtc                                                                                                <  2:1024126/101‑1 (MQ=255)
          gCCGGTGCCCGACACCGTCGCCGAGGCTTCGGTGTTGGTGGTCAGCGTCGCCGGATCGCGCTCATCGACGGTGGGCGTCGGATGCGTGGGCTCGGTCggat                                                                                            >  2:3217961/1‑101 (MQ=255)
          gCCGGTGCCCGACACCGTCGCCGAGGCTTCGGTGTTGGTGGTCAGCGTCGCCGGATCGCGCTCATCGACGGTGGGCGTCGGATGCGTGGGCTCGGTCggat                                                                                            <  1:2008377/101‑1 (MQ=255)
             ggTGCCCGACACCGTCGCCGAGGCTTCGGTGTTGGTGGTCAGCGTCGCCGGATCGCGCTCATCGACGGTGGGCGTCGGATGCGTGGGCTCGGTCggatcgg                                                                                         >  1:1831266/1‑101 (MQ=255)
             ggTGCCCGACACCGTCGCCGAGGCTTCGGTGTTGGTGGTCAGCGTCGCCGGATCGCGCTCATCGACGGTGGGCGTCGGATGCGTGGGCTCGGTCggatcgg                                                                                         >  2:4720356/1‑101 (MQ=255)
              gTGCCCGACACCGTCGCCGAGGCTTCGGTGTTGGTGGTCAGCGTCGCCGGATCGCGCTCATCGACGGTGGGCGTCGGATGCGTGGGCTCGGTCggatcggt                                                                                        >  1:470096/1‑101 (MQ=255)
                gCCCGACACCGTCGCCGAGGCTTCGGGGTTGGTGGTCAGCGTCGCCGGATCGCGCTCATCGACGGTGGGCGTCGGATGCGTGCGCCCGGTCggatcggtcg                                                                                      >  2:1895466/1‑101 (MQ=255)
                    gACACCGTCGCCGAGGCTTCGGTGTTGGTGGTCAGCGTCGCCGGATCGCGCTCATCGACGGTGGGCGTCGGATGCGTGGGCTCGGTCggatcggtcgggtc                                                                                  <  1:3367556/101‑1 (MQ=255)
                     acacCGTCGCCGAGGCTTCGGTGTTGGTGGTCAGCGTCGCCGGATCGCGCTCATCGACGGTGGGCGTCGGATGCGTGGGCTCGGTCggatcggtcgggtcc                                                                                 <  2:664049/101‑1 (MQ=255)
                          gTCGCCGAGGCTTCGGTGTTGGTGGTCAGCGTCGCCGGATCGCGCTCATCGACGGTGGGCGTCGGATGCGTGGGCTCGGTCggatcggtcgggtccgtggg                                                                            <  2:1463070/101‑1 (MQ=255)
                             gCCGAGGCTTCGGTGTTGGTGGTCAGCGTCGCCGGATCGCGCTCATCGACGGTGGGCGTCGGATGCGTGGGCTCGGTCggatcggtcgggtccgtgggatc                                                                         <  2:2754476/101‑1 (MQ=255)
                             gCCGAGGCTTCGGTGTTGGTGGTCAGCGTCGCCGGATCGCGCTCATCGACGGTGGGCGTCGGATGCGTGGGCTCGGTCggatcggtcgggtccgtgggatc                                                                         <  2:4545162/101‑1 (MQ=255)
                                gAGGCTTCGGTGTTGGTGGTCAGCGTCGCCGGATCGCGCTCATCGACGGTGGGCGTCGGATGCGTGGGCTCGGTCggatcggtcgggtccgtgggatcggt                                                                      >  1:2002773/1‑101 (MQ=255)
                                gAGGCTTCGGTGTTGGTGGTCAGCGTCGCCGGATCGCGCTCATCGACGGTGGGCGTCGGATGCGTGGGCTCGGTCggatcggtcgggtccgtgggatcggt                                                                      <  1:2106764/101‑1 (MQ=255)
                                gAGGCTTCGGTGTTGGTGGTCAGCGTCGCCGGATCGCGCTCATCGACGGTGGGCGTCGGATGCGTGGGCTCGGTCggatcggtcgggtccgtgggatcggt                                                                      <  2:2657097/101‑1 (MQ=255)
                                 aGGCTTCGGTGTTGGTGGTCAGCGTCGCCGGATCGCGCTCATCGACGGTGGGCGTCGGATGCGTGGGCTCGGTCggatcggtcgggtccgtgggatcggtc                                                                     <  1:3028577/101‑1 (MQ=255)
                                 aGGCTTCGGTGTTGGTGGTCAGCGTCGCCGGATCGCGCTCATCGACGGTGGGCGTCGGATGCGTGGGCTCGGTCggatcggtcgggtccgtgggatcggtc                                                                     <  2:4012253/101‑1 (MQ=255)
                                   gCTTCGGTGTTGGTGGTCAGCGTCGCCGGATCGCGCTCATCGACGGTGGGCGTCGGATGCGTGGGCTCGGTCggatcggtcgggtccgtgggatcggtcgg                                                                   <  1:88916/101‑1 (MQ=255)
                                     ttCGGTGTTGGTGGTCAGCGTCGCCGGATCGCGCTCATCGACGGTGGGCGTCGGATGCGTGGGCTCGGTCggatcggtcgggtccgtgggatcggtcgggt                                                                 <  1:3797770/101‑1 (MQ=255)
                                      tCGGTGTTGGGTGTCAGCGTCGCCGGATCGCGCTCATCGACGGTGGGCGTCGGATGCGTGGGCTCGGTCggatcgggcgggtccgtgggatcggtcgggtc                                                                <  2:156938/101‑1 (MQ=255)
                                       cGGTGTTGGTGGTCAGCGTCGCCGGATCGCGCTCATCGACGGTGGGCGTCGGATGCGTGGGCTCGGTCggatcggtcgggtccgtgggatcggtcgggtcc                                                               <  1:3016054/101‑1 (MQ=255)
                                       cGGTGTTGGTGGTCAGCGTCGCCGGATCGCGCTCATCGACGGTGGGCGTCGGATGCGTGGGCTCGGTCggatcggtcgggtccgtgggatcggtcgggtcc                                                               <  2:4503379/101‑1 (MQ=255)
                                       cGGTGTTGGTGGTCAGCGTCGCCGGATCGCGCTCATCGACGGTGGGCGTCGGATGCGTGGGCTCGGTCggatcggtcgggtccgtgggatcggtcgggtcc                                                               <  1:2445110/101‑1 (MQ=255)
                                       cGGTGTTGGTGGTCAGCGTCGCCGGATCGCGCTCATCGACGGTGGGCGTCGGATGCGTGGGCTCGGTCggatcggtcgggtccgtgggatcggtcgggtcc                                                               >  1:4428982/1‑101 (MQ=255)
                                        ggTGTTGGTGGTCAGCGTCGCCGGATCGCGCTCATCGACGGTGGGCGTCGGATGCGTGGGCTCGGTCggatcggtcgggtccgtgggatcggtcgggtccg                                                              <  2:2212775/101‑1 (MQ=255)
                                         gtgtTGGTGGTCAGCGTCGCCGGATCGCGCTCATCGACGGTGGGCGTCGGATGCGTGGGCTCGGTCggatcggtcgggtccgtgggatcggtcgggtccgt                                                             <  2:3196840/101‑1 (MQ=255)
                                          tgtTGGTGGTCAGCGTCGCCGGATCGCGCTCATCGACGGTGGGCGTCGGATGCGTGGGCTCGGTCggatcggtcggggccgtgggatcggtcggggccgtg                                                            >  2:4674213/1‑101 (MQ=255)
                                           gtTGGTGGTCAGCGTCGCCGGATCGCGCTCATCGACGGTGGGCGTCGGATGCGTGGGCTCGGTCGgatcggtcgggtccgtgggaccggtcggggccgggg                                                           >  1:691933/1‑101 (MQ=255)
                                           gtTGGTGGTCAGCGTCGCCGGATCGCGCTCATCGACGGTGGGCGTCGGACGCGTGGGCTCGGTCGgatcggtcgggtccgtgggatcggtcgggtccgtgg                                                           >  2:2082855/1‑101 (MQ=255)
                                           ggtggtggTCAGCGTCGCCGGAGCGGGCTCATCGACGGTGGGCGTCGGATGCGTGGGCTCGGTCGgatcggtcgggtccgtgggatcggtcgggtccgtgg                                                           <  2:4060823/99‑1 (MQ=255)
                                            ttGGTGGTCAGCGTCGCCGGATCGCGCTCATCGACGGTGGGCGTCGGATGCGTGGGCTCGGTCGGatcggtcgggtccgtgggatcggtcgggtccgtggg                                                          >  1:4638630/1‑101 (MQ=255)
                                               gtggtCAGCGTCGCCGGATCGCGCTCATCGACGGTGGGCGTCGGATGCGTGGGCTCGGTCGGATCggtcgggtccgtgggatcggtcgggtccgtgggatc                                                       <  2:3697826/101‑1 (MQ=255)
                                                tggtCAGCGTCGCCGGATCGCGCTCATCGACGGTGGGCGTCGGATGCGTGGGCTCGGTCGGATCGgtcgggtccgtgggatcggtcgggtccgtgggatcg                                                      >  2:2607995/1‑101 (MQ=255)
                                                   gcAGCGTCGCCGGATCGCGCTCATCGACGGTGGGCGTCGGATGCGTGGGCTCGGTCGGATCGGTCgggtccgtgggatcggtcgggtccgtgggatcggtc                                                   <  2:1656611/100‑1 (MQ=255)
                                                    cAGCGTCGCCGGATCGCGCTCATCGACGGTGGGCGTCGGATGCGTGGGCTCGGTCGGATCGGTCGggtccgtgggatcggtcggggccgtgggatcggtcg                                                  >  2:2802623/1‑101 (MQ=255)
                                                    cAGCGTCGCCGGATCGCGCTCATCGACGGTGGGCGTCGGATGCGTGGGCTCGGTCGGATCGGTCGggtccgtgggatcgggcgggtccgtggggtcggtcg                                                  >  2:2374426/1‑101 (MQ=255)
                                                    cAGCGTCCCCGGATCGCGCTCATCGACGGTGGGCGTCGGATGCGTGGGCTCGGTCGGATCGGTCGgggccgtgggatcggggggggccgtgggagcggtcg                                                  >  2:2900789/1‑101 (MQ=255)
                                                       cGTCGCCGGATCGCGCTCATCGACGGTGGGCGTCGGATGCGTGGGCTCGGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGAt                                               >  1:1030726/1‑101 (MQ=255)
                                                       cGTCGCCGGATCGCGCTCATCGACGGTGGGCGTCGGATGCGTGGGCTCGGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGAt                                               <  1:539137/101‑1 (MQ=255)
                                                        gTCGCCGGATCGCGCTCATCGACGGTGGGCGTCGGATGCGTGGGCTCGGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATc                                              <  2:819951/101‑1 (MQ=255)
                                                         tCGCCGGATCGCGCTCATCGACGGTGGGCGTCGGATGCGTGGGCTCGGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATcc                                             <  1:4146651/101‑1 (MQ=255)
                                                            ccGGATCGCGCTCATCGACGGTGGGCGTCGGATGCGTGGGCTCGGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGGTCGGGCGGAGCCggg                                          >  2:4345220/1‑99 (MQ=255)
                                                            ccGGATCGCGCTCATCGACGGTGGGCGTCGGATGCGTGGGCTCGGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCGGTg                                          >  2:2190263/1‑101 (MQ=255)
                                                            ccGGATCGCGCTCATCGACGGTGGGCGTCGGATGCGTGGGCTCGGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTg                                          >  1:2768716/1‑101 (MQ=255)
                                                            ccGGATCGCGCTCATCGACGGTGGGCGTCGGATGCGTGGGCTCGGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTg                                          >  2:169117/1‑101 (MQ=255)
                                                                 tCGCGCTCATCGACGGTGGGCGTCGGATGCGTGGGCTCGGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGAtc                                     <  1:4821728/101‑1 (MQ=255)
                                                                  cgcgCTCATCGACGGTGGGCGTCGGATGCGTGGGCTCGGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGAtcg                                    <  1:719765/101‑1 (MQ=255)
                                                                  cgcgCTCATCGACGGTGGGCGTCGGATGCGTGGGCTCGGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGAtcg                                    >  1:3199065/1‑101 (MQ=255)
                                                                    cgcTCATGGACGGTGGGCGTCGGATGCTTGGGCTCGGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGAtcggt                                  <  1:215715/101‑1 (MQ=255)
                                                                      cTCATCGACGGTGGGCGTCGGATGCGTGGGCTCGGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGAtcggtcg                                <  1:1433662/101‑1 (MQ=255)
                                                                      cTCATCGACGGTGGGCGTCGGATGCGTGGGCTCGGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGAtcggtcg                                <  1:4842160/101‑1 (MQ=255)
                                                                       tCATCGACGGTGGGCGTCGGATGCGTGGGCTCGGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGAtcggtcgg                               <  2:528046/101‑1 (MQ=255)
                                                                        cATCGACGGTGGGCGTCGGATGCGTGGGCTCGGGCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGGCggg                              >  1:265492/1‑101 (MQ=255)
                                                                           gggtggggggCGTCGGGGGGGTGGGGTCCGTGGGATCGGTCGGGTCCGTGGGGTCGGTGGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTcc                           <  1:3415430/94‑1 (MQ=1)
                                                                           ccgtgggatccGTCGGGTCCGGGGGATCCGTCGGATCTGTGGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTAGGATCGGTCGGGTcc                           <  2:3954526/91‑1 (MQ=11)
                                                                           cGACGGTGGGCGTCGGATGCGTGGGCTCGGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTcc                           >  2:4447330/1‑101 (MQ=255)
                                                                           cGACGGTGGGCGTCGGATGCGTGGGCTCGGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTcc                           >  1:2565137/1‑101 (MQ=255)
                                                                            gACGGTGGGCGTCGGATGCGTGGGCTCGGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCg                          <  2:3061190/101‑1 (MQ=255)
                                                                            gACGGTGGGCGTCGGATGCGTGGGCTCGGTCGGATCGGTCGGGTCCGTGGGATCGGACGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCg                          >  2:4204723/1‑101 (MQ=255)
                                                                             gcGGTGGGCGTCGGATGCGTGGGCTCGGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGt                         <  2:427864/100‑1 (MQ=255)
                                                                             gcGGTGGGCGTCGGATGCGTGGGCTCGGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGt                         <  1:4447649/100‑1 (MQ=255)
                                                                             aCGGTGGGCGTCGGATGCGTGGGCTCGGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGt                         >  2:3025185/1‑101 (MQ=255)
                                                                              tgggatcggTCGGGTCCGTGGGATCAGTCGGATCCGTGGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTAGGATCGGTCGGGTCCGtc                        <  2:880545/93‑2 (MQ=11)
                                                                              cGGTGGGCGTCGGATGCGTGGGCTCGGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTg                        <  1:288857/101‑1 (MQ=255)
                                                                               gggggcggTCGGTTCCGTGGGGTCGGGGGGGTCGGTGGGGTCGGTGGGATCGGTGGGGTCCGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGtcg                       <  2:1831133/94‑3 (MQ=2)
                                                                               gggatcggTCGGGTCCGTGGGATCGGTCGGATCCGTGGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTAGGATCGGTCGGGTCCGtcg                       >  2:3460898/8‑99 (MQ=11)
                                                                               gggatcggTCGGGTCCGTGGGATCGGTCGGATCCGTGGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTAGGATCGGTCGGGTCCGtcg                       >  1:3913257/8‑99 (MQ=11)
                                                                               ggTGGGCGTCGGATGCGTGGGCTCGGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGGCCGTgg                       >  2:2460259/1‑101 (MQ=255)
                                                                                    gCGTCGGATGCGTGGGCTGGGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGAtcg                  <  2:3410818/101‑1 (MQ=255)
                                                                                    cggTCGGGTCCGTGGGATGTGTCGGATCCGTGGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTAGGATCGGTCGGGTCCGTCGGGtcc                  <  1:4303634/99‑2 (MQ=11)
                                                                                    ccgtggggctcGTCGGATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGAtcg                  <  1:787424/91‑1 (MQ=11)
                                                                                    ccggctgATGCGTCTGATCCGTGGGATCCGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGAtcg                  <  2:4414395/95‑1 (MQ=11)
                                                                                    ccGTTGGATCCGTGGGGTCCGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGAtcg                  <  1:1554360/100‑1 (MQ=11)
                                                                                     cGTCGGATGCGTGGGCTCGGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGCCCGTgcatcaag                 <  2:1523748/101‑8 (MQ=255)
                                                                                      gttggcgccGTGGGATCCGTGGGATCGGTGGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTAGGATCGGTCGGGTCCGTCGGGtccgt                <  2:3599864/93‑1 (MQ=17)
                                                                                      gTCGGATGCGTGGGCTGGGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGAtcggt                <  2:3318434/101‑1 (MQ=255)
                                                                                      gTCGGATGCGTGGGCTCGGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGAtcggt                >  1:2181092/1‑101 (MQ=255)
                                                                                       tCGGATGCGTGGGCTCGGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGAtcggtc               >  1:3150941/1‑101 (MQ=255)
                                                                                         ggATGCGTGGGCTCGGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGGtcggtcgg             >  2:1120016/1‑101 (MQ=255)
                                                                                         gcagtgtcgggATCCGTCGCAACGGTCGGGGCCGTGGGATCGGGCGGATCGGTCGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGAtcggtcgg             <  2:2284463/93‑1 (MQ=11)
                                                                                          gATGCGTGGGCTCGGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGa            >  1:982663/1‑101 (MQ=255)
                                                                                            tcggtgtccTCGGTGGGCTCGGTCGTGTCCGTGGGAGCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTc          <  2:1047287/93‑1 (MQ=11)
                                                                                            tccGTGGGGTCGGTCGGGTCCGTGGGGTCCGGGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTAGGATCGGTCGGGTCCGTCGGGTCCGTGGGATc          <  2:4600624/99‑1 (MQ=11)
                                                                                            tccGTGGGATCGGTCGGGTCCGTCGGATCCGTGGGGTCCGTCGGGTCCGTGGGATCGGTCGGGTCCGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTc          <  1:3242147/99‑1 (MQ=11)
                                                                                             gCGTGGGCTCGGGCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATcc         <  2:1221517/101‑1 (MQ=255)
                                                                                             cggTGGGCTCGGTGGGATCCGTCGGGTCCGTTGGGTCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATcc         <  2:119583/99‑1 (MQ=14)
                                                                                              cctGGGCTCGGTCGGATCGGTCGGGTCCGGGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCTGTCGGGTCCGTGGGATCGGTCGGATCCg        <  1:627537/99‑1 (MQ=255)
                                                                                              cGTGGGCTCGGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCg        >  2:3279987/1‑101 (MQ=255)
                                                                                               gTGGGCTCGGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCgg       >  2:4804690/1‑100 (MQ=255)
                                                                                               gTGGGATCGGTCGGGTCCGTCGGATCCGTGGGGTCCGTCGGGTCCGTGGGATCGGTCGGGTCCGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGt       >  1:3708418/1‑101 (MQ=1)
                                                                                               gTGGGATCGGTCGGGTCCGTCGGATCCGTGGGGTCCGTCGGGTCCGTGGGATCGGTCGGGTCCGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGt       >  2:546907/1‑101 (MQ=11)
                                                                                                tGGGCTCGGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTg      >  1:185182/1‑101 (MQ=255)
                                                                                                tGGGATCGGTCGGGTCCGTCGGATCCGTGGGGTCCGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTcgggc      >  1:119599/1‑97 (MQ=11)
                                                                                                tGGGATCGGTCGGGTCCGTCGGATCCGTGGGGTCCGTCGGGTCCGTGGGATCGGTCGGGTCCGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGtc      >  2:2011910/1‑100 (MQ=11)
                                                                                                 gggGTCGGTGGGGTCGGTGGGGTCCGTGGGATCGGTCGGGTCCGTGGGGTCCGTCGGATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGATCCGTgg     <  1:3638579/101‑1 (MQ=11)
                                                                                                    cTCGGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGAt  <  1:169133/101‑1 (MQ=255)
                                                                                                    cTCGGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTCGGATCGGTCGGAACCGTGGGag  >  2:923653/1‑100 (MQ=255)
                                                                                                    |                                                                                                    
AGCCGGTTTCGCCGGTGCCCGACACCGTCGCCGAGGCTTCGGTGTTGGTGGTCAGCGTCGCCGGATCGCGCTCATCGACGGTGGGCGTCGGATGCGTGGGCTCGGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGAT  >  NC_014034/221849‑222049

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