Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_014034 586,916 +G 2.6% intergenic (‑144/+38) RCAP_rcc00548 ← / ← purL BolA family protein/phosphoribosylformylglycinamidine synthase II

Read alignment evidence...
  seq id position change freq score reads annotation genes product
*NC_014034586,9161.→G 2.6% 6.7 ‑4.2 0.0191intergenic (‑144/+38)RCAP_rcc00548/purLBolA family protein/phosphoribosylformylglycinamidine synthase II
Reads supporting (aligned to +/- strand):  new base (0/5):  ref base (108/78):  total (108/83)
Fisher's exact test strand distribution p-value = 1.44e-02
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 1.00e+00

Only 100 of 243 total aligned reads displayed.

AGGCGACAAGAGCCAAGCACCGGTTGCGTGTGCCAAGTGACGCAGCGGGCAGGGGCCGCGGCGAGGACGGGGGTCAGGACGGGGGCAAGTGTCGCCGAGC.T.C..GGGGGGGCGCTTCTGCCCGCCCCGGAGCTTGCGTTTTCAGCCTTCGACGGCGGCGGCGAAGGCGGTGCGGTAAAGCGTGCTGAGCTCTGCCAGCG  >  NC_014034/586816‑587012
                                                                                                      |                                                                                                  
agGCGACAAGAGCCAAGCACCGGTTGCGTGTGCCAAGTGACGCAGCGGGCAGGGGCCGCGGCGAGGACGGGGGTCAGGACGGGGGCAAGTGTCGCCGAGC.t                                                                                                     >  1:2815386/1‑101 (MQ=255)
agGCGACAAGAGCCAAGCACCGGTTGCGTGTGCCAAGTGACGCAGCGGGCAGGGGCCGCGGCGAGGACGGGGGTCAGGACGGGGGCAAGTGTCGCCGAGC.t                                                                                                     <  2:4986266/101‑1 (MQ=255)
    gACAAGAGCCAAGCACCGGTTGCGTGTGCCAAGTGACGCAGCGGGCAGGGGCCGCGGCGAGGACGGGGGTCAGGACGGGGGCAAGTGTCGCCGAGC.T.C..ggg                                                                                              >  1:4179436/1‑101 (MQ=255)
    gACAAGAGCCAAGCACCGGTTGCGTGTGCCAAGTGAAGCAGCGGGCAGGGGCCGCGGCGAGGACGGGGGTCAGGACGGGGGCAAGTGTCGCCGAGC.T.C..ggg                                                                                              >  2:5709917/1‑101 (MQ=255)
       aaGAGCTAAGCACCGGTTGCGTGTGCCAAGTGACGCAGCGGGCAGGGGCCGCGGCGAGGACGGGGGTCAGGACGGGGGCAAGGGTCGCCGAGC.T.C..gggggg                                                                                           >  2:1374048/1‑101 (MQ=255)
       aaGAGCCAAGCACCGGTTGCGTGTGCCAAGTGACGCAGCGGGCAGGGGCCGCGGCGAGGACGGGGGTCAGGACGGGGGCAAGTGTCGCCGAGC.T.C..gggggg                                                                                           >  2:3483185/1‑101 (MQ=255)
        agagCCAAGCACCGGTTGCGTGTGCCAAGTGACGCAGCGGGCAGGGGCCGCGGCGAGGACGGGGGTCAGGACGGGGGCAAGTGTCGCCGAGC.T.C..ggggggg                                                                                          <  1:3256159/101‑1 (MQ=255)
         gagCCAAGCACCGGTTGCGTGTGCCAAGTGACGCAGCGGGCAGGGGCCGCGGCGAGGACGGGGGTCAGGACGGGGGCAAGTGTCGCCGAGC.T.C..gggggggg                                                                                         >  1:1658573/1‑100 (MQ=255)
          atccAAGCACCGGTTGCGTGTGCCAAGTGACGCAGCGGGCAGGGGCCGCGGCGAGGACGGGGGTCAGGACGGGGGCAAGTGTCGCCGAGC.T.C..ggggggggc                                                                                        >  1:2454013/3‑99 (MQ=255)
          agCCAAGCGCGGGTTGCGTGGGCCAGGTGAGGCGGGGGGCAGGGGCCGCGGCGAGGGGGGGGGTCGGGAGGGGGGCAAGTGTCGCCGAGC.T.C..ggggggggc                                                                                        <  2:4785796/101‑3 (MQ=255)
          agCCAAGCACCGGTTGCGTGTGCCAAGTGACGCAGCGGGCAGGGGCCGCGGCGAGGACGGGGGTCAGGACGGGGGCAAGTGTCGCCGAGC.T.C..ggggggggg                                                                                        >  2:1728015/1‑99 (MQ=255)
          agCCAAGCACCGGTTGCGTGTGCCAAGTGACGCAGCGGGCAGGGGCCGCGGCGAGGACGGGGGTCAGGACGGGGGCAAGTGTCGCCGAGC.T.C..ggggggggc                                                                                        >  1:787942/1‑99 (MQ=255)
           gCCAAGCACCGGTTGCGTGTGCCAAGTGACGCAGCGGGCAGGGGCCGCGGCGAGGACGGGGGTCAGGACGGGGGCAAGTGTCGCCGAGC.T.C..ggggggggcg                                                                                       <  1:2797360/101‑4 (MQ=255)
             cAAGCACCGGTTGCGTGTGCCAAGTGACGCAGCGGGCAGGGGCCGCGGCGAGGACGGGGGTCAGGACGGGGGCAAGTGTCGCCGAGC.T.C..ggggggggcgct                                                                                     >  2:5384640/1‑96 (MQ=255)
             cAAGCACCGGTTGCGTGGGCCAAGTGACGCAGCGGGCAGGGGCCGCGGCGAGGACGGGGGTCAGGACGGGGGCAAGTGTCGCCGAGC.T.C..ggggggggcgct                                                                                     <  1:4816156/101‑6 (MQ=255)
               aGCACGTGTTGCGTGTGCCGAGTGACGCGGGGGGGGGGGGCCGCGGCGGGGAGGGGGGGGAGGAGGGGGGCAAGTGTCGCCGAGC.TGG..GGGGGGGCGCTTc                                                                                    <  2:4245189/101‑1 (MQ=255)
               aGCACCGGTTGCGTGTGCCAAGTGACGCAGCGGGCAGGGGCCGCGGCGAGGACGGGGGTCAGGACGGGGGCAAGTGTCGCCGAGC.T.CG.GGGGGGGCGCTTc                                                                                    >  1:648154/1‑101 (MQ=255)
               aGCACAGGTTGCGTGTGCCAAGTGACGCAGCGGGCAGGGGCAGCGGCGAGGATGGGGGTCAGGCCGGGGGCAAGTGTCCCCGAGCTT.G..GGGGGGGCGCTTc                                                                                    >  2:3876978/1‑101 (MQ=255)
                    cGGTTGCGTGTGCCAAGGGACGCAGCGGGCAGGGGCCGCGGCGAGGACGGGGGTCAGGACGGGGGCAAGTGTCGCCGAGC.T.CG.GGGGGGGCGCTTCTgccc                                                                               <  1:3649756/101‑1 (MQ=255)
                     ggTTGCGTGTGCCAAGTGACGCAGCGGGCAGGGGCCGCGGCGAGGACGGGGGTCAGGACGGGGGCAAGTGTCGCCGAGC.T.CG.GGGGGGGCGCTTCTgcccg                                                                              >  2:2860130/1‑101 (MQ=255)
                        tGCGTGTGCCAAGTGACGCAGCGGGCAGGGGCCGCGGCGAGGACGGGGGTCAGGACGGGGGCAAGTGTCGCCGAGC.T.CG.GGGGGGGCGCTTCTgcccgccc                                                                           >  1:1152283/1‑101 (MQ=255)
                        tGCGTGTGCCAAGTGACGCAGCGGGCAGGGGCCGCGGCGAGGACGGGGGGCAGGACGGGGGCAAGTGTCGCCGCGC.T.CG.GGGGGGGAGCTTCTgccagccc                                                                           >  1:3755552/1‑101 (MQ=255)
                         gCGTGTGCCAGGGGACGCAGCGGGGAGGGGCCGCGGGGAGGACGGGGGTCAGGACGGGGGCAAGTGTCGCCGAGC.T.CG.GGGGGGGCGCTTCTGCCCGcccc                                                                          <  2:2920811/101‑1 (MQ=255)
                         gCGTGTGCCAAGTGACGCAGCGGGCAGGGGCCGCGGCGAGGACGGGGGTCAGGACGGGGGCAAGTGTCGCCGAGC.T.CG.GGGGGGGCGCTTCTGCCCGcccc                                                                          <  1:1898606/101‑1 (MQ=255)
                         gCGTGGGCCAAGTGACGCAGCGGGGGGGGGCCGCGGCGAGGACGGGGGTCAGGACGGGGGCAAGTGTCGCCGAGC.T.CG.GGGGGGGCGCTTCTGCCCGcccc                                                                          <  2:4009693/101‑1 (MQ=255)
                              tgCCAAGTGACGCAGCGGGCAGGGGCCGCGGCGAGGACGGGGGTCAGGACGGGGGCAAGTGTCGCCGAGC.T.CG.GGGGGGGCGCTTCTGCCCGCCCCGGAGc                                                                     <  2:3359885/101‑1 (MQ=255)
                               gCCAAGGGACGCGGGGGGGGGGGGCGGGGGGGAGGGGGGGGGTCAGGACGGGGGCAAGTGTCGCCGGGC.T.CG.GGGGGGGCGCTTCTGCCCGCCCCGGAGCt                                                                    <  1:4437571/101‑1 (MQ=255)
                                cgaaGGGACGCAGCGGGCGGGGGGCGCGGCGGGGCCGGGGGTCAGGCCGGGGGCAAGTGTCGGAGAGC.TGG..GGGGGGGCGCTTCTGCCCGCCCCGGAGCtt                                                                   <  1:1789442/99‑1 (MQ=255)
                                   aGTGACGCAGCGGGCAGGGGCCGCGGCGAGGACGGGGGTCAGGACGGGGGCAAGTGTCGCCGAGC.T.CG.GGGGGGGAGCTTCTGCCCGCCCCGGAGCTTGCg                                                                >  2:2601319/1‑101 (MQ=255)
                                    gTGACGCAGCGGGCAGGGGCCGCGGCGAGGACGGGGGTCAGGACGGGGGCAAGTGTCGCCGAGC.T.CG.GGGGGGGCGCTTCTGCCCGCCCCGGAGCTTGCGt                                                               >  2:4823881/1‑101 (MQ=255)
                                      ggcGCAGGGGGGAGGGGACGCGGCGAGGAGGGGGGCCAGGGCGGGGAGAAATGTCGGGGAACGG.G..GGGGGGGCGCTTCTGCCCGCCCCGGAGCTTGCGttt                                                             <  1:297045/99‑1 (MQ=255)
                                      gACGCAGCGGGCAGGGGCCGCGGCGAGGACGGGGGTCAGGACGGGGGCAAGTGTCGCCGAGC.T.CG.GGGGGGGCGCTTCTGCCCGCCCCGGAGCTTGCGttt                                                             <  2:4262333/101‑1 (MQ=255)
                                      gACGCAGCGGGCAGGGGCCGCGGCGAGGACGGGGGTCAGGACGGGGGCAAGTGTCGCCGAGC.T.CG.GGGGGGGCGCTTCTGCCCGCCCCGGAGCTTGCGttt                                                             >  2:1542520/1‑101 (MQ=255)
                                      gACGCAGCGGGCAGGGGCCGCGGCGAGGACGGGGGCCAGGACGGGGGCAAGTGTCGCCGAGC.G.CG.GGGGGGGCGCTTCTGCCCGCCCCGGAGCTTGCGttt                                                             <  2:2529088/101‑1 (MQ=255)
                                       gcGCAGCGGGCAGGGGCCGCGGCGAGGACGGGGGTCAGGACGGGGGCAAGTGTCGCCGAGC.T.CG.GGGGGGGCGCTTCTGCCCGCCCCGGAGCTTGCGtttt                                                            <  2:1627560/100‑1 (MQ=255)
                                       aCGCAGCGGGCAGGGGCCGCGGCGAGGACGGGGGTCAGGACGGGGGCAAGTGTCGCCGAGC.T.CG.GGGGGGGCGCTTCTGCCCGCCCCGGAGCTTGCGtttt                                                            <  2:3736341/101‑1 (MQ=255)
                                        ctcAGCGGGCAGGGGCCGCGGCGAGGACGGGGGTCAGGACGGGGGCAAGTGTCGCCGAGC.T.CG.GGGGGGGCGCTTCTGCCCGCCCCGGAGCTTGCGTTTTc                                                           <  2:1424492/99‑1 (MQ=255)
                                        cGCAGCGGGCAGGGGCCGCGGCGAGGACGGGGGTCAGGACGGGGGCAAGTGTCGCCGAGC.T.CG.GGGGGGGCGCTTCTGCCCGCCCCGGAGCTTGCGTTTTc                                                           <  1:1379850/101‑1 (MQ=255)
                                         gCAGCGGGCAGGGGCCGCGGCGAGGACGGGGGTCAGGACGGGGGCAAGTGTCGCCGAGC.T.CG.GGGGGGGCGCTTCTGCCCGCCCCGGAGCTTGCGTTTTCa                                                          <  1:5704236/101‑1 (MQ=255)
                                          cAGCGGGCAGGGGCCGCGGCGAGGACGGGGGTCAGGACGGGGGCAAGTGTCGCCGAGC.T.CG.GGGGGGGCGCTTCTGCCCGCCCCGGAGCTTGCGTTTTCAg                                                         <  2:4491751/101‑1 (MQ=255)
                                          cAGCGGGCAGGGGCCGCGGCGAGGACGGGGGTCAGGACGGGGGCAAGTGTCGCCGAGC.T.CG.GGGGGGGCGCTTCTGCCCGCCCCGGAGCTTGCGTTTTCAg                                                         >  1:147775/1‑101 (MQ=255)
                                           aGCGGGCGGGGGCGGCGGCGAGGGCGGGGGACAGGACGGGGGCAAGTGTCGCCGGGC.T.CG.GGGGGGGCGCTTCTGCCCGCCCCGGAGCTTGCGTTTTCAGc                                                        <  1:5741019/101‑1 (MQ=255)
                                           aGCGGGCAGGGGCCGCGGCGAGGACGGGGGTCAGGACGGGGGCAAGTGTGGCCGAGC.T.CG.GGGGGGGCGCTTCTGCCCGCCCCGGAGCTCGCGTCTTCAGc                                                        >  1:952186/1‑101 (MQ=255)
                                             ggggCAGGGGCCGCGGCGAGGACGGGGGTCAGGACGGGGGCAAGTGTCGCCGAGC.TGG..GGGGGGGCGCTTCTGCCCGCCCCGGAGCTTGCGTTTTCAGCCt                                                      <  2:2489778/100‑1 (MQ=255)
                                             cGGGCAGGGGCCGCGGCGAGGACGGGGGTCAGGACGGGGGCAAGTGTCGCCGAGC.T.CG.GGGGGGGCGCTTCTGCCCGCCCCGGAGCTTGCGTTTTCAGCCt                                                      >  1:3872469/1‑101 (MQ=255)
                                              gggCAGGGGCGGCGGCGAGGACGGGGGTCAGGACGGGGGCAAGTGTCGCCGAGC.T.CG.GGGGGGGCGCTTCTGCCCGCCCCGGAGCTTGCGTTTTCAGCCtt                                                     <  2:237502/101‑1 (MQ=255)
                                              gggCAGGGGCCGCGGCGAGGACGGGGGTCAGGACGGGGGCAAGTGTCGCCGAGC.T.CG.GGGGGGGCGCTTCTGCCCGCCCCGGAGCTTGCGTTTTCAGCCtt                                                     <  1:517749/101‑1 (MQ=255)
                                              gggCAGGGGCCGCGGCGAGGACGGGGGTCAGGACGGGGGCAAGTGTCGCCGAGC.T.CG.GGGGGGGCGCTTCTGCCCGCCCCGGAGCTTGCGTTTTCAGCCtt                                                     >  2:712231/1‑101 (MQ=255)
                                               ggCGGGGGCCGGGGCGAGGACGGGGGGCAGGGGGGGGGCAAGTGTCGCCGAGCGG.G..GGGGGGGCGCTTCTGCCCGCCCCGGAGCTTGCGTTTTCAGCCTTc                                                    <  1:682627/101‑1 (MQ=255)
                                                gggggggCGGCGGCGAGGAGGGGGGTCGGGACGGGGGCAAGGGCCGCCGGGC.T.CG.GGGGGGGCGCTTCTGCCCGCCCCGGAGCTTGCGTTTTCAGCCTTCg                                                   <  2:1947656/98‑1 (MQ=255)
                                                gcgggggCCGGGGCGGGGCGGGGGGTGAGGAGGGGGGCAAGTGTCGCCGAGC.T.CGGGGGGGGGCGCTTCTGCCCGCCCCGGAGCTTGCGTTTTCAGCCTTc                                                    <  2:3431550/98‑1 (MQ=255)
                                                gCAGGGGCCGCGGCGAGGACGGGGGTCAGGACGGGGGCAAGTGTCGCCGAGC.T.CG.GGGGGGGCGCTTCTGCCCGCCCCGGAGCTTGCGTTTTCAGCCTTCg                                                   >  2:1814453/1‑101 (MQ=255)
                                                  aGGGGCCGCGGCGAGGAGGGGGGTCAGGACGGGGGCAAGTGTCGCCGAGC.T.CG.GGGGGGGCGCTTCTGCCCGCCCCGGAGCTTGCGTTTTCAGCCTTCGAc                                                 <  2:3104277/101‑1 (MQ=255)
                                                  aGGGGCCGCGGCGAGGACGGGGGTCAGGACGGGGGCAAGTGTAGCCGAGC.T.CG.GGGGGGGCGCTTCTGCCCGCCCCGGAGCTTGCGTTTTCAGCCTTCGAc                                                 >  2:2959827/1‑101 (MQ=255)
                                                   ggggCCGCGGCGAGGACGGGGGTCAGGACGGGGGCAAGTGTCGCCGAGC.T.CG.GGGGGGGCGCTTCTGCCCGCCCCGGAGCTTGCGTTTTCAGCCTTCGAcg                                                <  2:3551614/101‑1 (MQ=255)
                                                     ggCCGGGGCGAGGACGGGGGTCAGGACGGGGGCAAGTGTCGCCGAGC.T.CG.GGGGGGGCGCTTCTGCCCGCCCCGGAGCTTGCGTTTTCAGCCTTCGAcggc                                              <  2:4211948/101‑1 (MQ=255)
                                                      gCCGCGGCGGGGACGGGGGTCAGGCCGGGGGCAAGGGTCGCCGAGC.T.CG.GGGGGGGCGCTTCTACCCGCCCCGGAGCTTGCGTTTTCAGCCTTCGAcggcg                                             <  1:1005167/101‑1 (MQ=255)
                                                        cgcgGCGGGGGCGGGGTTCAGGAGGGGGGCAAGTCTCGCCGGGC.G.CG.GGGGGGGCGCTTCTGCCCGCCCCGGAGCTTGCGTTTTCAGCCTTCGAcggcggc                                           <  2:5153068/101‑1 (MQ=255)
                                                        cgcgGCGAGGACGGGGGTCAGGACGGGGGCAAGTGTCGCCGAGC.T.CG.GGGGGGGCGCTTCTGCCCGCCCCGGAGCTTGCGTTTTCAGCCTTCGAcggcggc                                           >  1:3060314/1‑101 (MQ=255)
                                                           ggcgAGGACGGGGGTCAGGACGGGGGCAAGTGTCGCCGAGC.T.CG.GGGGGGGCGCTTCTGCCCGCCCCGGAGCTTGCGTTTTCAGCCTTCGACggcggcggc                                        >  1:2124228/1‑101 (MQ=255)
                                                            gcgAGGACGGGGGTCAGGACGGGGGCAAGTGTCGCCGAGC.T.CG.GGGGGGGCGCTTCTGCCCGCCCCGGAGCTTGCGTTTTCAGCCTTCGACGgcggcggcg                                       >  1:5774950/1‑101 (MQ=255)
                                                               aGGACGGGGGTCAGGACGGGGGCAAGTGTCGCCGAGC.T.CG.GGGGGGGCGCTTCTGCCCGCCCCGGAGCTTGCGTTTTCAGCCTTCGACGGCGGCGGCGAAg                                    <  2:2746864/101‑1 (MQ=255)
                                                                ggACGGGGGTCAGGACGGGGGCAAGTGTCGCCGAGC.T.CG.GGGGGGGCGCTTCTGCCCGCCCCGGAGCTTGCGTTTTCAGCCTTCGACGGCGGCGGCGAAgg                                   >  2:5126352/1‑101 (MQ=255)
                                                                ggACGGGGGTCAGGACGGGGGCAAGTGTCGCCGAGC.T.CG.GGGGGGGCGCTTCTGCCCGCCCCGGAGCTTGCGTTTTCAGCCTTCGACGGCGGCGGCGAAgg                                   <  2:5607869/101‑1 (MQ=255)
                                                                ggACGGGGGTCAGGACGGGGGCAAGTGTCGCCGAGC.T.CG.GGGGGGGCGCTTCTGCCCGCCCCGGAGCTTGCGTTTTCAGCCTTCGACGGCGGCGGCGAAgg                                   >  2:870438/1‑101 (MQ=255)
                                                                 gACGGGGGTCAGGACGGGGGCAAGTGTCGCCGAGC.T.CG.GGGGGGGCGCTTCTGCCCGCCCCGGAGCTTGCGTTTTCAGCCTTCGACGGCGGCGGCGAAGgc                                  <  1:5450686/101‑1 (MQ=255)
                                                                  gcGGGGGTCAGGACGGGGGCAAGTGTCGCCGAGC.T.CG.GGGGGGGCGCTTCTGCCCGCCCCGGAGCTTGCGTTTTCAGCCTTCGACGGCGGCGGCGAAGgcg                                 <  1:4092286/100‑1 (MQ=255)
                                                                  aCGGGGGTCAGGACGGGGGCAAGTGTCGCCGAGC.T.CG.GGGGGGGCGCTTCTGCCCGCCCCGGAGCTTGCGTTTTCAGCCTTCGACGGCGGCGGCGAAGgcg                                 >  1:2675846/1‑101 (MQ=255)
                                                                    gggggTCGGGGGGGGGGCAGGTCGCGCCGAGC.TGG..GGGGGGGCGCTTCGGCCCGCCCCGGAGCTTGCGTTTTCAGCCTTCGACGGCGGCGGCGAAGtcggt                               <  2:3755324/101‑1 (MQ=255)
                                                                    gggggTCAGGACGGGGGCAAGTGTCGCCGAGC.T.CG.GGGGGGGCGCTTCTGCCCGCCCCGGAGCTTGCGTTTTCAGCCTTCGACGGCGGCGGCGAAGgcggg                               >  1:4670753/1‑100 (MQ=255)
                                                                      gggTCGGGACGGGGGCAAGTGTCGCCGAGC.T.CG.GGGGGGGCGCTTCTGCCCGCCCCGGAGCTTGCGTTTTCAGCCTTCGACGGCGGCGGCGAAGgcggtgc                             <  2:4165878/101‑1 (MQ=255)
                                                                      gggTCAGGACGGGGGCAAGTGTCGCCGAGC.T.CG.GGGGGGGCGCTTCTGCCCGCCCCGGAGCTTGCGTTTTCAGCCTTCGACGGCGGCGGCGAAGgcggtgc                             >  2:4558235/1‑101 (MQ=255)
                                                                          gaGGGCGGGGGAAGGTGTCCCCGAGC.G.CG.GGGGGGGCGCTTTTGCCCGCCCCGGGGCTTGCGTTTTCAGCCTTCGACGGCGGCGGCGAAGGCGGTGCGGTa                         <  1:1542614/100‑1 (MQ=255)
                                                                           aGGACGGGGGCAAGTGTCGCCGAGC.T.CG.GGGGGGGCGCTTCTGCCCGCCCCGGAGCTTGCGTTTTCAGCCTTCGACGGCGGCGGCGAAGGCGGGGCGGTaa                        >  2:1262727/1‑101 (MQ=255)
                                                                            ggAGGGGGGGAGGTGTGGCGGGGC.T.GG.GGGGGGGCGCTTCTGCCCGCCCCGGAGCTTGCGTTTTCAGCCTTCGACGGCGGCGGCGAAGGCGGTGCGGTaaa                       <  2:4945237/101‑1 (MQ=255)
                                                                            ggACGGGGGCAAGTGTCGCCGAGC.T.CG.GGGGGGGCGCTTCTGCCCGCCCCGGAGCTTGCGTTTTCAGCCTTCGACGGCGGCGGCGAAGGCGGTGCGGTaaa                       <  2:4384283/101‑1 (MQ=255)
                                                                             taCGGGGGCAAGTGTCGCCGAGC.T.CG.GGGGGGGCGCTTCTGCCCGCCCCGGAGCTTGCGTTTTCAGCCTTCGACGGCGGCGGCGAAGGCGGTGCGGTAAAg                      <  2:2239239/100‑1 (MQ=255)
                                                                             ggcGGGGGCAAGTGCCGCCGAGCGG.G..GGGGGGGCGCTTCTGCCCGCCCCGGAGCTTGCGTTTTCAGCCTTCGACGGCGGCGGCGAAGGCGGTGCGGTAAAg                      <  1:274785/99‑1 (MQ=255)
                                                                             gACGGGGGCAAGTGTCGCCGAGC.T.CG.GGGGGGGCGCTTCTGCCCGCCCCGGAGCTTGCGTTTTCAGCCTTCGACGGCGGCGGCGAAGGCGGTGCGGTAAAg                      >  1:2007352/1‑101 (MQ=255)
                                                                              aCGGGGGCAAGTGTCGCCGAGC.T.CG.GGGGGGGCGCTTCTGCCCGCCCCGGAGCTTGCGTTTTCAGCCTTCGACGGCGGCGGCGAAGGCGGTGCGGTAAAga                     >  1:2239365/1‑100 (MQ=255)
                                                                              aCGGGGGCAAGTGTCGCCGAGC.T.CG.GGGGGGGCGCTTCTGCCCGCCCCGGAGCTTGCGTTTTCAGCCTTCGACGGCGGCGGCGAAGGCGGTGCGGTAAAGc                     >  2:5269493/1‑101 (MQ=255)
                                                                               cGGGGGCAAGTGTCGCCGAGC.T.CG.GGGGGGGCGCTTCTGCCCGCCCCGGAGCTTGCGTTTTCAGCCTTCGACGGCGGCGGCGAAGGCGGTGCGGTAAAGCg                    >  2:4808987/1‑101 (MQ=255)
                                                                               cGGGGGCAAGTGTCGCCGAGC.T.CG.GGGGGGGCGCTTCTGCCCGCCCCGGAGCTTGCGTTTTCAGCCTTCGACGGCGGCGGCGAAGGCGGGGCGGTAAAGCg                    >  2:5482893/1‑101 (MQ=255)
                                                                                gggggCAAGTGTCGCCGAGC.T.CG.GGGGGGGCGCTTCTGCCCGCCCCGGAGCTTGCGTTTTCAGCCTTCGACGGCGGCGGCGAAGGCGGTGCGGTAAAGCGt                   >  1:866965/1‑101 (MQ=255)
                                                                                   ggCAAGTGTCGCCGGGC.T.CG.GGGGGGGCGCTTCTGCCCGCCCCGGAGCTTGCGTTTTCAGCCTTCGACGGCGGCGGCGAAGGCGGTGCGGTAAAGCGTGCt                <  2:552884/101‑1 (MQ=255)
                                                                                     cAAGTGTCGCCGAGC.T.CG.GGGGGGGCGCTTCTGCCCGCCCCGGAGCTTGCGTTTTCAGCCTTCGACGGCGGCGGCGAAGGCGGTGCGGTAAAGCGTGCTGa              <  2:3086946/101‑1 (MQ=255)
                                                                                     cAAGTGGCGGGGAGC.T.CG.GGGGGGGCGCTTCTGCCCGCCCCGGAGCTTGCGTTTTCAGCCTTCGACGGCGGCGGCGAAGGCGGTGCGGTAAAGCGTGCTGa              <  1:3189398/101‑1 (MQ=255)
                                                                                      aaGTGTCGCCGAGC.T.CG.GGGGGGGCGCTTCTGCCCGCCCCGGAGCTTGCGTTTTCAGCCTTCGACGGCGGCGGCGAAGGCGGTGCGGTAAAGCGTGCTGAg             >  1:4504102/1‑101 (MQ=255)
                                                                                       aGTGTCGCCGAGC.T.CG.GGGGGGGCGCTTCTGCCCGCCCCGGAGCTTGCGTTTTCAGCCTTCGACGGCGGCGGCGAAGGCGGTGCGGTAAAGCGTGCTGAGc            <  1:3426541/101‑1 (MQ=255)
                                                                                       aGTGTCGCCGAGC.T.CG.GGGGGGGCGCTTCTGCCCGCCCCGGAGCTTGCGTTTTCAGCCTTCGACGGCGGCGGCGAAGGCGGTGCGGTAAAGCGTGCTGAGc            >  1:3787613/1‑101 (MQ=255)
                                                                                       aGTGTCGCCGAGC.T.CG.GGGGGGGCGCTTCTGCCCGCCCCGGAGCTTGCGTTTTCAGCCTTCGACGGCGGCGGCGAAGGCGGGGCGGTAAAGCGTGCTGAGc            >  2:2385400/1‑101 (MQ=255)
                                                                                         tgtCGCCGAGC.T.CG.GGGGGGGCGCTTCTGCCCGCCCCGGAGCTTGCGTTTTCAGCCTTCGACGGCGGCGGCGAAGGCGGTGCGGTAAAGCGTGCTGAGctc          <  2:1152202/101‑1 (MQ=255)
                                                                                           tggCCGGGC.T.CG.GGGGGGGCGCTTCTGCCCGCCCCGGAGCTTGCGTTTTCAGCCTTCGACGGCGGCGGCGAAGGCGGTGCGGTAAAGCGTGCTGAGCTCTg        <  2:2296089/99‑1 (MQ=255)
                                                                                              ccGAGC.T.CG.GGGGGGGCGCTTCTGCCCGCCCCGGAGCTTGCGTTTTCAGCCTTCGACGGCGGCGGCGAAGGCGGTGCGGTAAAGCGTGCTGAGCTCTGCCa     <  2:3247963/101‑1 (MQ=255)
                                                                                              ccGAGC.T.CG.GGGGGGGCGCTTCTGCCCGCCCCGGAGCTTGCGTTTTCAGCCTTCGACGGCGGCGGCGAAGGCGGTGCGGTAAAGCGTGCTGAGCTCTGCCa     >  2:2046188/1‑101 (MQ=255)
                                                                                              ccGAGC.T.CG.GGGGGGGCGCTTCTGCCCGCCCCGGAGCTTGCGTTTTCAGCCTTCGACGGCGGCGGCGAAGGCGGTGCGGTAAAGCGTGCTGAGCTCTGCCa     <  2:2088034/101‑1 (MQ=255)
                                                                                              ccGAGC.T.CG.GGGGGGGCCCTTCTGCCCGCCCCGGAGCTTGCGGTTTCAGCCTTCGACGGCGGCGGCGAAGGCGGTGCGGTAAAGCGTGCTGAGCTCTGCCa     <  1:350480/101‑1 (MQ=255)
                                                                                               cGAGC.T.CG.GGGGGGGCGCTTCTGCCCGCCCCGGAGCTTGCGTTTTCAGCCTTCGACGGCGGCGGCGAAGGCGGTGCGGTAAAGCGTGCTGAGCTCTGCCAg    >  1:530910/1‑101 (MQ=255)
                                                                                                gAGC.T.CG.GGGGGGGCGCTTCTGCCCGCCCCGGAGCTTGCGTTTTCAGCCTTCGACGGCGGCGGCGAAGGCGGTGCGGTAAAGCNTGCTGAGCTCTGCCAgc   >  2:1570385/1‑101 (MQ=255)
                                                                                                 aGC.T.CG.GGGGGGGCGCTTCTGCCCGCCCCGGAGCTTGCGTTTTCAGCCTTCGACGGCGGCGGCGAAGGCGGTGCGGTAAAGCGTGCTGAGCTCTGCCAgcg  >  2:5283856/1‑101 (MQ=255)
                                                                                                      |                                                                                                  
AGGCGACAAGAGCCAAGCACCGGTTGCGTGTGCCAAGTGACGCAGCGGGCAGGGGCCGCGGCGAGGACGGGGGTCAGGACGGGGGCAAGTGTCGCCGAGC.T.C..GGGGGGGCGCTTCTGCCCGCCCCGGAGCTTGCGTTTTCAGCCTTCGACGGCGGCGGCGAAGGCGGTGCGGTAAAGCGTGCTGAGCTCTGCCAGCG  >  NC_014034/586816‑587012

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