Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_014034 2,494,023 C→A 6.3% intergenic (+71/‑171) RCAP_rcc02319 → / → RCAP_rcc02320 hypothetical protein/hypothetical protein

Read alignment evidence...
  seq id position change freq score reads annotation genes product
*NC_0140342,494,0230C→A 6.3% 19.5 ‑0.6 ‑0.2173intergenic (+71/‑171)RCAP_rcc02319/RCAP_rcc02320hypothetical protein/hypothetical protein
Reads supporting (aligned to +/- strand):  new base (7/4):  ref base (81/81):  total (88/85)
Fisher's exact test strand distribution p-value = 5.36e-01
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 6.69e-01

Only 100 of 188 total aligned reads displayed.

CGGGCGCGCGGCACGACAGGCCGAGT.AA.GGGGGGGCGCTGCCCCCGTTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGC..G.A..G.CAG..AAGGGGGGCAAGCG.C.CCCCTGCCCCGGCGGTCCGTGGCGGCGCCTTGACGGGCGGG  >  NC_014034/2493925‑2494122
                                                                                                    |                                                                                                              
 gggcgCGCGGCACGACAGGCCGAGT.AAGGGGGGGGCGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCCGGATCTTTGCAGCAAGATGAAAACCGAGGCCtt                                                                                                              >  1:1629492/1‑101 (MQ=17)
   gcgcgcGGCACGACAGGCCGAGT.AAGGGGGGGGCGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTc                                                                                                            <  2:455064/101‑1 (MQ=21)
   gcgcgcGGCACGACAGGCCGAGT.AAGGGGGGGGCGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTc                                                                                                            <  1:1545937/101‑1 (MQ=21)
    cgcgcgGCACGACAGGCCGAGT.AAGGGGGGGGCGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCa                                                                                                           <  2:887464/101‑1 (MQ=21)
     gcgcgGCACGGCGGGCAGAGTGAG.GGGGGGGCGCTGCCGCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAACCCGGATGCCTTTCAt                                                                                                          <  2:4382182/101‑1 (MQ=17)
     gcgcgGCACGACAGGCCGAGT.AAGGGGGGGGCGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCAt                                                                                                          >  1:4421786/1‑101 (MQ=21)
     gcgcgGCACGACAGGCCGAGT.AAGGGGGGGGCGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCAt                                                                                                          <  2:4753806/101‑1 (MQ=21)
     gcgcgGCACGACAGGCCGAGT.AAGGGGGGGGCGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCAt                                                                                                          >  2:4475525/1‑101 (MQ=21)
     gcgcgGCACGACAGGCCGAGT.AAGGGGGGGGCGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTAAt                                                                                                          >  1:3671805/1‑101 (MQ=18)
      cgcgGCACGACAGGCCGAGT.AAGGGGGGGGCGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATc                                                                                                         <  1:1456007/101‑1 (MQ=21)
       gcgGCACGACAGGCCGAGT.AAGGGGGGGGCGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCt                                                                                                        <  2:1119874/101‑1 (MQ=21)
        cgGCACGACAGGCCGAGT.AAGGGGGGGGCGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCtt                                                                                                       >  2:2651574/1‑101 (MQ=21)
        cgGCACGACAGGCCGAGT.AAGGGGGGGGCGCTGCCCCCCGTCCCCAAGCGGACTCCCCCCCGGATATTTGCAGCCAGATGAAACCGGGTGCCTTTCATCtt                                                                                                       >  2:3335491/1‑101 (MQ=18)
         ggCACGACAGGCCGAGT.AAGGGGGGGGCGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTg                                                                                                      <  2:4007606/101‑1 (MQ=21)
             cGACAGGCCGAGT.AAGGGGGGGGCGCTGCCCCCTGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCg                                                                                                  <  2:4202202/101‑1 (MQ=21)
             cGACAGGCCGAGT.AAGGGGGGGGCGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGCAATTTGCAGCAAGATGAAACCGGATGACTTTCATCTTGCCCg                                                                                                  >  2:2847826/1‑101 (MQ=18)
                cAGGCCGAGT.AAGGGGGGGGCGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGaaa                                                                                               >  2:2623504/1‑101 (MQ=21)
                 aGGCCGAGT.AAGGGGGGGGCGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAat                                                                                              >  2:1270726/1‑101 (MQ=21)
                 aGGCCGAGT.AAGGGGGGGGCGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAat                                                                                              >  1:2901878/1‑101 (MQ=21)
                  ggCCGAGT.AAGGGGGGGGCGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAata                                                                                             <  2:1003329/101‑1 (MQ=21)
                  ggCCGAGT.AAGGGGGGGGCGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAata                                                                                             <  1:4661547/101‑1 (MQ=21)
                  ggCCGAGT.AAGGGGGGGGCGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCCGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCCTCTTGCCCGAAata                                                                                             >  1:848932/1‑101 (MQ=18)
                   gCCGAGG.AAGGGGGGGGCGCTGCCCCCCGTNCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAatat                                                                                            <  1:962273/101‑1 (MQ=18)
                      gAGT.AAGGGGGGGGCGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATccc                                                                                         >  2:3439104/1‑101 (MQ=18)
                      gAGT.AAGGGGGGGGCGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATccc                                                                                         <  1:1960456/101‑1 (MQ=18)
                      gAGT.AAGGGGGGGGCGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATccc                                                                                         >  2:3097438/1‑101 (MQ=18)
                       gaggaaggggggggCGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATccc                                                                                         <  1:305139/94‑1 (MQ=18)
                       aGT.AAGGGGGGGGCGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAGTATCCCg                                                                                        >  2:4149881/1‑101 (MQ=18)
                       aGT.AAGGGGGGGGCGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCg                                                                                        <  2:911869/101‑1 (MQ=18)
                         gtaaggggggggCGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCgg                                                                                       <  2:1688279/96‑1 (MQ=255)
                           aaggggggggCGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCgggg                                                                                     >  2:3737533/4‑101 (MQ=255)
                           aaggggggggCGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCgggg                                                                                     >  2:3228300/4‑101 (MQ=255)
                            aggggggggCGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCggggg                                                                                    <  1:4354082/99‑1 (MQ=255)
                            aggggggggCGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCggggg                                                                                    >  1:1229022/3‑101 (MQ=255)
                             ggggggggCGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCgggggg                                                                                   >  1:1430865/2‑101 (MQ=255)
                              gggggggCGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGt                                                                                  <  1:2651728/101‑1 (MQ=255)
                              gggggggCGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGCGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGt                                                                                  >  2:4551779/1‑101 (MQ=255)
                               ggggggCGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTg                                                                                 >  2:3922708/1‑101 (MQ=255)
                                gggggCGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGa                                                                                >  1:4641123/1‑101 (MQ=255)
                                 ggggCGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGaa                                                                               <  1:4475775/101‑1 (MQ=255)
                                 ggggCGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGGGaa                                                                               >  1:4682235/1‑101 (MQ=255)
                                  gggCGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.t                                                                             >  2:788864/1‑101 (MQ=255)
                                   ggCGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.Tg                                                                            <  2:1976251/101‑1 (MQ=255)
                                   ggCGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.Tg                                                                            <  2:3759132/101‑1 (MQ=255)
                                     cgcTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTAAA.TGGc                                                                          >  1:789771/1‑101 (MQ=255)
                                       cTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGC..G.a                                                                     >  2:4291802/1‑101 (MQ=255)
                                          cccccG.TCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGC..G.aagcc                                                                 <  1:4605223/101‑5 (MQ=255)
                                          cccccG.TCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGC..G.aagcc                                                                 >  2:3871363/1‑97 (MQ=255)
                                              ccgtCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGC..G.aagccaga                                                              <  1:2848002/98‑8 (MQ=255)
                                               cgtCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGC..G.aagccagag                                                             <  2:2067283/99‑9 (MQ=255)
                                                 tCCGCAAGCGGACTCCCCCCGGGAGATTTGCAGCAAGATGAAACCGGATACCTTTCATCTTGCCCGAAATATCCCGGGGGGGGAA.TGGC..G.aagccagaggg                                                           >  1:43063/1‑91 (MQ=18)
                                                    gCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGC..G.AAGG.CAG..AGggggg                                                     >  1:1601862/1‑101 (MQ=18)
                                                      aaGCGGGCTCCCGCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGC..GAA..G.CCA..GAGGGGGGCa                                                  <  1:2260694/101‑1 (MQ=18)
                                                       aGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGC..GAA..G.CCA..GAGGGGGGCag                                                 <  1:4753551/101‑2 (MQ=18)
                                                        gCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGC..GAA..G.CCA..GAGGGGGGCagc                                                >  1:2176995/1‑99 (MQ=18)
                                                          ggACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGACCGAAATATCCCGGGGGGTGAA.TGGC..GAA..G.CCA..GAGGGGGGCagcgc                                              >  1:2769520/1‑97 (MQ=17)
                                                           gACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGGGAA.TGGC..GAA..G.CCA..GAGGGGGGC.AGCG.c.c                                          >  1:2725543/1‑101 (MQ=17)
                                                            aCTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGC..GAA..G.CCA..GAGGGGGGCagcggcc                                            >  2:2348947/1‑95 (MQ=17)
                                                            aCTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGC..GAA..G.CCA..GAGGGGGGC.AGCG.c.cc                                         >  2:2414556/1‑101 (MQ=18)
                                                             cTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGC..GAA..GCCAG...AGGGGGGC.AGCG.c.ccc                                        <  1:402027/101‑1 (MQ=17)
                                                             cTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGC..GAA..GCCAG...AGGGGGGC.AGCG.c.ccc                                        <  2:2954777/101‑1 (MQ=17)
                                                                cccccGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGC..GAA..G.CCA..GAGGGGGGC.AGCG.c.ccccct                                     <  2:848858/101‑3 (MQ=17)
                                                                   ccGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGC..GAA..G.CCA..GAGGGGGGC.AGCG.c.ccccctgcc                                  <  2:4018385/101‑6 (MQ=17)
                                                                     gggATATTTGCGGCAAGATGAAGGGGCATG.ATTTCATCTTGCCCGAAATATCCC.AGGGGTGAATTGGC..C.A..A.AGGCCAAGAGGGGGCAGCG.C.CCCCTGcccc                                 <  1:4206111/101‑1 (MQ=0)
                                                                     gggATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGC..GAA..G.CCA..GAGGGGGGC.AGCGCC.CCCCTGcccc                                 >  1:3846964/1‑101 (MQ=17)
                                                                     gggATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGC..GAA..G.CCA..GAGGGGGGC.AGCGCC.CCCCTGcccc                                 >  1:3588006/1‑101 (MQ=17)
                                                                     gggATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGC..GAA..G.CCA..GAGGGGGGC.AGCGCC.CCCCTGcccc                                 <  2:191057/101‑1 (MQ=17)
                                                                      ggATATTTGCGGCAAGATGAAGGGGCATG.ATTTCATCTTGCCCGAAATATCCC.AGGGGTGAATTGGCCAA.A..G.GCC..AAGAGGGGGCAGCG.C.CCCCTGCCCcg                                <  2:1362966/101‑1 (MQ=0)
                                                                      ggATATTTGCGGCAAGATGAAGGGGCATG.ATTTCATCTTGCCCGAAATATCCC.AGGGGTGAATTGGCCAA.A..G.GCC..AAGAGGGGGCAGCG.C.CCCCTGCCCcg                                >  1:659612/1‑101 (MQ=2)
                                                                      ggATATTTGCGGCAAGATGAAGGGGCATG.ATTTCATCTTGCCCGAAATATCCC.AGGGGTGAATTGGCCAA.A..G.GCC..AAGAGGGGGCAGCG.C.CCCCTGCCCcg                                <  1:4015443/101‑1 (MQ=0)
                                                                       gATATTTGCGGCAAGATGAAGGGGCATG.ATTTCATCTTGCCCGAAATATCCC.AGGGGTGAATTGGCCAA.A..G.GCC..AAGAGGGGGCAGCG.C.CCCCTGCCCcgc                               >  1:4404457/1‑100 (MQ=2)
                                                                        atatTTGCGGCAAGATGAAGGGGCATG.ATTTCATCTTGCCCGAAATATCCC.AGGGGTGAATTGGCCAA.A..G.GCC..AAGAGGGGGCAGCG.C.CCCCTGCCCcgcc                              >  1:2418385/1‑99 (MQ=2)
                                                                         tatTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGC..GAA..G.CCA..GAGGGGGGC.AGCGCC.CCCCTGCCCcggcg                             >  1:1636751/1‑101 (MQ=18)
                                                                          atTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGC..GAA..G.CCA..GAGGGGGGC.AGCG.CACCCCTGCCCcggcgg                            >  2:2029738/1‑101 (MQ=18)
                                                                           tttGCGGCAAGATGAAGGGGCATG.ATTTCATCTTGCCCGAAATATCCC.AGGGGTGAATTGGCCAA.A..G.GCC..AAGAGGGGGCAGCG.C.CCCCTGCCCcgccgga                           >  2:2605858/1‑100 (MQ=2)
                                                                           tttGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGC..GAA..G.CCA..GAGGGGGGC.AGCGCC.CCCCTGCCCcggcggc                           >  2:977702/1‑100 (MQ=17)
                                                                             tGCAGCAAGATGAAACCGGATGCCTTTCATCTGGCCCGAAATATCCCGGGGGGTGAA.TGGC..GAA..G.CCA..GAGGGGGGC.AGCGCC.CCCCTGCCCCGGCGGTcc                         <  1:3190927/101‑1 (MQ=17)
                                                                              gcagcaAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGC..GAA..G.CCA..GAGGGGGGC.AGCGCC.CCCCTGCCCCGGCGGTCCg                        >  1:4463377/1‑101 (MQ=17)
                                                                                  caAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGC..GAA..G.CCA..GAGGGGGGC.AGCGCC.CCCCTGCCCCGGCGGTCCGTggc                    <  1:2524071/101‑1 (MQ=17)
                                                                                  caAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGC..GAA..G.CCA..GAGGGGGG.AAGCGCC.CCCCTGCCCCGGCGGTCCGTggc                    >  1:345412/1‑101 (MQ=18)
                                                                                   aaGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGC..GAA..G.CCA..GAGGGGGGC.AGCGCC.CCCCTGCCCCGGCGGCCCGTggcg                   >  1:34950/1‑101 (MQ=17)
                                                                                    aGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGC..GAA..G.CCA..GAGGGGGGC.AGCGCC.CCCCTGCCCCGGCGGTCCGTggcgg                  >  2:1161814/1‑101 (MQ=18)
                                                                                     gATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGC..GAA..G.CCA..GAGGGGGGC.AGCGCC.CCCCTGCCCCGGCGGTCCGTggcggc                 >  2:4232837/1‑101 (MQ=17)
                                                                                        gAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGC..GAA..G.CCA..GAGGGGGGC.AGCGCC.CCCCTGCCCCGGCGGTCCGTGGCGGCGcc              >  1:2626861/1‑101 (MQ=17)
                                                                                         aaaCCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGGGGA.TGGC..GAA..G.CCA..GAGGGGGGCAGGGC.C.CCCCTGCCCCCGCGGTCAGAGGCGGCGccg             >  1:4463056/1‑100 (MQ=12)
                                                                                         aaaCCGGATGCCTTTCATCTTGCCCGAAATATCCCGGAGGGTGAA.TGGC..GAA..G.CCA..GAGGGGGGC.AGCGCC.CCCCTGCCCCGGCGGTCCGTGGCGGCGCCt             >  1:2831091/1‑101 (MQ=18)
                                                                                           aCCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGC..GAA..G.CCA..GAGGGGGGC.AGCGCC.CCCCTGCCCCGGCGGTCCGTGGCGGCGCCTTg           >  1:812992/1‑101 (MQ=17)
                                                                                           aCCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGC..GAA..G.CCA..GAGGGGGGC.AGCGCC.CCCCTGCCCCGGCGGTCCGTGGCGGCGCCTTg           <  1:3335693/101‑1 (MQ=17)
                                                                                            ccGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGC..GAA..G.CCA..GAGGGGGGC.AGCGCC.CCCCTGCCCCGGCGGTCCGTGGCGGCGCCTTGa          <  1:2623656/101‑1 (MQ=17)
                                                                                             ccgATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGC..GAA..G.CCA..GAGGGGGGC.AGCGCC.CCCCTGCCCCGGCGGTCCGTGGCGGCGCCTTGAc         <  2:1774723/99‑1 (MQ=17)
                                                                                             cGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGC..GAA..G.CCA..GAGGGGGGC.AGCGCC.CCCCTGCCCCGGCGGTCCGTGGCGGCGCCTTGAc         >  1:769470/1‑101 (MQ=18)
                                                                                             cGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGC..GAA..G.CCA..GAGGGGGGC.AGCGCC.CCCCTGCCCCGGCGGTCCGTGGCGGCGCCTTGAc         <  2:4640867/101‑1 (MQ=17)
                                                                                              ggATGCCTTTCATCTTGCCCGAAATGTCCCGGGGGTTGAA.TGGC..GAA..G.CCA..GAGGGGGGC.AGCGCC.CCCCTGCCCCGGCGGTCCGTGGCGGCGCCTTGAcg        <  1:4552032/101‑1 (MQ=17)
                                                                                              ggATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGC..GAA..G.CCA..GAGGGGGGC.AGCGCC.CCCCTGCCCCGGCGGTCCGTGGCGGCGCCTTGCcg        >  1:1223598/1‑101 (MQ=17)
                                                                                               gATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGC..GAA..G.CCA..GAGGGGGGC.AGCGCC.CCCCTGCCCCGGCGGTCCGTGGCGGCGCCTTGAcgg       <  1:2311185/101‑1 (MQ=17)
                                                                                                aTGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGC..GAA..G.CCA..GAGGGGGGC.AGCGCC.CCCCTGCCCCGGCGGTCCGTGGCGGCGCCTTGAcggg      <  2:1430771/101‑1 (MQ=17)
                                                                                                aTGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGC..GAA..G.CCA..GAGGGGGGC.AGCGCC.CCCCTGCCCCGGCGGTCCGTGGCGGCGCCTTGAcggg      <  1:1172984/101‑1 (MQ=17)
                                                                                                 tGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGC..GAA..G.CCA..GAGGGGGGC.AGCGCC.CCCCTGCCCCGGCGGTCCGTGGCGGCGCCTTGAcgggc     <  1:4037518/101‑1 (MQ=17)
                                                                                                   ccTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGC..GAA..G.CCA..GAGGGGGGC.AGCGCC.CCCCTGCCCCGGCGGTCCGTGGCGGCGCCTTGAcgggcgg   <  1:1254947/101‑1 (MQ=17)
                                                                                                    cTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGC..GAA..G.CCA..GAGGGGGGC.AGCGCC.CCCCTGCCCCGGCGGTCCGTGGCGGCGCCTTGAcgggcggg  >  1:3228256/1‑101 (MQ=18)
                                                                                                    |                                                                                                              
CGGGCGCGCGGCACGACAGGCCGAGT.AA.GGGGGGGCGCTGCCCCCGTTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGC..G.A..G.CAG..AAGGGGGGCAAGCG.C.CCCCTGCCCCGGCGGTCCGTGGCGGCGCCTTGACGGGCGGG  >  NC_014034/2493925‑2494122

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