Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_014034 2,494,070 G→A 7.1% intergenic (+118/‑124) RCAP_rcc02319 → / → RCAP_rcc02320 hypothetical protein/hypothetical protein

Read alignment evidence...
  seq id position change freq score reads annotation genes product
*NC_0140342,494,0700G→A 7.1% 21.7 ‑0.3 ‑0.3155intergenic (+118/‑124)RCAP_rcc02319/RCAP_rcc02320hypothetical protein/hypothetical protein
Reads supporting (aligned to +/- strand):  new base (7/4):  ref base (78/66):  total (85/70)
Fisher's exact test strand distribution p-value = 7.55e-01
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 5.59e-01

Only 100 of 175 total aligned reads displayed.

GCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGC..G.A.G.CAGA.AGGGGGGCAAGCG.CCCCCTGCCCCGGCGGTCCGTGGCGGCGCCTTGACGGGCGGGGGCGCGGGCATTAACCGAAATTGCAACAGATGCTGCGATCCTCGCCCA  >  NC_014034/2493975‑2494170
                                                                                                      |                                                                                                     
  aaGCGGGCTCCCGCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGC..GAA.G.CCAG.AGGGGGGCa                                                                                                 <  1:2260694/101‑1 (MQ=18)
  aaGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGC..GAA.G.CCAG.AGGGGGGCa                                                                                                 <  1:1431276/101‑1 (MQ=18)
   aGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGC..GAA.G.CCAG.AGGGGGGCag                                                                                                <  1:4753551/101‑2 (MQ=18)
    gCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGC..GAA.G.CCAG.AGGGGGGCagc                                                                                               >  1:2176995/1‑99 (MQ=18)
      ggACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAC.TGGC..GAA.G.CCAG.AGGGGGGCagcgc                                                                                             >  1:4638215/1‑97 (MQ=17)
      ggACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGC..GAA.G.CCAG.AGGGGGGCagcgc                                                                                             >  2:172563/1‑97 (MQ=18)
      ggACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGC..GAA.G.CCAG.AGGGGGGCagcgc                                                                                             >  1:4676075/1‑97 (MQ=18)
        aCTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGC..GAA.G.CCAG.AGGGGGGCagcggcc                                                                                           >  2:2348947/1‑95 (MQ=17)
        aCTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGC..GAA.G.CCAG.AGGGGGGC.AGCG.ccc                                                                                         >  1:283443/1‑101 (MQ=18)
        aCTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGC..GAA.G.CCAG.AGGGGGGC.AGCG.ccc                                                                                         >  2:3025943/1‑101 (MQ=18)
          tCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGC..GAA.G.CCAG.AGGGGGGC.AGCG.ccccc                                                                                       >  1:1565968/1‑101 (MQ=18)
          tCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGC..GAA.G.CCAG.AGGGGGGC.AGCG.ccccc                                                                                       >  1:806307/1‑101 (MQ=17)
            cccccGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGC..GAA.G.CCAG.AGGGGGGC.AGCG.cccccct                                                                                     <  2:848858/101‑3 (MQ=17)
            cccccGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGC..GAA.G.CCAG.AGGGGGGC.AGCG.cccccct                                                                                     <  1:4471390/101‑3 (MQ=17)
              cccGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGC..GAA.G.CCAG.AGGGGGGC.AGCGACCCCCTGc                                                                                    >  2:1978451/1‑101 (MQ=17)
               ccGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGC..GAA.G.CCAG.AGGGGGGC.AGCG.cccccctgcc                                                                                  <  2:4018385/101‑6 (MQ=17)
                cGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGC..GAA.G.CCAT.AGGGGGGC.AGCGCCCCCCTGccc                                                                                  >  1:2925682/1‑101 (MQ=18)
                 gggATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGC..GAA.G.CCAG.AGGGGGGC.AGCGCCCCCCTGcccc                                                                                 >  1:3846964/1‑101 (MQ=17)
                 gggATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGC..GAA.G.CCAG.AGGGGGGC.AGCGCCCCCCTGcccc                                                                                 <  2:191057/101‑1 (MQ=17)
                  ggATATTTGCGGCAAGATGAAGGGGCATG.ATTTCATCTTGCCCGAAATATCCC.AGGGGTGAATTGGCCAA.A.G.GCCA.AGAGGGGGCAGCG.CCCCCTGCCCcg                                                                                <  1:4015443/101‑1 (MQ=0)
                  ggATATTTGCGGCAAGATGAAGGGGCATG.ATTTCATCTTGCCCGAAATATCCC.AGGGGTGAATTGGCCAA.A.G.GCCA.AGAGGGGGCAGCG.CCCCCTGCCCcg                                                                                >  1:659612/1‑101 (MQ=2)
                  ggATATTTGCGGCAAGATGAAGGGGCATG.ATTTCATCTTGCCCGAAATATCCC.AGGGGTGAATTGGCCAA.A.G.GCCA.AGAGGGGGCAGCG.CCCCCTGCCCcg                                                                                <  2:1362966/101‑1 (MQ=0)
                  ggATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGC..GAA.G.CCAG.AGGGGGGC.AGCGCCCCCCTGCCCcg                                                                                <  1:4417932/101‑1 (MQ=17)
                    atatTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGC..GAA.G.CCAG.AGGGGGGC.AGCGCCCCCCTGCCCcggg                                                                              >  1:3337824/1‑100 (MQ=17)
                     tatTTGCGGCAAGATGAAGGGGCATG.ATTTCATCTTGCCCGAAATATCCC.AGGGGTGAATTGGCCAA.A.G.GCCA.AGAGGGGGCAGCG.CCCCCTGCCCcgccg                                                                             >  2:1139753/1‑101 (MQ=2)
                     tatTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGC..GAA.G.CCAG.AGGGGGGC.AGCGCCCCCCTGCCCcggcg                                                                             >  1:1636751/1‑101 (MQ=18)
                       tttGCGGCAAGATGAAGGGGCATG.ATTTCATCTTGCCCGAAATATCCC.AGGGGTGAATTGGCCAA.A.G.GCCA.AGAGGGGGCAGCG.CCCCCTGCCCcgccgga                                                                           >  2:2878051/1‑100 (MQ=2)
                       tttGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGC..GAA.G.CCAG.AGGGGGGC.AGCGCCCCCCTGCCCcggcggc                                                                           >  2:977702/1‑100 (MQ=17)
                         tGCAGCAATATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGC..GAA.G.CCAG.AGGGGGGC.AGCGCCCCCCTGCCCCGGCGGTcc                                                                         >  1:3303403/1‑101 (MQ=17)
                         tGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGC..GAA.G.CCAG.AGGGGGGC.AGCGCCCCCCTGCCCCGGCGGTcc                                                                         >  2:263673/1‑101 (MQ=17)
                           cagcaAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGC..GAA.G.CCAG.AGGGGGGC.AGCGCCCCCCTGCCCCGGCGGTCCGt                                                                       <  1:4346247/101‑1 (MQ=17)
                            agcaAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGC..GAA.G.CCAG.AGGGGGGC.AGCGCCCCCCTGCCCCGGCGGTCCGTg                                                                      <  1:998655/101‑1 (MQ=17)
                            agcaAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGC..GAA.G.CCAG.AGGGGGGC.AGCGCCCCCCTGCCCCGGCGGTCCGTg                                                                      <  1:2414691/101‑1 (MQ=17)
                              caAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGC..GAA.G.CCAG.AGGGGGGC.AGCGCCCCCCTGCCCCGGCGGTCCGTggc                                                                    <  1:2524071/101‑1 (MQ=17)
                               aaGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGC..GAA.G.CCAG.AGGGGGGC.AGCGCCCCCCTGCCCCGGCGGCCCGTggcg                                                                   >  1:34950/1‑101 (MQ=17)
                               aaGATCAAACCGGATGCCTTTCATCTTGCCAGAAATCGCCAGGGGGGTGAA.TGGC..GAA.G.CCAG.AGGGGGGC.AGCGCCCCCCTGCCCCGGCGGTCCGTggcg                                                                   <  1:836859/101‑1 (MQ=9)
                                    gAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGC..GAA.G.CCAG.AGGGGGGC.AGCGCCCCCCTGCCCCGGCGGTCCGTGGCGGCGcc                                                              >  1:2626861/1‑101 (MQ=17)
                                    gAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGC..GAA.G.CCAG.AGGGGGGC.AGCGCCCCCCTGCCCCGGCGGTCCGTGGCGGCGcc                                                              <  2:4421536/101‑1 (MQ=17)
                                    gAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGC..GAA.G.CCAG.AGGGGGG.AAGCGCCCCCCTGCCCCGGCGGTCCGGGGCGGGGcc                                                              >  2:3683987/1‑101 (MQ=17)
                                     aaaCCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGGGGA.TGGC..GAA.G.CCAG.AGGGGGGCAGGGC.CCCCCTGCCCCCGCGGTCAGAGGCGGCGccg                                                             >  1:4463056/1‑100 (MQ=12)
                                     aaaCCGGATGCCTTTCATCTTGCCCGAAATATCCCGGAGGGTGAA.TGGC..GAA.G.CCAG.AGGGGGGC.AGCGCCCCCCTGCCCCGGCGGTCCGTGGCGGCGCCt                                                             >  1:2831091/1‑101 (MQ=18)
                                         ccgATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGC..GAA.G.CCAG.AGGGGGGC.AGCGCCCCCCTGCCCCGGCGGTCCGTGGCGGCGCCTTGAc                                                         <  2:1774723/99‑1 (MQ=17)
                                         cGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGC..GAA.G.CCAG.AGGGGGGC.AGCGCCCCCCTGCCCCGGCGGTCCGTGGCGGCGCCTTGAc                                                         >  1:769470/1‑101 (MQ=18)
                                         cGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGC..GAA.G.CCAG.AGGGGGGC.AGCGCCCCCCTGCCCCGGCGGTCCGTGGCGGCGCCTTGAc                                                         <  2:4640867/101‑1 (MQ=17)
                                          ggATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGC..GAA.G.CCAG.AGGGGGGC.AGCGCCCCCCTGCCCCGGCGGTCCGTGGCGGCGCCTTGAcg                                                        >  1:3101667/1‑101 (MQ=17)
                                           gATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGC..GAA.G.CCAG.AGGGGGGC.AGCGCCCCCCTGCCCCGGCGGTCCGTGGCGGCGCCTTGAcgg                                                       <  1:2311185/101‑1 (MQ=17)
                                            aTGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGC..GAA.G.CCAG.AGGGGGGC.AGCGCCCCCCTGCCCCGGCGGTCCGTGGCGGCGCCTTGAcggg                                                      <  2:1430771/101‑1 (MQ=17)
                                             tGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGC..GAA.G.CCAG.AGGGGGGC.AGCGCCCCCCTGCCCCGGCGGTCCGTGGCGGCGCCTTGAcgggc                                                     <  1:4037518/101‑1 (MQ=17)
                                               ccTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGC..GAA.G.CCAG.AGGGGGGC.AGCGCCCCCCTGCCCCGGCGGTCCGTGGCGGCGCCTTGAcgggcgg                                                   <  1:1254947/101‑1 (MQ=17)
                                                cTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGC..GAA.G.CCAG.AGGGGGGC.AGCGCCCCCCTGCCCCGGCGGTCCGTGGCGGCGCCTTGAcgggcggg                                                  >  2:84451/1‑101 (MQ=17)
                                                cTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGC..GAA.G.CCAG.AGGGGGGC.AGCGCCCCCCTGCCCCGGCGGTCCGTGGCGGCGCCTTGAcgggcggg                                                  >  1:3228256/1‑101 (MQ=18)
                                                 tttCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGC..GAA.G.CCAG.AGGGGGGC.AGCGCCCCCCTGCCCCGGCGGTCCGTGGCGGCGCCTTGACGGGCgggg                                                 <  1:4605986/101‑1 (MQ=17)
                                                 tttCATCTTGCCCGAAATATCCCGGGGGGGGAA.TGCG..A.A.A.CCAA.AGGGGGGCCAGGGCCCCCCGGCCCCGGCGG.CCGGGGCCGGGCCCTTGCCGGCggggg                                                >  1:1221501/1‑101 (MQ=9)
                                                   tCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGC..GAA.G.CCAG.AGGGGGGC.AGCGCCCCCCTGCCCCGGCGGTCCGTGGCGGCGCCTTGACGGGCGGGGgc                                               <  2:2564012/101‑1 (MQ=18)
                                                   tCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGC..GAA.G.CCAG.AGGGGGGC.AGCGCCCCCCTGCCCCGGCGGTCCGTGGCGGCGCCTTGACGGGCGGGGgc                                               <  2:4272800/101‑1 (MQ=18)
                                                   tCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGC..GAA.G.CCAG.AGGGGGGC.AGCGCCCCCCTGCCCCGGCGGTCCGAGGCGGCGCCTTGACGGGCGGGGgc                                               >  2:2012970/1‑101 (MQ=17)
                                                    cATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGC..GAA.G.CCAG.AGGGGGGC.AGCGCCCCCCTGCCCCGGCGGTCCGTGGCGGCGCCTTGACGGGCGGGGgcg                                              >  1:3806231/1‑101 (MQ=18)
                                                      tCTTTCCCGAAATATCCCGGGGGGTGAA.TGGC..GAA.G.CCAG.AGGGGGGC.AGCGCCCCCCTGCCCCGGCGGTCCGTGGCGGCGCCTTGACGGGCGGGGGcgcg                                            <  2:2176871/101‑1 (MQ=18)
                                                      tCTTGCCCGAAATATCCCGGGGGGTGAA.TGGC..GAA.G.CCAG.AGGGGGGC.AGCGCCCCCCTGCCCCGGCGGTCCGTGGCGGCGCCTTGACGGGCGGGGGcgcg                                            >  2:4376916/1‑101 (MQ=18)
                                                       cTTGCCCGAAATATCCCGGGGGGTGAA.TGGC..GAA.G.CCAG.AGGGGGGC.AGCGCCCCCCTGCCCCGGCGGTCCGTGGCGGCGCCTTGACGGGCGGGGGCGCgg                                           >  2:1633663/1‑101 (MQ=18)
                                                         tGCCCGAAATATCCCGGGGGGTGAA.TGGC..GAA.G.CCAG.AGGGGGGC.AGCGCCCCCCTGCCCCGGCGGTCCGTGGCGGCGCCTTGACGGGCGGGGGCGCGGGc                                         <  1:617765/101‑1 (MQ=18)
                                                          gCCCGAAATATCCCGGGGGGTGAA.TGGC..GAA.G.CCAG.AGGGGGGC.AGCGCCCCCCTGCCCCGGCGGTCCGTGGCGGCGCCTTGACGGGCGGGGGCGCGGGCa                                        >  2:4636090/1‑101 (MQ=18)
                                                          gCCCGAAATATCCCGGGGGGTGAA.TGGC..GAA.G.CCAG.AGGGGGGC.AGCGCCCCCCTGCCCCGGCGGTCCGTGGCGGCGCCTTGACGGGCGGGGGCGCGGGCa                                        <  2:2094526/101‑1 (MQ=18)
                                                           cccGAAATATCCCGGGGGGTGAA.TGGC..GAA.G.CCAG.AGGGGGGC.AGCGCCCCCCTGCCCCGGCGGTCCGTGGCGGCGCCTTGACGGGCGGGGGCGCGGGCAt                                       <  1:172579/101‑1 (MQ=18)
                                                           cccGAAATATCCCGGGGGGTGAA.TGGC..GAA.G.CCAG.AGGGGGGC.AGCGCCCCCCTGCCCCGGCGGTCCGTGGCGGCGCCTTGACGGGCGGGGGCGCGGGCAt                                       <  2:3538976/101‑1 (MQ=18)
                                                            ccGAAATATCCCGGGGGGTGAA.TGGC..GAA.G.CCAG.AGGGGGGC.AGCGCCCCCCTGCCCCGGCGGTCCGTGGCGGGGCCTTGACGGGCGGGGGCGCGGGCAtt                                      >  1:2628535/1‑101 (MQ=18)
                                                            ccGAAATATCCCGGGGGGTGAA.TGGC..GAA.G.CCAG.AGGGGGGC.AGCGCCCCCCTGCCCCGGCGGTCCGTGGCGGCGCCTTGACGGGCGGGGGCGCGGGCAtt                                      <  2:2304365/101‑1 (MQ=18)
                                                            ccGAAATATCCCGGGGGGTGAA.TGGC..GAA.G.CCAG.AGGGGGGC.AGCGCCCCCCTGCCCCGGCGGTCCGTGGCGGCGCCTTGACGGGCGGGGGCGCGGGCAtt                                      <  2:3593736/101‑1 (MQ=18)
                                                            ccGAAATACCCCGGGGGGTGAA.TGGC..GAA.G.CCAG.AGGGGGGC.AGCGCCCCCCTGCCCCGGCGGTCCGTGGCGGCGCCTTGACGGGCGGGGGCGCGGGCAtt                                      <  1:2122/101‑1 (MQ=18)
                                                              gAAATATCCCGGGGGGTGAA.TGGC..GAA.G.CCAG.AGGGGGGC.AGCGCCCCCCTGCCCCGGCGGTCCGTGGCGGCGCCTTGACGGGCGGGGGCGCGGGCATTaa                                    <  1:98413/101‑1 (MQ=18)
                                                               aaaTATCCCGGGGGGTGAA.TGGC..GAA.G.CCAG.AGGGGGGC.AGCGCCCCCCTGCCCCGGCGGTCCGTGGCGGCGCCTTGACGGGCGGGGGCGCGGGCATTAAc                                   <  2:4734183/101‑1 (MQ=18)
                                                               aaaTATCCCGGGGGGTGAA.TGGC..GAA.G.CCAG.AGGGGGGC.AGCGCCCCCCTGCCCCGGCGGTCCGTGGCGGCGCCTTGACGGGCGGGGGCGCGGGCATTAAc                                   <  1:459589/101‑1 (MQ=18)
                                                                aaTATCCCGGGGGGTGAA.TGGC..GAA.G.CCAG.AGGGGGGC.AGCGCCCCCCTGCCCCGGCGGTCCGTGGCGGCGCCTTGACGGGCGGGGGCGCGGGCATTAAcc                                  >  2:90304/1‑101 (MQ=21)
                                                                aaTATCCCGGGGGGTGAA.TGGC..GAA.G.CCAG.AGGGGGGC.AGCGCCCCCCTGCCCCGGCGGTCCGTGGCGGCGCCTTGACGGGCGGGGGCGCGGGCATTAAcc                                  >  1:1127494/1‑101 (MQ=21)
                                                                aaTATCCCGGGGGGTGAA.TGGC..GAA.G.CCAG.AGGGGGGC.AGCGCCCCCCTGCCCCGGCGGTCCGTGGCGGCGCCTTGACGGGCGGGGGCGCGGGCATTAAcc                                  <  1:2740830/101‑1 (MQ=21)
                                                                 atatCCCGGGGGGTGAA.TGGC..GAA.G.CCAG.AGGGGGGC.AGCGCCCCCCTGCCCCGGCGGTCCGTGGCGGCGCCTTGACGGGCGGGGGCGCGGGCATTAACCg                                 <  2:4398025/101‑1 (MQ=21)
                                                                  taccccGGGGGGTGAA.TGGC..G.A.A.GCCA.GAGGGGGCCAGCGCCCCCCTGCCCCGGCGGTCCGTGGCGGCGCCTTGACGGGCGGGGGCGCGGGCATTAACCGa                                <  2:1329142/98‑1 (MQ=21)
                                                                   atCCCGGGGGGTGAA.TGGC..GAA.G.CCAG.AGGGGGGC.AGCGCCCCCCTGCCCCGGCGGTCCGTGGCGGCGCCTTGACGGGCGGGGGCGCGGGCATTAACCGaa                               >  2:1228239/1‑101 (MQ=21)
                                                                   atCCCGGGGGGTGAA.TGGC..GAA.G.CCAG.AGGGGGG.GAGCGCCCCCCTGCCCCGGCGGTCCGTGGCGGGCGCTTGACGGGCGGGGGCGCGGGCATTAACCGaa                               >  1:538066/1‑101 (MQ=16)
                                                                    tCCCGGGGGGTGAA.TGGC..GAA.G.CCAG.AGGGGGGC.AGCGCCCCCCTGCCCCGGCGGTCCGTGGCGGCGCCTTGACGGGCGGGGGCGCGGGCATTAACCGaaa                              <  1:4150127/101‑1 (MQ=21)
                                                                        ggggggTGAA.TGGC..GAA.G.CCAG.AGGGCGGC.AGCGCCCCCCTGCCCCGGCGGTCCGTGGCGGCGCCTTGACGGGCGGGGGCGCGGGCATTAACCGAAATTGc                          >  1:4254827/1‑101 (MQ=25)
                                                                         gggggTGAA.TGGC..G.AGG.CCAG.AGGGGGGC.AGCGCCCCCCTGCCCCGGCGGTCCGTGGCGGCGCCTTGACGGGCGGGGGCGCGGGCATTAACCGAAATTGCa                         <  2:417052/101‑1 (MQ=25)
                                                                           gggTGAA.TGGC..GAA.GCCAG..AGGGGGGC.AGCGCCCCCCTGCCCCGGCGGTCCGTGGCGGCGCCTTGACGGGCGGGGGCGCGGGCATTAACCGAAATTGCAAc                       >  1:3130257/1‑101 (MQ=25)
                                                                           gggTGAA.TGGC..GAA.G.CCAG.AGGGGGGC.AGCGCCCCCCTGCCCCGGCGGTCCGTGGCGGCGCCTTGACGGGCGGGGGCGCGGGCATTAACCGAAATTGCAAc                       <  1:2590566/101‑1 (MQ=25)
                                                                           gggTGAA.TGGC..G.A.A.GAGAGAGGGGGGC.AGCGCCCCCCTGCCCCGGCGGTCCGTGGCGGCGCCTTGACGGGCGGGGGCGCGGGCATTAACCGAAATTGCAAc                       <  2:71212/101‑1 (MQ=37)
                                                                            ggTGAA.TGGC..GAA.GCCAG..AGGGGGGC.AGCGCCCCCCTGCCCCGGCGGTCCGTGGCGGCGCCTTGACGGGCGGGGGCGCGGGCATTAACCGAAATTGCAACa                      >  2:377256/1‑101 (MQ=25)
                                                                             gTGAA.TGGC..GAA.GCCAG..AGGGGGGC.AGCGCCCCCCTGCCCCGGCGGTCCGTGGCGGCGCCTTGACGGGCGGGGGCGCGGGCATTAACCGAAATTGCAACAg                     >  1:3642685/1‑101 (MQ=25)
                                                                             gTGAA.TGGC..GAA.G.CCAG.AGGGGGGC.AGCGCCCCCCTGCCCCGGCGGTCCGTGGCGGCGCCTTGACGGGCGGGGGCGCGGGCATTAACCGAAATTGCAACAg                     <  1:3047534/101‑1 (MQ=25)
                                                                                          gcgaagccagaGGGGGGC.AGCGCCCCCCTGCCCCGGCGGTCCGTGGCGGCGCCTTGACGGGCGGGGGCGCGGGCATTAAACGAAATTGCAACAGAtgctgc              >  2:547187/11‑101 (MQ=37)
                                                                                           cgaagccagaGGGGGGC.AGCGCCCCCCTGCCCCGGCGGTCCGTGGCGGCGCCTTGACGGGCGGGGGCGCGGGCATTAACCGAAATTGCAACAGATGCTGCg             >  2:3657279/10‑101 (MQ=38)
                                                                                            gaagccagaGGGGGGC.AGCGCCCCCCTGCCCCGGCGGTCCGTGGCGGCGCCTTGACGGGCGGGGGCGCGGGCATTAACCGAAATTGCAACAGATGCTGCGa            >  1:1992962/9‑101 (MQ=38)
                                                                                             aagccagaGGGGGGC.AGCGCCCCCCTGCCCCGGCGGTCCGTGGCGGCGCCTTGACGGGCGGGGGCGCGGGCATTAACCGAAATTGCAACAGATGCTGCGAt           >  1:4579513/8‑101 (MQ=38)
                                                                                                ccagaGGGGGGC.AGCGCCCCCCTGCCCCGGCGGTCCGTGGCGGCGCCTTGACGGGCGGGGGCGCGGGGATTAACCGGAATTGCAACAGCTGCTGCGATCCt        >  2:2927109/5‑101 (MQ=37)
                                                                                                 cagaGGGGGGC.AGCGCCCCCCTGCCCCGGCGGTCCGTGGCGGCGCCTTGACGGGCGGGGGCGCGGGCATTAACCGAAATTGCAACAGATGCTGCGATCCTc       >  2:1961035/4‑101 (MQ=38)
                                                                                                 cagaGGGGGGC.AGCGCCCCCCTGCCCCGGCGGTCCGTGGCGGCGCCTTGACGGGCGGGGGCGCGGGCATTAACCGAAATTGCAACAGATGCTGCGATCCTc       >  1:81612/4‑101 (MQ=38)
                                                                                                   agggggggGCAGCGCCCCCCTGCCCCGGCGGTCCGTGGCGGCGCCTTGACGGGCGGGGGCGCGGGCATTAACCGAAATTGCAACAGATGCTGCGATCCTCg      >  1:2711651/3‑101 (MQ=255)
                                                                                                     ggggggC.AGCGCCCCCCTGCCCCGGCGGTCCGTGGCGGCGCCTTGACGGGCGGGGGCGCGGGCATTAACCGAAATTGCAACAGATGCTGCGATCCTCGccc   <  1:3169857/101‑1 (MQ=255)
                                                                                                     ggggggC.AGCGCCCCCCTGCCCCGGCGGTCCGTGGCGGCGCCTTGACGGGCGGGGGCGCGGGCATTAACCGAAATTGCAACAGATGCTGCGATCCTCGccc   <  2:806234/101‑1 (MQ=255)
                                                                                                      gggggC.AGCGCCCCCCTGCCCCGGCGGTCCGTGGCGGCGCCTTGACGGGCGGGGGCGCGGGCATTAACCGAAATTGCAACAGATGCTGCGATCCTCGCCCa  >  1:3266503/1‑101 (MQ=255)
                                                                                                      gggggC.AGCGCCCCCCTGCCCCGGCGGTCCGTGGCGGCGCCTTGACGGGCGGGGGCGCGGGCATTAACCGAAATTGCAACAGATGCTGCGAACCTCGCCCa  <  1:1079081/101‑1 (MQ=255)
                                                                                                      |                                                                                                     
GCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGC..G.A.G.CAGA.AGGGGGGCAAGCG.CCCCCTGCCCCGGCGGTCCGTGGCGGCGCCTTGACGGGCGGGGGCGCGGGCATTAACCGAAATTGCAACAGATGCTGCGATCCTCGCCCA  >  NC_014034/2493975‑2494170

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