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

Read alignment evidence...
  seq id position change freq score reads annotation genes product
*NC_0140342,494,0160C→G 3.4% 17.2 ‑0.3 ‑0.2325intergenic (+64/‑178)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 = 6.44e-01

Only 100 of 356 total aligned reads displayed.

CGCCGCGGCGGGCGCGCGGCACGACAGGCCGAGTAA.GGGGGGG..CGCTGCCCCCGTTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGC..G.AG.CAGA.AGGGGGGCAAGCG.CCCCCTGCCCCGGCGGTCCGTGGCGGCGCCTTGAC  >  NC_014034/2493917‑2494115
                                                                                                      |                                                                                                          
cgccgcGGCGGGCGCGCGGCACGACAGGCCGAGTAAGGGGGGGG..CGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAAcc                                                                                                            <  1:6363469/101‑1 (MQ=21)
cgccgcGGCGGGCCCGCGGCACGACAGGCCGAGTAAGGGGGGGG..CGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAAcc                                                                                                            <  2:5317057/101‑1 (MQ=18)
 gccgcGGCGGGCGCGCGGCACGACAGGCCGAGTAAGGGGGGGG..CGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCg                                                                                                           <  1:6885361/101‑1 (MQ=21)
   cgcgGCGGGCGCGCGGCACGACAGGCCGAGTAAGGGGGGGG..CGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGa                                                                                                         <  1:154616/101‑1 (MQ=21)
    gcggcgGGGGCGCGGCACGACAGGCCGGGTAAGGGGGGGG..CGCTGCGCCCCCTCCCGAAGCGGCCTTCCCCCGGGGAAGTTGCAGGACGATCAAACCGGag                                                                                                        >  2:162524/1‑100 (MQ=16)
    gcggcgGGCGCGCGGCACGACAGGCCGAGTAAGGGGGGGG..CGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGAt                                                                                                        <  1:8366429/101‑1 (MQ=21)
      ggcgggcgCGCGGCACGACAGGCCGAGTAAGGGGGGGG..CGCTGCCCCCCGTCCGCAAGGGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGc                                                                                                      >  1:5295295/1‑101 (MQ=21)
      ggcgggcgCGCGGCACGACAGGCCGAGTAAGGGGGGGG..CGCTGCCCCCCGTCCGCAAGAGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGc                                                                                                      <  1:4318653/101‑1 (MQ=21)
       gcgggcgCACGGCACGACAGGCCGAGTAAGGGGGGGG..CGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGcc                                                                                                     <  1:8601610/101‑1 (MQ=18)
         gggcgCGCGGCACGACAGGCCGAGTAAGGGGGGGG..CGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCtt                                                                                                   <  2:1045837/101‑1 (MQ=21)
         gggcgCGCGGCACGACAGGCCGAGTAAGGGGGGGG..CGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCtt                                                                                                   <  2:2294198/101‑1 (MQ=21)
         gggcgCGCGGCACGACAGGCCGAGTAAGGGGGGGG..CGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCtt                                                                                                   <  2:4257301/101‑1 (MQ=21)
          ggcgCGCGGCACGACAGGCCGAGTAAGGGGGGGG..CGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCttt                                                                                                  <  2:92579/101‑1 (MQ=21)
           gcgcgcGGCACGACAGGCCGAGTAAGGGGGGGG..CGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGGTATTTGCAGCAAGATGAAACCGGCTGCCTTTc                                                                                                 >  1:1802607/1‑101 (MQ=18)
           gcgcgcGGCACGACAGGCCGAGTAAGGGGGGGG..CGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTc                                                                                                 <  1:5943665/101‑1 (MQ=21)
           gcgcgcGGCACGACAGGCCGAGTAAGGGGGGGG..CGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTc                                                                                                 <  1:3715740/101‑1 (MQ=21)
           gcgcgcGGCACGACAGGCCGAGTAAGGGGGGGG..CGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTc                                                                                                 <  1:786711/101‑1 (MQ=21)
            cgcgcgGCACGACAGGCCGAGTAAGGGGGGGG..CGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCa                                                                                                <  2:3501035/101‑1 (MQ=21)
             gcgcgGCACGACAGGCCGAGTAAGGGGGGGG..CGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCAt                                                                                               <  2:2861940/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)
                 ggCACGACAGGCCGAGTAAGGGGGGGG..CGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTg                                                                                           >  1:8177107/1‑101 (MQ=21)
                 ggCACGACAGGCCGAGTAAGGGGGGGG..CGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTg                                                                                           <  2:5643495/101‑1 (MQ=21)
                 ggCACGACAGGCCGAGTAAGGGGGGGG..CGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCCTTCATTTTg                                                                                           >  1:5643702/1‑101 (MQ=18)
                 ggCACAACAGGCCGAGTAAGGGGGGGG..CGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTg                                                                                           <  1:6489526/101‑1 (MQ=18)
                   cACGACAGGCCGAGTAAGGGGGGGG..CGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGcc                                                                                         <  2:5947047/101‑1 (MQ=21)
                   cACGACAGGCCGAGTAAGGGGGGGG..CGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGcc                                                                                         >  1:6997411/1‑101 (MQ=21)
                   cACGACAGGCCGAGTAAGGGGGGGG..CGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGcc                                                                                         >  2:7296832/1‑101 (MQ=21)
                    aCGACAGGCCGAGTAAGGGGGGGG..CGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCCTCTTGccc                                                                                        >  2:1895580/1‑101 (MQ=18)
                    aCGACAGGCCGAGTAAGGGGGGGG..CGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGccc                                                                                        >  1:3225756/1‑101 (MQ=21)
                     cGACAGGCCGAGTAAGGGGGGGG..CGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCg                                                                                       >  1:1210863/1‑101 (MQ=21)
                     cGACAGGCCGAGTAAGGGGGGGG..CGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCg                                                                                       <  1:6124724/101‑1 (MQ=21)
                     cGACAGGCCGAGTAAGGGGGGGG..CGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCg                                                                                       >  1:5141554/1‑101 (MQ=21)
                      gACAGGCCGAGTAAGGGGGGGG..CGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGa                                                                                      >  2:7073967/1‑101 (MQ=21)
                       aCAGGCCGAGTAAGGGGGGGG..CGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGaa                                                                                     >  2:6218555/1‑101 (MQ=21)
                        cAGGCCGAGTAAGGGGGGGG..CGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGaaa                                                                                    >  2:7452255/1‑101 (MQ=21)
                         aGGCCGAGTAAGGGGGGGG..CGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAat                                                                                   >  2:3357068/1‑101 (MQ=21)
                          ggCCGAGTAAGGGGGGGG..CGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAata                                                                                  <  2:6865874/101‑1 (MQ=21)
                           gCCGGGTAA.GGGGGGGGCCCCTCCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAata                                                                                  <  1:8946637/101‑1 (MQ=18)
                           gCCGAGTAAGGGGGGGG..CGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAatat                                                                                 <  2:2527562/101‑1 (MQ=21)
                           gCCGAGTAAGGGGGGGG..CGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAatat                                                                                 <  1:934760/101‑1 (MQ=21)
                           gCCGAGTAAGGGGGGGG..CGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAatat                                                                                 <  2:2140873/101‑1 (MQ=21)
                             cGAGTAAGGGGGGGG..CGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATcc                                                                               <  2:812952/101‑1 (MQ=21)
                              gAGTAAGGGGGGGG..CGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATccc                                                                              >  1:8974454/1‑101 (MQ=18)
                               aGTAAGGGGGGGG..CGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCg                                                                             <  1:6236828/101‑1 (MQ=18)
                                gtaagggggggg..CGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCgg                                                                            >  2:7827681/6‑101 (MQ=255)
                                  aagggggggg..CGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCgggg                                                                          <  2:5049963/98‑1 (MQ=255)
                                    tcttccgatccgcTGCCCCCGGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCcggg                                                                          <  1:7232375/91‑1 (MQ=21)
                                    gggggggg..CGCTGCCCCCCGTCCGCAAGCGGACTCCCCACGGGCTATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCgggggg                                                                        >  1:4038957/2‑101 (MQ=18)
                                       ggggg..CGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGa                                                                     <  1:1052772/101‑1 (MQ=255)
                                         ggg..CGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGaag                                                                   >  2:1114671/1‑100 (MQ=255)
                                          gg..CGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.Tg                                                                 <  2:3220547/101‑1 (MQ=255)
                                           g..cgcTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.Tgg                                                                <  2:2595400/101‑1 (MQ=255)
                                               gcTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGC..g                                                            <  1:3322414/101‑1 (MQ=255)
                                                 tGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGC..G.aa                                                         <  1:7350725/101‑2 (MQ=255)
                                                  gCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGC..G.aag                                                        <  1:984402/101‑3 (MQ=255)
                                                   cccccG.TCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGC..G.aagcc                                                      <  1:772453/101‑5 (MQ=255)
                                                      cccgtCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGGGAA.TGGC..G.aagccag                                                    >  2:1972654/5‑95 (MQ=18)
                                                        cgtCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGC..G.aagccagag                                                  <  1:8117465/99‑9 (MQ=255)
                                                          tCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.GGGG..G.aagccagaggg                                                >  1:2373765/1‑91 (MQ=18)
                                                                    gACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGC..GAAG.CCAG.AGGGGGGC.AGCG.cc                                   >  2:7126503/1‑101 (MQ=18)
                                                                    gACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGGGAA.TGGC..GAAG.CCAG.AGGGGGGC.AGCG.cc                                   >  2:4464824/1‑101 (MQ=17)
                                                                     aCTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGC..GAAG.CCAG.AGGGGGGC.AGCG.ccc                                  >  2:8319079/1‑101 (MQ=18)
                                                                     aCTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGC..GAAG.CCAG.AGGGGGGC.AGCG.ccc                                  >  1:489283/1‑101 (MQ=18)
                                                                      cTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGC..GAAGCCAG..AGGGGGGC.AGCG.cccc                                 <  2:882595/101‑1 (MQ=17)
                                                                      cTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGC..GAAGCCAG..AGGGGGGC.AGCG.cccc                                 <  1:3447213/101‑1 (MQ=17)
                                                                      cTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGC..GAAG.CCAG.AGGGGGGC.AGCG.cccc                                 >  2:759813/1‑101 (MQ=18)
                                                                         cccccGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGC..GAAG.CCAG.AGGGGGGC.AGCG.cccccct                              <  1:3132562/101‑3 (MQ=17)
                                                                           cccGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGC..GAAG.CCAG.AGGGGGGC.AGCG.cccccctgc                            <  2:7704352/101‑5 (MQ=17)
                                                                                gATATTTGCGGCAAGATGAAGGGGCATG.ATTTCATCTTGCCCGAAATATCCC.AGGGGTGAATTGGCCAA.AG.GCCA.AGAGGGGGCAGCG.CCCCCTGCCCcgc                        <  2:2808375/101‑2 (MQ=0)
                                                                                 atatTTGCGGCAAGATGAAGGGTCATG.ATTTCATCTTGCCCGAAATATCCC.AGGGGTGAATTGGCCAA.AG.GCCA.AGAGGGGGCAGCG.CCCCCTGCCCcgcc                       <  2:8953323/101‑3 (MQ=2)
                                                                                 atatTTGCGGCAAGATGAAGGGGCATG.ATTTCATCTTGCCCGAAATATCCC.AGGGGTGAATTGGCCAA.AG.GCCA.AGAGGGGGCAGCG.CCCCCTGCCCcgcc                       >  2:5005103/1‑99 (MQ=2)
                                                                                 atatTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGC..GAAG.CCAG.AGGGGGGC.AGCGCCCCCCTGCCCcggc                       >  2:2661722/1‑101 (MQ=17)
                                                                                   atTTGCGGCAAGATGAAGGGGCATG.ATTTCATCTTGCCCGAAATATCCC.AGGGGTGAATTGGCCAA.AG.GCCA.AGAGGGGGCAGCG.CCCCCTGCCCcgccgg                     <  2:3069142/101‑1 (MQ=2)
                                                                                   atTTGCGGCAAGATGAAGGGGCATG.ATTTCATCTTGCCCGAAATATCCC.AGGGGTGAATTGGCCAA.AG.GCCA.AGAGGGGGCAGCG.CCCCCTGCCCcgccgg                     <  1:7555424/101‑1 (MQ=2)
                                                                                   atTTGCGGCAAGATGAAGGGGCATG.ATTTCACCTTGCCCGAAATATCCC.AGGGGTGAATTGGCCAA.AG.GCCA.AGAGGGGGCAGCG.CCCCCTGCCCcgccgg                     <  1:6583533/101‑1 (MQ=2)
                                                                                   atTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGC..G.AA.GACAGAGGGGGGC.AGCGCCCCCCTGCCCcggcgg                     >  1:4170280/1‑101 (MQ=18)
                                                                                     ttGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGAGAA.TGGC..GAAG.CCAG.AGGGGGGC.AGCGCCCCCCTGCCCcggcggcc                   >  1:8007750/1‑99 (MQ=17)
                                                                                      tGCAGCAAGATGAAACCGGATGCTTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGC..GAAG.CCAG.AGGGGGGC.AGCGCCCCCCTGCCCCGGCGGTcc                  <  1:4666843/101‑1 (MQ=17)
                                                                                      tGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGC..GAAG.CCAG.AGGGGGGC.CGCGCCCCCCTGCCCcggcggccg                  >  2:3906303/1‑98 (MQ=17)
                                                                                      tGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGC..GAAG.CCAG.AGGGGGGC.AGCGCCCCCCTGCCCCGGCGGTcc                  <  2:6445306/101‑1 (MQ=17)
                                                                                       gcagcaAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGC..GAAG.CCAG.AGGGGGGC.AGCGCCCCCCTGCCCCGGCGGTCCg                 >  2:664137/1‑101 (MQ=17)
                                                                                       gcagcaAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGC..GAAG.CCAG.AGGGGGGC.AGCGCCCCCCTGCCCCGGCGGGCCg                 >  2:6821784/1‑101 (MQ=17)
                                                                                        cagcaAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGC..GAAG.CCAG.AGGGGGGC.AGCGCCCCCCTGCCCCGGCGGTCCGt                <  1:1876870/101‑1 (MQ=17)
                                                                                          gcaAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATAGCCCGGGGGGTGAA.TGGC..GAAG.CCAG.AGGGGGGC.AGCGCCCCCCTGCCCCGGCGGTCCGTgg              <  1:1972714/101‑1 (MQ=17)
                                                                                           caAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGC..GAAG.CCAG.AGGGGGGC.AGCGCCCCCCTGCCCCGGCGGTCCGTggc             >  1:132532/1‑101 (MQ=18)
                                                                                            aaGATGAAACCGGATGCCTTTCATGTTGCCCGAAATATCCCGGGGGGTGAA.TGGC..GAAG.CCAG.AGGGGGGC.AGCGCCCCCCTGCCCCGGCGGTCCGTggcg            >  2:4385111/1‑101 (MQ=17)
                                                                                              gATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGC..GAAG.CCAG.AGGGGGGC.AGCGCCCCCCTGACCCAGCGGTCCGTggcggc          >  2:600238/1‑101 (MQ=17)
                                                                                                 gAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGC..GAAG.CCAG.AGGGGGGCAGCTC.CCCCCTGCCCCGGCGGTACGTGgccgcggc       >  2:1494372/1‑99 (MQ=14)
                                                                                                 gAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGC..GAAG.CCAG.AGGGGGGC.AGCGCCCCCCTGCCCCGGCGGTCCGTGGCGGCGcc       <  2:5453497/101‑1 (MQ=17)
                                                                                                  aaaCCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGC..GAAG.CCAG.AGGGGGGC.AGCGCCCCCCTGCCCCGGCGGTCCGTGGCGGCGCCt      >  2:3606184/1‑101 (MQ=18)
                                                                                                  aaaCCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGC..GAAG.CCAG.AGGGGGGC.AGCGCCCCCCTGCCCCGGCGGTCCGTGGCGGCGCCt      <  2:1258087/101‑1 (MQ=17)
                                                                                                  aaaCCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGGGAA.TGGC..GAAG.CCAG.AGGGGGGCAGCCC.CCCCCTGCCCCGCCGGCCGGTGGGGGCGCCt      >  1:1127422/1‑101 (MQ=12)
                                                                                                   aaCCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGC..GAAG.CCAG.AGGGGGGC.AGCGCCCCCCTGCCCCGGCGGTCCGTGGCGGCGCCtt     <  2:5378781/101‑1 (MQ=17)
                                                                                                   aaCCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGC..GAAG.CCAG.AGGGGGGC.AGCGCCCCCCTGCCCCGGCGGTCCGTGGCGGCGCCtt     <  2:4600707/101‑1 (MQ=17)
                                                                                                     ccGGATGCCTTTCATCTTGCCCGAAATGTCCCGGGGGGTGAA.TGGC..GAAG.CCAG.AGGGGGGC.AGCGCCCCCCTGCCCCGGCGGTCCGTGGCGGCGCCTTGa   <  1:382543/101‑1 (MQ=17)
                                                                                                      cGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGC..GAAG.CCAG.AGGGGGGC.AGCGCCCCCCTGCCCCGGCGGTCCGTGGCGGCGCCTTGAc  <  2:3641976/101‑1 (MQ=17)
                                                                                                      cGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGC..GAAG.CCAG.AGGGGGGC.AGCGACCCCCTGCCCCGGCGGTCCGTGGCGGCGCCTTGAc  >  2:4739024/1‑101 (MQ=18)
                                                                                                      |                                                                                                          
CGCCGCGGCGGGCGCGCGGCACGACAGGCCGAGTAA.GGGGGGG..CGCTGCCCCCGTTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGC..G.AG.CAGA.AGGGGGGCAAGCG.CCCCCTGCCCCGGCGGTCCGTGGCGGCGCCTTGAC  >  NC_014034/2493917‑2494115

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