Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_014034 3,717,668 +G 1.8% coding (268/699 nt) hemD → uroporphyrinogen‑III synthase

Read alignment evidence...
  seq id position change freq score reads annotation genes product
*NC_0140343,717,6681.→G 1.8% 4.2 ‑4.3 ‑0.1227coding (268/699 nt)hemDuroporphyrinogen‑III synthase
Reads supporting (aligned to +/- strand):  new base (4/0):  ref base (74/149):  total (78/149)
Fisher's exact test strand distribution p-value = 1.32e-02
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 8.76e-01

Only 100 of 303 total aligned reads displayed.

CATCTTCACCTCGGAACAGGCGGTCAAACCGCTGGCCGAAAGCGCCCTTTCGGGCCGTCCGGCCTTTTGCGTCGGCGAGCGGACGGCGGCGGTGGTGCGGC.GGGC.GGGTTTCCCGGTGCT.GGCCGTCGCCCCCGATGCCGAGCATCTGTTCGAAGTGATCCTGTCCGCCGCCGATCCGGGGCCGATCCTGCATGCCCGCGG  >  NC_014034/3717568‑3717768
                                                                                                     |                                                                                                      
cATCTTCACCTCGGAACAGGCGGTCAAACCGCTGGCCGAAAGCGCCCTTTCGGGCCGTCCGGCCTTTTGCGTCGGCGAGCGGACGGCGGCGGTGGTgcggc                                                                                                         <  2:7079540/101‑1 (MQ=255)
  tCTTCACCTCGGAACAGGCGGTCAAACCGCTGGCCGAAAGCGCCCTTTCGGGCCGTCCGGCCTTTTGCGTCGGCGAGCGGACGGCGGCGGTGGTGcggc.gg                                                                                                      >  1:6752935/1‑101 (MQ=255)
   cTTCACCTCGGAACAGGCGGTCAAACCGCTGGCCGAAAGCGCCCTTTCGGGCCGTCCGGCCTTTTGCGTCGGCGAGCGGACGGCGGCGGTGGTGCggc.ggg                                                                                                     >  1:3685729/1‑101 (MQ=255)
    tgcACCTCGGAACAGGCGGTCAAACCGCTGGCCGAAAGCGCCCTTTCGGGCCGTCCGGCCTTTTGCGTCGGCGAGCGGACGGCGGCGGTGGTGCggc.gggc                                                                                                    <  2:6177767/99‑1 (MQ=255)
     tCACCTCGGAACAGGCGGTCAAACCGCTGGCCGAAAGCGCCCTTTCGGGCCGTCCGGCCTTTTGCGTCGGCGAGCGGACGGCGGCGGTGGTGCggc.ggggg                                                                                                   >  2:20315/1‑99 (MQ=255)
     tCACCTCGGAACAGGCGGTCAAACCGCTGGCCGAAAGCGCCCTTTCGGGCCGTCCGGCCTTTTGCGTCGGCGAGCGGACGGCGGCGGTGGTGCggc.gggc.g                                                                                                  <  2:1217714/101‑1 (MQ=255)
       aCCTCGGAACAGGCGGTCAAACCGCTGGCCGAAAGCGCCCTTTCGGGCCGTCCGGCCTTTTGCGTCGGCGAGCGGACGGCGGCGGTGGTGCGGc.gggg.ggg                                                                                                >  2:6254827/1‑101 (MQ=255)
         cTCGGAACAGGCGGTCAAACCGCTGGCCGAAAGCGCCCTTTCGGGCCGTCCGGCCTTTTGCGTCGGCGAGCGGACGGCGGCGGTGGTGCGGc.gggg.ggggg                                                                                              >  2:177220/1‑99 (MQ=255)
         cTCGGAACAGGCGGTCAAACCGCTGGCCGAAAGCGCCCTTTCGGGCCGTCCGGCCTTTTGCGTCGGCGAGCGGACGGCGGCGGGGGGGCGGg.gggg.ggggg                                                                                              >  1:3817982/1‑99 (MQ=255)
          tCGGAACAGGCGGTCAAACCGCTGGCCGAAAGCGCCCTTTCGGGCCGTCCGGCCTTTTGCGTCGGCGAGCGGACGGCGGCGGTGGTGCGGC.GGGC.GGGttt                                                                                             <  1:5886246/101‑1 (MQ=255)
          tCGGAACAGGCGGTCAAACCGCTGGCCGAAAGCGCCCTTTCGGGCCGTCCGGCCTTTTGCGTCGGCGAGCGGACGGCGGCGGTGGTGCGGC.GGGC.GGGttt                                                                                             <  1:6113791/101‑1 (MQ=255)
          tCGGAACAGGCGGTCAAACCGCTGGCCGAAAGCGCCCCTTCGGGCCGTCCGGCCTTTTGCGTCGGCGAGCGGACGGCGGCGGTGGTGCGGC.GGGC.GGGttt                                                                                             <  1:1360243/101‑1 (MQ=255)
            ggAACAGGCGGTCAAACCGCTGGCCGAAAGCGCCCTTTCGGGCCGTCCGGCCTTTTGCGTCGGCGAGCGGACGGCGGCGGTGGTGCGGC.GGGC.GGGTTTcc                                                                                           <  2:1424255/101‑1 (MQ=255)
            ggAACAGGCGGTCAAACCGCTGGCCGAAAGCGCCCTTTCGGGCCGTCCGGCCTTTTGCGTCGGCGAGCGGACGGCGGCGGGGGGGCGGc.gggg.ggggggtc                                                                                           >  2:2302502/1‑96 (MQ=255)
             gAACAGGCGGTCAAACCGCTGGCCGAAAGCGCCCTTTCGGGCCGTCCGGCCTTTTGCGTCGGCGAGCGGACGGCGGCGGTGGTGCGGC.GGGC.GGGGGTccc                                                                                          >  2:6746230/1‑101 (MQ=255)
              aaCAGGCGGTCAAACCGCTGGCCGAAAGCGCCCTTTCGGGCCGTCCGGCCTTTTGCGTCGGCGAGCGGACGGCGGCGGTGGTGCGGC.GGGC.GGGGTTctcc                                                                                         >  1:5536627/1‑98 (MQ=255)
              aaCAGGCGGTCAAACCGCTGGCCGAAAGCGCCCTTTCGGGCCGTCCGGCCTTTTGCGTCGGCGAGCGGACGGCGGCGGGGGTGCGGC.GGGG.GGGGTTTCCg                                                                                         >  2:5107199/1‑101 (MQ=255)
               aCAGGCGGTCAAACCGCTGGCCGAAAGCGCCCTTTGGGGCCGTCCGGCCTTTTGCGTCGGCGAGCGGACGGCGGCGGTGGTGCGGc.gggg.gggggtgcccg                                                                                        >  1:6412160/1‑93 (MQ=255)
                cAGGCNGTCAAACCGCTGGCCGAAAGCGCCCTTTCGGGCCGTCCGGCCTTTTGCGTCGGCGAGCGGACGGCGGCGGTGGTGCGGC.GGGC.GGGTTTCCCGGt                                                                                       <  1:2010590/101‑1 (MQ=255)
                 aGGCGGTCAAACCGCTGGCCGAAAGCGCCCTTTCGGGCCGTCCGGCCTTTTGCGTCGGCGAGCGGACGGCGGCGGTGGTGCGGCGGGGG.GGGTTTCCCgtg                                                                                       >  2:1623566/1‑99 (MQ=255)
                  ggCGGTCAAACCGCTGGCCGAAAGCGCCCTTTCGGGCCGTCCGGCCTTTTGCGTCGGCGAGCGGACGGCGGCGGTGGTGCGGC.GGGC.GGGTTTCCCGGTGc                                                                                     <  2:4205948/101‑1 (MQ=255)
                    cGGTCAAACCGCTGGCCGAAAGCGCCCTTTCGGGCCGTCCGGCCTTTTGCGTCGGCGAGCGGACGGCGGCGGTGGTGCGGCGGGGG.GGGGTTCCCgggggg                                                                                    >  1:5811693/1‑97 (MQ=255)
                    cGGTCAAACCGCTGGCCGAAAGCGCCCTTTCGGGCCGTCCGGCCTTTTGCGTCGGCGAGCGGACGGCGGCGGTGGTGCGGC.GGGC.GGGTTTCCCGGTGCT.g                                                                                  <  1:3624123/101‑1 (MQ=255)
                     ggTCAAACCGCTGGCCGAAAGCGCCCTTTCGGGCCGTCCGGCCTTTTGCGTCGGCGAGCGGACGGCGGCGGTGGTGCGGC.GGGG.GGGGGTCCCCGGGCT.gg                                                                                 >  1:5121053/1‑101 (MQ=255)
                         aaaCCGCTGGCCGAAAGCGCCCTTTCGGGCCGTCCGGCCTTTTGCGTCGGCGAGCGGACGGCGGCGGTGGTGCGGC.GGGG.GGGGGTTCCGGGGCT.GGGCgg                                                                             >  1:7115695/1‑100 (MQ=255)
                               cTGGCCGAAAGCGCCCTTTCGGGCCGTCCGGCCTTTTGCGTCGGCGAGCGGACGGCGGCGGTGGTGCGGC.GGGG.GGGGTTCCCCGGGGT.GGGCGGCGcccc                                                                       >  2:3375477/1‑101 (MQ=255)
                               cTGGCCGAAAGCGCCCTTTCGGGCCGTCCGGCCTTTTGCGTCGGCGAGCGGACGGCGGCGGTGGTGCGGC.GGGCGGGGTTTCCCGG.GCT.GGGCGGGGcccc                                                                       >  1:4946578/1‑101 (MQ=255)
                                tGGCCGAAAGCGCCCTTTCGGGCCGTCCGGCCTTTTGCGTCGGCGAGCGGACGGCGGCGGTGGTGCGGC.GGGG.GGGGGTTCCGGGGGT.GGGCGTCGcccga                                                                      >  2:2514236/1‑99 (MQ=255)
                                   ccGAAAGCGCCCTTTCGGGCCGTCCGGCCTTTTGCGTCGGCGAGCGGACGGCGGCGGTGGTGCGGC.GGGC.GGGTTTCCCTGTGCT.GGCCGTCGCCCCCGAt                                                                   <  2:2108719/101‑1 (MQ=255)
                                    cGAAAGCGCCCTTTCGGGCCGTCCGGCCTTTTGCGTCGGCGAGCGGACGGCGGCGGTGGTGCGGC.GGGG.GGGGTTTGCCGTGGT.GGGCGGTGGCCCCGagg                                                                  >  2:428364/1‑99 (MQ=255)
                                     gAAAGCGCCCTTTCGGGCCGTCCGGCCTTTTGCGTCGGCGAGCGGACGGCGGCGGTGGTGCGGC.GGGC.GGGGGTCCCGGGGCC.CGCCGTCGcccccaccgc                                                                 >  2:6981931/1‑96 (MQ=255)
                                      aaaGCGCCCTTTCGGGCCGTCCGGCCTTTTGCGTCGGCGAGCGGACGGCGGCGGTGGTGCGGC.GGGG.GGGGTTTCCGGGGGT.GGGCGTCGCCCCCCAGGcc                                                                >  1:3450504/1‑101 (MQ=255)
                                      aaaGCGCCCTTTCGGGCCGTCCGGCCTTTTGCGTCGGCGAGCGGACGGCGGCGGTGGTGCGGC.GGGC.GGGGTTCCCCGGGCC.GGCCGGCGCCCCCGATGcc                                                                >  1:4549342/1‑101 (MQ=255)
                                       aaGCGCCCTTTCGGGCCGTCCGGCCTTTTGCGTCGGCGAGCGGACGGCGGCGGTGGTGCGGCGGGGG.GGGTTTCCCGG.GCG.GGGCGGCGCCCCCCCTGccc                                                               >  1:5629574/1‑100 (MQ=255)
                                       aaGCGCCCTTTCGGGCCGTCCGGCCTTTTGCGTCGGCGAGCGGACGGCGGCGGTGGTGCGGC.GGGC.GGGTTTCCCGGTGCT.GGCCGTCGCCCCCGATGCCg                                                               <  2:1895396/101‑1 (MQ=255)
                                        aGCGCCCTTTCGGGCCGTCCGGCCTTTTGCGTCGGCGAGCGGACGGCGGCGGTGGTGCGGC.GGGC.GGGTTTCCCGGTGCT.GGCCGTCGCCCCCGATGTCGa                                                              <  1:2934972/101‑1 (MQ=255)
                                        aGCGCCCTTTCGGGCCGTCCGGCCTTTTGCGTCGGCGAGCGGACGGCGGCGGTGGTGCGGC.GGGC.GGGTTTCCCGGTGCT.GGCCGTCGCCCCCGATGCCGa                                                              <  1:2302597/101‑1 (MQ=255)
                                        aGCGCCCTTTCGGGCCGCCCGCCCTTTTGCGCCGGCGAGCGGCCGGCGGCGGTGGGGCGGC.GGGG.GGGGTTCCCGGGGCT.GGCCGCCCCCCCCCCTGCCgg                                                              >  2:1000555/1‑100 (MQ=255)
                                          cgcCCTTTCGGGCCGTCCGGCCTTTTGCGTCGGCGAGCGGACGGCGGCGGTGGTGCGGC.GGGC.GGGTTTCCCGGTGCT.GGCCGTCGCCCCCGATGCCGAGc                                                            <  1:1755165/101‑1 (MQ=255)
                                          cgcCCTTTCGGGCCGTCCGGCCTTTTGCGTCGGCGAGCGGACGGCGGCGGTGGTGCGGC.GGGC.GGGTTTCCCGGTGCT.GGCCGTCGCCCCCGATGCCGAGc                                                            <  2:491682/101‑1 (MQ=255)
                                           gcCCTTTCGGGCCGTCCGGCCTTTTGCGTCGGCGAGCGGACGGCGGCGGTGGTGCGGC.GGGG.GGGGTTTCCGGGGCT.GGGCGGCGCCCCCCCTGcccggca                                                           >  2:1492783/1‑96 (MQ=255)
                                             ccTTTCGGGCCGTCCGGCCTTTTGCGTCGGCGAGCGGACGGCGGCGGTGGTGCGGC.GGGC.GGGGTTCCCCGGGGT.GGGCGTCGCCCCCGCTGccccgcacc                                                         >  1:828646/1‑94 (MQ=255)
                                              cTTTCGGGCCGTCCGGCCTTTTGCGTCGGCGAGCGGACGGCGGCGGTGGTGCGGC.GGGG.GGGGTTTCCCGTGCT.GGCCGTCGCCCCCCCTGCCggggctct                                                        >  1:2838957/1‑94 (MQ=255)
                                               tttCGGGCCGTCCGGCCTTTTGCGTCGGCGAGCGGACGGCGGCGGTGGTGCGGC.GGGG.GGGGTTTCCGGGGGT.GGGCGGCGGCCCCGCTGCCCAggctccg                                                       >  2:6908510/1‑95 (MQ=255)
                                                 tCGGGCCGTCCGGCCTTTTGCGTCGGCGAGCGGACGGCGGCGGTGGTGCGGC.GGGC.GGGGTTCCCGGGGCT.GGCCGGCGCCCCCCCGGCCGGGCCTCTGtt                                                     >  1:3880493/1‑101 (MQ=255)
                                                  cGGGCCGTCCGGCCTTTTGCGTCGGCGAGCGGACGGCGGCGGTGGTGCGGC.GGGC.GGGTTTCCCGGTGCT.GGCCGTCGCCCCCGATGCCGAGCATCTGTTc                                                    <  2:270949/101‑1 (MQ=255)
                                                  cGGGCCGTCCGGCCTTTTGCGTCGGCGAGCGGACGGCGGCGGTGGTGCGGC.GGGC.GGGTTTCCCGGTGCT.GGCCGTCGCCCCCGATGCCGAGCATCTGTTc                                                    <  2:3817819/101‑1 (MQ=255)
                                                   gggCCGTCCGGCCTTTTGCGTCGGCGAGCGGACGGCGGCGGTGGTGCGGC.GGGG.GGGGTTTCCGGGGGT.GGGCGTCGCCCCCGATGCCGCGCATCTGTTCg                                                   >  1:4680507/1‑101 (MQ=255)
                                                   gggCCGTCCGGCCTTTTGCGTCGGCGAGCGGACGGCGGCGGTGGTGCGGC.GGGC.GGGGTTTCCGGGGCT.GGGCGGCGCCCCCGATGCCGAGCATCTGTTCg                                                   >  1:4518781/1‑101 (MQ=255)
                                                      ccGTCCGGCCTTTTGCGTCGGCGAGCGGACGGCGGCGGTGGTGCGGC.GGGC.GGGTTTCCCGGTGCT.GGCCGTCGCCCCCGATGCCGAGCATCTGTTCGAAg                                                <  2:3048617/101‑1 (MQ=255)
                                                        gTCCGGCCTTTTGCGTCGGCGAGCGGACGGCGGCGGTGGTGCGGC.GGGG.GGGGTTTCCGGGGGT.GGCCGGCGCCCCCGATGCCGGGCCTCTGTTCGGAggg                                              >  1:2123340/1‑99 (MQ=255)
                                                        gTCCGGCCTTTTGCGTCGGCGAGCGGACGGCGGCGGTGGTGCGGC.GGGC.GGGTTTCCCGGTGCT.GGCCGTCGCCCCCGATGCCGAGCATCTGTTCGAAGTg                                              <  1:5238195/101‑1 (MQ=255)
                                                          ccGGCCTTTTGCGTCGGCGAGCGGACGGCGGCGGTGGTGCGGC.GGGC.GGGTTTCCCGGTGCT.GGCCGTCGCCCCCGATGCCGAGCATCTGTTCGAAGTGAt                                            <  1:2530227/101‑1 (MQ=255)
                                                          ccGGCCTTTTGCGTCGGCGAGCGGACGGCGGCGGTGGTGCGGC.GGGC.GGGGTTTCCGGGGCT.GGGCGGCGCCCCCGCTGCCGAGCATCTGTTTGAAGGGAt                                            >  1:62697/1‑101 (MQ=255)
                                                           cGGCCTTTTGCGTCGGCGAGCGGACGGCGGCGGTGGTGCGGC.GGGC.GGGTTTCCCGGTGCT.GGCCGTCGCCCCCGATGCCGAGCATCTGTTCGAAGTGATc                                           <  1:5733690/101‑1 (MQ=255)
                                                            ggCCTTTTGCGTCGGCGAGCGGACGGCGGCGGTGGTGCGGC.GGGG.GGGGTTCCCCGTGCT.GGGCGGCGCCCCCCCTGCCGGGCACCTGTTCGTAGGGAttc                                          >  1:2507315/1‑99 (MQ=255)
                                                            ggCCTTTTGCGTCGGCGAGCGGACGGCGGCGGTGGTGCGGC.GGGC.GGGTTTCCCGGTGCT.GGCCGTCGCCCCCGATGCCGAGCATCTGTTCGAAGTGATcc                                          <  2:3944274/101‑1 (MQ=255)
                                                              ccgatctcGTCGGCGAGCGGACGGCGGCGGTGGTGCGGC.GGGC.GGGTTTCCCGGTGCT.GGCCGTCGCCCCCGATGCCGAGCATCTGTTCGAAGTGATCCTg                                        <  2:1114701/94‑1 (MQ=255)
                                                               cTTTTGCGTCGGCGAGCGGACGGCGGCGGTGGTGCGGCGGGGG.GGGTTTCCCGG.GCT.GGCCGGCGCCCCCGCTGCCGGGCCTCTGTTCGAAGTGATCCTgg                                       >  2:6670478/1‑100 (MQ=255)
                                                               cTTTTGCGTCGGCGAGCGGACGGCGGCGGTGGTGCGGC.GGGC.GGGTTTCCCGGTGCT.GGCCGTCGCCCCCGATGCCGAGCATCTGTTCGAAGTGATCCTGt                                       <  1:1826243/101‑1 (MQ=255)
                                                               cTTTTGCGTCGGCGAGCGGACGGCGGCGGTGGTGCGGC.GGGC.GGGTTTCCCGGTGCT.GGCCGTCGCCCCCGATGCCGAGCATCTGTTCGAAGTGATCCTGt                                       <  1:6529671/101‑1 (MQ=255)
                                                                ttttGCGTCGGCGAGCGGACGGCGGCGGTGGTGCGGC.GGGC.GGGCTTCCCCGGGGT.GGCCTTGGCCCCCCATGCCGAGGCTCTGCTCGAAGTGGTTCTGTc                                      >  2:4746156/1‑101 (MQ=255)
                                                                 tttGCGTCGGCGAGCGGACGGCGGCGGTGGTGCGGC.GGGC.GGGTTTCCCGGTGCT.GGCCGTCGCCCCCGATGCCGAGCATCTGTTCGAAGTGATCCTGTcc                                     <  1:1000590/101‑1 (MQ=255)
                                                                 tttGCGTCGGCGAGCGGACGGCGGCGGTGGTGCGGC.GGGC.GGGTTTCCCGGTGCT.GGCCGTCGCCCCCGATGCCGAGCATCTGTTCGAAGTGATCCTGTcc                                     <  1:1623627/101‑1 (MQ=255)
                                                                  ttGCGTCGGCGAGCGGACGGCGGCGGTGGTGCGGC.GGGC.GGGTTTCCCGGTGCT.GGCCGTCGCCCCCGATGCCGAGCATCTGTTCGAAGTGATCCTGTccg                                    <  2:214206/101‑1 (MQ=255)
                                                                  ttGCGTCGGCGAGCGGACGGCGGCGGTGGTGCGGC.GGGC.GGGGTTCCCGGGGCT.GGGCGGCGCCCCCGATGCCGAGCATCTGTTCGGAGTGATCCTGTccg                                    >  2:5959019/1‑101 (MQ=255)
                                                                    gCGTCGGCGAGCGGACGGCGGCGGTGGTGCGGC.GGGC.GGGTTTCCCGGTGCT.GGCCGTCGCCCCCGATGCCGAGCATCTGTTCGAAGTGATCCTGTccgcc                                  <  2:1202434/101‑1 (MQ=255)
                                                                     cGTCGGCGAGCGGACGGCGGCGGTGGTGCGGC.GGGC.GGGGTTCCCGGGGCT.GGCCGGCGCCCCCCCTGCCCAGCCTCTGTTCGGAGTGGTcccgaggccgg                                 >  1:1114736/1‑92 (MQ=255)
                                                                      gTCGGCGAGCGGACGGCGGCGGTGGTGCGGC.GGGC.GGGGTTTCCGGTGCT.GGCCGGCGCCCCCCATGCCCAGCATCTGTTCGAAGTGATCCTGTccgccgc                                >  1:4581102/1‑101 (MQ=255)
                                                                        cGGGGAGCGGACGGCGGCGGTGGTGCGGC.GGGC.GGGTTTCCCGGTGCT.GGCCGTCGCCCCCGATGCCGAGCATCTGTTCGAAGTGATCCTGTccgccgccg                              <  2:636078/101‑1 (MQ=255)
                                                                        cGGCGAGCGGACGGCGGCGGTGGTGCGGC.GGGC.GGGTTTCCCGGTGCT.GGCCGTCGCCCCCGATGCCGAGCATCTGTTCGAAGTGATCCTGTccgccgccg                              <  2:5165058/101‑1 (MQ=255)
                                                                        cGGCGAGCGGACGGCGGCGGTGGTGCGGC.GGGC.GGGGTTTCCGGTGCT.GGCCGTCGCCCCCCCTGCCGGGCCTCTGTTCGAAGGGATCCCGTccgcccccc                              >  1:403463/1‑100 (MQ=255)
                                                                         ggCGAGCGGACGGCGGCGGTGGTGCGGC.GGGC.GGGGTTCCCGGGGCT.GGGCGTCGCCCCCCATGCCGGGCATCCGTTCGAAGGGATCCCGTccccccccca                             >  1:6818837/1‑99 (MQ=255)
                                                                          gCGAGCGGACGGCGGCGGTGGTGCGGC.GGGC.GGGTTTCCCGGTGCT.GGCCGTCGCCCCCGATGCCGAGCATCTGTTCGAAGTGATCCTGTCCGCCGCCGAt                            <  2:6467358/101‑1 (MQ=255)
                                                                            gAGCGGACGGCGGCGGTGGTGCGGC.GGGC.GGGGTTCCCGGGGCT.GGCCGGCGCCCCCCATGCCGAGCATCTGTTCGAAGTGATCCTGTCCGCCGCCGATcc                          >  1:768648/1‑101 (MQ=255)
                                                                            gAGCGGACGGCGGCGGTGGTGCGGC.GGGC.GGGGTTCCCGGGGCG.GGCCGTCGCCCCCGGTGCCGAGGATCCGGTCGGAGTGAGCCCGTCCGCCGCCGGTcc                          >  1:2055307/1‑101 (MQ=255)
                                                                                ggACGGCGGCGGTGGTGCGGC.GGGC.GGGTTTCCCGGTGCT.GGCCGTCGCCCCCGATGCCGAGCATCTGTTCGAAGTGATCCTGTCCGCCGCCGATCCgggg                      <  2:5830776/101‑1 (MQ=255)
                                                                                   cggcggcggTGGTGCGGC.GGGC.GGGGTTCCCGGGGCT.GGCCGTCGCCCCCGATGCCGAGCATCTGTTCGAAGTGATCCTGGCCGCCGCCGGTCCCGGGCCg                   >  1:4373526/1‑101 (MQ=255)
                                                                                    ggcggcggTGGTGCGGC.GGGC.GGGGTTCCCGGGGCT.GGCCGTCGCCCCCGATGCCGAGCATCTGTTCGAAGTGATCCTGTCCGCCGCCGATCCGGGGCCGa                  >  1:1428131/1‑101 (MQ=255)
                                                                                       ggcggTGGTGCGGC.GGGC.GGGTTTCCCGGTGCT.GGCCGTCGCCCCCGATGCCGAGCATCTGTTCGAAGTGATCCTGTCCGCCGCCGATCCGGGGCCGATcc               <  2:2221538/101‑1 (MQ=255)
                                                                                       ggcggTGGTGCGGC.GGGC.GGGTTTCCCGGTGCT.GGCCGTCGCCCCCGATGCCGAGCATCTGTTCGAAGTGATCCTGTCCGCCGCCGATCCGGGGCCGATcc               <  1:5350478/101‑1 (MQ=255)
                                                                                       ggcggTGGTGCGGC.GGGC.GGGTTTCCCGGGGCT.GGGCGTCGCCCCCGATGCCGAGCATCTGTTCGAAGTGATCCTGTCCGCCGCCGATCCCGGGCCGAGcc               >  2:5365618/1‑101 (MQ=255)
                                                                                        gcggTGGTGCGGC.GGGC.GGGTTTCCCGGTGCT.GGCCGTCGCCCCCGATGCCGAGCATCTGTTCGAAGTGATCCTGTCCGCCGCCGATCTGGGGCCGATCCt              <  1:5097270/101‑1 (MQ=255)
                                                                                        gcggTGGTGCGGC.GGGC.GGGTTTCCCGGTGCT.GGCCGTCGCCCCCGATGCCGAGCATCTGTTCGAAGTGATCCTGTCCGCCGCCGATCCGGGGCCGATCCt              <  2:3582009/101‑1 (MQ=255)
                                                                                        gcggTGGTGCGGC.GGGC.GGGTTTCCCGGTGCT.GGCCGTCGCCCCCGATGCCGAGCATCTGTTCGAAGTGATCCTGTCCGCCGCCGATCCGGGGCCGATCCt              <  1:5823403/101‑1 (MQ=255)
                                                                                         cggTGGTGCGGC.GGGC.GGGGTTCCCGGGGCTGGGCCGCGGCCCCCG.TGCCGAGCATCTGTTCGAAGGGATCCTGTCCCCCGCCGCTCCGGGGCCGATCCtc             >  2:5712615/1‑100 (MQ=255)
                                                                                         cggTGGTGCGGC.GGGC.GGGGTTCCCGGGGCT.GGCCGTCGCCCCCGATGCCGAGCATCTGTTCGAAGTGATCCTGTCCGCCGCCGATCCGGGGCCGATCCtc             >  2:3752979/1‑100 (MQ=255)
                                                                                          ggtggtGCGGC.GGGC.GGGGTTCCCGGTGCT.GGCCGTCGCCCCCGATGCCGAGCATCTGTTCGAAGTGATCCTGTCCGCCGCCGATCCGGGGCCGATCCTGc            >  2:3467990/1‑101 (MQ=255)
                                                                                            tggcgCGGC.GGGC.GGGTTTCCCGGTGCT.GGCCGTCGCCCCCGATTCCGAGCATCTGTTCGAAGTGATCCTGTCCGCCGCCGATCCGGGGCCGATCCTGCAt          <  2:5275660/101‑1 (MQ=255)
                                                                                             ggtgCGGC.GGGC.GGGTTTTCCCGTGCC.TGCCGCCGCCCCCCATGCCGAGCGTCTGTTTGAAGTGCTCCTGTTCGCCGCCGACCCCGCGGCCATCCTGCacc         >  2:414722/1‑99 (MQ=255)
                                                                                             ggtgCGGC.GGGC.GGGGTTCCCGGGGCT.GGGCGTCGCCCCCGCTGCCGAGCATCTGTTCGAAGTGATCCTGTCCGCCGCCGATCCCGGGGCGATCCTGCATg         >  2:5531941/1‑101 (MQ=255)
                                                                                               tgCGGC.GGGC.GGGTTTCCCGGTCCT.GGCCGTCGCCCCCGATGCCGAGCATCTGTTCGAAGTGATCCTGTCCGCCGCCGATCCGGGGCCGATCCTGCATGcc       <  1:3294304/101‑1 (MQ=255)
                                                                                               gtgcgc.gggc.ggGTTTCCCGGTGCT.GGCCGTCGCCTCCGATGCCGAGCATCTGTTCGAAGTGATCCTGTCCGCCGCCGATCCGGGGCCGATCCTGCATGcc       <  2:5309767/97‑1 (MQ=255)
                                                                                                gcggc.gGGC.GGGTTTCCCGGTGCT.GGCCGTCGCCCTCGATGCCGAGCATCTGTTCGAAGTGATCCTGTCCGCCGCCGATCCGGGGCCGATCCTGCATGccc      <  1:3375600/101‑1 (MQ=255)
                                                                                                gcggc.gGGC.GGGTTTCCCGGTGCT.GGCCGTCGCCCCCGATGCCGAGCATCTGTTCGAAGTGATCCTGTCCGCCGCCGATCCGGGGCCGATCCTGCATGccc      >  2:6274628/1‑101 (MQ=255)
                                                                                                gcggc.gGGC.GGGTTTCCCGGTGCT.GGCCGTCGCCCCCGATGCCGAGCATCTGTTCGAAGTGATCCTGTCCGCCGCCGATCCGGGGCCGATCCTGCATGccc      >  1:4318/1‑101 (MQ=255)
                                                                                                 cggc.ggGC.GGGTTTCTCGGTGCT.GGCCGTCGCCCCCGATGCCGAGCATCTGTTCGAAGTGATCCTGTCCGCCGCCGATCCGGGGCCGATCCTGCATGCCcg     >  1:6854021/1‑101 (MQ=255)
                                                                                                  ggc.gggc.gGGTTTCCCGGTGCT.GGCCGTCGCCCCCGATGCCGAGCATCTGTTCGAAGTGATCCTGTCCGCCGCCGATCCGGGGCCGATCCTGCATGCCcgc    >  1:4197699/1‑101 (MQ=255)
                                                                                                  ggc.gggc.gGGTTTCCCGGTGCT.GGCCGTCGCCCCCGATGCCGAGCATCTGTTCGAAGTGATCCTGTCCGCCGCCGATCCGGGGCCGATCCTGCATGCCcgc    <  2:4580863/101‑1 (MQ=255)
                                                                                                    c.gggc.gggTTTCCCGGTGCT.GGCCGTCGCCCCCGATGCCGAGCATCTGTTCGAAGTGATCCTGTCCGCCGCCGATCCGGGGCCGATCCTGCATGCCCGCgg  >  2:780120/1‑101 (MQ=255)
                                                                                                     |                                                                                                      
CATCTTCACCTCGGAACAGGCGGTCAAACCGCTGGCCGAAAGCGCCCTTTCGGGCCGTCCGGCCTTTTGCGTCGGCGAGCGGACGGCGGCGGTGGTGCGGC.GGGC.GGGTTTCCCGGTGCT.GGCCGTCGCCCCCGATGCCGAGCATCTGTTCGAAGTGATCCTGTCCGCCGCCGATCCGGGGCCGATCCTGCATGCCCGCGG  >  NC_014034/3717568‑3717768

Base quality scores:  ATCG  < 3 ≤  ATCG  < 18 ≤  ATCG  < 26 ≤  ATCG  < 33 ≤  ATCG  < 40 ≤  ATCG