Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_014034 2,494,015 C→G 3.4% intergenic (+63/‑179) RCAP_rcc02319 → / → RCAP_rcc02320 hypothetical protein/hypothetical protein

Read alignment evidence...
  seq id position change freq score reads annotation genes product
*NC_0140342,494,0150C→G 3.4% 16.4 ‑0.3 ‑0.2325intergenic (+63/‑179)RCAP_rcc02319/RCAP_rcc02320hypothetical protein/hypothetical protein
Reads supporting (aligned to +/- strand):  new base (6/5):  ref base (146/168):  total (152/173)
Fisher's exact test strand distribution p-value = 7.61e-01
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 5.81e-01

Only 100 of 356 total aligned reads displayed.

GCGCCGCGGCGGGCGCGCGGCACGACAGGCCGAGTAA.GGGGGGG.CGCTGCCCCCGTTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGC..G.AG.CAGAAGGGGGGCAAGCG.CCCCCTGCCCCGGCGGTCCGTGGCGGCGCCTTGA  >  NC_014034/2493916‑2494114
                                                                                                     |                                                                                                         
gcgcCGCGGCGGGCGCGCGGCACGACAGGCCGAGTAAGGGGGGGG.CGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAAc                                                                                                           >  2:7116961/1‑101 (MQ=21)
 cgccgcGGCGGGCGCGCGGCATGACAGGCCGAGTAAGGGGGGGG.CGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAAcc                                                                                                          <  1:5745430/101‑1 (MQ=18)
  gccgcGGCGGGCGCGCGGCACGACAGGCCGAGTAAGGGGGGGG.CGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCg                                                                                                         <  2:7460995/101‑1 (MQ=21)
   ccgcGGCGGGCGCGCGGCACGACAGGCCGAGTAAGGGGGGGG.CGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCgg                                                                                                        <  1:935108/101‑1 (MQ=21)
    cgcgGCGGGCGCGCGGCACGACAGGCCGAGTAAGGGGGGGG.CGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGa                                                                                                       <  1:154616/101‑1 (MQ=21)
      cggcggGCGCGCGGCACGACAGGCCGGGTAAGGGGGGGG.CGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGGAAGATGAAACCGGATg                                                                                                     >  2:5961040/1‑101 (MQ=18)
      cggcggGCGCGCGGCACGACAGGCCGAGTAAGGGGGGGG.CGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATg                                                                                                     <  2:2095322/101‑1 (MQ=21)
      cggcggGCGCGCGGCACGACAGGCCGAGTAAGGGGGGGG.CGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATg                                                                                                     <  1:311754/101‑1 (MQ=21)
          gggcgCGCGGCACGACAGGCCGAGTAAGGGGGGGG.CGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCtt                                                                                                 >  1:5995945/1‑101 (MQ=21)
          gggcgCGCGGCACGACAGGCCGAGTAAGGGGGGGG.CGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCtt                                                                                                 <  2:4257301/101‑1 (MQ=21)
           ggcgCGCGGCACGACAGGCCGAGTAAGGGGGGGG.CGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCttt                                                                                                <  2:2465632/101‑1 (MQ=21)
           ggcgCGCGGCACGACAGGCCGAGTAAGGGGGGGG.CGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCttt                                                                                                <  2:92579/101‑1 (MQ=21)
            gcgcgcGGCACGACAGGCCGAGTAAGGGGGGGG.CGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGGTATTTGCAGCAAGATGAAACCGGCTGCCTTTc                                                                                               >  1:1802607/1‑101 (MQ=18)
             cgcgcgGCACGACAGGCCGAGTAAGGGGGGGG.CGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCa                                                                                              >  1:4839243/1‑101 (MQ=21)
               cgcgGCACGACAGGCCGAGTAAGGGGGGGG.CGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATc                                                                                            <  1:7733596/101‑1 (MQ=21)
               cgcgGCACGACAGGCCGAGTAAGGGGGGGG.CGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATc                                                                                            <  1:2578386/101‑1 (MQ=21)
                gcgGCACGACAGGCCGAGTAAGGGGGGGG.CGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCt                                                                                           >  1:2829069/1‑101 (MQ=21)
                 cgGCACGACAGGCCGAGTAAGGGGGGGG.CGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCtt                                                                                          <  2:3788717/101‑1 (MQ=21)
                 cgGCACGACAGGCCGAGTAAGGGGGGGG.CGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCtt                                                                                          <  2:4001678/101‑1 (MQ=21)
                 cgGCACGACAGGCCGAGTAAGGGGGGGG.CGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCtt                                                                                          <  2:6857187/101‑1 (MQ=21)
                      cggcAGGCCGAGTAAGGGGGGGG.CGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCg                                                                                     >  2:2774401/4‑101 (MQ=18)
                      cGACAGGCCGAGTAAGGGGGGGG.CGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCg                                                                                     >  1:5141554/1‑101 (MQ=21)
                      cGACAGGCCGAGTAAGGGGGGGG.CGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCg                                                                                     >  1:1210863/1‑101 (MQ=21)
                      cGACAGGCCGAGTAAGGGGGGGG.CGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCg                                                                                     >  1:8030789/1‑101 (MQ=21)
                      cGACAGGCCGAGTAAGGGGGGGG.CGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCg                                                                                     >  1:6788868/1‑101 (MQ=21)
                       gACAGGCCGAGTAAGGGGGGGG.CGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGa                                                                                    <  2:1527735/101‑1 (MQ=21)
                        aCAGGCCGAGTAAGGGGGGGG.CGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGaa                                                                                   <  1:8498056/101‑1 (MQ=21)
                         cAGGCCGAGTAAGGGGGGGG.CGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGaaa                                                                                  >  2:6323861/1‑101 (MQ=21)
                         cAGGCCGAGTAA.GGGGGGGTCGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGaaa                                                                                  <  1:3587970/101‑1 (MQ=21)
                          aGGCCGAGTAAGGGGGGGG.CGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTAATCTTGCCCGAAat                                                                                 >  2:3586766/1‑101 (MQ=18)
                           ggCGGAGTAAGGGGGGGG.CGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAata                                                                                <  1:6291246/101‑1 (MQ=18)
                           ggCCGAGTAAGGGGGGGG.CGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAata                                                                                >  2:6446711/1‑101 (MQ=21)
                           ggCCGAGTAAGGGGGGGG.CGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAata                                                                                >  2:7981973/1‑101 (MQ=21)
                                 gtaagggggggg.CGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGGTATTTGCAGCAAAATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCgg                                                                          >  2:5664204/6‑101 (MQ=17)
                                 gtaagggggggg.CGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCGTTCATCTTGCCCGAAATATCCCgg                                                                          >  2:5025270/6‑101 (MQ=18)
                                  taagggggggg.CGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCggg                                                                         >  2:7402294/5‑101 (MQ=255)
                                   aagggggggg.CGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCgggg                                                                        <  1:3192936/98‑1 (MQ=255)
                                   aagggggggg.CGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCgggg                                                                        >  2:3051073/4‑101 (MQ=255)
                                   aagggggggc.CGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCcggg                                                                        <  1:8801286/98‑1 (MQ=18)
                                     gggggggg.CGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGAGGAAACCGGATGCCTTTCATCTTGCCCGAAATATCACggcggg                                                                      >  1:8301384/2‑101 (MQ=18)
                                      ggggggg.CGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCggggggg                                                                     >  1:5378965/1‑100 (MQ=255)
                                       gggggg.CGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTg                                                                    <  2:1946823/101‑1 (MQ=255)
                                       gggggg.CGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTg                                                                    <  1:3931415/101‑1 (MQ=255)
                                        ggggg.CGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGa                                                                   <  2:6127768/101‑1 (MQ=255)
                                        ggggg.CGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGa                                                                   <  1:1052772/101‑1 (MQ=255)
                                        ggggg.CGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGa                                                                   <  1:7696743/101‑1 (MQ=255)
                                        ggggg.CGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATAGCCCGGGGGGTGa                                                                   <  1:8155471/101‑1 (MQ=255)
                                          ggg.CGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.t                                                                <  1:99830/101‑1 (MQ=255)
                                          ggg.CGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.t                                                                <  2:8176758/101‑1 (MQ=255)
                                          ggg.CGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.t                                                                <  2:7193470/101‑1 (MQ=255)
                                          ggg.CGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.t                                                                >  2:5189108/1‑101 (MQ=255)
                                           gg.CGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.Tg                                                               <  1:7297136/101‑1 (MQ=255)
                                              cgcTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGc                                                             >  1:4154401/1‑101 (MQ=255)
                                              cgcTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGGGAA.TGGc                                                             >  1:6159433/1‑101 (MQ=255)
                                               gcTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCGGCGGGTGAA.TGGC..g                                                          >  2:700788/1‑101 (MQ=255)
                                                   cccccG.TCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGC..G.aagcc                                                    <  1:9107171/101‑5 (MQ=255)
                                                       ccgtCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGC..G.aagccaga                                                 >  2:3327175/4‑94 (MQ=255)
                                                       ccgtCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGC..G.aagccaga                                                 <  2:2574335/98‑8 (MQ=255)
                                                        cgtCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGCATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGC..G.aagccagag                                                <  1:8955443/99‑9 (MQ=21)
                                                         gtCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGC..G.aagccagagg                                               <  2:1152574/100‑10 (MQ=255)
                                                          tCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.GGGG..G.aagccagaggg                                              >  1:2373765/1‑91 (MQ=18)
                                                                   ggACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGC..GAAG.CCAGAGGGGGGCagcgc                                     >  1:3768029/1‑97 (MQ=17)
                                                                    gACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGC..GAAG.CCAGAGGGGGGC.AGCG.cc                                  >  2:3447109/1‑101 (MQ=18)
                                                                    gACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGC..GAAG.CCAGAGGGGGG.AAGCG.cc                                  >  2:264577/1‑101 (MQ=18)
                                                                    gACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGGGAA.TGGC..GAAG.CCAGAGGGGGGC.AGCG.cc                                  >  2:4464824/1‑101 (MQ=17)
                                                                      cTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGC..GAAGCCAG.AGGGGGGC.AGCG.cccc                                <  2:882595/101‑1 (MQ=17)
                                                                          ccccGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGGGAA.TGGC..GAAG.CCAGAGGGGGGC.AGCG.cccccccg                            >  1:6917471/1‑98 (MQ=17)
                                                                           cccGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGC..GAAG.CCAGAGGGGGGC.AGCG.cccccctgc                           <  1:3318694/101‑5 (MQ=17)
                                                                            ccGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGC..GAAG.CCAGAGGGGGGCACCGC.CCCCCTGcc                           >  1:4531221/1‑101 (MQ=17)
                                                                              gggATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGC..GAAG.CCAGAGGGGGGC.AGCGCCCCCCTGcccc                         >  2:1083124/1‑101 (MQ=18)
                                                                              gggATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGC..GAAG.CCAGAGGGGGGC.AGCGCCCCCCTGcccc                         <  2:2158618/101‑1 (MQ=17)
                                                                                gATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGC..GAAG.CCAGAGGGGGGC.AGCGCCCCCCTGCCCcgg                       <  2:314687/101‑1 (MQ=17)
                                                                                 atatTTGCGGCAAGATGAAGGGTCATG.ATTTCATCTTGCCCGAAATATCCC.AGGGGTGAATTGGCCAA.AG.GCCAAGAGGGGGCAGCG.CCCCCTGCCCcgcc                      <  2:8953323/101‑3 (MQ=2)
                                                                                 atatTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGC..GAAG.CCAGAGGGGGGC.AGCGCCCCCCTGCCCcggc                      >  2:4744540/1‑101 (MQ=17)
                                                                                  tatTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGC..GAAG.CCAGAGGGGGGC.AGCGCCCCCCTGCCCcggcg                     >  1:6499480/1‑101 (MQ=18)
                                                                                  tatTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGC..GAAG.CCAGAGGGGGGC.AGCGCCCCCCTGCCCcggcg                     >  2:7766288/1‑101 (MQ=18)
                                                                                  tatTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGCGAA.TGGC..GAAG.CCAGAGGGGGGC.AGCGCCCCCCTGCCCcggcg                     >  2:5374617/1‑101 (MQ=17)
                                                                                   atTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGC..GAAG.CCAGAGGGGGGCAGCCC.CCCCCTGCCCcggcgg                    >  2:5454249/1‑101 (MQ=18)
                                                                                     ttGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGC..GAAG.CCAGAGGGAGGC.AGCGCCCCCCTGCCCCGGCGGTc                  >  1:5013433/1‑101 (MQ=17)
                                                                                      tGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGC..GAAG.CCAGAGGGGGGC.CGCGCCCCCCTGCCCcggcggccg                 >  2:3906303/1‑98 (MQ=17)
                                                                                      tGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGC..GAAG.CCAGAGCGGGGC.AGCGCCCCCCTGCCCCCGCGGTcc                 >  2:1769121/1‑101 (MQ=17)
                                                                                       gcagcaAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGC..GAAG.CCAGAGGGGGGC.AGCGCCCCCCTGCCCCGGCGGTCCg                >  2:1361329/1‑101 (MQ=18)
                                                                                       gcagcaAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGC..GAAG.CCAGAGGGGGGC.AGCGCCCCCCTGCCCCGGCGGTCCg                >  2:5384668/1‑101 (MQ=18)
                                                                                       gcagcaAGATGAAACCGGATGCCTTTCACCTTGCCCGAAATATCCCGGGGGGTGAA.TGGC..GAAG.CCAGAGGGGGGC.AGCGCCCCCCTGCCCCGGCGGTCCg                <  1:6395668/101‑1 (MQ=17)
                                                                                        cagcaAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGC..GAAG.CCAGAGGGGGGC.AGCGCCCCCCTGCCCCGGCGGTCCGt               >  2:853078/1‑101 (MQ=17)
                                                                                        cagcaAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGC..GAAG.CCAGAGGGGGGC.AGCGCCCCCCTGCCCCGGCGGTCCGt               <  1:1876870/101‑1 (MQ=17)
                                                                                        cagcaAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGC..GAAG.CCAGAGGGGGGC.AGCGCCCCCCTGCCCCGGCGGTCCGt               <  1:3359253/101‑1 (MQ=17)
                                                                                          gcaAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATAGCCCGGGGGGTGAA.TGGC..GAAG.CCAGAGGGGGGC.AGCGCCCCCCTGCCCCGGCGGTCCGTgg             <  1:1972714/101‑1 (MQ=17)
                                                                                           caAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGC..GAAG.CCAGAGGGGGGC.AGCGCCCCCCTGCCCCGGCGGTCCGTggc            <  1:7074252/101‑1 (MQ=17)
                                                                                           caAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGC..GAAG.CCAGAGGGGGGC.AGCGCCCCCCTGCCCCGGCGGTCCGTggc            >  1:132532/1‑101 (MQ=18)
                                                                                            aaGATGAAACCGGATGCCTTTCATGTTGCCCGAAATATCCCGGGGGGTGAA.TGGC..GAAG.CCAGAGGGGGGC.AGCGCCCCCCTGCCCCGGCGGTCCGTggcg           >  2:4385111/1‑101 (MQ=17)
                                                                                             aGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGC..GAAG.CCAGAGGGGGGC.AGCGCCCCCCTGCCCCGGCGGTCCGTggcgg          >  2:3636282/1‑101 (MQ=18)
                                                                                                 gAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGC..GAAG.CCAGAGGGGGGC.AGCGCCCCCCTGCCCCGGCGGTCCGTGGCGGCGcc      >  2:6677205/1‑101 (MQ=18)
                                                                                                  aaaCCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGC..GAAG.CCAGAGGGGGGC.AGCGCCCCCCTGCCCCGGCGGTCCGTGGCGGCGCCt     <  2:2828971/101‑1 (MQ=17)
                                                                                                  aaaCCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGGGAA.TGGC..GAAG.CCAGAGGGGGGCAGCCC.CCCCCTGCCCCGCCGGCCGGTGGGGGCGCCt     >  1:1127422/1‑101 (MQ=12)
                                                                                                   aaCCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGC..GAAG.CCAGAGGGGGGC.AGCGCCCCCCTGCCCCGGCGGTCCGTGGCGGCGCCtt    >  2:766081/1‑101 (MQ=18)
                                                                                                   aaCCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGC..GAAG.CCAGAGGGGGGC.AGCGCCCCCCTGCCCCGGCGGTCCGTGGCGGCGCCtt    >  1:4189766/1‑101 (MQ=18)
                                                                                                   aaCCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGC..GAAG.CCAGAGGGGGGC.AGCGCCCCCCTGCCCCGGCGGTCCGTGGCGGCGCCtt    <  2:4600707/101‑1 (MQ=17)
                                                                                                    aCCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGC..GAAG.CCAGAGGGGGGC.AGCGCCCCCCTGCCCCGGCGGTCCGTGGCGGCGCCTTg   <  1:7982313/101‑1 (MQ=17)
                                                                                                     ccGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGC..GAAG.CCAGAGGGGGGC.AGCGCCCCCCTGCCCCGGCGGTCCGTGGCGGCGCCTTGa  <  1:7402618/101‑1 (MQ=17)
                                                                                                     |                                                                                                         
GCGCCGCGGCGGGCGCGCGGCACGACAGGCCGAGTAA.GGGGGGG.CGCTGCCCCCGTTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGC..G.AG.CAGAAGGGGGGCAAGCG.CCCCCTGCCCCGGCGGTCCGTGGCGGCGCCTTGA  >  NC_014034/2493916‑2494114

Base quality scores:  ATCG  < 3 ≤  ATCG  < 16 ≤  ATCG  < 25 ≤  ATCG  < 33 ≤  ATCG  < 41 ≤  ATCG