Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_014034 221,940 A→G 2.7% H239H (CAT→CAC RCAP_rcc00185 ← hypothetical protein

Read alignment evidence...
  seq id position change freq score reads annotation genes product
*NC_014034221,9400A→G 2.7% 5.1 ‑0.9 ‑0.7224H239H (CAT→CACRCAP_rcc00185hypothetical protein
Reads supporting (aligned to +/- strand):  new base (4/2):  ref base (101/117):  total (105/119)
Fisher's exact test strand distribution p-value = 4.23e-01
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 2.05e-01

Only 100 of 260 total aligned reads displayed.

CCACGGTCGAGCCGGTTTCGCCGGTGCCCGACACCGTCGCCGAGGCTTCGGTGTTGGTGGTCAGCGTCGCCGGATCGCGCTCATCGACGGTGGGCGTCGGATGCGTGGGCTCGGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGAT  >  NC_014034/221840‑222040
                                                                                                    |                                                                                                    
  acGGTCGAGCCGGTTTCGCCGGTGCCCGACACCGTCGCCGAGGCTTCGGTGTTGGTGGTCAGCGTCGCCGGATCGCGCTCATCGACGGTGGGCGTCGGATg                                                                                                    >  2:527837/1‑101 (MQ=255)
   cGGTGGAGCCGGTTTCGCCGGTGCCCGACACCGTCGCCGAGGCTTCGGTGTTGGTGGTCAGCGTCGCCGGATCGCGCTCATCGACGGTGGGCGTCGGATGc                                                                                                   <  1:5631756/101‑1 (MQ=255)
    ggTCGAGCCGGTTTCGCCGGTGCCCGACACCGTCGCCGAGGCTTCGGTGTTGGTGGTCAGCGTCGCCGGATCGCGCTCATCGACGGTGGGCGTCGGATGCg                                                                                                  >  1:1453540/1‑101 (MQ=255)
     gTCGAGCCGGTTTCGCCGGTGCCCGACACCGTCGCCGAGGCTTCGGTGTTGGTGGTCAGCGTCGCCGGATCGCGCTCATCGACGGTGGGCGTCGGATGCGt                                                                                                 <  2:3943978/101‑1 (MQ=255)
       cGAGCCGGTTTCGCCGGTGCCCGACACCGTCGCCGAGGCTTCGGTGTTGGTGGTCAGCGTCGCCGGATCGCGCTCATCGACGGTGGGCGTCGGATGCGTgg                                                                                               >  1:3047003/1‑101 (MQ=255)
           ccGGTTTCGCCGGTGCCCGACACCGTCGCCGAGGCTTCGGTGTTGGTGGTCAGCGTCGCCGGATCGCGCTCATCGACGGTGGGCGTCGGATGCGTGGGCtc                                                                                           <  2:773042/101‑1 (MQ=255)
            cGGTTTCGCCGGTGCCCGACACCGTCGCCGAGGCTTCGGTGTTGGTGGTCAGCGTCGCCGGATCGCGCTCATCGACGGTGGGGGTCGGATGCGTGGGCtcg                                                                                          >  1:3183029/1‑101 (MQ=255)
            cGGTTTCGCCGGTGCCCGACACCGTCGCCGAGGCTTCGGTGTTGGTGGTCAGCGTCGCCGGATCGCGCTCATCGACGGTGGGCGTCGGATGCGTGGGCtcg                                                                                          >  2:1202657/1‑101 (MQ=255)
             ggTTTCGCCGGTGCCCGACACCGTCGCCGAGGCTTCGGTGTTGGTGGTCAGCGTCGCCGGATCGCGCTCATCGACGGTGGGCGTCGGATGCGTGGGCtcgg                                                                                         >  1:934315/1‑101 (MQ=255)
                ttCGCCGGTGCCCGACACCGTCGCCGAGGCTTCGGTGTTGGTGGTCAGCGTCGCCGGATCGCGCTCATCGACGGTGGGCGTCGGATGCGTGGGCtcggtcg                                                                                      >  2:3212172/1‑101 (MQ=255)
                 tCGCCGGTGCCCGACACCGTCGCCGAGGCTTCGGTGTTGGTGGTCAGCGTCGCCGGATCGCGCTCATCGACGGTGGGCGTCGGATGCGTGGGCtcggtcgg                                                                                     <  1:3907545/101‑1 (MQ=255)
                 tCGCCGGTGCCCGACACCGTCGCCGAGGCTTCGGTGTTGGTGGTCAGCGTCGCCGGATCGCGCTCATCGACGGTGGGCGTCGGATGCGTGGGCtcgggcgg                                                                                     >  1:4452585/1‑101 (MQ=255)
                  cGCCGGTGCCCGACACCGTCGCCGAGGCTTCGGTGTTGGTGGTCAGCGTCGCCGGATCGCGCTCATCGACGGTGGGCGTCGGATGCGTGGGCTCGGTCgga                                                                                    <  1:4079656/101‑1 (MQ=255)
                  cGCCGGTGCCCGACACCGTCGCCGAGGCTTCGGTGTTGGTGGTCAGCGTCGCCGGATCGCGCTCATCGACGGTGGGCGTCGGATGCGTGGGCTCGGTCgga                                                                                    <  1:2372898/101‑1 (MQ=255)
                   gCCGGTGCCCGACACCGTCGCCGAGGCTTCGGTGTTGGTGGTCAGCGTCGCCGGATCGCGCTCATCGACGGTGGGCGTCGGATGCGTGGGCTCGGTCggat                                                                                   >  2:1519357/1‑101 (MQ=255)
                   gCCGGTGCCCGACACCGTCGCCGAGGCTTCGGTGTTGGTGGTCAGCGTCGCCGGATCGCGCTCATCGACGGTGGGCGTCGGATGCGTGGGCTCGGTCggat                                                                                   >  1:4026728/1‑101 (MQ=255)
                   gCCGGTGCCCGACACCGTCGCCGAGGCTTCGGTGTTGGTGGTCAGCGTCGCCGGATCGCGCTCATCGACGGTGGGCGTCGGATGCGTGGGCTCGGTCggat                                                                                   <  1:2186886/101‑1 (MQ=255)
                    ccGGTGCCCGACACCGTCGCCGAGGCTTCGGTGTTGGTGGTCAGCGTCGCCGGATCGCGCTCATCGACGGTGGGCGTCGGATGCGTGGGCTCGGTCggatc                                                                                  <  1:1989272/101‑1 (MQ=255)
                     cGGTGCCCGACACCGTCGCCGAGGCTTCGGTGTTGGTGGTCAGCGTCGCCGGATCGCGCTCATCGACGGTGGGCGTCGGATGCGTGGGCTCGGTCggatcg                                                                                 <  1:2025273/101‑1 (MQ=255)
                     cGGTGCCCGACACCGTCGCCGAGGCTTCGGTGTTGGTGGTCAGCGTCGCCGGATCGCGCTCATCGACGGTGGGCGTCGGATGCGTGGGCTCGGTCggatcg                                                                                 >  2:867918/1‑101 (MQ=255)
                        tGCCCGACACCGTCGCCGAGGCTTCGGTGTTGGTGGTCAGCGTCGCCGGATCGCGCTCATCGACGGTGGGCGTCGGATGCGTGGGCTCGGTCggatcggtc                                                                              <  2:2629724/101‑1 (MQ=255)
                        tGCCCGACACCGTCGCCGAGGCTTCGGTGTTGGTGGTCAGCGTCGCCGGATCGCGCTCATCGACGGTGGGCGTCGGATGCGTGGGCTCGGTCggatcggtc                                                                              <  1:2231480/101‑1 (MQ=255)
                        tGCCCGACACCGTCGCCGAGGCTTCGGTGTTGGTGGTCAGCGTCGCCGGATCGCGCTCATCGACGGTGGGCGTCGGATGCGTGGGCTCGGTCggatcggtc                                                                              <  1:1076952/101‑1 (MQ=255)
                             gACACCGTCGCCGAGGCTTCGGTGTTGGTGGTCAGCGTCGCCGGATCGCGCTCATCGACGGTGGGCGTCGGATGCGTGGGCTCGGTCggatcggtcgggtc                                                                         >  2:2404455/1‑101 (MQ=255)
                              gcacCGTCGCCGGGGCTTTGGTGTTGGTGGTCAGCGTCGCCGGATCGCGCTCATCGACGGTGGGCGTCGGATGCGTGGGCTCGGTCggatcggtcgggtcc                                                                        <  2:2311967/100‑1 (MQ=255)
                               ccccGTCGCCGAGGCTTCGGTGTTGGTGGTCAGCGTCGCCGGATCGCGCTCATCGACGGTGGGCGTCGGATGCGTGGGCTCGGTCggatcggtcgggtccg                                                                       <  2:1650911/99‑1 (MQ=255)
                               cacCGTCGCCGAGGCTTCGGTGTTGGTGGTCAGCGTCGCCGGATCGCGCTCATCGACGGTGGGCGTCGGATGCGTGGGCTCGGTCggatcggtcgggtccg                                                                       <  2:4628118/101‑1 (MQ=255)
                                acCGTCGCCGAGGCTTCGGTGTTGGTGGTCAGCGTCGCCGGATCGCGCTCATCGACGGTGGGCGTCGGATGCGTGGGCTCGGTCggatcggtcgggtccgt                                                                      <  2:2782169/101‑1 (MQ=255)
                                 ccGTCGCCGAGGCTTCGGTGTTGGTGGTCAGCGTCGCCGGATCGCGCTCATCGACGGTGGGCGTCGGATGCGTGGGCTCGGTCggatcggtcgggtccgtg                                                                     <  2:934239/101‑1 (MQ=255)
                                   gTCGCCGAGGCTTCGGTGTTGGTGGTCAGCGTCGCCGGATCGCGCTCATCGACGGTGGGCGTCGGATGCGTGGGCTCGGTCggatcggtcgggtccgtggg                                                                   <  1:2528413/101‑1 (MQ=255)
                                      gCCGAGGCTTCGGTGTTGGTGGTCAGCGTCGCCGGATCGCGCTCATCGACGGTGGGCGTCGGATGCGTGGGCTCGGTCggatcggtcgggtccgtgggatc                                                                <  2:142847/101‑1 (MQ=255)
                                       ccGAGGCTTCGGTGTTGGTGGTCAGCGTCGCCGGATCGCGCTCATCGACGGTGGGCGTCGGATGCGTGGGCTCGGTCggatcggtcgggtccgtgggatcg                                                               <  1:4652181/101‑1 (MQ=255)
                                        cGAGGCTTCGGTGTTGGTGGTCAGCGTCGCCGGATCGCGCTCATCGACGGTGGGCGTCGGATGCGTGGGCTCGGTCggatcggtcgggtccgtgggatcgg                                                              <  1:3302589/101‑1 (MQ=255)
                                          aGGCTTCGGTGTTGGTGGTCAGCGTCGCCGGATCGCGCTCATCGACGGTGGGCGTCGGATGCGTGGGCTCGGTCggatcggtcgggtccgtgggatcggtc                                                            <  2:2974997/101‑1 (MQ=255)
                                          aGGCTTCGGTGTTGGTGGTCAGCGTCGCCGGATCGCGCTCATCGACGGTGGGCGTCGGATGCGTGGGCTCGGTCggatcggtcgggtccgtgggatcggtc                                                            <  2:5374910/101‑1 (MQ=255)
                                           ggCTTCGGTGTTGGTGGTCAGCGTCGCCGGATCGCGCTCATCGACGGTGGGCGTCGGATGCGTGGGCTCGGTCggatcggtcgggtccgtgggatcggtcg                                                           <  1:2408704/101‑1 (MQ=255)
                                            gCTTCGGTGTTGGTGGTCAGCGTCGCCGGATCGCGCTCATCGACGGTGGGCGTCGGATGCGTGGGCTCGGTCggatcggtcgggtccgtgggatcggtcgg                                                          <  2:715774/101‑1 (MQ=255)
                                            gCTTCGGTGTTGGTGGTCAGCGTCGCCGGATCGCGCTCATCGACGGTGGGCGTCGGATGCGTGGGCTCGGTCggatcggtcgggtccgtgggatcggtcgg                                                          >  1:5300326/1‑101 (MQ=255)
                                             cTTCGGTGTTGGTGGTCAGCGTCGCGGGTCCGCCCTCATCGCGGGGGGGGGCGGGAGGCGGGGCCTCGGTGggatcggtggggcccgtggggcggggcggg                                                         >  2:4997613/1‑91 (MQ=255)
                                             cTTCGGTGTTGGTGGTCAGCGTCGCCGGATCGCGCTCATCGACGGTGGGCGTCGGATGCGTGGGCTCGGTCggatcggtcgggtccgtgggatcggtcggg                                                         >  2:2706925/1‑101 (MQ=255)
                                              ttCGGTGTTGGTGGTCAGCGTCGCCGGATCGCGCTCATCGACGGTGGGCGTCGGATGCGTGGGCTCGGTCggatcggtcgggtccgtgggatcggtcgggt                                                        <  1:4316893/101‑1 (MQ=255)
                                               tCGGTGTTGGTGGTCAGCGTCGCCGGATCGCGCTCATCGACGGTGGGCGTCGGATGCGTGGGCTCGGTCggatcggccgggtccgtgggatcggtcgggtc                                                       >  2:1554566/1‑101 (MQ=255)
                                                   tgtTGGTGGTCAGCGTCGCCGGATCGCGCTCATCGACGGTGGGCGTCGGATGCGTGGGCTCGGTCggatcggtcgggtccgtgggatcggtcgggtccgtg                                                   >  2:100369/1‑101 (MQ=255)
                                                   tgtTGGTGGTCAGCGTCGCCGGATCGCGCTCATCGACGGTGGGCGTCGGATGCGTGGGCTCGGTCggatcggtcgggtccgtgggatcggtcgggtccgtg                                                   <  2:2151345/101‑1 (MQ=255)
                                                     ttGGTGGTCAGCGTCGCCGGATCGCGCTCATCGACGGTGGGCGTCGGATGCGTGGGCTCGGTCGGatcggtcgggtccgtgggatcggtcgggtccgtggg                                                 >  1:5511237/1‑101 (MQ=255)
                                                       ggtggtCAGCGTCGCCGGATCGCGCTCATCGACGGTGGGCGTCGGATGCGTGGGCTCGGTCGGATcggtcgggtccgtgggatcggtcgggcccgtgggat                                               >  1:4800205/1‑101 (MQ=255)
                                                       ggtggtCAGCGTCGCCGGATCGCGCTCATCGACGGTGGGCGTCGGATGCGTGGGCTCGGTCGGATcggtcgggtccgtgggatcggccgggtccgtgggat                                               >  1:816576/1‑101 (MQ=255)
                                                         gggtCAGCGTCGCCGGATCGCGCTCATCGACGGTGGGCGTCGGATGCGTGGGCTCGGTCGGATCGgtcgggtccgtgggatcggtcgggtccgtgggatcg                                             <  2:4328694/100‑1 (MQ=255)
                                                          ggtCAGCGTCGCCGGATCGCGCTCATCGACGGTGGGCGTCGGATGCGTGGGCTCGGTCGGATCGGtcgggtccgtgggatcggtcgggtccgtggggtcgg                                            >  2:1305788/1‑101 (MQ=255)
                                                             cAGCGTCGCCGGATCGCGCTCATCGACGGTGGGCGTCGGATGCGTGGGCTCGGTCGGATCGGTCGggtccgtgggatcggtcgggtccgtgggatcggtcg                                         >  1:4111094/1‑101 (MQ=255)
                                                              aGCGTCGCCGGATCGCGCTCATCGACGGTGGGCGTCGGATGCGTGGGCTCGGTCGGATCGGTCGGgtccgtgggatcggtcgggtccgtgggatcggtcgg                                        <  2:2991851/101‑1 (MQ=255)
                                                              aGCGTCGCCGGATCGCGCTCATCGACGGTGGGCGTCGGATGCGTGGGCTCGGTCGGATCGGTCGGgtccgtgggatcggtcgggtccgtgggatcggtcgg                                        >  2:4449078/1‑101 (MQ=255)
                                                              aGCGTCGCCGGATCGCGCTCATCGACGGTGGGCGTCGGATGCGTGGGCTCGGTCGGATCGGTCGGgtccgtgggatcggtcgggtccgtgggatcggtcgg                                        <  1:4209440/101‑1 (MQ=255)
                                                                cGTCGCCGGATCGCGCTCATCGACGGTGGGCGTCGGATGCGTGGGCTCGGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGAt                                      <  2:4799924/101‑1 (MQ=255)
                                                                 gTCGCCGGATCGCGCTCATCGACGGTGGGCGTCGGATGCGTGGGCTCGGTCGGATCGGTCGGGTCCGTGGGATCGGTCCGGTCCGTGGGGTCGGGCGGacc                                     >  2:5237227/1‑99 (MQ=255)
                                                                    gCCGGATCGCGCTCATCGACGGTGGGCGTCGGATGCGTGGGCTCGGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGt                                  >  1:5029945/1‑101 (MQ=255)
                                                                     ccGGATCGCGCTCATCGACGGTGGGCGTCGGATGCGTGGGCTCGGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTg                                 >  1:1636627/1‑101 (MQ=255)
                                                                      cGGATCGCGCTCATCGACGGTGGGCGTCGGATGCGTGGGCTCGGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTgg                                <  1:4869751/101‑1 (MQ=255)
                                                                      cGGATCGCGCTCATCGACGGTGGGCGTCGGATGCGTGGGCTCGGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTgg                                >  1:3384451/1‑101 (MQ=255)
                                                                      cGGATCGCGCTCATCGACGGTGGGCGTCGGATGCGTGGGCTCGGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGGCCGTGGGGTCGGTCGGATCCGTgg                                >  1:1154642/1‑101 (MQ=255)
                                                                       ggATCGCGCTCATCGACGGTGGGCGTCGGATGCGTGGGCTCGGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTggg                               >  2:2900764/1‑101 (MQ=255)
                                                                         aTCGCGCTCATCGACGGTGGGCGTCGGATGCGTGGGCTCGGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGAt                             <  1:3627000/101‑1 (MQ=255)
                                                                          tCGCGCTCATCGACGGTGGGCGTCGGATGCGTGGGCTCGGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGAtc                            >  2:1914091/1‑101 (MQ=255)
                                                                            gcgcTCATCGACGGTGGGCGTCGGATGCGTGGGCTCGGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGAtcgg                          >  1:4060305/1‑101 (MQ=255)
                                                                            gcgcTCATCGACGGTGGGCGTCGGATGCGTGGGCTCGGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGAtcgg                          <  1:5045283/101‑1 (MQ=255)
                                                                             cgcTCATCGACGGTGGGCGTCGGATGCGTGGGCTCGGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGAtcggt                         <  1:118875/101‑1 (MQ=255)
                                                                              gcTCATCGACGGTGGGCGTCGGATGCGTGGGCTCGGTCGGATCGGTCGGGTCCGTGGGTTCGGGCGGGTCCGGGGGATCGGTCGGATCCGTGGGAttggtc                        >  2:2276944/1‑101 (MQ=25)
                                                                              gcTCATCGACGGTGGGCGTCGGATGCGTGGGCTCGGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGAtcggtc                        <  2:306366/101‑1 (MQ=255)
                                                                              gcTCATCGACGGTGGGCGTCGGATGCGTGGGCTCGGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGAtcggtc                        >  2:3144217/1‑101 (MQ=255)
                                                                              gcTCATCGACGGTGGGCGTCGGATGCGTGGGCTCGGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGAtcggtc                        <  1:1575881/101‑1 (MQ=255)
                                                                                tCATCGACGGTGGGCGTCGGATGCGTGGGCTCGGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGAtcggtcgg                      <  1:648204/101‑1 (MQ=255)
                                                                                  aTCGACGGTGGGCGTCGGATGCGTGGGCTCGGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGGTCGGTCGGATCCGTGGGATCGGTCgggg                    >  1:539305/1‑100 (MQ=255)
                                                                                  aTCGACGGTGGGCGTCGGATGCGTGGGCTCGGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGt                    >  1:3722078/1‑101 (MQ=255)
                                                                                   tCGACGGTGGGCGTCGGATGCGTGGGCTCGGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATTGGTCGGGTc                   <  2:1973044/101‑1 (MQ=255)
                                                                                    cGACGGTGGGCGTCGGATGCGTGGGCTCGGTCGGATCTGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTcc                  <  1:4942114/101‑1 (MQ=255)
                                                                                    cGACGGTGGGCGTCGGATGCGTGGGCTCGGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTcc                  >  2:3460952/1‑101 (MQ=255)
                                                                                      ggggaatctttCGGGTCCGGGGGATCGGTCGTATCCTTGGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTAGGATCGGTCGGGTCCGt                <  2:4109081/91‑1 (MQ=2)
                                                                                      gcGGTGGGCGTCGGATGCGTGGGCTCGGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGt                <  1:1386912/100‑1 (MQ=255)
                                                                                        gggttcgatCGGGTCCGTGGGATCGGTCCGATCCGTGGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTAGGATCGGTCGGGTCCGtcg              <  2:3578538/93‑3 (MQ=11)
                                                                                         ggggtgggtggATGCGTTGGGTCCGGTGGATCGGTCGGATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTggg             <  1:5695518/92‑1 (MQ=11)
                                                                                         gTGGGCGTCGGATGGGTGGGCTCGGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTggg             <  1:2277070/101‑1 (MQ=255)
                                                                                          tGGGCGTCGGATGCGTGGGCTCGGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTgggg            >  2:4227518/1‑100 (MQ=255)
                                                                                          tGGGCGTCGGATGCGTGGGCTCGGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGa            <  1:2797046/101‑1 (MQ=255)
                                                                                          ggccggcgggATGCGGGGGGTCGGGGGGGTCGGTCGGGTCTGTGCGATCGGGCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGa            <  2:2586579/93‑1 (MQ=14)
                                                                                           gggCGTCGGATGCGTGGGCTCGGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGAt           <  2:4066746/101‑1 (MQ=255)
                                                                                           gggCGTCGGATGCGTGGGCTCGGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGAt           <  2:484663/101‑1 (MQ=255)
                                                                                           gggCGTCGGATGCGTGGGCTCGGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTAGGATCCGTGGGATCGGTCGGGTCCGTGGGAt           >  2:133832/1‑101 (MQ=255)
                                                                                           gggCGGCGAATGCGTGGGCTCGGTCGGATCGGGCGGGTCCGTGGGATCGGTGGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGAt           <  2:2033743/101‑1 (MQ=18)
                                                                                           gccgtgcGGATGCGTGGGCTCGGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGAt           <  2:4708678/95‑1 (MQ=255)
                                                                                            tcggTCGGGTCCGGGGGGTCGGTCGGATCCGTGGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTAGGATCGGTCGGGTCCGTCGGGtc          <  1:5099218/98‑1 (MQ=1)
                                                                                              ggTCGGGTCCGTGGGATCGGTCGGATCCGTGGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTAGGATCGGTCGGGTCCGTCGGGtccg        >  1:3538633/2‑99 (MQ=1)
                                                                                              cGTCGGATGCGTGGGCTCGGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGAtcgg        <  2:816508/101‑1 (MQ=255)
                                                                                               gTCGGATGCGTGGGCTCGGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGAtcggt       >  1:899938/1‑101 (MQ=255)
                                                                                               gTCGGATGCGTGGGCTCGGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGGCGGGTCCGTGGGAtcggg       >  2:3353016/1‑100 (MQ=255)
                                                                                                tggtaTGAGTTGGTTCAGTGGGAGCGGTCGCATCCGTGTGATCGGTCGGGTCCGTGGGATCCGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGAtcggtc      <  2:5723133/97‑1 (MQ=11)
                                                                                                tCGGATGCGTGGGCTCGGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGAtcggtc      <  1:4592732/101‑1 (MQ=255)
                                                                                                 cGGATGCGTGGGCTCGGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGAtcggtcg     >  2:3667348/1‑101 (MQ=255)
                                                                                                 cGGATGCGTGGGCTCGGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGAtcggtcg     >  2:5494938/1‑101 (MQ=255)
                                                                                                  ggATGCGTGGGCTCGGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGAtcggtcgg    >  1:583036/1‑101 (MQ=255)
                                                                                                    aTGCGTGGGCTCGGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGAt  >  1:2877781/1‑101 (MQ=255)
                                                                                                    |                                                                                                    
CCACGGTCGAGCCGGTTTCGCCGGTGCCCGACACCGTCGCCGAGGCTTCGGTGTTGGTGGTCAGCGTCGCCGGATCGCGCTCATCGACGGTGGGCGTCGGATGCGTGGGCTCGGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGAT  >  NC_014034/221840‑222040

Base quality scores:  ATCG  < 3 ≤  ATCG  < 17 ≤  ATCG  < 24 ≤  ATCG  < 33 ≤  ATCG  < 40 ≤  ATCG