Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_014034 2,493,952 +G 98.0% coding (252/252 nt) RCAP_rcc02319 → hypothetical protein

Read alignment evidence...
  seq id position change freq score reads annotation genes product
*NC_0140342,493,9521.→G98.0% 5.4 ‑0.5 ‑1.0247coding (252/252 nt)RCAP_rcc02319hypothetical protein
Reads supporting (aligned to +/- strand):  new base (64/178):  ref base (2/3):  total (66/181)
Fisher's exact test strand distribution p-value = 6.12e-01
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 9.66e-02

Only 100 of 291 total aligned reads displayed.

TGAGCCCCGAAGAGCGCCGCGCGGTGATGGCCGAGCTGCCCACCCGCGCGCCGCGTCTGGCCCAGCGCCGCGGCGGGCGCGCGGCACGACAGGCCGAGTAA..GG.GG.GGGCGCTGCCCCCGTTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCG  >  NC_014034/2493852‑2494047
                                                                                                     |                                                                                                  
tGAGCCCCGAAGAGCGCCGCGCGGTGATGGCCGAGCTGCCCACCCGCGCGCCGCGTCTGGCCCAGCGCCGCGGCGGGCGCGCGGCACGACAGGCCGAGTaa                                                                                                     <  1:5329467/101‑1 (MQ=255)
 gAGCCCCGAAGAGCGCCGTGCGGTGATGGCCGAGCTGCCCACCCGCGCGCCGCGTCTGGCCCAGCGCCGCGGCGGGCGCGCGGCACGACAGGCCGAGCAA..g                                                                                                  >  2:6778567/1‑101 (MQ=255)
 gAGCCCCGAAGAGCGCCGCGCGGTGATGGCCGAGCTGCCCACCCGCGCGCCGCGTCTGGCCCAGCGCCGCGGCGGGCGCGCGGCACGACAGGCCGAGTAA..g                                                                                                  <  1:2621665/101‑1 (MQ=255)
 gAGCCCCGAAGAGCGCCGCGCGGTGATGGCCGAGCTGCCCACCCGCGCGCCGCGTCTGGCCCAGCGCCGCGGCGGGCGCGCGGCACGACAGGCCGAGGAA..g                                                                                                  >  1:1123659/1‑101 (MQ=255)
  aGCCCCGAAGAGCGCCGCGCGGTGATGGCCGAGCTGCCCACCCGCGCGCCGCGTCTGGCCCAGCGCCGCGGCGGGCGCGCGGCACGACAGGCCGAGCAA..gg                                                                                                 >  2:6198389/1‑101 (MQ=255)
   gCCCCGAAGAGCGCCGCGCGGTGATGGCCGAGCTGCCCACCCGCGCGCCGCGTCTGGCCCAGCGCCGCGGCGGGCGCGCGGCACGACAGGCCGAGTAA..gg.g                                                                                               <  1:932905/101‑1 (MQ=255)
   cccccGAAGAGCGCCGCGCGGTGATGGCCGAGCTGCCCACCCGCGCGCCGCGCCTGGCCCAGCGCCGCGGCGGGCGCGCGGCACGACAGGCCGAGTAA..gg.g                                                                                               <  1:2797822/100‑1 (MQ=255)
     cccGAAGAGCGCCGCGCGGTGATGGCCGAGCTGCCCACCCGCGCGCCGCGTCTGGCCCAGCGCCGCGGCGGGCGCGCGGGACGGCCGGCCGCGTAA..gg.gg.g                                                                                            >  2:6940032/1‑101 (MQ=255)
         aaGAGCGCCGCGCGGTGATGGCCGAGCTGCCCACCCGCGCGCCGCGTCTGGCCCAGCGCCGCGGCGGGCGCGCGGCACGACAGGCCGAGTAA..gg.gg.ggggg                                                                                        >  1:1803594/1‑99 (MQ=255)
            ggcgcCGCGCGGTGATGGCCGAGCTGCCCACCCGCGCGCCGCGTCTGGCCCAGCGCCGCGGCGGGCGCGCGGCACGACAGGCCGAGTAA..gg.gg.ggggcgct                                                                                     <  2:4558375/100‑6 (MQ=255)
            agCGCCGCGCGGTGATGGCCGAGCTGCCCACCCGCGCGCCGCGTCTGGCCCAGCGCCGCGGCGGGCGCGCGGCACGACAGGCCGAGTAA..gg.gg.ggggcgct                                                                                     <  1:4677198/101‑6 (MQ=255)
                  gcgcgGTGATGGCCGAGCTGCCCACCCGCGCGCCGCGTCTGGCCCAGCGCCGCGGCGGGCGCGCGGCACGACAGGCCGAGTAAG.GG.GG.GGGCGCTGccccc                                                                                <  2:877258/101‑1 (MQ=255)
                  gcgcgGTGATGGCCGAGCTGCCCACCCGCGCGCCGCGTCTGGCCCAGCGCCGCGGCGGGCGCGCGGCACGACAGGCCGAGTAAG.GG.GG.GGGCGCTGccccc                                                                                <  2:2096441/101‑1 (MQ=255)
                   cgcgGTGATGGCCGAGCTGCCCCCCCGCGCGCCGCGTCTGGCCCAGCGCCGCGGCGGGCGCGCTGCACGACAGGCCGAGTAAG.GG.GG.GGGCGCTGcccccc                                                                               <  1:6322406/101‑2 (MQ=255)
                   cgcgGTGATGGCCGAGCTGCCCACCCGCGCGCCGCGTCTGGCCCAGCGCCGCGGCGGGCGCGCGGCACGACAGGCCGAGTAAG.GG.GG.GGGCGCTGcccccc                                                                               <  1:4329320/101‑2 (MQ=255)
                   cgcgGTGATGGCCGAGCTGCCCACCCGCGCGCCGCGTCTGGCCCAGCGCCGCGGCGGGCGCGCGGCACGACAGGCCGAGTAAG.GG.GG.GGGCGCTGcccccc                                                                               <  2:4392178/101‑2 (MQ=255)
                    gcgGTGATGGCCGAGCTGCCCACCCGCGCGCCGCGTCTGGCCCAGCGCCGCGGCGGGCGCGCGGCACGACAGGCCGAGTAAG.GG.GG.GGGCGCTGccccccg                                                                              <  1:4727421/101‑3 (MQ=255)
                     cgGTGATGGCCGAGCTGCCCACCCGCGCGCCGCGTCTGGCCCAGCGCCGCGGCGGGCGCGCGGCACGACAGGCCGAGTAAG.GG.GG.GGGCGCTGccccccgt                                                                             <  1:1697260/101‑4 (MQ=255)
                       gTGATGGCCGAGCTGCCCGCCCGCGCGCCGCGTCTGGCCCAGCGCCGCGGCGGGCGCGCGGCACGACAGGCCGAGTAAG.GG.GG.GGGCGCTGccccccgtcc                                                                           <  2:6968981/101‑6 (MQ=255)
                       gTGATGGCCGAGCTGCCCACCCGCGCGCCGCGTCTGGCCCAGCGCCGCGGCGGGCGCGCGGCACGACAGGCCGAGTAAG.GG.GG.GGGCGCTGccccccgtcc                                                                           <  2:6917592/101‑6 (MQ=255)
                        tGATGGCCGAGCTGCCCACCCGCGCGCCGCGTCTGGCCCAGCGCCGCGGCGGGCGCGCGGCACGACAGGCCGAGTAAG.GG.GG.GGGGGCTGCCCCCCGTCCg                                                                          >  2:3299322/1‑101 (MQ=255)
                         gATGGCCGAGCTGCCCACCCGCGCGCCGCGTCTGGCCCAGCGCCGCGGCGGGCGCGCGGCACGACAGGCCGAGTAAG.GG.GG.GGGCGCTGCCCCCCGTCCGc                                                                         <  2:2562432/101‑1 (MQ=39)
                          aTGGCCGAGCTGCCCACCCGCGCGCCGCGTCTGGCCCAGCGCCGCGGCGGGCGCGCGGCACGACAGGCCGAGTAAG.GG.GG.GGGCGCTGgcccccggccgca                                                                        >  2:3169678/1‑93 (MQ=255)
                          aTGGCCGAGCTGCCCACCCGCGCGCCGCGTCTGGCCCAGCGCCGCGGCGGGCGCGCGGCACGACAGGCCGAGTAAG.GG.GG.GGGCGCTGCCCCCCGTCCGCa                                                                        <  2:2959443/101‑1 (MQ=39)
                           tGGCCTAGCTGCCCACCCGCGCGCCGCGTCTGGCCCAGCGCCGCGGCGGGCGCGCGGCACGACAGGCCGAGTAAG.GG.GG.GGGCGCTGCCCCCCGTCCGCaa                                                                       <  2:6545544/101‑1 (MQ=38)
                           tGGCCGAGCTGCCCACCCGCGCGCCGCGTCTGGCCCAGCGCCGCGGCGGGCGCGCGGCACGACAGGCCGAGTAAG.GG.GG.GGGCGCTGCCCCCCGTCCGCaa                                                                       <  1:7097902/101‑1 (MQ=38)
                           gggCCGAGCGCCCCCCCCGCGCGCCGCGTCTGGCCCAGCGCCGCGGCGGGCGCGCGGCACGCCGGGCCGAGTAAG.GG.GG.GGGCGCTGCCCCCCGTCCGCaa                                                                       <  2:5905489/100‑1 (MQ=25)
                            ggCCGAGCTGCCCCCCCGCGCGCCGCGTCTGGCCCAGCGCCGCGGCGGGCGCGCGGCACGACAGGCCGAGTAAG.GG.GG.GGGCGCTGCCCCCCGTCCGCAAg                                                                      <  1:3650248/101‑1 (MQ=38)
                            ggCCGAGCTGCCCCCCCGCGCGCCGCGTCTGGCCCAGCGCCGCGGCGGGCGCGCGGCACGACAGGCCGAGTAAG.GG.GG.GGGCGCTGCCCCCCGTCCGCAAg                                                                      <  1:1848139/101‑1 (MQ=38)
                             gCCGAGCTGCCCCCCCGCGCGCCGCGTCTGGCCCAGCGCCGCGGCGGGCGCGCGGCACGACAGGCCGAGTAAG.GG.GG.GGGCGCTGCCCCCCGTCCGCAAGc                                                                     <  1:6739299/101‑1 (MQ=38)
                               cGAGCTGCCCACCCGCGCGCCGCGTCTGGCCCAGCGCCGCGGCGGGCGCGCGGCACGACAGGCCGAGTAAG.GG.GG.GGGCGCTGCCCCCCGTCCGCAAGCgg                                                                   <  1:5980542/101‑1 (MQ=38)
                               cGAGCTGCCCACCCGCGCGCCGCGTCTGGCCCAGCGCCGCGGCGGGCGCGCGGCACGACAGGCCGAGTAAG.GG.GG.GGGCGCTGCCCCCCGTCCGCAAGCgg                                                                   <  1:4563872/101‑1 (MQ=38)
                                gAGCTGCCCACCCGCGCGCCGCGTCTGGCCCAGCGCCGCGGCGGGCGCGCGGCACGACAGGCCGAGTAAG.GG.GG.GGGGGCTGCCCCCCGGCCGCCAGCGGa                                                                  >  1:5653440/1‑101 (MQ=255)
                                 aGCTGCCCCCCCGCGCGCCGCGTCTGGCTCAGCGCCGCGGCGGGCGCGCGGCACGACAGGCCGAGTAAG.GG.GG.GGGCGCTGCCCCCCGTCCGCAAGCGGAc                                                                 <  1:6186800/101‑1 (MQ=37)
                                 aGCTGCCCACCCGCGCGCCGCGTCTGGCCCAGCGCCGCGGCGGGCGCGCGGCACGACAGGCCGAGTAAG.GG.GG.GGGCGCTGCCCCCCGTCCGCAAGCgggc                                                                 >  2:3644341/1‑99 (MQ=38)
                                 aGCTGCCCACCCGCGCGCCGCGTCTGGCCCAGCGCCGCGGCGGGCGCGCGGCACGACAGGCCGAGTAAG.GG.GG.GGGCGCTGCCCCCCGTCCGCAAGCGGAc                                                                 <  2:6034580/101‑1 (MQ=37)
                                   cTGCCCACCCGCGCGCCGCGTCTGGCCCAGCGCCGCGGCGGGCGCGCGGCACGACAGGCCGAGTAAG.GG.GG.GGGGGCTGCCCCCCGTCCGCAAGCGGACTc                                                               >  1:2416325/1‑101 (MQ=37)
                                    tGCCCACCCGCGCGCCGCGTCTGGCCCAGCGCCGCGGCGGGCGCGCGGCACGACAGGCCGAGTAAG.GG.GG.GGGCGTTGGCCCCCGGCCGCAAGCGGACGcc                                                              >  2:1009549/1‑101 (MQ=255)
                                    tGCCCACCCGCGCGCCGCGTCTGGCCCAGCGCCGCGGCGGGCGCGCGGCACGACAGGCCGAGTAAG.GG.GG.GGGCGCTGCCCCCCGTCCGCAAGCGGACTcc                                                              <  1:5226157/101‑1 (MQ=37)
                                      cccACCCGCGCGCCGCGTCTGGCCCAGCGCCGCGGCGGGCGCGCGGCACTACAGGCCGAGTAAG.GG.GG.GGGCGCTGCCCCCCGTCCGCAAGCGGACTcccc                                                            <  1:1037378/101‑1 (MQ=37)
                                      cccACCCGCGCGCCGCGTCTGGCCCAGCGCCGCGGCGGGCGCGCGGCACGACAGGCCGAGTAAG.GG.GG.GGGCGCTGCCCCCCGTCCGCAAGCGGACTcccc                                                            >  1:2499712/1‑101 (MQ=37)
                                        cACCCGCGCGCCGCGTCTGGCCCAGCGCCGCGGCGGGCGCGCGGCACGACAGGCCGAGTAAG.GG.GG.GGGCGCTGCCCCCCGTCCGCAAGCGGACTcccccc                                                          <  1:1302290/101‑1 (MQ=37)
                                         ccccGCGCGCCGCGTCTGGCCCAGCGCCGCGGCGGGCGCGCGGCACGACAGGCCGAGTAAG.GG.GG.GGGCGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCg                                                         <  2:4711275/100‑1 (MQ=25)
                                           ccGCGCGCCGCGTCTGGCCCAGCGCCGCGGCGGGCGCGCGGCACGACAGGCCGAGTAAG.GG.GG.GGGCGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCggg                                                       >  1:1636112/1‑101 (MQ=25)
                                           ccGCGCGCCGCGTCTGGCCCAGCGCCGCGGCGGGCGCGCGGCACGACAGGCCGAGTAAG.GG.GG.GGGCGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCggg                                                       <  2:1644715/101‑1 (MQ=25)
                                            cgcgccCCGCGTCCGCCCCCGCCCCGCGGCTTGCGCGCGGCACGACAGGCCGAGTAA..GGTTG.GGGCGCTTCCCCCCGTCCGCAAGCGGACTCCCCCCGGGa                                                      <  2:2284499/101‑1 (MQ=16)
                                             gcgcgcCGCGTCTGGCCCAGCGCCGCGGCGGGCGCGCGGCACGACAGGCCGAGTAAG.GG.GG.GGGCGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGat                                                     <  2:1074776/101‑1 (MQ=25)
                                             gcgcgcCGCGTCTGGCCCAGCGCCGCGGCGGGCGCGCGGCACGACAGGCCGAGGAAG.GG.GG.GGGCGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGat                                                     >  2:1534831/1‑101 (MQ=25)
                                                 cccgcGTCTGGCCCAGCGCCGCGGCGGGCGCGCGGCACGACAGGCCGAGTAAG.GG.GG.GGGCGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATAttt                                                 <  2:5546794/100‑1 (MQ=25)
                                                  ccgcGTCTGGCCCAGCGCCGCGGCGGGCGCGCGGCACGACAGGCCGAGTAAG.GG.GG.GGGCGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTg                                                <  1:6661287/101‑1 (MQ=25)
                                                  ccgcGTCTGGCCCAGCGCCGCGGCGGGCGCGCGGCACGACAGGCCGAGTAAG.GG.GG.GGGCGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGAGATTTg                                                >  1:761342/1‑101 (MQ=25)
                                                      gTCTGGCCCAGCGCCGCGGCGGGCGCGCGGCACGACAGGCCGAGTAAG.GG.GG.GGGCGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTgcagc                                            <  1:827700/101‑1 (MQ=25)
                                                        cTGGCCCAGCGCCGCGGCGGGCGCGCGGCACGACAGGCCGAGTAAG.GG.GG.GGGCGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCaa                                          >  1:3550791/1‑101 (MQ=25)
                                                        cTGGCCCAGCGCCGCGGCGGGCGCGCGGCACGACAGGCCGAGTAAG.GG.GG.GGGCGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCaa                                          <  1:5399307/101‑1 (MQ=25)
                                                         tGGCCCAGCGCCGCGGCGGGCGCGCGGCACGACAGGCCGAGTAAG.GG.GG.GGGCGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAg                                         <  1:3345514/101‑1 (MQ=21)
                                                         tGGCCCAGCGCCGCGGCGGGCGCGCGGCACGACAGGCCGAGTAAG.GG.GG.GGGCGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAg                                         <  2:3350789/101‑1 (MQ=21)
                                                         tGGCCCAGCGCCGCGGCGGGCGCGCGGCACGACAGGCCGAGTAAG.GG.GG.GGGCGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAg                                         <  2:1348650/101‑1 (MQ=21)
                                                         tGGCCCAGCGCCGCGGCGGGCGCGCGGCACGACAGGCCGAGTAAG.GG.GG.GGGCGCTGCCCCCCGTCCCCAAGCGGACTCCCCCCGGGATATTTGCAGCAAg                                         <  1:5022385/101‑1 (MQ=21)
                                                           gCCCAGCGCCGCGGCGGGCGCGCGGCACGAGAGGCCGAGTAAG.GG.GG.GGGCGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGAt                                       <  1:5541737/101‑1 (MQ=21)
                                                           gCCCAGCGCCGCGGCGGGCGCGCGGCACGACAGGCCGAGTAAG.GG.GG.GGGCGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGAt                                       <  1:954978/101‑1 (MQ=21)
                                                            cccAGCGCCGCGGCGGGCGCGCGGCACGACAGGCCGAGTAAG.GG.GG.GGGCGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATg                                      <  1:585173/101‑1 (MQ=21)
                                                            cccAGCGCCGCGGCGGGCGCGCGGCACGACAGGCCGAGTAAG.GG.GG.GGGCGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATg                                      <  1:6930124/101‑1 (MQ=21)
                                                              cAGCGCCGCGGCGGGCGCGCGGCACGACAGGCCGAGTAAG.GG.GG.GGGCGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGaa                                    <  1:228495/101‑1 (MQ=21)
                                                              cAGCGCCGCGGCGGGCGCGCGGCACGACAGGCCGAGTAAG.GG.GG.GGGCGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGaa                                    <  2:525570/101‑1 (MQ=21)
                                                              cAGCGCCGCGGCGGGCGCGCGGCACGACAGGCCGAGTAAG.GG.GG.GGGCGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGaa                                    <  2:1803530/101‑1 (MQ=21)
                                                               aGCGCCGCGGCGGGCGCGCGGCACGACAGGCCGAGTAAG.GG.GG.GGGCGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGaaa                                   <  2:6685182/101‑1 (MQ=21)
                                                               aGCGCCGCGGCGGGCGCGCGGCACGACAGGCCGAGTAAG.GG.GG.GGGCGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGaaa                                   <  1:3984449/101‑1 (MQ=21)
                                                                 ggccgcGGCGGGCGCGCGGCACGACAGGCCGAGTAAG.GG.GG.GGGCGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAAcc                                 <  2:299424/100‑1 (MQ=21)
                                                                 cgcggcGGCGGGCGCGCGGCACGACAGGCCGAGTAAG.GG.GG.GGGCGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAAcc                                 <  2:4941203/101‑1 (MQ=18)
                                                                 cgccgcGGCGGGCGCGCGGCACGACAGGCCGAGTAAG.GG.GG.GGGCGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAAcc                                 <  1:5471404/101‑1 (MQ=21)
                                                                  gccgcGGCGGGCGCGCGGCACGACAGGCCGAGTAAG.GG.GG.GGGCGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCg                                <  1:6574276/101‑1 (MQ=21)
                                                                    cgggcggggcgCGCGGC.CGGCCGGCCGGGTAAGGGG.GG.GGGCGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGa                              <  2:4891201/95‑1 (MQ=14)
                                                                     gcggcgGGCGCGCGGCACGACAGGCCGAGTAAG.GG.GG.GGGCGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGAt                             <  2:4510851/101‑1 (MQ=21)
                                                                       ggcgggcgCGCGGCACGACAGGCCGAGTAAG.GG.GG.GGGCGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGc                           >  1:2034279/1‑101 (MQ=21)
                                                                       ggcgggcgCGCGGCACGACAGGCCGAGTAAG.GG.GG.GGGCGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGc                           <  1:3737208/101‑1 (MQ=21)
                                                                          gggcgCGCGGCACGACAGGCCGAGTAAG.GG.GG.GGGCGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCtt                        <  1:3236547/101‑1 (MQ=21)
                                                                            gcgcgcGGCACGACAGGCCGAGTAAG.GG.GG.GGGCGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTc                      <  2:1128414/101‑1 (MQ=21)
                                                                             cgcgcgGCACGACAGGCCGAGTAAG.GG.GG.GGGCGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCa                     <  2:1249204/101‑1 (MQ=21)
                                                                              gcgggGCCCGACAGGCCGAGTAAG.GG.GG.GGGCGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCAt                    <  1:6050482/101‑1 (MQ=18)
                                                                               cgcgGCACGACAGGCCGAGTAAG.GG.GG.GGGCGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATc                   <  1:4354939/101‑1 (MQ=21)
                                                                                 cgGCACGACAGGCCGAGTAAG.GG.GG.GGGCGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCtt                 <  2:72494/101‑1 (MQ=21)
                                                                                 cgGCACGACAGGCCGAGTAAG.GG.GG.GGGCGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCtt                 >  2:4226979/1‑101 (MQ=21)
                                                                                  ggCACGACAGGCCGAGTAAG.GG.GG.GGGCGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTg                <  2:6596797/101‑1 (MQ=21)
                                                                                    cACGACAGGCCGAGTAAG.GG.GG.GGGCGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGcc              <  2:6470756/101‑1 (MQ=21)
                                                                                     ccGACAGGCCGAGTAAG.GG.GG.GGGCGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGccc             <  2:1995956/100‑1 (MQ=21)
                                                                                      cGACAGGCCGAGTAAG.GG.GG.GGGCGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCg            <  2:4467361/101‑1 (MQ=21)
                                                                                       gACAGGCCGAGTAAG.GG.GG.GGGCGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGa           >  2:6135284/1‑101 (MQ=21)
                                                                                          aGGCCGAGTAAG.GG.GG.GGGCGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAat        >  1:316536/1‑101 (MQ=21)
                                                                                          aGGCCGAGTAAG.GG.GG.GGGCGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAat        <  1:4587345/101‑1 (MQ=21)
                                                                                           ggCCGAGTAAG.GG.GG.GGGCGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAata       <  2:2762090/101‑1 (MQ=21)
                                                                                           ggCCGAGTAAG.GG.GG.GGGCGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAata       <  1:3481457/101‑1 (MQ=21)
                                                                                           ggCCGAGTAAG.GG.GG.GGGCGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAata       <  2:3722321/101‑1 (MQ=21)
                                                                                            gCCGAGTAAG.GG.GG.GGGCGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAatat      <  1:4924078/101‑1 (MQ=21)
                                                                                            gCCGAGTAAG.GG.GG.GGGCGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAatat      <  1:2094458/101‑1 (MQ=21)
                                                                                             ccGATTAAG.GG.GG.GGGCGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATc     <  2:40144/101‑1 (MQ=18)
                                                                                             ccGAGTAAG.GG.GG.GGGCGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATc     <  2:5613774/101‑1 (MQ=21)
                                                                                             ccGAGTAAG.GG.GG.GGGCGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATc     <  2:535315/101‑1 (MQ=21)
                                                                                             ccGAGTAA..GG.GGTGAGCGCTGCCCCCCGTCCGCAAGCGTACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATc     <  2:2836554/101‑1 (MQ=21)
                                                                                              cGAGTAAG.GG.GG.GGGCGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATcc    <  1:4164654/101‑1 (MQ=21)
                                                                                                aGTAAG.GG.GG.GGGCGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCg  <  1:1425872/101‑1 (MQ=18)
                                                                                                     |                                                                                                  
TGAGCCCCGAAGAGCGCCGCGCGGTGATGGCCGAGCTGCCCACCCGCGCGCCGCGTCTGGCCCAGCGCCGCGGCGGGCGCGCGGCACGACAGGCCGAGTAA..GG.GG.GGGCGCTGCCCCCGTTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCG  >  NC_014034/2493852‑2494047

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