Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_014034 586,917 +G 100% intergenic (‑145/+37) RCAP_rcc00548 ← / ← purL BolA family protein/phosphoribosylformylglycinamidine synthase II

Read alignment evidence...
  seq id position change freq score reads annotation genes product
*NC_014034586,9171.→G100.0% 37.3 201intergenic (‑145/+37)RCAP_rcc00548/purLBolA family protein/phosphoribosylformylglycinamidine synthase II
Reads supporting (aligned to +/- strand):  new base (91/96):  ref base (0/14):  total (91/110)
Fisher's exact test strand distribution p-value = 1.53e-04
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 6.74e-01

Only 100 of 240 total aligned reads displayed.

GCGACAAGAGCCAAGCACCGGTTGCGTGTGCCAAGTGACGCAGCGGGCAGGGGCCGCGGCGAGGACGGGGGTCAGGACGGGGGCAAGTGTC.GCCGAGC.T.C..GGGGGGGCGCTTCTGCCCGCCCCGGAGCTTGCGTTTTCAGCCTTCGACGGCGGCGGCGAAGGCGGTGCGGTAAAGCGTGCTGAGCTCTGCCAGCG  >  NC_014034/586818‑587012
                                                                                                       |                                                                                                
 cGACAAGAGCCAAGCACCGGTTGCGTGTGCCAAGTGACGCAGCGGGCAGGGGCCGCGGCGAGGACGGGGGTCAGGACGGGGGCAAGTGTC.GCCGAGC.T.C..gg                                                                                               >  2:3049990/1‑101 (MQ=255)
  gACAAGAGCCAAGCACCGGTTGCGTGTGCCAAGTGACGCAGCGGGCAGGGGCCGCGGCGAGGACGGGGGTCAGGACGGGGGCAAGTGTC.GCCGAGC.T.C..ggg                                                                                              >  1:4179436/1‑101 (MQ=255)
  gACAAGAGCCAAGCACCGGTTGCGTGTGCCAAGTGAAGCAGCGGGCAGGGGCCGCGGCGAGGACGGGGGTCAGGACGGGGGCAAGTGTC.GCCGAGC.T.C..ggg                                                                                              >  2:5709917/1‑101 (MQ=255)
     aaGAGCTAAGCACCGGTTGCGTGTGCCAAGTGACGCAGCGGGCAGGGGCCGCGGCGAGGACGGGGGTCAGGACGGGGGCAAGGGTC.GCCGAGC.T.C..gggggg                                                                                           >  2:1374048/1‑101 (MQ=255)
     aaGAGCCAAGCACCGGTTGCGTGTGCCAAGTGACGCAGCGGGCAGGGGCCGCGGCGAGGACGGGGGTCAGGACGGGGGCAAGTGTC.GCCGAGC.T.C..gggggg                                                                                           >  2:3546371/1‑101 (MQ=255)
      agagCCAAGCACCGGTTGCGTGTGCCAAGTGACGCAGCGGGCAGGGGCCGCGGCGAGGACGGGGGTCAGGACGGGGGCAAGTGTC.GCCGAGC.T.C..ggggggg                                                                                          <  1:3256159/101‑1 (MQ=255)
       gagCCAAGCACCGGTTGCGTGTGCCAAGTGACGCAGCGGGCAGGGGCCGCGGCGAGGACGGGGGTCAGGACGGGGGCAAGTGTC.GCCGAGC.T.C..gggggggg                                                                                         >  1:1658573/1‑100 (MQ=255)
        atccAAGCACCGGTTGCGTGTGCCAAGTGACGCAGCGGGCAGGGGCCGCGGCGAGGACGGGGGTCAGGACGGGGGCAAGTGTC.GCCGAGC.T.C..ggggggggc                                                                                        >  1:2454013/3‑99 (MQ=255)
        agCCAAGCACCGGTTGCGTGTGCCAAGTGACGCAGCGGGCAGGGGCCGCGGCGAGGACGGGGGTCAGGACGGGGGCAAGTGTC.GCCGAGC.T.C..ggggggggg                                                                                        >  2:1728015/1‑99 (MQ=255)
        agCCAAGCACCGGTTGCGTGTGCCAAGTGACGCAGCGGGCAGGGGCCGCGGCGAGGACGGGGGTCAGGACGGGGGCAAGTGTC.GCCGAGC.T.C..ggggggggc                                                                                        >  2:3426336/1‑99 (MQ=255)
        agCCAAGCACCGGTTGCGTGTGCCAAGTGACGCAGCGGGCAGGGGCCGCGGCGAGGACGGGGGTCAGGACGGGGGCAAGTGTC.GCCGAGC.T.C..ggggggggc                                                                                        >  1:787942/1‑99 (MQ=255)
           cAAGCACCGGTTGCGTGGGCCAAGTGACGCAGCGGGCAGGGGCCGCGGCGAGGACGGGGGTCAGGACGGGGGCAAGTGTC.GCCGAGC.T.C..ggggggggcgct                                                                                     <  1:4816156/101‑6 (MQ=255)
             aGCACGTGTTGCGTGTGCCGAGTGACGCGGGGGGGGGGGGCCGCGGCGGGGAGGGGGGGGAGGAGGGGGGCAAGTGTC.GCCGAGC.TGG..GGGGGGGCGCTTc                                                                                    <  2:4245189/101‑1 (MQ=255)
             aGCACCGGTTGCGTGTGCCAAGTGACGCAGCGGGCAGGGGCCGCGGCGAGGACGGGGGTCAGGACGGGGGCAAGTGTC.GCCGAGC.T.CG.GGGGGGGCGCTTc                                                                                    >  1:648154/1‑101 (MQ=255)
             aGCACAGGTTGCGTGTGCCAAGTGACGCAGCGGGCAGGGGCAGCGGCGAGGATGGGGGTCAGGCCGGGGGCAAGTGTC.CCCGAGCTT.G..GGGGGGGCGCTTc                                                                                    >  2:3876978/1‑101 (MQ=255)
              gccccGGTTGCGTGTGCCAAGTGACGCGGGGGGCAGGGGCCGCGGCGAGGACGGGGGTCAGGACGGGGGCAAGTGTC.GCCGAGC.T.CG.GGGGGGGCGCTTCt                                                                                   <  2:1462959/98‑1 (MQ=255)
                aCCGGTTGCGTGTGCCAAGTGACGCAGCGGGCAGGGGCCGCGGCGAGGACGGGGGTCAGGACGGGGGCAAGTGTC.GCCGAGC.T.CG.GGGGGGGCGCTTCTgc                                                                                 >  2:517699/1‑101 (MQ=255)
                  cGGTTGCGTGTGCCAAGGGACGCAGCGGGCAGGGGCCGCGGCGAGGACGGGGGTCAGGACGGGGGCAAGTGTC.GCCGAGC.T.CG.GGGGGGGCGCTTCTgccc                                                                               <  1:3649756/101‑1 (MQ=255)
                      tGCGTGTGCCAAGTGGCGCAGCGGGCAGGGGCCGCGGCGAGGACGGGGGTCAGGACGGGGGCAAGTGTC.GCCGAGC.T.CG.GGGGGGGCGCTTCTgcccgccc                                                                           >  1:937061/1‑101 (MQ=255)
                      tGCGTGTGCCAAGTGACGCAGCGGGCAGGGGCCGCGGCGAGGACGGGGGTCAGGACGGGGGCAAGTGTC.GCCGAGC.T.CG.GGGGGGGCGCTTCTgcccgccc                                                                           >  2:5591094/1‑101 (MQ=255)
                      tGCGTGTGCCAAGTGACGCAGCGGGCAGGGGCCGCGGCGAGGACGGGGGTCAGGACGGGGGCAAGTGTC.GCCGAGC.T.CG.GGGGGGGCGCTTCTgcccgccc                                                                           >  1:1152283/1‑101 (MQ=255)
                       gCGTGGGCCAAGTGACGCAGCGGGGGGGGGCCGCGGCGAGGACGGGGGTCAGGACGGGGGCAAGTGTC.GCCGAGC.T.CG.GGGGGGGCGCTTCTGCCCGcccc                                                                          <  2:4009693/101‑1 (MQ=255)
                         gtgtGCCAAGTGACGCAGCGGGCAGGGGCCGCGGCGAGGACGGGGGTCAGGACGGGGGCAAGTGTC.GCCGAGC.T.CG.GGGGGGGCGCTTCTGCCCGCCCCgg                                                                        <  2:1936219/101‑1 (MQ=255)
                            tgCCAAGGGCCGCAGCGGGCAGGGGCCGCGGCGAGGACGGGGGTCAGGACGGGGGCAAGTGTC.GCCGAGC.T.CG.GGGGGGGCGCTTCTGCCCGCCCCGGAGc                                                                     <  2:3103584/101‑1 (MQ=255)
                             gCCAAGGGACGCGGGGGGGGGGGGCGGGGGGGAGGGGGGGGGTCAGGACGGGGGCAAGTGTC.GCCGGGC.T.CG.GGGGGGGCGCTTCTGCCCGCCCCGGAGCt                                                                    <  1:4437571/101‑1 (MQ=255)
                              cggagggACGCAGCGGGCGGGGGCCGCGGCGAGGGCGGGGGTCGGGACGGGGGCAAGGGTCGGCGGCCC.G.G..GGGGGGGCGCTTCTGCCCGCCCCGGAGCtt                                                                   <  1:2785808/95‑1 (MQ=255)
                              cgaaGGGACGCAGCGGGCGGGGGGCGCGGCGGGGCCGGGGGTCAGGCCGGGGGCAAGTGTC.GGAGAGC.TGG..GGGGGGGCGCTTCTGCCCGCCCCGGAGCtt                                                                   <  1:1789442/99‑1 (MQ=255)
                               cAAGGGACGCAGGGGGCGGGGGCCGCGGCGAGGACGGGGGTCAGGACGGGGGCAAGTGTC.GCCGAGC.T.CG.GGGGGGGCGCTTCTGCCCGCCCCGGAGCTTg                                                                  <  2:2342422/101‑1 (MQ=255)
                                aaGTGACGCAGCGGGCAGGGGCCGCGGCGAGGACGGGGGTCAGGACGGGGGCAAGTGTC.GCCGAGC.T.CG.GGGGGGGCGCTTCTGCCCGCCCCGGAGCTTGc                                                                 >  2:2596608/1‑101 (MQ=255)
                                 aGTGACGCAGCGGGCAGGGGCCGCGGCGAGGACGGGGGTCAGGACGGGGGCAAGTGTC.GCCGAGC.T.CG.GGGGGGGCGCTTCTGCCCGCCCCGGAGCTTGCg                                                                >  1:3193285/1‑101 (MQ=255)
                                  gTGACGCAGCGGGCAGGGGCCGCGGCGAGGACGGGGGTCAGGACGGGGGCAAGTGTC.GCCGAGC.T.CG.GGGGGGGCGCTTCTGCCCGCCCCGGAGCTTGCGt                                                               >  2:4823881/1‑101 (MQ=255)
                                  gTGACGCAGCGGGCAGGGGCCGCGGCGAGGACGGGGGTCAGGACGGGGGCAAGTGTC.GCCGAGC.T.CG.GGGGGGGCGCTTCTGCCCGCCCCGGAGCTTGCGt                                                               >  2:2882767/1‑101 (MQ=255)
                                    gACGCAGCGGGCAGGGGCCGCGGCGAGGACGGGGGTCAGGACGGGGGCAAGTGTC.GCCGAGC.T.CG.GGGGGGGCGCTTCTGCCCGCCCCGGAGCTTGCGttt                                                             <  2:4262333/101‑1 (MQ=255)
                                    gACGCAGCGGGCAGGGGCCGCGGCGAGGACGGGGGCCAGGACGGGGGCAAGTGTC.GCCGAGC.G.CG.GGGGGGGCGCTTCTGCCCGCCCCGGAGCTTGCGttt                                                             <  2:2529088/101‑1 (MQ=255)
                                     aCGCAGCGGGCAGGGGCCGCGGCGAGGACGGGGGTCAGGACGGGGGCAAGTGTC.GCCGAGC.T.CG.GGGGGGGCGCTTCTGCCCGCCCCGGAGCTTGCGtttt                                                            <  2:5400544/101‑1 (MQ=255)
                                      ctcAGCGGGCAGGGGCCGCGGCGAGGACGGGGGTCAGGACGGGGGCAAGTGTC.GCCGAGC.T.CG.GGGGGGGCGCTTCTGCCCGCCCCGGAGCTTGCGTTTTc                                                           <  2:1424492/99‑1 (MQ=255)
                                      cGCAGCGGGCAGGGGCCGCGGCGAGGACGGGGGTCAGGACGGGGGCAAGTGTC.GCCGAGC.T.CG.GGGGGGGCGCTTCTGCCCGCCCCGGAGCTTGCGTTTTc                                                           <  1:1379850/101‑1 (MQ=255)
                                       gCAGCGGGCAGGGGCCGCGGCGAGGACGGGGGTCAGGACGGGGGCAAGTGTC.GCCGAGC.T.CG.GGGGGGGCGCTTCTGCCCGCCCCGGAGCTTGCGTTTTCa                                                          >  2:3280033/1‑101 (MQ=255)
                                       gCAGCGGGCAGGGGCCGCGGCGAGGACGGGGGTCAGGACGGGGGCAAGTGTC.GCCGAGC.T.CG.GGGGGGGCGCTTCTGCCCGCCCCGGAGCTTGCGTTTTCa                                                          <  1:5704236/101‑1 (MQ=255)
                                        cAGCGGGCAGGGGCCGCGGCGAGGACGGGGGTCAGGACGGGGGCAAGTGTC.GCCGAGC.T.CG.GGGGGGGCGCTTCTGCCCGCCCCGGAGCTTGCGTTTTCAg                                                         >  1:147775/1‑101 (MQ=255)
                                         aGCGGGCGGGGGCGGCGGCGAGGGCGGGGGACAGGACGGGGGCAAGTGTC.GCCGGGC.T.CG.GGGGGGGCGCTTCTGCCCGCCCCGGAGCTTGCGTTTTCAGc                                                        <  1:5741019/101‑1 (MQ=255)
                                         aGCGGGCAGGGGCCGCGGCGAGGACGGGGGTCAGGACGGGGGCAAGTGTC.GCCGAGC.T.CG.GGGGGGGCGCTTCTGCCCGCCCCGGAGCTTGCGTTTTCAGc                                                        >  2:4916119/1‑101 (MQ=255)
                                           ggggCAGGGGCCGGGGCGAGGACGGGGGTCAGGACGGGGGCAAGTGTC.GCCGAGC.T.CG.GGGGGGGCGCTTCTGCCCGCCCCGGAGCTTGCGTTTTCAGCCt                                                      <  1:2596750/100‑1 (MQ=255)
                                            gggggggggCCGGGGGGAGGGCGGGGGTGAGGGCGGGGGCAAGTGTC.GCCGGGC.G.CG.GGGGGGGCGCTTCTGCCCGCCCCGGAGCTTGCGTTTTCAGCCtt                                                     <  1:2958389/96‑1 (MQ=255)
                                            gggCGGGGGCCGCGGCGAGGACGGGGGTCAGGACGGGGGCAAGTGTC.GCCGAGC.T.CG.GGGGGGGCGCTTCTGCCCGCCCCGGAGCTTGCGTTTTCAGCCtt                                                     <  2:5006348/101‑1 (MQ=255)
                                            gggCAGGGGCCGCGGCGAGGACGGGGGTCAGGACGGGGGCAAGTGTC.GCCGAGC.T.CG.GGGGGGGCGCTTCTGCCCGCCCCGGAGCTTGCGTTTTCAGCCtt                                                     <  1:517749/101‑1 (MQ=255)
                                            gggCAGGGGCCGCGGCGAGGACGGGGGTCAGGACGGGGGCAAGTGTC.GCCGAGC.T.CG.GGGGGGGCGCTTCTGCCCGCCCCGGAGCTTGCGTTTTCAGCCtt                                                     >  2:712231/1‑101 (MQ=255)
                                             ggCGGGGGCCGGGGCGAGGACGGGGGGCAGGGGGGGGGCAAGTGTC.GCCGAGCGG.G..GGGGGGGCGCTTCTGCCCGCCCCGGAGCTTGCGTTTTCAGCCTTc                                                    <  1:682627/101‑1 (MQ=255)
                                              gCAGGGGCCGCGGCGAGGACGGGGGTCAGGACGGGGGCAAGTGTC.GCCGAGC.T.CG.GGGGGGGCGCTTCTGCCCGCCCCGGAGCTTGCGTTTTCAGCCTTCg                                                   >  2:1814453/1‑101 (MQ=255)
                                               cAGGGGCCGCGGCGGGGAGGGGGGTGAGGGCGGGGGCAAGTGTC.GCCGAGC.G.CG.GGGGGGGCGCTTCTGCCCGCCCCGGAGCTTGCGTTTTCAGCCTTCGa                                                  <  1:1799735/101‑1 (MQ=255)
                                               cAGGGGCCGCGGCGAGGACGGGGGTCAGGACGGGGGCAAGTGTC.GCCGAGC.T.CG.GGGGGGGCGCTTCTGCCCGCCCCGGAGCTTGCGTTTTCAGCCTTCgc                                                  >  2:4092039/1‑100 (MQ=255)
                                                gggggCCGCGGCGGGCACGGGGGTCGGGGAGGGGGCAAGTCCC.GCCGAGC.T.CG.GGGGGGGCGCTTTTGCCCGCCCCGGAGCTTGCGTTTTCAGCCTTCGAc                                                 <  2:2741364/100‑1 (MQ=255)
                                                aGGGGCCGCGGCGAGGACGGGGGTCAGGACGGGGGCAAGTGTA.GCCGAGC.T.CG.GGGGGGGCGCTTCTGCCCGCCCCGGAGCTTGCGTTTTCAGCCTTCGAc                                                 >  2:2959827/1‑101 (MQ=255)
                                                   ggCCGGGGCGAGGACGGGGGTCAGGACGGGGGCAAGTGTC.GCCGAGC.T.CG.GGGGGGGCGCTTCTGCCCGCCCCGGAGCTTGCGTTTTCAGCCTTCGAcggc                                              <  2:4211948/101‑1 (MQ=255)
                                                    gCCGCGGCGGGGACGGGGGTCAGGCCGGGGGCAAGGGTC.GCCGAGC.T.CG.GGGGGGGCGCTTCTACCCGCCCCGGAGCTTGCGTTTTCAGCCTTCGAcggcg                                             <  1:1005167/101‑1 (MQ=255)
                                                      cgcgGCGAGGACGGGGGTCAGGACGGGGGCAAGTGTC.GCCGAGC.T.CG.GGGGGGGCGCTTCTGCCCGCCCCGGAGCTTGCGTTTTCAGCCTTCGAcggcggc                                           >  2:3444231/1‑101 (MQ=255)
                                                      cgcgGCGAGGACGGGGGTCAGGACGGGGGCAAGTGTC.GCCGAGC.T.CG.GGGGGGGCGCTTCTGCCCGCCCCGGAGCTTGCGTTTTCAGCCTTCGAcggcggc                                           >  1:3060314/1‑101 (MQ=255)
                                                        cggcgAGGACGGGGGTCAGGACGGGGGCAAGTGTC.GCCGAGC.T.CG.GGGGGGGCGCTTCTGCCCGCCCCGGAGCTTGCGTTTTCAGCCTTCGAcggcggcgg                                         >  1:3887886/1‑101 (MQ=255)
                                                         ggcgAGGACGGGGGTCAGGACGGGGGCAAGTGTC.GCCGAGC.T.CG.GGGGGGGCGCTTCTGCCCGCCCCGGAGCTTGCGTTTTCAGCCTTCGACggcggcggc                                        >  1:2124228/1‑101 (MQ=255)
                                                          gggggcgcGGGGGGCGGGACGGGGGCAAGTGTC.GCCGAGC.T.CG.GGGGGGGCGCTTCTGCCCGCCCCGGAGCTTGCGTTTTCAGCCTTCGACGgcggcggcg                                       <  2:227436/94‑1 (MQ=255)
                                                          gcgAGGACGGGGGTCAGGACGGGGGCAAGTGTC.GCCGAGC.T.CG.GGGGGGGCGCTTCTGCCCGCCCCGGAGCTTGCGTTTTCAGCCTTCGACGgcggcggcg                                       >  1:5774950/1‑101 (MQ=255)
                                                           cgAGGACGGGGGTCAGGACGGGGGCAAGTGTC.GCCGAGC.T.CGGGGGGGGGGGCTTCTCCCCCCCCCGGGGCTGCGGTTTTACCCCTTCGCcggcggggcgg                                       >  1:2983587/1‑96 (MQ=255)
                                                            gAGGACGGGGGTCAGGACGGGGGCAAGTGTC.GCCGAGC.T.CG.GGGGGGGCGCTTCTGCCCGCCCCGGAGCTTGCGTTTTCAGCCTTCGACGGCGGCGGCGaa                                     >  1:4707016/1‑101 (MQ=255)
                                                            gAGGACGGGGGTCAGGACGGGGGCAAGTGTC.GCCGAGC.T.CG.GGGGGGGCGCTTCTGCCCGCCCCGGAGCTTGCGTTTTCAGCCTTCGACGGCGGCGCCGaa                                     >  2:3755145/1‑101 (MQ=255)
                                                              ggACGGGGGTCAGGACGGGGGCAAGTGTC.GCCGAGC.T.CG.GGGGGGGCGCTTCTGCCCGCCCCGGAGCTTGCGTTTTCAGCCTTCGACGGCGGCGGCGAAgg                                   >  2:5126352/1‑101 (MQ=255)
                                                              ggACGGGGGTCAGGACGGGGGCAAGTGTC.GCCGAGC.T.CG.GGGGGGGCGCTTCTGCCCGCCCCGGAGCTTGCGTTTTCAGCCTTCGACGGCGGCGGCGAAgg                                   <  2:5607869/101‑1 (MQ=255)
                                                              ggACGGGGGTCAGGACGGGGGCAAGTGTC.GCCGAGC.T.CG.GGGGGGGCGCTTCTGCCCGCCCCGGAGCTTGCGTTTTCAGCCTTCGACGGCGGCGGCGAAgg                                   >  2:870438/1‑101 (MQ=255)
                                                               gACGGGGGTCAGGACGGGGGCAAGTGTC.GCCGAGC.T.CG.GGGGGGGCGCTTCTGCCCGCCCCGGAGCTTGCGTTTTCAGCCTTCGACGGCGGCGGCGAAGgc                                  <  1:5450686/101‑1 (MQ=255)
                                                                gcGGGGGTCAGGACGGGGGCAAGTGTC.GCCGAGC.T.CG.GGGGGGGCGCTTCTGCCCGCCCCGGAGCTTGCGTTTTCAGCCTTCGACGGCGGCGGCGAAGgcg                                 <  1:4092286/100‑1 (MQ=255)
                                                                    gggTCAGGACGGGGGCAAGTGTC.GCCGAGC.T.CG.GGGGGGGCGCTTCTGCCCGCCCCGGAGCTTGCGTTTTCAGCCTTCGACGGCGGCGGCGAAGgcggtgc                             >  2:4558235/1‑101 (MQ=255)
                                                                    gggGCAGGGCGGGGGCAAGTGTC.TCCGAGC.T.CG.GGGGGGGCGCTTCTGCCCGCCCCGGAGCTTGCGTTTTCAGCCTTCGACGGCGGCGGCGAAGgcggtgc                             <  2:1598977/101‑1 (MQ=255)
                                                                      gTCAGGACGGGGGCAAGTGTC.GCCGAGC.T.CG.GGGGGGGCGCTTCTGCCCGCCCCGGAGCTTGCGTTTTCAGCCTTCGACGGCGGCGGCGAAGgcggtgcgg                           <  1:3440048/101‑1 (MQ=255)
                                                                        gaGGGCGGGGGAAGGTGTC.CCCGAGC.G.CG.GGGGGGGCGCTTTTGCCCGCCCCGGGGCTTGCGTTTTCAGCCTTCGACGGCGGCGGCGAAGGCGGTGCGGTa                         <  1:1542614/100‑1 (MQ=255)
                                                                         aGGACGGGGGCAAGTGTC.GCCGAGC.T.CG.GGGGGGGCGCTTCTGCCCGCCCCGGAGCTTGCGTTTTCAGCCTTCGACGGCGGCGGCGAAGGCGGGGCGGTaa                        >  2:1262727/1‑101 (MQ=255)
                                                                          ggAGGGGGGTAAGTTTC.GCCGGGG.T.CG.GGGGGGGCGCTTCTGCCCGCCCCGGAGCTTGCGTTTTCAGCCTTCGACGGCGGCGGCGAAGgcggtgcagatcg                       <  2:2076027/101‑7 (MQ=255)
                                                                          ggACGGGGGCAAGTGTC.GCCGAGC.T.CG.GGGGGGGCGCTTCTGCCCGCCCCGGAGCTTGCGTTTTCAGCCTTCGACGGCGGCGGCGAAGGCGGTGCGGTaaa                       <  2:4384283/101‑1 (MQ=255)
                                                                           taCGGGGGCAAGTGTC.GCCGAGC.T.CG.GGGGGGGCGCTTCTGCCCGCCCCGGAGCTTGCGTTTTCAGCCTTCGACGGCGGCGGCGAAGGCGGTGCGGTAAAg                      <  2:2239239/100‑1 (MQ=255)
                                                                           ggcGGGGGCAAGTGCC.GCCGAGCGG.G..GGGGGGGCGCTTCTGCCCGCCCCGGAGCTTGCGTTTTCAGCCTTCGACGGCGGCGGCGAAGGCGGTGCGGTAAAg                      <  1:274785/99‑1 (MQ=255)
                                                                           gACGGGGGCAAGTGTC.GCCGAGC.T.CG.GGGGGGGCGCTTCTGCCCGCCCCGGAGCTTGCGTTTTCAGCCTTCGACGGCGGCGGCGAAGGCGGTGCGGTAAAg                      >  1:2007352/1‑101 (MQ=255)
                                                                            aCGGGGGCAAGTGTC.GCCGAGC.T.CG.GGGGGGGCGCTTCTGCCCGCCCCGGAGCTTGCGTTTTCAGCCTTCGACGGCGGCGGCGAAGGCGGTGCGGTAAAga                     >  1:2239365/1‑100 (MQ=255)
                                                                            aCGGGGGCAAGTGTC.GCCGAGC.T.CG.GGGGGGGCGCTTCTGCCCGCCCCGGAGCTTGCGTTTTCAGCCTTCGACGGCGGCGGCGAAGGCGGTGCGGTAAAGc                     >  2:5269493/1‑101 (MQ=255)
                                                                             cGGGGGCAAGTGTC.GCCGAGC.T.CG.GGGGGGGCGCTTCTGCCCGCCCCGGAGCTTGCGTTTTCAGCCTTCGACGGCGGCGGCGAAGGCGGTGCGGTAAAGCg                    >  2:4808987/1‑101 (MQ=255)
                                                                             cGGGGGCAAGTGTC.GCCGAGC.T.CG.GGGGGGGCGCTTCTGCCCGCCCCGGAGCTTGCGTTTTCAGCCTTCGACGGCGGCGGCGAAGGCGGGGCGGTAAAGCg                    >  2:5482893/1‑101 (MQ=255)
                                                                              gggggCAAGTGTC.GCCGAGC.T.CG.GGGGGGGCGCTTCTGCCCGCCCCGGAGCTTGCGTTTTCAGCCTTCGACGGCGGCGGCGAAGGCGGTGCGGTAAAGCGt                   >  1:866965/1‑101 (MQ=255)
                                                                               ggggCAAGTGTC.CCCGAGC.T.CG.GGGGGGGCGCTTCTGCCCGCCCCGGAGCTTGCGTTTTCAGCCTTCGACGGCGGCGGCGAAGGCGGTGCGGTAAAGCGTg                  <  2:2817462/101‑1 (MQ=255)
                                                                                gggCAAGTGTC.GCCGAGC.T.CG.GGGGGGGCGCTTCTGCCCGCCCCGGAGCTTGCGTTTTCAGCCTTCGACGGCGGCGGCGAAGGCGGTGCGGTAAAGCGTGc                 <  2:3226014/101‑1 (MQ=255)
                                                                                    aaGTGTC.GCCGAGC.T.CG.GGGGGGGCGCTTCTGCCCGCCCCGGAGCTTGCGTTTTCAGCCTTCGACGGCGGCGGCGAAGGCGGTGCGGTAAAGCGTGCTGAg             >  1:4504102/1‑101 (MQ=255)
                                                                                    aaGTGTC.GCCGAGC.T.CG.GGGGGGGCGCTTCTGCCCGCCCCGGAGCTTGCGTTTTCAGCCTTCGACGGCGGCGGCGAAGGCGGTGCGGTAAAGCGTGCTGAg             <  2:3706318/101‑1 (MQ=255)
                                                                                     aGTGTC.GCCGAGC.T.CG.GGGGGGGCGCTTCTGCCCGCCCCGGAGCTTGCGTTTTCAGCCTTCGACGGCGGCGGCGAAGGCGGTGCGGTAAAGCGTGCTGAGc            >  1:3787613/1‑101 (MQ=255)
                                                                                     aGTGTC.GCCGAGC.T.CG.GGGGGGGCGCTTCTGCCCGCCCCGGAGCTTGCGTTTTCAGCCTTCGACGGCGGCGGCGAAGGCGGGGCGGTAAAGCGTGCTGAGc            >  2:2385400/1‑101 (MQ=255)
                                                                                      gtgtC.GCCGAGC.T.CG.GGGGGGGCGCTTCTGCCCGCCCCGGAGCTTGCGTTTTCAGCCTTCGACGGCGGCGGCGAAGGCGGTGCGGTAAAGCGTGCTGAGct           <  1:3759581/101‑1 (MQ=255)
                                                                                       tgtC.GCCGAGC.T.CG.GGGGGGGCGCTTCTGCCCGCCCCGGAGCTTGCGTTTTCAGCCTTCGACGGCGGCGGCGAAGGCGGTGCGGTAAAGCGTGCTGAGctc          <  2:1152202/101‑1 (MQ=255)
                                                                                         gc.GCCGGGC.T.CG.GGGGGGGCGCTTCTGCCCGCCCCGGAGCTTGCGTTTTCAGCCTTCGACGGCGGCGGCGAAGGCGGTGCGGTAAAGCGTGCTGAGCTCTg        <  2:248298/100‑1 (MQ=255)
                                                                                             ccGAGC.T.CG.GGGGGGGCGCTTCTGCCCGCCCCGGAGCTTGCGTTTTCAGCCTTCGACGGCGGCGGCGAAGGCGGTGCGGTAAAGCGTGCTGAGCTCTGCCa     >  2:2046188/1‑101 (MQ=255)
                                                                                             ccGAGC.T.CG.GGGGGGGCCCTTCTGCCCGCCCCGGAGCTTGCGGTTTCAGCCTTCGACGGCGGCGGCGAAGGCGGTGCGGTAAAGCGTGCTGAGCTCTGCCa     <  1:350480/101‑1 (MQ=255)
                                                                                              cGAGC.T.CG.GGGGGGGCGCTTCTGCCCGCCCCGGAGCTTGCGTTTTCAGCCTTCGACGGCGGCGGCGAAGGCGGTGCGGTAAAGCGTGCTGAGCTCTGCCAg    >  2:461923/1‑101 (MQ=255)
                                                                                              cGAGC.T.CG.GGGGGGGCGCTTCTGCCCGCCCCGGAGCTTGCGTTTTCAGCCTTCGACGGCGGCGGCGAAGGCGGTGCGGTAAAGCGTGCTGAGCTCTGCCAg    >  1:530910/1‑101 (MQ=255)
                                                                                               gAGC.T.CG.GGGGGGGCGCTTCTGCCCGCCCCGGAGCTTGCGTTTTCAGCCTTCGACGGCGGCGGCGAAGGCGGTGCGGTAAAGCNTGCTGAGCTCTGCCAgc   >  2:1570385/1‑101 (MQ=255)
                                                                                               gAGC.T.CG.GGGGGGGCGCTTCTGCCCGCCCCGGAGCTTGCGTTTTCAGCCTTCGACGGCGGCGGCGAAGGCGGTGCGGTAAAGCGTGCTGAGCTCTGCCAgc   >  2:4430812/1‑101 (MQ=255)
                                                                                                aGC.T.CG.GGGGGGGCGCTTCTGCCCGCCCCGGAGCTTGCGTTTTCAGCCTTCGACGGCGGCGGCGAAGGCGGTGCGGTAAAGCGTGCTGAGCTCTGCCAgcg  >  2:5283856/1‑101 (MQ=255)
                                                                                                       |                                                                                                
GCGACAAGAGCCAAGCACCGGTTGCGTGTGCCAAGTGACGCAGCGGGCAGGGGCCGCGGCGAGGACGGGGGTCAGGACGGGGGCAAGTGTC.GCCGAGC.T.C..GGGGGGGCGCTTCTGCCCGCCCCGGAGCTTGCGTTTTCAGCCTTCGACGGCGGCGGCGAAGGCGGTGCGGTAAAGCGTGCTGAGCTCTGCCAGCG  >  NC_014034/586818‑587012

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