Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_014034 586,916 +G 3.2% 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 3.2% 8.2 ‑3.9 0.0154intergenic (‑144/+38)RCAP_rcc00548/purLBolA family protein/phosphoribosylformylglycinamidine synthase II
Reads supporting (aligned to +/- strand):  new base (0/5):  ref base (82/67):  total (82/72)
Fisher's exact test strand distribution p-value = 2.07e-02
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 1.00e+00

Only 100 of 192 total aligned reads displayed.

AGGCGACAAGAGCCAAGCACCGGTTGCGTGTGCCAAGTGACGCAGCGGGCAGGGGCCGCGGCGAGGACGGGGGTCAGGACGGGGGCAAGTGTCGCCGAGCT.C.GGGGGG.GCGCTTCTGCCCGCCCCGGAGCTTGCGTTTTCAGCCTTCGACGGCGGCGGCGAAGGCGGTGCGGTAAAGCGTGCTGAGCTCTGCCAGCG  >  NC_014034/586816‑587012
                                                                                                     |                                                                                                  
agGCGACAAGAGCCAAGCCCCGGTTGCGTGTGCCAAGTGACGCAGCGGGCAGGGGCCGCGGCGAGGACGGGGGTCAGGACGGGGGCAAGTGTCGCCGAGCt                                                                                                     <  2:2222189/101‑1 (MQ=255)
agGCGACAAGAGCCAAGCACCGGTTGCGTGTGCCAAGTGACGCAGCGGGCAGGGGCCGCGGCGAGGACGGGGGTCAGGACGGGGGCAAGTGTCGCCGAGCt                                                                                                     >  2:3080516/1‑101 (MQ=255)
 ggCGACAAGAGCCAAGCACCGGTTGCGTGTGCCAAGTGACGCAGCGGGCAGGGGCCGCGGCGAGGACGGGGGTCAGGACGGGGGCAAGTGTCGCCGAGCT.c                                                                                                   >  1:4089661/1‑101 (MQ=255)
 ggCGACAAGAGCCAAGCACCGGTTGCGTGTGCCAAGTGACGCAGCGGGCAGGGGCCGCGGCGAGGACGGGGGTCAGGACGGGGGCAAGTGTCGCCGAGCT.c                                                                                                   <  2:774262/101‑1 (MQ=255)
  gCGACAAGAGCCAAGCACCGGTTGCGTGTGCCAAGTGACGCAGCGGGCAGGGGCCGCGGCGAGGACGGGGGTCAGGACGGGGGCAAGTGTCGCCGAGCT.C.g                                                                                                 >  1:1859001/1‑101 (MQ=255)
   cGACAAGAGCCAAGCACCGGTTGCGTGTGCCAAGTGACGCAGCGGGCAGGGGCCGCGGCGAGGACGGGGGTCAGGACGGGGGCAAGTGTCGCCGAGCT.C.gg                                                                                                <  1:1021562/101‑1 (MQ=255)
   cGACAAGAGCCAAGCACCGGTTGCGTGTGCCAAGTGACGCAGCGGGCAGGGGCCGCGGCGAGGACGGGGGTCAGGACGGGGGCAAGTGTCGCCGAGCT.C.gg                                                                                                >  2:470785/1‑101 (MQ=255)
      cAAGAGCCAAGCACCGGTTGCGTGTGCCAAGTGACGCAGCGGGCAGGGGCCGCGGCGAGGACGGGGGTCAGGACGGGGGCAAGTGTCGCCGAGCT.C.ggggg                                                                                             <  2:1051645/101‑1 (MQ=255)
      cAAGAGCCAAGCACCGGTTGCGTGTGCCAAGTGACGCAGCGGGCAGGGGCCGCGGCGAGGACGGGGGTCAGGACGGGGGCAAGTGTCGCCGAGCT.C.ggggg                                                                                             >  1:1449963/1‑101 (MQ=255)
       aaGAGCCAAGCACCGGTTGCGTGTGCCAAGTGACGCAGCGGGCAGGGGCCGCGGCGAGGACGGGGGTCAGGACGGGGGCAAGTGTCGCCGAGCT.C.gggggg                                                                                            >  1:979381/1‑101 (MQ=255)
       aaGAGCCAAGCACCGGTTGCGTGTGCCAAGTGACGCAGCGGGCAGGGGCCGCGGCGAGGACGGGGGTCAGGACGGGGGCAAGTGTCGCCGAGCT.C.gggggg                                                                                            <  1:2486138/101‑1 (MQ=255)
        agagCCAAGCACCGGTTGCGTGTGCCAAGTGACGCAGCGGGCAGGGGCCGCGGCGAGGACGGGGGTCAGGACGGGGGCAAGTGTCGCCGAGCT.C.gggggg.g                                                                                          >  2:3680588/1‑101 (MQ=255)
         gagCCAAGCACCGGTTGCGTGTGCCAAGTGACGCAGCGGGCAGGGGCCGCGGCGAGGACGGGGGTCAGGACGGGGGCAAGTGTCGCCGAGCT.C.gggggg.gg                                                                                         >  2:3882430/1‑100 (MQ=255)
             cAAGCACCGGTTGCGTGTGCCAAGTGACGCAGCGGGCAGGGGCCGCGGCGAGGACGGGGGTCAGGACGGGGGCAAGTGTCGCCGAGCT.C.gggggg.ggcgct                                                                                     <  2:1826133/101‑6 (MQ=255)
              aaGCACCGGTTGCGTGTGCCAAGTGACGCAGCGGGCAGGGGCCGCGGCGAGGACGGGGGTCAGGACGGGGGCAAGTGTCGCCGAGCT.CGGGGGGG.GCGCtt                                                                                     >  1:475308/1‑101 (MQ=255)
               aGCCCCGGTTGCGTGTGCCAAGGGACGCGGCGGGGGGGGGCCGGGGCGGGGACGGGGGTCAGGAGGGGGGCAAGTGTCGCCGAGCT.CGGGGGGG.GCGCTTc                                                                                    <  1:2391333/101‑1 (MQ=255)
               aGCACCGGTTGCGTGTGCCAAGTGACGCAGCGGGCAGGGGCCGCGGCGAGGACGGGGGTCAGGACGGGGGCAAGTGTCGCCGAGCT.CGGGGGGG.GCGCTTc                                                                                    >  2:1149646/1‑101 (MQ=255)
               aGCACCGGTTGCGTGTGCCAAGTGACGCAGCGGGCAGGGGCCGCGGCGAGGACGGGGGTCAGGACGGGGGCAAGTGTCGCCGAGCT.CGGGGGGG.GCGCTTc                                                                                    >  2:2154433/1‑101 (MQ=255)
                gccccGGTTGCGTGGGCCAAGTGACGCAGCGGGCGGGGGCCGCGGCGGGGGCGGGGGTGAGGACGGGGGCAAGTGTCGCCGAGCT.CGGGGGGG.GCGCTTCt                                                                                   <  1:3664928/98‑1 (MQ=255)
                    cGGTTGCGTGTGCCAAGTGACGCAGCGGGCAGGGGCCGCGGCGAGGACGGGGGTCAGGACGGGGGCAAGTGTCGCCGAGCT.CGGGGGGG.GCGCTTCTgccc                                                                               <  2:1858877/101‑1 (MQ=255)
                      gTTGCGTGTGCCAAGTGACGCAGCGGGCAGGGGCCGCGGCGAGGACGGGGGTCAGGACGGGGGCAAGTGTCGCCGAGCT.CGGGGGGG.GCGCTTCTgcccgc                                                                             <  1:2782000/101‑1 (MQ=255)
                       ttGCGTGTGCCAAGGGACGCAGCGGGCAGGGGCCGCGGCGAGGACGGGGGTCAGGACGGGGGCAAGTGTCGCCGAGCT.CGGGGGGG.GCGCTTCTgcccgcc                                                                            <  2:361865/101‑1 (MQ=255)
                        tGCGTGTGCCAAGTGACGCAGCGGGCAGGGGCCGCGGCGAGGGCGGGGGTCAGGACGGGGGCAAGTGTCGCCGAGCT.CGGGGGGG.GCGCTTCTgcccgccc                                                                           >  1:431900/1‑101 (MQ=255)
                        tGCGTGTGCCAAGGGACGCGGCGGGCGGGGGCCGCGGCGAGGACGGGGGTCAGGACGGGGGCAAGTGTCGCCGAGCT.CGGGGGGG.GCGCTTCTgcccgccc                                                                           <  2:2290340/101‑1 (MQ=255)
                         gCGTGGGCCAAGTGACGCAGCGGGCAGGGGCCGCGGCGAGGACGGGGGTCAGGACGGGGGCAAGTGTCGCCGAGCT.CGGGGGGG.GCGCTTCTGCCCGcccc                                                                          <  2:149408/101‑1 (MQ=255)
                          ggtggGCCAAGGGGCGCAGGGGGCAGGGGCCGCGGCGAGGAGGGGGGGCGGGAGGGGGGCAAGTGTCGCCGAGCT.CGGGGGGG.GCGCTTCTGCCCGCCCCg                                                                         <  2:2474131/100‑1 (MQ=255)
                              tgCCAAGTGAGGCAGCGGGCAGGGGCCGCGGCGAGGACGGGGGCCAGGAGGGGGGCAAGTGTGGCCGAGCT.CGGGGGGG.GCGGTTCTGCCCGCCCCGGAGc                                                                     >  1:614429/1‑101 (MQ=255)
                               gCCAAGTGACGCAGCGGGCAGGGGCCGCGGCGAGGACGGGGGTCAGGACGGGGGCAAGTGTCGCCGAGCT.CGGGGGGG.GCGCTTCTGCCCGCCCCGGAGCt                                                                    >  1:2671885/1‑101 (MQ=255)
                                 cAAGTGACGCAGCGGGCAGGGGCCGCGGCGAGGACGGGGGTCAGGACGGGGGCAAGTGTCGCCGAGCT.CGGGGGGG.GCGCTTCTGCCCGCCCCGGAGCTTg                                                                  <  2:245353/101‑1 (MQ=255)
                                   aGTGACGCAGCGGGCAGGGGCCGCGGCGAGGACGGGGGTCAGGACGGGGGCAAGTGTCGCCGAGCT.CGGGGGGG.GCGCTTCTGCCCGCCCCGGAGCTTGCg                                                                >  1:377894/1‑101 (MQ=255)
                                   aGTGACGCAGCGGGCAGGGGCCGCGGCGAGGACGGGGGTCAGGACGGGGGCAAGTGTCGCCGAGCT.CGGGGGGG.GCGCTTCTGCCCGCCCCGGAGCTTGCg                                                                <  2:2431689/101‑1 (MQ=255)
                                      ggcGCAGCGGGCAGGGGCCGCGGCGAGGACGGGGGGCAGGACGGGGGCAAGTGTCGCCGAGCT.CGGGGGGG.GCGCTTCTGCCCGCCCCGGAGCTTGCGttt                                                             <  1:4443134/99‑1 (MQ=255)
                                      gACGGAGGGGGCGGGGGCCGCGGCGAGGACGGGGGTCGGGACGGGGGCAAGTGTCGCCGAGCT.CGGGGGGG.GCGCTTCTGCCCGCCCCGGAGCTTGCGttt                                                             <  1:4426258/101‑1 (MQ=255)
                                      gACGCAGCGGGCAGGGGCCGCGGCGAGGACGGGGGTCAGGACGGGGGCAAGTGTCGCCGAGCT.CGGGGGGG.GCGCTTCTGCCCGCCCCGGAGCTTGCGttt                                                             >  1:1192827/1‑101 (MQ=255)
                                       gcGCAGCGGGCAGGGGCCGCGGGGAGGACGGGGGTCAGGACGGGGGCAAGTGTCGCCGAGCT.CGGGGGGG.GCGCTTCTGCCCGCCCCGGAGCTTGCGtttt                                                            <  1:2209513/100‑1 (MQ=255)
                                       gcGCAGCGGGCAGGGGCCGCGGCGAGGACGGGGGTCAGGGCGGGGGCAAGTGTCGCCGAGCT.CGGGGGGG.GCGCTTCTGCCCGCCCCGGAGCTTGCGtttt                                                            <  1:766551/100‑1 (MQ=255)
                                        ggCAGCGGGCGGGGGCCGCGGCGAGGACGGGGGTCAGGACGGGGGCAAGTGTCGCCGAGCT.CGGGGGGG.GCGCTTCTGCCCGCCCCGGAGCTTGCGTTTTc                                                           <  2:1194092/100‑1 (MQ=255)
                                        cGCAGCGGGCAGGGGCCGCGGCGAGGACGGGGGTCAGGACGGGGGCAAGTGTCGCCGAGCT.CGGGGGGG.GCGCTTCTGCCCGCCCCGGAGCTTGCGTTTTc                                                           <  2:380038/101‑1 (MQ=255)
                                         gCAGCGGGCAGGGCCCGCGGCGAGGACGGGGGTCAGGACGGGGGCAAGTGTCGCCGAGCT.CGGGGGGG.GCGCTTCTGCCCGCCCCGGAGCTTGCGTTTTCa                                                          <  2:3574860/101‑1 (MQ=255)
                                           aGCGGGCGGGGGCCGCGGGGAGGACGGGGGTCAGGACGGGGGCAAGTGTCGCCGAGCT.CGGGGGGG.GCGCTTCTGCCCGCCCCGGAGCTTGCGTTTTCAGc                                                        <  1:2701423/101‑1 (MQ=255)
                                           aGCGGGCAGGGGCCGCGGCGAGGACGGGGGTCAGGACGGGGGCAAGTGTCGCCGAGCT.CGGGGGGG.GCGCTTCTGCCCGCCCCGGAGCTTGCGTTTTCAGc                                                        >  2:1605031/1‑101 (MQ=255)
                                            gCGGGCGGGGGGCGGGGCGAGGGCGGGGCTCAGGACGGGGGCAAGTGCCGTCGAGCT.CGGGGGGG.GCGCTTTTGCCCGCCCCGGAGCTTGCGTTTTCAGcc                                                       <  1:3585862/101‑1 (MQ=255)
                                             ggggGAGGGGCCGCGGCGAGGACGGGGGTCAGGACGGGGGCAGGTGGCGCCGAGGTGG.GGGGGG.GCGCTTCTGCCCGCCCCGGAGCTTGCGTTTTCAGCCt                                                      <  2:645255/100‑1 (MQ=255)
                                              gggggggggCCGGGGGGGGGAGGGGGGTCGGGACGGGGGCAAGTGTCGCCGGGCT.CGGGGGGG.GCGCTTCTGCCCGCCCCGGAGCTTGCGTTTTCAGCCtt                                                     <  1:1178561/96‑1 (MQ=255)
                                                gggggggCCGGGGCGAGGAGGGGGGTCGGGCCGGGGGCAAGTGTCGCCGGGCG.CGGGGGGG.GCGCTTCTGCCCGCCCCGGAGCTTGCGTTTTCAGCCTTCg                                                   <  1:3386677/98‑1 (MQ=255)
                                                gCAGGGGCCGCGGCGAGGACGGGGGTCAGGACGGGGGCAAGTGTCGCCGAGCT.CGGGGGGG.GCGCTTCTGCCCGCCCCGGAGCTTGCGTTTTCAGCCTTCg                                                   >  2:4022597/1‑101 (MQ=255)
                                                gCAGGGGCCGCGGCGAGGACGGGGGTCAGGACGGGGGCAAGTGTCGCCGAGCT.CGGGGGGG.GCGCTTCTGCCCGCCCCGGAGCTTGCGTTTTCAGCCTTCg                                                   >  2:3133742/1‑101 (MQ=255)
                                                 gaGGGGCGGGGGCGGGGGCGGGGGGCAGGACGGGGGGAAGTGCCGCGGGGCG.CGGGGGGG.GCGCTTCTGCCCGCCCCGGAGCTTGCGTTTTCAGCCTTCGa                                                  <  2:2502270/100‑1 (MQ=255)
                                                 cAGGGGCCGCGGCGAGGACGGGGGTCAGGACGGGGGCAAGTGTCGCCGAGCT.CGGGGGGG.GCGCTTCTGCCCGCCCCGGAGCTTGCGTTTTCAGCCTTCGa                                                  >  2:287643/1‑101 (MQ=255)
                                                  gggggCCGCGGCGAGGACGGGGGTCAGGAGGGGGGCAAGTGTCGCCGAGCT.CGGGGGGG.GCGCTTCTGCCCGCCCCGGAGCTTGCGTTTTCAGCCTTCGAc                                                 <  2:530791/100‑1 (MQ=255)
                                                  gggcgCCGGGGCGAGGGGGGGGGTCGGGCCGGGGGAAAGTGTCGCCGAGCT.CGGGGGGG.GCGCTTCTGCCCGCCCCGGAGCTTGCGTTTTCAGCCTTCGAc                                                 <  2:3635721/97‑1 (MQ=255)
                                                  aGGGGCCGCGGCGAGGACGGGGGTCAGGACGGGGGCAAGTGTCGCCGAGCT.CGGGGGGG.GCGCTTCTGCCCGCCCCGGAGCTTGCGTTTTCAGCCTTCGAc                                                 <  2:1213157/101‑1 (MQ=255)
                                                  aGGGGCCGCGGCGAGGACGGGGGTCAGGACGGGGGCAAGTGTCGCCGAGCT.CGGGGGGG.GCGCTTCTGCCCGCCCCGGAGCTTGCGTTTTCAGCCTTCGAc                                                 >  2:1239560/1‑101 (MQ=255)
                                                  aGGGGCCGCGGCGAGGACGGGGGTCAGGACGGGGGCAAGTGTCGCCGAGCT.CGGGGGGG.GCGCTTCTGCCCGCCCCGGAGCTTGCGTTTTCAGCCTTCGAc                                                 >  1:3600507/1‑101 (MQ=255)
                                                  aGGGGCCGCGGCGAGGACGGGGGTCAGGACGGGGGCAAGTGTCGCCGAGCT.CGGGGGGG.GCGCTTCTGCCCGCCCCGGAGCTTGCGTTTTCAGCCTTCGAc                                                 >  2:3234528/1‑101 (MQ=255)
                                                   ggggCCGCGGCGAGGACGGGGGTCAGGACGGGGGCAAGTGTCGCCGAGCT.CGGGGGGG.GCGCTTCTGCCCGCCCCGGAGCTTGCGTTTTCAGCCTTCGAcg                                                >  1:3331472/1‑101 (MQ=255)
                                                   ggggCCGCGGCGAGGACGGGGGTCAGGACGGGGGCAAGTGTCGCCGAGCT.CGGGGGGG.GCGCTTCTGCCCGCCCCGGAGCTTGCGTTTTCAGCCTTCGAcg                                                >  1:3979211/1‑101 (MQ=255)
                                                       ccGCGGCGAGGACGGGGGTCAGGACGGGGGCAAGTGTCGCCGAGCT.CGGGGGGG.GCGCTTCTGCCCGCCCCGGAGCTTGCGTTTTCAGCCTTCGAcggcgg                                            >  2:1895434/1‑101 (MQ=255)
                                                       ccGCGGCGAGGACGGGGGTCAGGACGGGGGCAAGTGTCGCCGAGCT.CGGGGGGG.GCGCTTCTGCCCGCCCCGGAGCTTGCGTTTTCAGCCTTCGAcggcgg                                            >  1:2446120/1‑101 (MQ=255)
                                                           ggcgGGGGCGGGGGTGAGGAGGGGGGCAGGTGTCGCCGGGCTGG.GGGGGG.GCGCTTCTGCCCGCCCCGGAGCTTGCGTTTTCAGCCTTCGACggcggcggc                                        <  2:3407412/101‑1 (MQ=255)
                                                           ggcgAGGACGGGGGTCAGGACGGGGGCAAGTGTCGCCGAGCT.CGGGGGGG.GCGCTTCTGCCCGCCCCGGAGCTTGCGTTTTCAGCCTTCGACggcggcggc                                        <  2:3981832/101‑1 (MQ=255)
                                                           ggcgAAGACGGGGGTCAGGACGGGGGCAAGTGTCGCCCAGCT.CGGGGGGG.GCGCTTCTGCCCGCCCCGGAGCTTGCGTTTTCAGCCTTCGACggcggcggc                                        <  1:3905038/101‑1 (MQ=255)
                                                            gcgAGGACGGGGGTCAGGACGGGGGCAAGTGTCGCCGAGCT.CGGGGGGG.GCGCTTCTGCCCGCCCCGGAGCTTGCGGTTTCAGCCTTCGACGgcggcggcg                                       >  1:2566742/1‑101 (MQ=255)
                                                             cgAGGACGGGGGTCAGGCCGGGGGCAATTGTCGCCGAGCT.CGGGGGGG.GCGCTTCTGCCCGCCCCGGAGCTTGCGTTTTCAGCCTTCGAcggcggcgggtg                                      <  2:1442273/101‑4 (MQ=255)
                                                              ggggAGGGGGGTCAGGACGGGGGCAAGTGTCGCCGAGCT.CGGGGGGG.GCGCTTCTGCCCGCCCCGGAGCTTGCGTTTTCAGCCTTCGACGGCGGCGGCGaa                                     <  2:4293822/99‑1 (MQ=255)
                                                              gAGGACGGGGGTCAGGACGGGGGCAAGTGTCGCCGAGCT.CGGGGGGG.GCGCTTCTGCCCGCCCCGGAGCTTGCGTTTTCAGCCTTCGACGGCGGCGGCGaa                                     >  1:4112662/1‑101 (MQ=255)
                                                               aGGACGGGGGTCAGGACGGGGGCAAGTGTCGCCGAGCT.CGGGGGGG.GCGCTTCTGCCCGCCCCGGAGCTTGCGTTTTCAGCCTTCGACGGCGGCGGCGAAg                                    <  1:1400017/101‑1 (MQ=255)
                                                               aGGACGGGGGTCAGGACGGGGGCAAGTGTCGCCGAGCT.CGGGGGGG.GCGCTTCTGCCCGCCCCGGAGCTTGCGTTTTCAGCCTTCGACGGCGGCGGCGAAg                                    >  2:831615/1‑101 (MQ=255)
                                                               aGGACGGGGGTCAGGACGGGGGCAAGTGTCGCCGAGCT.CGGGGGGG.GCGCTTCTGCCCGCCCCGGAGCTTGCGTTTTCAGCCTTCGACGGCGGCGGCGAAg                                    >  1:2274683/1‑101 (MQ=255)
                                                                 gACGGGGGTCAGGACGGGGGCAAGTGTCGCCGAGCT.CGGGGGGG.GCGCTTCTGCCCGCCCCGGAGCTTGCGTTTTCAGCCTTCGACGGCGGCGGCGAAGgc                                  >  2:2184567/1‑101 (MQ=255)
                                                                 gACGGGGGTCAGGACGGGGGCAAGTGTCGCCGAGCT.CGGGGGGG.GCGCTTCTGCCCGCCCCGGAGCTTGCGTTTTCAGCCTTCGACGGCGGCGGCGAAGgc                                  >  2:1073115/1‑101 (MQ=255)
                                                                   cGGGGGTCAGGACGGGGGCAAGTGTCGCCGAGCT.CGGGGGGG.GCGCTTCTGCCCGCCCCGGAGCTTGCGTTTTCAGCCTTCGACGGCGGCGGCGAAGgcgg                                >  2:944234/1‑101 (MQ=255)
                                                                    gggggTCAGGACGGGGGCAAGTGTCGCCGAGCT.CGGGGGGG.GCGCTTCTGCCCGCCCCGGAGCTTGCGTTTTCAGCCTTCGACGGCGGCGGCGAAGgcggt                               >  1:2052771/1‑101 (MQ=255)
                                                                      ggggcgggACGGGGCCAAGCGCCCCCCACCC.CGGGGGGG.GCGCTTCTGCCCGCCCCGGAGCTTGCGTTTTCAGCCTTCGACGGCGGCGGCGAAGgcggtgc                             <  1:1781561/95‑1 (MQ=255)
                                                                      gggACAGGACGGGGGCAAGTGTCGCCGAGCT.CGGGGGGG.GCGCTTCTGCCCGCCCCGGAGCTTGCGCTTTCAGCCTTAGACGGCGGCGGCGAAGgcggtgc                             >  2:2792380/1‑101 (MQ=255)
                                                                       ggTCAGGACGGGGGCAAGTGTCGCCGAGCT.CGGGGGGG.GCGCTTCTGCCCGCCCCGGAGCTTGCGTTTTCAGCCTTCGACGGCGGCGGCGAAGgcggtgcg                            >  2:4072760/1‑101 (MQ=255)
                                                                          cggggggggggGAAGTGGCGCCGAGCTGG.GGGGGG.GCGCTTCTGCCCGCCCCGGAGCTTGCGTTTTCAGCCTTCGACGGCGGCGGCGAAGGCGGTGCGGTa                         <  2:940138/95‑1 (MQ=255)
                                                                           aGGACGGGGGCAAGTGTCGCCGAGCT.CGGGGGGG.GCGCTTCTGCCCGCCCCGGAGCTTGCGTTTTCAGCCTTCGACGGCGGCGGCGAAGGCGGTGCGGTaa                        >  2:4323437/1‑101 (MQ=255)
                                                                           aGGACGGGGCCAAGTGTCGCCGAGCT.CGGGGGGG.GCGCTTCTGCCCGCCCCGGAGCTTGCGTTTTCAGCCTTCGACGGCGGCGGCGAAGGCGGTGCGGTaa                        <  2:2671713/101‑1 (MQ=255)
                                                                            ggAGGGGGGCAAGTGGCGCCGAGGTGG.GGGGGG.GCGCTTCTGCCCGCCCCGGAGCTTGCGTTTTCAGCCTTCGACGGCGGCGGCGAAGGCGGTGCGGTaaa                       <  2:3216569/101‑1 (MQ=255)
                                                                            ggACGGGGGCAAGTGTCGCCGAGCT.CGGGGGGG.GCGCTTCTGCCCGCCCCGGAGCTTGCGTTTTCAGCCTTCGACGGCGGCGGCGAAGGCGGTGCGGTaaa                       <  2:394608/101‑1 (MQ=255)
                                                                              gggggggCAAGTGTCGCGGAGCTGG.GGGGGG.GCGCTTCTGCCCGCCCCGGAGCTTGCGTTTTCAGCCTTCGACGGCGGCGGCGAAGGCGGTGCGGTAAAGc                     <  1:544387/99‑1 (MQ=255)
                                                                               cGGGGGCAAGTGTCGCCGAGCT.CGGGGGGG.GCGCTTCTGCCCGCCCCGGAGCTTGCGTTTTCAGCCTTCGACGGCGGCGGCGAAGGGGGGGCGGTAAAGCg                    >  2:1220724/1‑101 (MQ=255)
                                                                                gggggCAAGTGTCGCCGAGCT.CGGGGGGG.GCGCTTCTGCCCGCCCCGGAGCTTGCGTTTTCAGCCTTCGACGGCGGCGGCGAAGGCGGTGCGTTAAAGCGt                   <  2:1669345/101‑1 (MQ=255)
                                                                                 ggggCAAGTGTCGCCGAGCT.CGGGGGGG.GCGCTTCTGCCCGCCCCGGAGCTTGCGTTTTCAGCCTTCGACGGCGGCGGCGAAGGCGGTGCGGTAAAGCGTg                  >  2:1536615/1‑101 (MQ=255)
                                                                                 ggggCAAGTGTCGCCGAGCT.CGGGGGGG.GCGCTTCTGCCCGCCCCGGAGCTTGCGTTTTCAGCCTTCGACGGCGGCGGCGAAGGCGGTGCGGTAAAGCGTg                  >  2:3136330/1‑101 (MQ=255)
                                                                                  gggCAAGTGTCGCCGAGCT.CGGGGGGG.GCGCTTCTGCCCGCCCCGGAGCTTGCGTTTTCAGCCTTCGACGGCGGCGGCGAAGGCGGTGCGGTAAAGCGTGc                 <  2:3331235/101‑1 (MQ=255)
                                                                                  gggCAAGTGTCGCCGAGCT.CGGGGGGG.GCGCTTCTGCCCGCCCCGGAGCTTCCGTTTTCAGCCTTCGACGGCGGCGGCGAAGGCGGTGCGGTAAACCGTGc                 >  1:4116673/1‑101 (MQ=255)
                                                                                    gCAAGTGTCGCCGAGCT.CGGGGGGG.GCGCTTCTGCCCGCCCCGGAGCTTGCGTTTTCAGCCTTCGACGGCGGCGGCGAAGGCGGTGCGGTAAAGCGTGCTg               <  1:1239653/101‑1 (MQ=255)
                                                                                     cAAGTGTCGCCGAGCT.CGGGGGGG.GCGCTTCTGCCCGCCCCGGAGCTTGCGTTTTCAGCCTTCGACGGCGGCGGCGAAGGCGGGGCGGTAAAGCGTGCTGa              >  1:3084385/1‑101 (MQ=255)
                                                                                      aaGTGTCGCCGAGCTGG.GGGGGG.GCGCTTCTGCCCGCCCCGGAGCTTGCGTTTTCAGCCTTCGACGGCGGCGGCGAAGGCGGTGCGGTAAAGCGTGCTGAg             <  2:4114118/101‑1 (MQ=255)
                                                                                      aaGTGTCGCCGAGCT.CGGGGGGG.GCGCTTCTGCCCGCCCCGGAGCTTGCGTTTTCAGCCTTCGACGGCGGCGGCGAAGGCGGTGCGGTAAAGCGTGCTGAg             >  1:1066607/1‑101 (MQ=255)
                                                                                      aaGTGTCGCCGAGCT.CGGGGGGG.GCGCTTCTGCCCGCCCCGGAGCTTGCGTTTTCAGCCTTCGACGGCGGCGGCGAAGGCGGTGCGGTAAAGCGTGCTGAg             >  1:2726473/1‑101 (MQ=255)
                                                                                      aaGTGTCGCCGAGCT.CGGGGGGG.GCGCTTCTGCCCGCCCCGGAGCTTGCGTTTTCAGCCTTCGACGGCGGCGGCGAAGGCGGTGCGGTAAAGCGTGCTGAg             <  2:138716/101‑1 (MQ=255)
                                                                                             gCCGAGCT.CGGGGGGG.GCGCTGCTGCCCGCCCCGGAGCTTGCGTTTTCAGCCTTCGACGGCGGCGGCGAAGGCGGTGCGGTAAAGCGTGCTGAGCTCTGcc      >  1:3292351/1‑101 (MQ=255)
                                                                                              ccGAGCT.CGGGGGGG.GCGCTTCTGCCCGCCCCGGAGCTTGCGTTTTCAGCCTTCGACGGCGGCGGCGAAGGCGGTGCGGTAAAGAGTGCTGAGCTCTGCCa     <  1:4022854/101‑1 (MQ=255)
                                                                                               ggggCT.C.GGGGGGCGCACTTCTGCCCGCCCCGGAGCTTGCGTTTTCAGCCTTCGACGGCGGCGGCGAAGGCGGTGCGGTAAAGCGTGCTGAGCTCTGCCAg    <  1:2962342/98‑1 (MQ=255)
                                                                                                gAGCT.CGGGGGGG.GCGCTTCTGCCCGCCCCGGAGCTTGCGTTTTCAGCCTTCGACGGCGGCGGCGAAGGCGGTGCGGTAAAGCGTGCTGAGCTCTGCCAgc   >  1:3879468/1‑101 (MQ=255)
                                                                                                gAGCT.CGGGGGGG.GCGCTTCTGCCCGCCCCGGAGCTTGCGTTTTCAGCCTTCGACGGCGGCGGCGAAGGCGGTGCGGTAAAGCGTGCTGAGCTCTGCCAgc   >  1:677200/1‑101 (MQ=255)
                                                                                                 aGCT.CGGGGGGG.GCGCTTCTGCCCGCCCCGGAGCTTGCGTTTTCAGCCTTCGACGGCGGCGGCGAAGGCGGTGCGGTAAAGCGTGCTGAGCTCTGCCAgcg  >  2:583757/1‑101 (MQ=255)
                                                                                                     |                                                                                                  
AGGCGACAAGAGCCAAGCACCGGTTGCGTGTGCCAAGTGACGCAGCGGGCAGGGGCCGCGGCGAGGACGGGGGTCAGGACGGGGGCAAGTGTCGCCGAGCT.C.GGGGGG.GCGCTTCTGCCCGCCCCGGAGCTTGCGTTTTCAGCCTTCGACGGCGGCGGCGAAGGCGGTGCGGTAAAGCGTGCTGAGCTCTGCCAGCG  >  NC_014034/586816‑587012

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