Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_014034 2,566,536 +G 97.9% intergenic (+8/+186) purD → / ← sdaA phosphoribosylamine‑‑glycine ligase/L‑serine ammonia‑lyase

Read alignment evidence...
  seq id position change freq score reads annotation genes product
*NC_0140342,566,5361.→G97.9% 5.4 ‑1.2 ‑1.5238intergenic (+8/+186)purD/sdaAphosphoribosylamine‑‑glycine ligase/L‑serine ammonia‑lyase
Reads supporting (aligned to +/- strand):  new base (49/184):  ref base (2/3):  total (51/187)
Fisher's exact test strand distribution p-value = 2.92e-01
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 3.42e-02

Only 100 of 298 total aligned reads displayed.

GCCGAGGCGCAACAGCGCGCCTATGCGATGGTCGATGCGATCGACTGGCCCGAGGGGTTCTGCCGCCGCGACATCGGCTGGCGCGCGCTCTGAGACG.AAA.C.GG.GG.G.CGCTGCCCCCCTCC.T..G..CTTCGGCCA.ATTCACCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAA  >  NC_014034/2566436‑2566629
                                                                                                       |                                                                                                      
gCCGAGGCGCAACAGCGCGCCTATGCGATGGTCGATGCGATCGACTGGCCCGAGGGGTTCTGCCGCCGCGACATCGGCTGGCGCGCGCTCTGAGACG.AAA.c                                                                                                         <  1:1240990/101‑1 (MQ=255)
        gcAACAGCGCGCCTATGCGATGGTCGATGCGATCGACTGGCCCGAGGGGTTCTGCCGCCGCGACATCGGCTGGCGCGCGCCCTGAGACG.AAA.C.ggccagga                                                                                                >  1:4513506/1‑95 (MQ=255)
          aaCAGCGCGCCTATGCGATGGTCGATGCGATCGACTGGCCCGAGGGGTTCTGCCGCCGCGACATCGGCTGGCGCGCGCTCTGCGACG.AAA.ccgaggggggg                                                                                               >  2:1885107/1‑91 (MQ=255)
          aaCAGCGCGCCTATGCGATGGTCGATGCGATCGACTGGCCCGAGGGGTTCTGCCGCCGCGACATCGGCTGGCGCGCGCTCTGAGACG.AAA.C.gg.gg.ggcgct                                                                                            <  1:1835373/101‑6 (MQ=255)
           aCAGCGCGCCTATGCGATGGTCGATGCGATCGACTGGCCCGAGGGGTTCTGCCGCCGCGACATCGGCTGGGGGGCGCCCCGGGGCG.AAA.C.gg.gg.ggcgggg                                                                                           >  1:2263395/1‑95 (MQ=255)
            cAGCGCGCCTATGCGATGGTCGATGCGATCGACTGGCCCGAGGGGTTCTGCCGCCGCGACATCGGCTGGCGCGCGCTCTGAGACG.AAA.C.cg.gg.gggggggg                                                                                          >  2:4306342/1‑94 (MQ=255)
               cgcgcCTATGCGGGGGTCGATGCGATCGACTGGCCCGAGGGGTTCTGCCGCCGCGACATCGGCTGGCGCGCGCTCTGAGACG.AAA.CGGG.GG.G.CGCTGcccc                                                                                       <  1:4101395/101‑1 (MQ=255)
               cgcgcCTATGCGATGGTCGATGCGATCGACTGGCCCGAGGGGTTCTGCCGCCGCGACATCGGCTGGCGCGCGCTCTGAGACG.AAA.CGGG.GG.G.CGCTTcccc                                                                                       <  2:2563068/101‑1 (MQ=255)
                gcgcCTATGCGATGGTCGATGCGATCGACTGGCCCGAGGGGTTCTGCCGCCGCGACATCGGCTGGCGCGCGCTCTGAGACG.AAA.CGGG.GG.G.CGCTGccccc                                                                                      <  1:1037177/101‑1 (MQ=255)
                 cgcCTATGCGATGGTCGATGCGATCGACTGGCCCGAGGGGTTCTGCCGCCGCGACATCGGCTGGCGCGCGCTCTGAGACGAAAA.G.GG.GG.G.GGCTGgccccc                                                                                     >  1:315395/1‑101 (MQ=255)
                    cTATGCGATGGTCGATGCGATCGACTGGCCCGAGGGGTTCTGCCGCCGCGACATCGGCTGGCGCGCGCTCTGAGACG.AAA.CGGG.GG.G.CGCTGcccccccct                                                                                  >  2:6861942/1‑98 (MQ=255)
                    cTATGCGATGGTCGATGCGATCGACTGGCCCGAGGGGTTCTGCCGCCGCGACATCGGCTGGCGCGCGCTCTGAGACG.AAA.CGGG.GG.G.CGCTGCCCCcctct                                                                                  <  2:2443510/101‑2 (MQ=255)
                    cTATGCGATGGTCGATGCGATCGACTGGCCCGAGGGGTTCTGCCGCCGCGACATCGGCTGGCGCGCGATCTGAGACGAAAA.G.GG.GG.G.CGCAGcccccccct                                                                                  >  1:2023703/1‑98 (MQ=255)
                     tATGCGATGGTCGATGCGATCGACTGGCCCGAGGGGTTCTGCCGCCGCGACATCGGCTGGCGCGCGCTCTGAGACG.AAA.CGGG.GG.G.GGCTGGCCCcctctt                                                                                 >  2:4927446/1‑99 (MQ=255)
                     tATGCGATGGTCGATGCGATCGACTGGCCCGAGGGGTTCTGCCGCCGCGACATCGGCTGGCGCGCGCTCTGAGACG.AAA.CGGG.GG.G.CGCTGCCCCcctctt                                                                                 <  1:6188491/101‑3 (MQ=255)
                      aTGCGATGGCCGATGCGATCGACTGGCCCGAGGGGTTCTGCCGCCGCGACATCGGCTGGCGCGCGCTCTGAGACG.AAA.CGGG.GG.G.CGCTGCCCCCCTCT.T..g                                                                             <  1:6537275/101‑1 (MQ=255)
                          gATGGTCGATGCGATCGACTGGCCCGAGGGGTTCTGCCGCCGCGACATCGGCTGGCGCGCGCTCTGAGACG.AAACC.GG.GG.G.GGGCGccccccccttggcct                                                                            >  1:2578133/1‑92 (MQ=255)
                          gATGGTCGATGCGATCGACTGGCCCGAGGGGTTCTGCCGCCGCGACATCGGCTGGCGCGCGCTCTGAGACG.AAA.CGGG.GG.G.GGCTGccccccccttggcct                                                                            >  1:3698280/1‑92 (MQ=255)
                          gATGGTCGATGCGATCGACTGGCCCGAGGGGTTCTGCCGCCGCGACATCGGCTGGCGCGCGCTCTGAGACG.AAA.CGGG.GG.G.CGCTGCCCCCCTCT.T..ggcct                                                                         <  2:744056/101‑5 (MQ=39)
                           atcgTCGATGCGATCGACTGGCCCGAGGGGTTCTGCCGCCGCGACATCGGCTGGCGCGCGCTCTGAGACG.AAA.CGGG.GG.G.CGCTGCCCCCCTCT.T..ggcctt                                                                        <  2:1325268/98‑6 (MQ=38)
                           aTGGTCGATGCGATCGACTGGCCCGAGGGGTTCTGCCGCCGCGACATCGGCTGGCGGGCGCTCTGAGACGAAAA.G.GG.GG.G.CGCTGccccccccttggccct                                                                           >  2:2547494/1‑91 (MQ=255)
                           aTGGTCGATGCGATCGACTGGCCCGAGGGGTTCTGCCGCCGCGACATCGGCTGGCGCGCGCTCTGAGACC.AAA.CGGG.GG.G.GGCTGCCCCCCTCT.T..gggctt                                                                        >  2:2980044/1‑96 (MQ=255)
                            tGGTCGATGCGATCGACTGGCCCGAGGGGTTCTGCCGCCGCGACATCGGCTGGCGCGCGCTCTGAGACG.AAA.CGGG.GG.G.GGCTGCCCCCCTCT.T..ggccttc                                                                       >  1:2651264/1‑95 (MQ=255)
                            tGGTCGATGCGATCGACTGGCCCGAGGGGTTCTGCCGCCGCGACATCGGCTGGCGCGCGCTCTGAGACG.AAA.CGGG.GG.G.CGCTGCCCCCCTCT.T..ggccttc                                                                       <  1:2877183/101‑7 (MQ=38)
                            tGGTCGATGCGATCGACTGGCCCGAGGGGTTCTGCCGCCGCGACATCGGCTGGCGCGCGCTCTGAGACG.AAA.CGGG.GG.G.CGCTGCCCCCCTCT.T..ggccttc                                                                       <  2:3970458/101‑7 (MQ=38)
                             ggTCGATGCGATCGACTGGCCCGAGGGGTTCTGCCGCCGCGACATCGGCTGGCGCGCGCTCTGAGACG.AAA.CGGG.GG.G.CGCTGCCCCCCTCT.T..ggccttcg                                                                      <  2:3567438/101‑8 (MQ=38)
                             ggTCGATGCGATCGACTGGCCCGAGGGGTTCTGCCGCCGCGACATCGGCTGGCGCGCGCTCTGAGACG.AAA.CGGG.GG.G.CGCTGCCCCCCTCT.T..ggcctccg                                                                      <  1:1939216/101‑8 (MQ=38)
                                cGATGCGATCGAGTGGCCCGAGGAGTTCTGCCGCCGCGACATCGGCTGGCGCGCGCTCTGAGACG.AAA.CGGG.GG.G.CGCTGCCCCCCTCT.T..GGCCTTCGGcc                                                                   <  1:7075642/101‑1 (MQ=25)
                                   tGCGATCGACTGGCCCGAGGGGTTCTGCCGCCGCGACATCGGCTGGCGCGCGCTCTGAGACG.AAA.CGGG.GG.G.CGCTGCCCCCCTCT.T..GGCCTTCGGCCA.At                                                               <  2:2926537/101‑1 (MQ=37)
                                     cGATCGACTGGCCCGAGGGGTTCTGCCGCCGCGACATCGGCTGGCGCGCGCTCTGAGACG.AAA.CGGG.GG.G.GGCTGCCCCCCTCC.TTGG..CCTCGGCCC.ATTc                                                             >  2:5970852/1‑101 (MQ=255)
                                     cGATCGACTGGCCCGAGGGGTTCTGCCGCCGCGACATCGGCTGGCGCGCGCTCTGAGACG.AAA.CGGG.GG.G.CGCTGCCCCCCTCT.T..GGCCTTCGGGCA.ATTc                                                             >  2:1172561/1‑101 (MQ=37)
                                     cGATCGACTGGCCCGAGGGGTTCTGCCGCCGCGACATCGGCTGGCGCGCGCTCTGAGACG.AAA.CGGG.GG.G.CGCTGCCCCCCTCT.T..GGCCTTCGGCCA.ATTc                                                             <  2:5290542/101‑1 (MQ=37)
                                         cGACTGGCCCGAGGGGTTCTGCCGCCGCGACATCGGCTGGCGCGCGCTCTGAGACG.AAA.CGGG.GG.G.CGCTGCCCCCCTCT.T..GGCCTTCGGCCA.ATTCAccc                                                         <  2:4793399/101‑1 (MQ=37)
                                          gACTGGCCCGAGGGGTTCTGCCGCCGCGACATCGGCTGGCGCGCGCCCTGAGACG.AAA.CGGG.GG.G.GGCTGCCCCCCCCCTT..GC.CTTCGGCCACATTCCccc                                                         >  2:83510/1‑101 (MQ=255)
                                           gcTGGCCCGAGGGGTTCTGCCGCCGCGACATCGGCTGGCGCGCGCTCTGAGACG.AAA.C.GGTGG.G.CGCTGCCCCCCTCT.T..GGCCTTCGGCCA.ATTCAccccc                                                       <  2:3046775/100‑1 (MQ=37)
                                           aCTGGCCCGAGGGGTTCTGCCGCCGCGACATCGGCTGGCGCGCGCTCTGAGACG.AAA.CGGG.GG.G.CGCTGCCCCCCTCT.T..GGCCTTCGGCCA.ATTCAccccc                                                       <  1:1752097/101‑1 (MQ=25)
                                           aCTGGCCCGAGGGGTTCTGCCGCCGCGACATCGGCTGGCGCGCGCTCTGAGACG.AAA.CGGG.GG.G.CGCTGCCCCCCTCT.T..GGCCTTCGGCCA.ATTCAccccc                                                       <  1:3988735/101‑1 (MQ=25)
                                               gCCCGAGGGGTTCTGCCGCCGCGACATCGGCTGGCGCGCGCTCTGAGACG.AAA.CGGG.GG.G.CGCTGCCCCCCTCT.T..GGCCTTCGGCCA.ATTCACCCCCCggg                                                   <  2:197939/101‑1 (MQ=25)
                                               gCCCGAGGGGTTCTGCCGCCGCGACATCGGCTGGCGCGCGCTCTGAGACG.AAA.CGGG.GG.G.CGCTGCCCCCCTCT.T..GGCCTTCGGCCA.ATTCACCCCCCggg                                                   <  2:2631863/101‑1 (MQ=25)
                                                 ccGAGGGGTTCTGCCGCCGCGACATCGGCTGGCGCGCGCTCTGAGACG.AAA.CGGG.GG.G.CGCTGCCCCCCTCT.T..GGCCTTCGGCCA.ATTCACCCCCCGGGat                                                 <  1:5537372/101‑1 (MQ=25)
                                                 ccGAGGGGTTCTGCCGCCGCGACATCGGCTGGCGCGCGCTCTGAGACG.AAA.CGGG.GG.G.CGCTGCCCCCCCCT.T..GGCCTTCGGCCA.ATTCACCCCCCCGGat                                                 >  2:7002147/1‑101 (MQ=25)
                                                    gggggTTCTGCCGCCGCGACATCGGCTGGCGCGCGCTCTGAGACG.AAA.CGGG.GG.G.CGCTGCCCCCCTCT.T..GGCCTTCGGCCA.ATTCACCCCCCGGGATAtt                                              <  2:4684788/100‑1 (MQ=21)
                                                    aGGGGTTCTGCCGCCGCGACATCGGCTGGCGCGCGCTCTGAGACG.AAA.CGGG.GG.G.CGCTGCCCCCCTCT.T..GGCCTTCGGCCC.ATTCACCCCCCCGGCTAtt                                              >  1:6354303/1‑101 (MQ=21)
                                                     ggtgTTCTGCCGCCGCGACATCGGCTGGCGCGCGCTCTGAGACG.AAA.CGGG.GG.G.CGCTGCCCCCCTCT.T..GGCCTTCGGCCA.ATTCACCCCCCGGGATAttt                                             <  2:6633905/98‑1 (MQ=21)
                                                      gggTTGTGCCGCCGCGACATCGGCTGGCGCGCGCTCTGAGACG.AAA.CGGG.GG.G.CGCTGCCCCCCTCT.T..GGCCTTCGGCCA.ATTCACCCCCCGGGATATTTc                                            <  1:5731627/101‑1 (MQ=21)
                                                      gggTTCTGCCGCCGCGACATCGGCTGGCGCGCGCTCTGAGACG.AAA.CGGG.GG.G.CGCTGCCCCCCTCT.TTGG..CCTCGGCCA.ATTCACCCCCCCGGATCTTTc                                            >  2:1133395/1‑101 (MQ=25)
                                                       ggTTCTGCCGCCGCGACATCGGCTGGGGGGCGCTCTGAGACGAAAA.G.GG.GG.G.GGCTGCCCCCCCCT.T..GGCCTTCGGCCC.ATTCCCCCCCCCGCGTATTTCg                                           >  1:4713374/1‑101 (MQ=255)
                                                       ggTTCTGCCGCCGCGACATCGGCTGGCGCGCGCTCTGAGACG.AAA.CGGG.GG.G.CGCTGCCCCCCTCT.T..GGCCTTCGGCCA.ATTCACCCCCCGGGATATTTCg                                           <  1:6697804/101‑1 (MQ=21)
                                                           cTGCCGCCGCGACATCGGCTGGCGCGCGCTCTGAGACG.AAA.CGGG.GG.G.CGCTGCCCCCCTCT.T..GGCCTTCGGCCA.ATTCACCCCCCGGGATATTTCGGGCa                                       <  2:3394892/101‑1 (MQ=21)
                                                           cTGCCGCCGCGACATCGGCTGGCGCGCGCTCTGAGACG.AAA.CGGG.GG.G.CGCTGCCCCCCTCT.T..GGCCTTCGGCCA.ATTCACCCCCCGGGATATTTCGGGCa                                       <  1:463114/101‑1 (MQ=21)
                                                           cTGCCGCCGCGACATCGGCTGGCGCGCGCTCTGAGACG.AAA.CGGG.GG.G.CGCTGCCCCCCTCT.T..GGCCTTCGGCCA.ATTCACCCCCCGGGATATTTCGGGCa                                       <  2:6299472/101‑1 (MQ=21)
                                                             gccgccGCGACATCGGCTGGCGCGCGCTCTGAGACG.AAA.CGGG.GG.G.CGCTGCCCCCCTCT.T..GGCCTTCGGCCA.ATTCACCCCCCGGGATATTTCGGGCAAg                                     <  1:2904986/101‑1 (MQ=21)
                                                             gccgccGCGACATCGGCTGGCGCGCGCTCTGAGACG.AAA.CGGG.GG.G.CGCTGCCCCCCTCT.T..GGCCTTCGGCCA.ATTCACCCCCCGGGATATTTCGGGCAAg                                     <  1:5214826/101‑1 (MQ=21)
                                                              ccgccgCGACATCGGCTGGCGCGCGCTCTGAGACG.AAA.CGGG.GG.G.CGCTGCCCCCCTCT.T..GGCCTTCGGCCA.ATTCACCCCCCGGGATATTTCGGGCAAGa                                    <  1:1342072/101‑1 (MQ=21)
                                                               cgccgcGACATCGGCTGGCGCGCGCTCTGAGACG.AAA.CGGG.GG.G.CGCTGCCCCCCTCT.T..GGCCTTCGGCCA.ATTCACCCCCCGGGATATTTCGGGCAAGAt                                   <  2:5073261/101‑1 (MQ=21)
                                                                gccgcGACATCGGCTGGCGCGCGCTCTGAGACG.AAA.CGGG.GG.G.CGCTGCCCCCCTCT.T..GGCCTTCGGCCA.ATTCACCCCCCGGGATATTTCGGGCAAGATa                                  <  1:3848319/101‑1 (MQ=21)
                                                                gccgcGACATCGGCTGGCGCGCGCTCTGAGACG.AAA.CGGG.GG.G.CGCTGCCCCCCTCT.T..GGCCTTCGGCCA.ATTCACCCCCCGGGATATTTCGGGCAAGATa                                  <  2:3678331/101‑1 (MQ=21)
                                                                gccgcGACATCGGCTGGCGCGCGCTCTGAGACG.AAA.CGGG.GG.G.CGCTGCCCCCCTCT.T..GGCCTTCGGCCA.ATTCACCCCCCGGGATATTTCGGGCAAGATa                                  <  2:606707/101‑1 (MQ=21)
                                                                gccgcGACATCGGCTGGCGCGCGCTCTGAGACG.AAA.CGGG.GG.G.CGCTGCCCCCCTCT.T..GGCCTTCGGCCA.ATTCACCCCCCGGGATATTTCGGGCAAGATa                                  <  1:5047668/101‑1 (MQ=21)
                                                                 ccgcGACATCGGCTGGCGCGCGCTCTGAGACG.AAA.CGGG.GG.G.CGCTGCCCCCCTCT.T..GGCCTTCGGCCA.ATTCACCCCCCGGGATATTTCGGGCAAGATaa                                 <  2:5601198/101‑1 (MQ=21)
                                                                 ccacgACATCGGCTGGGGCACGCTCTGAGACG.AAA.C.GG.GG.GACGCTGCCCCCCTCT.T..GGCCTTCGGCCA.ATTCACCCCCCGGGATATTTCGGGCAAGATaa                                 <  2:7078549/98‑1 (MQ=18)
                                                                  cgcgACATCGGCTGGCGCGCGCTCTGAGACG.AAA.CGGG.GG.G.CGCTGCCCCCCTCT.T..GGCCTTCGGCCA.ATTCACCCCCCGGGATATTTCGGGCAAGATaaa                                <  1:1607458/101‑1 (MQ=21)
                                                                  caggACAGCGGCTGGCGCGCGCTCTGAGACG.AAA.C.GG.GGTG.CGCTGCCCCCCTCT.T..GTCCTTCGGCCA.ATTCACCCCCCGGGATATTTCGGGCAAGATaaa                                <  2:2147433/98‑1 (MQ=21)
                                                                     gACATCGGCTGGCGCGCGCTCTGAGACG.AAA.CGGG.GG.G.CGCTGCCCCCCTCT.T..GGCCTTCGGCCA.ATTCACCCCCCGGGATATTTCGGGCAAGATAaagaa                             >  1:6463733/1‑101 (MQ=21)
                                                                     gACATCGGCTGGCGCGCGCTCTGAGACG.AAA.CGGG.GG.G.CGCTGCCCCCCTCT.T..GGCCTTCGGCCA.ATTCACCCCCCGGGATATTTCGGGCAAGATAaagaa                             <  1:5441870/101‑1 (MQ=21)
                                                                      aCATCGGCTGGCGCGCGCTCTGAGACG.AAA.CGGG.GG.G.CGCTGCCCCCCTCT.T..GGCCTTCGGCCA.ATTCACCCCCCGGGATATTTCGGGCGAGATAaagaag                            <  2:317249/101‑1 (MQ=21)
                                                                      aCATCGGCTGGCGCGCGCTCTGAGACG.AAA.CGGG.GG.G.CGCTGCCCCCCTCT.T..GGCCTTCGGCCA.ATTCACCCCCCGGGATATTTCGGGCAAGATAaagaag                            <  1:5076772/101‑1 (MQ=21)
                                                                      aCATCGGCTGGCGCGCGCTCTGAGACG.AAA.CGGG.GG.G.CGCTGCCCCCCTCT.T..GGCCTTCGGCCA.ATTCACCCCCCGGGATATTTCGGGCAAGATAaagaag                            <  2:1003257/101‑1 (MQ=21)
                                                                        aTCGGCTGGCGCGCGCTCTGAGACG.AAA.CGGG.GG.G.CGCTGCCCCCCTCT.T..GGCCTTCGGCCA.ATTCACCCCCCGGGATATTTCGGGCAAGATAAAGAAggg                          >  1:3054090/1‑100 (MQ=21)
                                                                        aTCGGCTGGCGCGCGCTCTGAGACG.AAA.CGGG.GG.G.CGCTGCCCCCCTCT.T..GGCCTTCGGCCA.ATTCACCCCCCGGGATATTTCGGGCAAGATAAAGAAGGt                          <  1:4425903/101‑1 (MQ=21)
                                                                          cGGCTGGCGCGCGCTCTGAGACG.AAA.CGGG.GG.G.CGCTGCCCCCCTCT.T..GGCCTTCGGCCA.ATTCACCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCt                        <  1:3665733/101‑1 (MQ=21)
                                                                          cGGCTGGCGCGCGCTCTGAGACG.AAA.CGGG.GG.G.CGCTGCCCCCCTCT.T..GGCCTTCGGCCA.ATTCACCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCt                        <  2:3860472/101‑1 (MQ=21)
                                                                           ggCTGGCGCGCGCTCTGAGACG.AAA.CGGG.GG.G.CGCTGCCCCCCTCT.T..GGCCTTCGGCCA.ATTCACCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCtt                       <  1:4252317/101‑1 (MQ=21)
                                                                           ggCTGGCGCGCGCTCTGAGACG.AAA.CGGG.GG.G.CGCTGCCCCCCTCT.T..GGCCTTCGGCCA.ATTCACCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCtt                       <  2:1449346/101‑1 (MQ=21)
                                                                             cttgcgcgctcTCTGAGACG.AAA.CGGG.GG.G.CGCTGCCCCCCTCT.T..GTCCTTCGGCCA.ATTCACCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCtttt                     <  1:894438/98‑1 (MQ=21)
                                                                             cTGGCGCGCGCTCTGAGACG.AAA.CGGG.GG.G.CGCTGCCCCCCTCT.T..GGCCTTCGGCCA.ATTCACCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCtttt                     <  2:5679357/101‑1 (MQ=21)
                                                                              tGGCGCGCGCTCTGAGACG.AAG.CGGG.GG.G.CGCTGCCCCCCTCT.T..GGCCTTCGGCCA.ATTCACCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTc                    <  2:2734820/101‑1 (MQ=21)
                                                                              tGGCGCGCGCTCTGAGACG.AAA.CGGG.GG.G.CGGTGCCCCCCTCT.T..GGCCTTCGGGCA.ATTCACCCCCCGGGAAATTTCGGCGAAGATAAAGAAGGCCTTTTc                    >  1:6946152/1‑101 (MQ=18)
                                                                                gcgcgcgcTCTGAGACG.AAA.CGGG.GG.G.CGCTGCCCCCCTCT.T..GGCCTTCGGCCA.ATTCACCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCt                  >  2:2787478/1‑101 (MQ=21)
                                                                                 cgcgcgcTCTGAGACG.AAA.CGGG.GG.G.CGCTGCCCCCCTCT.T..GGCCTTCGGCCA.ATTCACCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCtt                 <  2:5420300/101‑1 (MQ=21)
                                                                                  gcgcgcTCTGAGACG.AAA.CGGG.GG.G.CGCTGCCCCCCTCT.T..GGCCTTCGGCCA.ATTCACCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCttt                <  2:3488378/101‑1 (MQ=21)
                                                                                  gcgcgcTCTGAGACG.AAA.CGGG.GG.G.CGCTGCCCCCCTCT.T..GGCCTTCGGCCA.ATTCACCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCttt                <  2:4426601/101‑1 (MQ=21)
                                                                                  gcgcgcTCTGAGACG.AAA.CGGG.GG.G.CGCTGCCCCCCTCT.T..GGCCTTCGGCCA.ATTCACCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCttt                <  1:2170609/101‑1 (MQ=21)
                                                                                   cgcgcTCTGAGACG.AAA.CGGG.GG.G.CGCTGCCCTCCTCT.T..GGCCTTCGGCCA.ATTCACCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTc               <  2:5247186/101‑1 (MQ=21)
                                                                                    gcgcTCTGAGACG.AAA.CGGG.GG.G.CGCTGCCCCCCTCT.T..GGCCTTCGGCCA.ATTCACCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCa              <  1:2482257/101‑1 (MQ=21)
                                                                                    gcgcTCTGAGACG.AAA.CGGG.GG.G.CGCTGCCCCCCTCT.T..GGCCTTCGGCCA.ATTCACCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCa              >  1:808102/1‑101 (MQ=21)
                                                                                    gcgcTCTGAGACG.AAA.CGGG.GG.G.CGCTGCCCCCCTCT.T..GGCCTTCGGCCA.ATTCACCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCa              <  2:6213452/101‑1 (MQ=21)
                                                                                     cgcTCTGAGACG.AAA.CGGG.GG.G.CGCTGCCCCCCTCT.T..GGCCTTCGGCCA.ATTCACCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCAt             <  1:3411622/101‑1 (MQ=21)
                                                                                     cgcTCTGAGACG.AAA.CGGG.GG.G.CGCTGCCCCCCTCT.T..GGCCTTCGGCCA.ATTCACCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCAt             <  1:533247/101‑1 (MQ=21)
                                                                                     cgcTCTGAGACG.AAA.CGGG.GG.G.CGCTGCCCCCCTCT.T..GGCCTTCGGCCA.ATTCACCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCAt             >  1:2331370/1‑101 (MQ=21)
                                                                                      gcTCTGAGACG.AAA.CGGG.GG.G.CGCTGCCCCCCTCT.T..GGCCTTCGGCCA.ATTCACCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATc            >  1:5698528/1‑101 (MQ=21)
                                                                                      gcTCTGAGACG.AAA.CGGG.GG.G.CGCTGCCCCCCTCT.T..GGCCTTCGGCCA.ATTCACCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATc            <  2:1218978/101‑1 (MQ=21)
                                                                                       ctctGAGACG.AAA.CGGG.GG.G.CGCTGCCCCCCTCT.T..GGCCTTCGGCCA.ATTCACCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCt           <  1:1014956/101‑1 (MQ=21)
                                                                                       ctctGAGACG.AAA.CGGG.GG.G.CGCTGCCCCCCTCT.T..GGCCTTCGGCCA.ATTCACCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCt           <  2:2333179/101‑1 (MQ=21)
                                                                                          tGAGACG.AAA.CGGG.GG.G.CGCTGCCCCCCTCT.T..GGCCTTCGGCCA.ATTCACCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGc        <  2:7017151/101‑1 (MQ=21)
                                                                                            agaCG.AAA.CGGG.GG.G.CGCTGCCCCCCTCT.T..GGCCTTCGGCCA.ATTCACCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGccc      <  2:3224992/101‑1 (MQ=21)
                                                                                             gaCG.AAA.CGGG.GG.G.CGCTGCCCCCCTCT.T..GGCCTTCGGCCA.ATTCACCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCg     <  2:4081874/101‑1 (MQ=21)
                                                                                              aCG.AAA.CGGG.GG.G.CGCTGCCCCCCTCT.T..GGCCTTCGGCCA.ATTCACCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGa    >  1:5882428/1‑101 (MQ=21)
                                                                                              aCG.AAA.CGGG.GG.G.CGCTGCCCCCCTCT.T..GGCCTTCGGCCA.ATTCACCCCCCGGGATAATTTCGGCAAGATAAAGAAGGGCTTTTTCCTTCACCTTTCCCgc    >  1:2427410/1‑100 (MQ=9)
                                                                                                g.AAA.CGGG.GG.G.CGCGGCCCCCCTCT.T..GGCCTTCGGCCA.ATTCACCCCCCGGGATATTTCGGGCAAGATAAAGAAGGGCTTTTCCTTTCATCTTGCCCGaaa  <  2:1775307/101‑1 (MQ=21)
                                                                                                       |                                                                                                      
GCCGAGGCGCAACAGCGCGCCTATGCGATGGTCGATGCGATCGACTGGCCCGAGGGGTTCTGCCGCCGCGACATCGGCTGGCGCGCGCTCTGAGACG.AAA.C.GG.GG.G.CGCTGCCCCCCTCC.T..G..CTTCGGCCA.ATTCACCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAA  >  NC_014034/2566436‑2566629

Base quality scores:  ATCG  < 3 ≤  ATCG  < 15 ≤  ATCG  < 25 ≤  ATCG  < 33 ≤  ATCG  < 41 ≤  ATCG