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% 26.3 196intergenic (‑145/+37)RCAP_rcc00548/purLBolA family protein/phosphoribosylformylglycinamidine synthase II
Reads supporting (aligned to +/- strand):  new base (99/86):  ref base (1/10):  total (100/96)
Fisher's exact test strand distribution p-value = 4.40e-03
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 2.77e-01

Only 100 of 235 total aligned reads displayed.

GGCGACAAGAGCCAAGCACCGGTTGCGTGTGCCAAGTGACGCAGCGGGCAGGGGCCGCGGCGAGGACGGGGGTCAGGACGGGGGCAAGTGTCGCCGAGCT.C..G.GGGGGGCGCTTCTGCCCGCCCCGGAGCTTGCGTTTTCAGCCTTCGACGGCGGCGGCGAAGGCGGTGCGGTAAAGCGTGCTGAGCTCTGCCAGCG  >  NC_014034/586817‑587012
                                                                                                      |                                                                                                 
ggCGACAAGAGCCAAGCACCGGTTGCGTGTGCCAAGTGACGCAGCGGGCAGGGGCCGCGGCGAGGACGGGGGTCAGGACGGGGGCAAGTGTCGCCGAGCT.c                                                                                                    <  1:1633932/101‑1 (MQ=255)
 gCGACAAGAGCCAAGCACCGGTTGCGTGTGCCAAGTGACGCAGCGGGCAGGGGCCGCGGCGAGGACGGGGGTCAGGACGGGGGCAAGTGTCGCCGAGCT.C..g                                                                                                 <  2:4101505/101‑1 (MQ=255)
  cGACAAGAGCCAAGCACCGGTTGCGTGTGCCAAGTGACGCAGCGGGCAGGGGCCGCGGCGAGGACGGGGGTCAGGACGGGGGCAAGGGTCGCCGGGGT.C..g.g                                                                                               >  2:3690418/1‑101 (MQ=255)
   gACAAGAGCCAAGCACCGGTTGCGTGTGCCAAGTGACGCAGCGGGCAGGGGCCGCGGCGAGGACGGGGGTCAGGACGGGGGCAAGTGTCGCCGAGCT.C..g.gg                                                                                              >  1:401593/1‑101 (MQ=255)
   gACAAGAGCCAAGCACCGGTTGCGTGTGCCAAGTGACGCAGCGGGCAGGGGCCGCGGCGAGGACGGGGGTCAGGACGGGGGCAAGTGTCGCCGAGCT.C..g.gg                                                                                              >  2:801743/1‑101 (MQ=255)
     cAAGAGCCAAGCACCGGTTGCGTGTGCCAAGTGACGCAGCGGGCAGGGGCCGCGGCGAGGACGGGGGTCAGGACGGGGGCAAGTGTCGCCGAGCT.C..g.gggg                                                                                            <  2:752081/101‑1 (MQ=255)
       agagCCAAGCACCGGTTGCGTGTGCCAAGTGACGCAGCGGGCAGGGGCCGCGGCGAGGACGGGGGTCAGGACGGGGGCAAGTGTCGCCGAGCT.C..g.gggggg                                                                                          <  1:1077191/101‑1 (MQ=255)
        gagCCAAGCACCGGTTGCGTGTGCCAAGTGACGCAGCGGGCAGGGGCCGCGGCGAGGACGGGGGTCAGGACGGGGGCAAGTGTCGCCGAGCT.C..g.ggggggg                                                                                         >  2:3232564/1‑100 (MQ=255)
         agCCAAGCACCGGTTGCGTGTGCCAAGTGACGCAGCGGGCAGGGGCCGCGGCGAGGACGGGGGTCAGGACGGGGGCAAGTGTCGCCGAGCT.C..g.gggggggc                                                                                        <  1:2191324/101‑3 (MQ=255)
         agCCAAGCACCGGTTGCGTGTGCCAAGTGACGCAGCGGGCAGGGGCCGCGGCGAGGACGGGGGTCAGGACGGGGGCAAGTGTCGCCGAGCT.C..g.gggggggc                                                                                        <  1:3086979/101‑3 (MQ=255)
          gCCAAGCGCGGGTTGCGTGTGCCAAGGGGCGCGGGGGGCGGGGGCCGCGGCGAGGGCGGGGGTCAGGACGGGGGCAAGTGTCGCCGAGCT.C..g.gggggggcg                                                                                       <  1:3694165/101‑4 (MQ=255)
          gCCAAGCACCGGTTGCGTGTGCCAAGTGACGCAGCGGGCAGGGGCCGCGGCGAGGACGGGGGTCAGGACGGGGGCAAGTGTCGCCGAGCT.C..g.gggggggcg                                                                                       <  1:3722137/101‑4 (MQ=255)
            cAAGCACCGGTTGCGTGTGCCAAGTGACGCAGCGGGCAGGGGCCGCGGCGAGGACGGGGGTCAGGACGGGGGCAAGTGTCGCCGAGCT.C..g.gggggggcgct                                                                                     >  1:4277136/1‑96 (MQ=255)
             aaGCGCCGGTGGCGTGTGCCAAGTGACGCAGCGGGCGGGGGCCGCGGGGAGGAGGGGGGTCAGGACGGGGGCAAGTGTCGCCGAGCT.CG.G.GGGGGGCGCtt                                                                                     <  2:389797/101‑1 (MQ=255)
             aaGCCCCGGTTGCGTGTGCCAAGTGACGCAGCGGGCGGGGGCCGGGGCGAGGACGGGGGTCAGGACGGGGGCAAGTGTCGCCGAGCT.CG.G.GGGGGGCGCtt                                                                                     <  2:4420114/101‑1 (MQ=255)
             aaGCACCGGTTGCGTGTGCCAAGTGACGCAGCGGGCAGGGGCCGCGGCGAGGACGGGGGTCAGGACGGGGGCAAGTGTCGCCGAGCT.CG.G.GGGGGGCGCtt                                                                                     >  2:2796089/1‑101 (MQ=255)
             aaGCACCGGTTGCGTGTGCCAAGTGACGCAGCGGGCAGGGGCCGCGGCGAGGACGGGGGTCAGGACGGGGGCAAGTGTCGCCGAGCT.CG.G.GGGGGGCGCtt                                                                                     >  1:2891081/1‑101 (MQ=255)
              aGCACCGGTTGCGTGTGCCAAGTGACGCAGCGGGCAGGGGCCGCGGCGAGGACGGGGGTCAGGACGGGGGCAAGTGTCGCCGAGCT.CG.G.GGGGGGCGCTTc                                                                                    >  1:52855/1‑101 (MQ=255)
                cgccGGTGGCGTGTGCCAAGGGACGCAGCGGGCAGGGGCCGCGGCGAGGACGGGGGTCAGGACGGGGGCAAGTGTCGCCGAGCT.CG.G.GGGGGGCGCTTCTg                                                                                  <  2:559466/99‑1 (MQ=255)
                cACCGGTTGCGTGTGCCAAGTGACGCAGCGGGCAGGGGCCGCGGCGAGGACGGGGGTCAGGACGGGGGCAAGTGTCGCCGAGCT.CG.G.GGGGGGCGCTTCTg                                                                                  >  1:3011563/1‑101 (MQ=255)
                 aCCGGTTGCGTGTGCCAAGTGACGCAGCGGGCAGGGGCCGCGGCGAGGACGGGGGTCAGGACGGGGGCAAGTGTCGCCGAGCT.CG.G.GGGGGGCGCTTCTgc                                                                                 >  2:4028940/1‑101 (MQ=255)
                 aCCGGTTGCGTGTGCCAAGTGACGCAGCGGGCAGGGGCCGCGGCGAGGACGGGGGTCAGGACGGGGGCAAGTGTCGCCGAGCT.CG.G.GGGGGGCGCTTCTgc                                                                                 >  2:1200533/1‑101 (MQ=255)
                  ccGGTTGCGTGTGCCAAGTGACGCAGCGGGCAGGGGCCGCGGCGGGGACGGGGGTCAGGACGGGGGCAAGTGTCGCCGAGCT.CG.G.GGGGGGCGCTTCTgcc                                                                                <  1:1200625/101‑1 (MQ=255)
                   cGGTTGCGTGTGCCAAGTGACGCAGCGGGCAGGGGCCGCGGCGAGGACGGGGGTCAGGACGGGGGCAAGTGTCGCCGAGCT.CG.G.GGGGGGCGCTTCTgccc                                                                               <  1:668117/101‑1 (MQ=255)
                   cGGTTGCGTGTGCCAAGTGACGCAGCGGGCAGGGGCCGCGGCGAGGACGGGGGTCAGGACGGGGGCAAGTGTCGCCGAGCT.CG.G.GGGGGGCGCTTCTgccc                                                                               >  1:837042/1‑101 (MQ=255)
                      ttGCGTGTGCCAAGTGACGCGGCGGGCAGGGGCCGCGGCGAGGACGGGGGTCAGGACGGGGGCAAGTGTCGCCGAGCT.CG.G.GGGGGGCGCTTCTgcccgcc                                                                            <  2:2881023/101‑1 (MQ=255)
                      ttGCGTGTGCCAAGTGACGCAGCGGGCAGGGGCCGCGGCGAGGACGGGGGTCAGGACGGGGGCAAGTGTCGCCGAGCT.CG.G.GGGGGGCGCTTCTgcccgcc                                                                            >  1:3385297/1‑101 (MQ=255)
                      tggCGTGGGCCAAGTGACGCAGCGGGCAGGGGCCGCGGCGAGGACGGGGGTCAGGACGGGGGCAAGTGGCGCCGAGCT.CG.G.GGGGGGCGCTTCTgcccgcc                                                                            <  2:4239592/99‑1 (MQ=255)
                       tGCGTGTGCCAAGTGACGCAGCGGGCGGGGGCCGCGGCGAGGACGGGGGTCAGGACGGGGGCAAGTGTCGCCGAGCT.CG.G.GGGGGGCGCTTCTgcccgccc                                                                           <  1:2959249/101‑1 (MQ=255)
                        gggggcgCCAGGTGACGCAGCGGGGGGGGGCCGCGGGGAGGAGGGGGGTCAGGGGGGGGGCAAGTGTCGCCGAGCT.CG.G.GGGGGGCGCTTCTGCCCGcccc                                                                          <  2:4013299/95‑1 (MQ=255)
                        gCGTGTGCCAAGTGACGCAGCGGGCAGGGGCCGCGGCGAGGACGGGGGTCAGGACGGGGGCAAGTGTCGCCGAGCT.CG.G.GGGGGGCGCTTCTGCCCGcccc                                                                          >  2:88690/1‑101 (MQ=255)
                        gCGTGTGCAAAGTGACGCAGCGGGCAGGGGCCGCGGCGAGGACGGGGGTCAGGACGGGGGCAAGTGTCGCCGAGCT.CG.G.GGGGGGCGCTTCTGCCCGcccc                                                                          <  2:2422027/101‑1 (MQ=255)
                         ggtgtGCCAAGTGACGCAGCGGGCAGGGGCCGCGGCGAGGACGGGGGTCAGGACGGGGGCAAGTGTCGCCGAGCT.CG.G.GGGGGGCGCTTCTGCCCGCCCCg                                                                         <  2:2137776/100‑1 (MQ=255)
                         ggggggCCAAGTGACGCAGCGGGCAGGGGCCGCGGGGAGGACGGGGGTCAGGACGGGGGCAAGTGTCGCCGAGCT.CG.G.GGGGGGCGCTTCTGCCCGCCCCg                                                                         <  1:1152748/96‑1 (MQ=255)
                         cGTGTGCCAAGTGACGCAGCGGGCAGGGGCCGCGGCGAGGACGGGGGTCAGGACGGGGGCAAGTGTCGCCGAGCT.CG.G.GGGGGGCGCTTCTGCCCGCCCCg                                                                         >  1:2728381/1‑101 (MQ=255)
                          gggtgCCAAGGGGCGGGGGGGGCAGGGGCCGCGGGGAGGGCGGGGGTCAGGACGGGGGCAAGGGTCGCCCAGCT.CG.G.GGGGGGCGCTTCTGCCCGCCCCgg                                                                        <  1:4132595/99‑1 (MQ=255)
                            gtgCCAAGTGACGCAGCGGGCAGGGGCCGCGGCGAGGACGGGGGTCAGGACGGGGGCAAGTGTCGCCGAGCT.CG.G.GGGGGGCGCTTCTGCCCGCCCCGGAg                                                                      >  2:3527567/1‑101 (MQ=255)
                             ggCCAAGTGACGCAGCGGGCAGGGGCCGCGGCGAGGACGGGGGTCAGGACGGGGGCAAGTGTCGCCGAGCT.CG.G.GGGGGGCGCTTCTGCCCGCCCCGGAGc                                                                     <  2:695069/100‑1 (MQ=255)
                              gCCAAGTGACGCAGCGGGCAGGGGCCGCGGCGAGGACGGGGGTCAGGACGGGGGCAAGTGTCGCCGAGCT.CG.G.GGGGGGCGCTTCTGCCCGCCCCGGAGCt                                                                    >  1:1779160/1‑101 (MQ=255)
                               ccAAGTGGCGCAGCGGGCAGGGGCCGCGGCGAGGAGGGGGGTGAGGGCGGGGGCAAGTGTCGCCGAGCT.CG.G.GGGGGGCGCTTCTGCCCGCCCCGGAGCtt                                                                   <  1:1901123/101‑1 (MQ=255)
                                cAAGTGACGCAGCGGGGGGGGGCCGCGGGGAGGACGGGGGTCAGGACGGGGGCAAGTGTCGCCGAGCT.CG.G.GGGGGGCGCTTCTGCCCGCCCCGGAGCTTg                                                                  <  1:4537585/101‑1 (MQ=255)
                                 aaGTGACGCAGCGGGCAGGGGCCGCGGCGAGGACGGGGGTCAGGACGGGGGCAAGTGTCGCCGAGCT.CG.G.GGGGGGCGCTTCTGCCCGCCCCGGAGCTTGc                                                                 <  2:2277009/101‑1 (MQ=255)
                                  aGTGACGCAGCGGGCAGGGGCCGCGGCGAGGACGGGGGTCAGGACGGGGGCAAGTGTCGCCGAGCT.CG.G.GGGGGGCGCTTCTGCCCGCCCCGGAGCTTGCg                                                                >  2:1712207/1‑101 (MQ=255)
                                     gACGCAGCGGGGGGGGGCGGGGGGGGGGACGGGGGTGAGGACGGGGGCAAGTGTCGCCGGGCTGG..G.GGGGGGCGCTTCTGCCCGCCCCGGAGCTTGCGttt                                                             <  2:4126607/101‑1 (MQ=255)
                                     gACGCAGCGGGCAGGGGCCGCGGCGAGGACGGGGGTCAGGACGGGGGCAAGTGTCGCCGAGCT.CG.G.GGGGGGCGCTTCTGCCCGCCCCGGAGCTTGCGttt                                                             >  1:1119721/1‑101 (MQ=255)
                                      aCGCAGCGGGCAGGGGCCGCGGCGAGGACGGGGGTCAGGACGGGGGCAAGTGTCGCCGAGCT.CG.G.GGGGGGCGCTTCTGCCCGCCCCGGAGCTTGCGtttt                                                            <  2:3049266/101‑1 (MQ=255)
                                        gCAGCGGGCAGGGGCCGCGGCGAGGACGGGGGTCAGGACGGGGGCAAGTGTCGCCGAGCT.CG.G.GGGGGGCGCTTCTGCCCGCCCCGGAGCTTGCGTTTTCa                                                          >  2:1518812/1‑101 (MQ=255)
                                           gCGGGCAGGGGCCGCGGCGAGGACGGGGGTCAGGACGGGGGCAAGTGTCGCCGAGCT.CG.G.GGGGGGCGCTTCTGCCCGCCCCGGAGCTTGCGTTTTCAGcc                                                       >  1:3884845/1‑101 (MQ=255)
                                             gggCGGGGGCGGCGGCGAGGAGGGGGGTGAGGTCGGGGGCAAGTGGCGCCGAGCT.CG.G.GGGGGGCGCTTCTGCCCGCCCCGGAGCTTGCGTTTTCAGCCtt                                                     <  2:2022423/101‑1 (MQ=255)
                                             gggCAGGGGCCGCGGCGAGGACGGGGGTCAGGACGGGGGCAAGTGTCGCCGAGCT.CG.G.GGGGGGCGCTTCTGCCCGCCCCGGAGCTTGCGTTTTCAGCCtt                                                     >  2:1097698/1‑101 (MQ=255)
                                             gggCAGGGGCCGCGGCGAGGACGGGGGTCAGGACGGGGGCAAGTGTCGCCGAGCT.CG.G.GGGGGGCGCTTCTGCCCGCCCCGGAGCTTGCGTTTTCAGCCtt                                                     <  2:1148291/101‑1 (MQ=255)
                                              ggCGGGGGCGGCGGGGGGGGCGGGGGTCAGGACGGGGGCAAGTGTCGCCGAGCT.CG.G.GGGGGGCGCTTCTGCCCGCCCCGGAGCTTGCGTTTTCAGCCTTc                                                    <  2:1657047/101‑1 (MQ=255)
                                                cAGGGGCCGCGGCGAGGACGGGGGTCAGGACGGGGGCAAGTGTCGCCGAGCT.CG.G.GGGGGGCGCTTCTGCCCGCCCCGGAGCTTGCGTTTTCAGCCTTCGa                                                  <  2:1579506/101‑1 (MQ=255)
                                                 gggggGCGGGGGGAGGGCGGGGGTCAGGACGGGGGCAAGTGTCGCCGAGCT.CG.G.GGGGGGCGCTTCTGCCCGCCCCGGAGCTTGCGTTTTCAGCCTTCGAc                                                 <  2:1336601/100‑1 (MQ=255)
                                                 aGGGGCCGCGGCGAGGACGGGGGTCAGGACGGGGGCAAGTGTCGCCGAGCT.CG.G.GGGGGGCGATTCTGCCCGCCCCGGAGCTTGCGTTTTCAGCCTTCGAc                                                 >  1:3769691/1‑101 (MQ=255)
                                                  ggggCCGCGGCGAGGACGGGGGTCAGGACGGGGGCAAGTGTCGCCGAGCT.CG.G.GGGGGGCGCTTCTGCCCGCCCCGGAGCTTGCGTTTTCAGCCTTCGAcg                                                >  1:3484324/1‑101 (MQ=255)
                                                    ggCGGCGGGGGGGACGGGGGTCAGGGCGGGGGCAAGTGTCGCCGAGCT.CG.G.GGGGGGCGCTTCTGCCCGCCCCGGAGCTTGCGTTTTCAGCCTTCGAcggc                                              <  1:2042911/101‑1 (MQ=255)
                                                    ggCCGCGGCGAGGACGGGGGTCAGGACGGGGGCAAGTGTCGCCGAGCT.CG.G.GGGGGGCGCTTCTGCCCGCCCCGGAGCTTGCGTTTTCAGCCTTCGAcggc                                              <  1:4610020/101‑1 (MQ=255)
                                                       ggcggccAGGTCGGGGGTCGGGAGGGGGGCAAGTGTCGCCGAGCT.CG.G.GGGGGGCGCTTCTGCCCGCCCCGGAGCTTGCGTTTTCAGCCTTCGAcggcggc                                           <  2:1970898/100‑1 (MQ=255)
                                                       cgcgGCGAGGACGGGGGTCAGGACGGGGGCAAGTGTCGCCGAGCT.CG.G.GGGGGGCGCTTCTGCCCGCCCCGGGGCATCCGTTTTCAGCCTTCGAcgggggc                                           >  2:2493010/1‑101 (MQ=255)
                                                        gcggcgAGGACGGGGGTCAGGACGGGGGCAAGTGTCGCCGAGCT.CG.G.GGGGGGCGCTTCTGCCCGCCCCGGAGCTTGCGTTTTCAGCCTTCGAcggcggcg                                          >  1:4475190/1‑101 (MQ=255)
                                                         cggcgAGGACGGGGGTCAGGACGGGGGCAAGTGTCGCCGAGCT.CG.G.GGGGGGCGCTTCTGCCCGCCCCGGAGCTTGCGTTTTCAGCCTTCGAcggcggcgg                                         >  1:1421382/1‑101 (MQ=255)
                                                          ggggAGGGCGGGGGTGAGGGCGGGGGCAAGTGGCCCCGAGCT.CG.G.GGGGGGCGCTTCTGCCCGCCCCGGAGCTTGCGTTTTCAGCCTTCGACggcggcggc                                        <  2:1044612/98‑1 (MQ=255)
                                                           gcgAGGACGGGGGTCAGGACGGGGGCAAGTGTCGCCGGGCT.CG.G.GGGGGGCGCTTCTGCCCGCCCCGGAGCTTGCGTTTTCAGCCTTCGACGgcggcggcg                                       <  2:3482369/101‑1 (MQ=255)
                                                           gcgAGGACGGGGGTCAGGACGGGGGCAAGTGTCGCCGAGCT.CGGG.GGGGGGGGCTTCGGCCCGCCCCGGGGTTGGGGTTTTCACCTTTCGCcggcggcggg                                        >  1:1541078/1‑100 (MQ=255)
                                                           gcgAGGACGGGGGTCAGGACGGGGGCAAGTGTCGCCGAGCT.CG.G.GGGGGGCGCTTCTGCCCGCCCCGGAGCTTGCGTTTTCAGCCTTCGACGgcggcggcg                                       >  2:4281045/1‑101 (MQ=255)
                                                            cgAGGACGGGGGTCAGGACGGGGGCAAGTGTCGCCGAGCT.CG.G.GGGGGGCGCTTCTGCCCGCCCCGGAGCTTGCGTTTTCAGCCTTCGACGGCGGCGGCGa                                      >  2:2526973/1‑101 (MQ=255)
                                                              aGGAGGGGGGTCAGGACGGGGGCAAGTGTCGCCGAGCT.CG.G.GGGGGGCGCTTCTGCCCGCCCCGGAGCTTGCGTTTTCAGCCTTCGACGGCGGCGGCGAAg                                    <  2:2947230/101‑1 (MQ=255)
                                                              aGGACGGGGGTCAGGACGGGGGCAAGGGTCGCCGAGCT.CG.G.GGGGGGCGCTTCTGCCCGCCCCGGAGCTTGCGTTTTCAGCCTTCGACGGCGGCGGCGAAg                                    <  2:3126232/101‑1 (MQ=255)
                                                               ggACGGGGGTCAGGACGGGGGCAAGTGTCGCCGAGCT.CG.G.GGGGGGCGCTTCTGCCCGCCCCGGAGCTTGCGTTTTCAGCCTTCGACGGCGGCGGCGAAgg                                   <  1:801817/101‑1 (MQ=255)
                                                                 ggggtggTAAGGAAGGGGGAAAGTGTCGCCGAGCT.CG.G.GGGGGGCGCTTCTGCCCGCCCCGGAGCTGGCGTTTTCAGCCTTCGACGGCGGCGGCGAAGgcg                                 <  2:210107/96‑1 (MQ=255)
                                                                   gggggTGAGGACGGGGGCAAGTGTCGCCGAGCT.CG.G.GGGGGGCGCTTCTGCCCGCCCCGGAGCTTGCGTTTTCAGCCTTCGACGGCGGCGGCGAAGgcggt                               <  2:3328428/101‑1 (MQ=255)
                                                                   gggggTCAGGACGGGGGCAAGTGTCGCCGAGCT.CG.G.GGGGGGCGCTTCTGCCCGCCCCGGAGCTTGCGTTTTCAGCCTTCGACGGCGGCGGCGAAGgcggt                               >  2:4635067/1‑101 (MQ=255)
                                                                    ggggTCGGGCCGGGGGCAAGTGTCTCCGAGCTGG..G.GGGGGGCGCTTCTGCCCGCCCCGGAGCTTGCGTTTTCAGCCTTCGACGGCGGCGGCGAAGgcggtg                              <  2:641992/101‑1 (MQ=255)
                                                                    ggggTCAGGACGGGGGCAAGTGTCGCCGAGCT.CG.G.GGGGGGCGCTTCTGCCCGCCCCGGAGCTTGCGTTTTCAGCCTTCGACGGCGGCGGCGAAGgcggtg                              >  1:1805180/1‑101 (MQ=255)
                                                                     gtgTCAGTGCGGGGGCAAGTGTCGCCGAGCT.CG.G.GGGGGGCGCTTCTGCCCGCCCCGGAGCTTGCGTTTTCAGCCTTCGACGGCGGCGGCGAAGgcggtgc                             <  2:3769416/99‑1 (MQ=255)
                                                                     gggTCAGGACGGGGGCAAGTGTCGCCGAGCTAG..G.GGGGGGCGCTTCTGCCCGCCCCGGAGCTTGCGTTTTCAGCCTTCGACGGCGGCGGCGAAGgcggtgc                             >  1:4153836/1‑101 (MQ=255)
                                                                     gggGCAGGAAGGGGGCAAGTGTGGAGGAGCT.CG.G.GGGGGGCGCTTCTGCCCGCCCCGGAGCTTGCGTTTTCAGCCTTCGACGGCGGCGGCGAAGgcggtgc                             <  2:2728597/101‑1 (MQ=255)
                                                                      gggcggggggggggCAAGTGTCGCCGGGCT.CG.G.GGGGGGCGCTTCTGCCCGCCCCGGAGCTTGCGTTTTCAGCCTTCGACGGCGGCGGCGAAGgcggtgcg                            <  1:1392572/92‑1 (MQ=255)
                                                                        tCAGGACGGGGGCAAGTGTCGCCGAGCT.CG.G.GGGGGGCGCTTCTGCCCGCCCCGGAGCTTGCGTTTTCAGCCTTCGACGGCGGCGGCGAAGgcggtgcggt                          >  1:3328451/1‑101 (MQ=255)
                                                                        gcgggAGGGGGGCAGGTGTCGCCGAGCG.CG.G.GGGGGGCGCTTCTGCCCGCCCCGGAGCTTGCGTTTTCAGCCTTCGACGTCGGCGGCGAAGgcggtgcggt                          <  1:2643359/98‑1 (MQ=255)
                                                                          aGGACGGGGGCAAGTGTCGCCGAGCT.CG.G.GGGGGGCGCTTCTGCCCGCCCCGGAGCTTGCGTTTTCAGCCTTCGACGGCGGCGGCGAAGGCGGTGCGGTaa                        >  1:3538931/1‑101 (MQ=255)
                                                                              cGGGGGCAAGTGTCGCCGAGCT.CG.G.GGGGGGCGCTTCTGCCCGCCCCGGAGCTTGCGTTTTCAGCCTTCGACGGCGGCGGCGAAGGCGGTGCGGTAAAGCg                    >  1:4059478/1‑101 (MQ=255)
                                                                              cGGGGGCAAGTGTCGCCGAGCT.CG.G.GGGGGGCGCTTCTGCCCGCCCCGGAGCTTGCGTTTTCAGCCTTCGACGGCGGCGGCGAAGGCGGTGCGGTAAAGCg                    <  2:4474895/101‑1 (MQ=255)
                                                                                ggggCAAGTGTCGCCGAGCT.CG.G.GGGGGGCGCTTCTGCCCGCCCCGGAGCTTGCGTTTTCAGCCTTCGACGGCGGCGGCGAAGGCGGTGCGGGAAAGCgcg                  >  1:2515827/1‑99 (MQ=255)
                                                                                 gggCAAGTGTCGCCGAGCTGG..G.GGGGGGCGCTTCTGCCCGCCCCGGAGCTTGCGTTTTCAGCCTTCGACGGCGGCGGCGAAGGCGGTGCGGTAAAGCGTGc                 <  2:423350/101‑1 (MQ=255)
                                                                                 gggCAAGTGTCGCCGAGCT.CG.G.GGGGGGCGCTTCTGCCCGCCCCGGAGCTTGCGTTTTCAGCCTTCGACGGCGGCGGCGAAGGCGGTGCGGTAAAGCGTGc                 >  2:2401198/1‑101 (MQ=255)
                                                                                 gggCAAGTGTCGCCGAGCT.CG.G.GGGGGGCGCTTCTGCCCGCCCCGGAGCTTGCGTTTTCAGCCTTCGACGGCGGCGGCGAAGGCGGTGCGGTAAAGCGTGc                 >  1:2059634/1‑101 (MQ=255)
                                                                                  ggCAAGTGTCGCCGAGCT.CG.G.GGGGGGCGCTTCTGCCCGCCCCGGAGCTTGCGTTTTCAGCCTTCGACGGCGGCGGCGAAGGCGGTGCGGTAAAGCGTGCt                <  2:2076671/101‑1 (MQ=255)
                                                                                     aggtgtCGCCGAGCTGG..G.GGGGGGCGCTTCTGCCCGCCCCGGAGCTTGCGTTTTCAGCCTTCGACGGCGGCGGCGAAGGCGGTGCGGTAAAGCGTGCTGAg             <  1:3205689/99‑1 (MQ=255)
                                                                                     aggtctCGCCGAGCTGG..G.GGGGGGCGCTTCTGCCCGCCCCGGAGCTTGCGTTTTCAGCCTTCGACGGCGGCGGCGAAGGCGGTGCGGTAAAGCGTGCTGAg             <  2:1827827/96‑1 (MQ=255)
                                                                                     aaGTGTCGCCGAGCT.C..GNGGGGGGCGCTTCTGCCCGCCCCGGAGCTTGCGTTTTCAGCCTTCGACGGCGGCGGCGAAGGCGGTGCGGTAAAGCGTGCTGAg             <  2:270391/101‑1 (MQ=255)
                                                                                      aGTGTCGCCGAGCT.CG.G.GGGGGGCGCTTCTGCCCGCCCCGGAGCTTGCGTTTTCAGCCTTCGACGGCGGCGGCGAAGGCGGTGCGGTAAAGCGTGCTGAGc            <  2:52849/101‑1 (MQ=255)
                                                                                      aGTGTCGCCGAGCT.CG.G.GGGGGGCGCTTCTGCCCGCCCCGGAGCTTGCGTTTTCAGCCTTCGACGGCGGCGGCGAAGGCGGGGCGGGAAAGCGTGCTGAGc            >  1:4249197/1‑101 (MQ=255)
                                                                                        tgtCGCCGAGCT.CG.G.GGGGGGCGCTTCTGCCCGCCCCGGAGCTTGCGTTTTCAGCCTTCGACGGCGGCGGCGAAGGCGGTGCGGTAAAGCGTGCTGAGctc          >  1:2484835/1‑101 (MQ=255)
                                                                                        tgtCGCCGAGCT.CG.G.GGGGGGCGCTTCTGCCCGCCCCGGAGCTTGCGTTTTCAGCCTTCGACGGCGGCGGCGAAGGCGGTGCGGTAAAGCGTGCTGAGctc          <  1:1510256/101‑1 (MQ=255)
                                                                                        tgtCGCCGAGCT.CG.G.GGGGGGCGCTTCTGCCCGCCCCGGAGCTTGCGTTTTCAGCCTTCGACGGCGGCGGCGAAGGCGGTGCGGTAAAGCGTGCTGAGctc          <  2:3840868/101‑1 (MQ=255)
                                                                                            gCCGAGCT.CG.G.GGGGGGCGCTTCTGCCCGCCCCGGAGCTTGCGTTTTCAGCCTTCGACGGCGGCGGCGAAGGCGGTGCGGTAAAGCGTGCTGAGCTCTGcc      <  2:3011683/101‑1 (MQ=255)
                                                                                             ccGAGCT.CG.G.GGGGGGCGCTTCTGCCCGCCCCGGAGCTTGCGTTTTCAGCCTTCGACGGCGGCGGCGAAGGCGGTGCGGTAAAGCGTGCTGAGCTCTGCCa     <  1:2371523/101‑1 (MQ=255)
                                                                                              cGAGCT.CG.G.GGGGGGCGCTTCTGCCCGCCCCGGAGCTTGCGTTTTCAGCCTTCGACGGCGGCGGCGAAGGCGGTGCGGTAAAGCGTGCTGAGCTCTGCCAg    >  2:3404947/1‑101 (MQ=255)
                                                                                                      |                                                                                                 
GGCGACAAGAGCCAAGCACCGGTTGCGTGTGCCAAGTGACGCAGCGGGCAGGGGCCGCGGCGAGGACGGGGGTCAGGACGGGGGCAAGTGTCGCCGAGCT.C..G.GGGGGGCGCTTCTGCCCGCCCCGGAGCTTGCGTTTTCAGCCTTCGACGGCGGCGGCGAAGGCGGTGCGGTAAAGCGTGCTGAGCTCTGCCAGCG  >  NC_014034/586817‑587012

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