Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_014034 2,566,551 +T 3.3% intergenic (+23/+171) 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,5511.→T 3.3% 8.7 0.0 0.0121intergenic (+23/+171)purD/sdaAphosphoribosylamine‑‑glycine ligase/L‑serine ammonia‑lyase
Reads supporting (aligned to +/- strand):  new base (1/3):  ref base (28/89):  total (29/92)
Fisher's exact test strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 1.00e+00

Only 100 of 137 total aligned reads displayed.

GCGCCTATGCGATGGTCGATGCGATCGACTGGCCCGAGGGGTTCTGCCGCCGCGACATCGGCTGGCGCGCGCTCTGAGACGAAAC.GGG.GGCGCTGCCCCC.CTC.CT.G..CTTC.GGCCAATTCA.CCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCCGGGGGAGCC  >  NC_014034/2566452‑2566644
                                                                                                      |                                                                                                   
  tccTAGGCGATGGTCGATGCGTTCGACGGGCCTGAGGGGTTCTGCCGCCGCGACATCGGCTGGCGCGCGCTCTGAGACGAAACGGGG.GGCGCTGCCCCc.ct                                                                                                   <  2:585315/100‑1 (MQ=255)
  gcCTATGCGATGGTCGATGCGATCGACTGGCCCGAGGGGTTCTGCCGCCGCGACATCGGCTGGCGCGCGCTCTGAGACGAAACGGGG.GGCGCTGCCCCc.ct                                                                                                   <  2:846157/101‑1 (MQ=255)
   ccTATGCGATGGTCGATGCGATCGACTGGCCCGAGGGGTTCTGCCGCCGCGACATCGGCTGGCGCGCGCTCTGAGACGAAACGGGG.GGCGCTGCCCCc.ctc                                                                                                  <  2:878730/101‑1 (MQ=255)
     tATGCGATGGTCGATGCGATCGACTGGCCCGAGGGGTTCTGCCGCCGCGACATCGGCTGGCGCGCGCTCTGAGACGAAACGGGG.GGGGCTGCCCCc.ccc.ct                                                                                               >  1:1789161/1‑101 (MQ=255)
      aTGCGATGGTCGATGCGATCGACTGGCCCGAGGGGTTCTGCCGCCGCGACATCGGCTGGCGCGCGCTCTGAGACGAAACGGGG.GGCGCTGCCCCC.CTC.TT.g                                                                                             <  1:2279289/101‑1 (MQ=255)
        gCGATGGTCGATGCGATCGACTGGCCCGAGGGGTTCTGCCGCCGCGACATCGGCTGGCGCGCGCTCTGAGACGAAACGGGG.GGCGCTGCCCCC.CTC.TT.ggc                                                                                           <  1:3437533/101‑3 (MQ=255)
         cGATGGTCGATGCCATCGACTGGCCCGAGGGGTTCTGCCGCCGCGACATCGGCTGGCGCGCGCTCTGAGACGAAACGGGG.GGCGCTGCCCCC.CTC.TT.ggcc                                                                                          <  2:2937542/101‑4 (MQ=38)
          gATGGTCGATGCGATCGACTGGCCCGAGGGGTTCTGCCGCCGCGACATCGGCTGGCGCGCGCTCTGAGACGAAACGGGG.GGCGCTGCCCCC.CTC.TT.ggcct                                                                                         <  1:1058341/101‑5 (MQ=39)
                 gATGCGCTCGGCTGGCCCGAGGGGTTCTGCCGCCGCGACATCGGCTGGCGCGCGCTCTGAGACGAAACGGGG.GGCGCTGCCCCC.CTC.TT.GGCCTTC.GGCCa                                                                                 <  1:705356/101‑1 (MQ=25)
                 gATGCGATCGACTGGCCCGAGGGGTTCTGCCGCCGCGACATCGGCTGGCGCGCGCTCTGAGACGAAACGGGG.GGCGCTGCCCCC.CTC.TT.GGCCTTC.GGCCa                                                                                 <  2:2743559/101‑1 (MQ=37)
                   tGCGATCGACTGGCCCGAGGGGTTCTGCCGCCGCGACATCGGCTGGCGCGCGCTCTGAGACGAAACGGGG.GGCGCTGCCCCC.CTC.TT.GGCCTTC.GGCCAAt                                                                               <  2:3129681/101‑1 (MQ=37)
                    gCGATCGACTGGCCCGAGGGGTTCTGCCGCCGCGACATCGGCTGGCGCGCGCTCTGAGACGAAACGGGG.GGCGCTGCCCCC.CTC.TT.GGCCTTC.GGCCAAtt                                                                              <  1:3672165/101‑1 (MQ=37)
                    gCGATCGACTGGCCCGAGGGGTTCTGCCGCCGCGACATCGGCTGGCGCGCGCTCTGAGACGAAACGGGG.GGCGCTGCCCCC.CTC.TA.GGCCTTC.GGCCAAtt                                                                              <  1:3566316/101‑1 (MQ=37)
                     cGATCGACTGGCCCGAGGGGTTCTGCCGCCGCGACATCGGCTGGCGCGCGCTCTGAGACGAAACGGGG.GGCGCTGCCCCC.CTC.TT.GGCCTTC.GGCCAATTc                                                                             <  1:2286059/101‑1 (MQ=37)
                        tCGACTGGCCCGAGGGGTTCTGCCGCCGCGACATCGGCTGGCGCGCGCTCTGAGACGAAACGGGG.GGCGCTGCCCCC.CTC.TT.GGCCTTC.GGCCAATTCA.cc                                                                         <  1:486815/101‑1 (MQ=37)
                        tCGACTGGCCCGAGGGGTTCTGCCGCCGCGACATCGGCTGGCGCGCGCTCTGAGACGAAACGGGG.GGCGCTGCCCCC.CCC.TT.GGCCCTC.GGCCAATTCC.cc                                                                         >  1:2553452/1‑101 (MQ=25)
                         cGACTGGCCCGAGGGGTTCTGCCGCCGCGACATCGGCTGGCGCGCGCTCTGAGACGAAACGGGG.GGCGCTGCCCCC.CTC.TT.GGCCTTC.GGCCAATTCA.ccc                                                                        <  2:659481/101‑1 (MQ=37)
                           gcgggCCCGAGGGGTTCTGCCGCCGCGACATCGGCTGGCGCGCGCTCTGAGACGAAACGGGG.GGCGCTGCCCCC.CTC.TT.GGCCTTC.GGCCAATTCA.ccccc                                                                      <  2:2391411/98‑1 (MQ=25)
                           ccTGGCCCGAGGGGTTCTGCCGCCGCGACATCGGCTGGCGCGCGCTCTGAGACGAAACGGGG.GGCGCTGCCCCC.CTC.TT.GGCCTTC.GGCCAATTCA.ccccc                                                                      <  2:203294/100‑1 (MQ=25)
                           aCTGGCCCGAGGGGTTCTGCCGCCGCGACATCGGCTGGCGCGCGCTCTGAGACGAAACGGGG.GGCGCTGCCCCC.CTC.TT.GGCCTTC.GGCCAATTCA.ccccc                                                                      <  1:1187826/101‑1 (MQ=25)
                            cTGGCCCGAGGGGTTCTGCCGCCGCGACATCGGCTGGCGCGCGCTCTGAGACGAAACGGGG.GGCGCTGCCCCC.CTC.TT.GGCCTTC.GGCCAATTCA.cccccc                                                                     <  2:2895798/101‑1 (MQ=25)
                             tGGCCCGAGGGGTTCTGCCGCCGCGACATCGGCTGGCGCGCGCTCTGAGACGAAACGGGG.GGCGCTGCCCCC.CTC.TT.GGCCTTC.GGCCAATTCA.CCCCCCg                                                                    <  1:2181231/101‑1 (MQ=25)
                              ggCCCGAGGGGTTCTGCCGCCGCGACATCGGCTGGCGCGCGCTCTGAGACGAAACGGGG.GGCGCTGCCCCC.CTC.TT.GGCCTTC.GGCCAATTCA.CCCCCCgg                                                                   <  1:2944444/101‑1 (MQ=25)
                              ggCCCGAGGGGTTCTGCCGCCGCGACATCGGCTGGCGCGCGCTCTGAGACGAAACGGGG.GGCGCTGCCCCC.CTC.TT.GGCCTTC.GGCCAATTCA.CCCCCCgg                                                                   <  1:503621/101‑1 (MQ=25)
                              ggCCCGAGGGGTTCTGCCGCCGCGACATCGGCTGGCGCGCGCTCTGAGACGAAACGGGG.GGCGCTGCCCCC.CTC.TT.GGCCTTC.GGCCAATTCA.CCCCCCgg                                                                   <  1:1099121/101‑1 (MQ=25)
                               gCCCGAGGGGTTCTGCCGCCGCGACATCGGCTGGCGCGCGCTCTGAGACGAAAC.GGGTGGGGCTGCCCCC.CTC.TT.GGCCCTC.CGCCAATTCA.CCCCCCcgg                                                                  >  2:3490303/1‑101 (MQ=25)
                                cccGAGGGGTTCTGCCGCCGCGACATCGGCTGGCGCGCGCTCTGAGACGAAACGGGG.GGCGCTGCCCCC.CTC.TT.GGCCTTC.GGCCAATTCA.CCCCCCGGGa                                                                 <  1:1829286/101‑1 (MQ=25)
                                 ccGAGGGGTTCTGCCGCCGCGACATCGGCTGGCGCGCGCTCTGAGACGAAACGGGG.GGGGCCGCCCCC.CCCTTT.GC.CTTC.GGCCAAATCA.CCCCCCcggga                                                                >  1:1978096/1‑99 (MQ=255)
                                 ccGAGGGGTTCTGCCGCCGCGACATCGGCTGGCGCGCGCTCTGAGACGAAACGGGG.GGCGCTGCCCCC.CTC.TT.GGCCTTC.GGCCAATTCA.CCCCCCGGGat                                                                <  2:2168907/101‑1 (MQ=25)
                                 ccGAGGGGTTCTGCCGCCGCGACATCGGCTGGCGCGCGCTCTGAGACGAAACGGGG.GGCGCTGCCCCC.CTC.TT.GGCCTTC.GGCCAATTCA.CCCCCCGGGat                                                                <  1:3194193/101‑1 (MQ=25)
                                 ccGAGGGGTTCTGCCGCCGCGACATCGGCTGGCGCGCGCTCTGAGACGAAACGGGG.GGCGCTGCCCCC.CTC.TT.GGCCTTC.GGCCAATTCA.CCCCCCGGGat                                                                <  2:4180170/101‑1 (MQ=25)
                                  cGAGGGGTTCTGCCGCCGCGACATCGGCTGGCGCGCGCTCTGAGACGAAACGGGG.GGGGCTGCCCCC.CCC.TT.GGCCTTC.CGCCAAATCC.CCCCCCCGGgta                                                               >  2:3168341/1‑101 (MQ=255)
                                  cGAGGGGTTCTGCCGCCGCGACATCGGCTGGCGCGCGCTCTGAGACGAAACGGGG.GGCGCTGCCCCC.CTC.TT.GGCCTTC.GGCCAATTCA.CCCCCCGGGata                                                               <  1:1708987/101‑1 (MQ=25)
                                   gAGGGGTTCTGCCGCCGCGACATCGGCTGGCGCGCGCTCTGAGACGAAACGGGG.GGCGCTGCCCCC.CTC.TT.GGCCTTC.GGCCAATTCA.CCCCCCCGGatat                                                              >  1:2119448/1‑101 (MQ=25)
                                    aGGGGTTCTGCCGCCGCGACATCGGCTGGCGCGCGCTCTGAGACGAAACGGGG.GGCGCTGCCCCC.CTC.TT.GGCCTTC.GGCCAATTCA.ACCCCCGGGATAtt                                                             >  2:1052333/1‑101 (MQ=21)
                                     ggggTTCTGCCGCCGCGACATCGGCTGGCGCGCGCTCTGAGACGAAACGGGG.GGCGCTGCCCCC.CTC.TT.GGCCTTC.GGCCAATTCA.CCCCCCGGGATAttt                                                            <  2:77585/101‑1 (MQ=21)
                                     ggggTTCTGCCGCCGCGACATCGGCTGGCGCGCGCTCTGAGACGAAACGGGG.GGCGCTGCCCCC.CTC.TT.GGCCTTC.GGCCAATTCA.CCCCCCGGGATAttt                                                            <  1:1697153/101‑1 (MQ=21)
                                          tCTGCCGCCGCGACATCGGCTGGCGCGCGCTCTGAGACGAAACGGGG.GGCGCTGCCCCC.CTC.TT.GGCCTTC.GGCCCATTAA.ACCCCCCGGATATTTTgggg                                                       >  2:1592017/1‑100 (MQ=21)
                                          tCTGCCGCCGCGACATCGGCTGGCGCGCGCTCTGAGACGAAACGGGG.GGCGCTGCCCCC.CTC.TT.GGCCTTC.GGCCAATTCA.CCCCCCGGGATATTTCGGGc                                                       <  1:1268244/101‑1 (MQ=21)
                                          tCTGCCGCCGCGACATCGGCTGGCGCGCGCTCTGAGACGAAACGGGG.GGCGCTGCCCCC.CTC.TT.GGCCTTC.GGCCAATTCA.CCCCCCGGGATATTTCGGGc                                                       <  2:229061/101‑1 (MQ=21)
                                           cTGCCGCCGCGACATCGGCTGGCGCGCGCTCTGAGACGAAACGGGG.GGCGCTGCCCCC.CTC.TT.GGCCTTC.GGCCAATTCA.CCCCCCCGGCTATTTCGGGCa                                                      >  2:861613/1‑101 (MQ=21)
                                            tGCCGCCGCGACATCGGCTGGCGCGCGCTCTGAGACGAAGCGGGG.GGCGCTGCCCCC.CTC.TT.GGCCTTC.GGCCAATTCA.CCCCCCGGGATATTTCGGGCaa                                                     <  1:473092/101‑1 (MQ=21)
                                            tGCCGCCGCGACATCGGCTGGCGCGCGCTCTGAGACGAAACGGGG.GGCGCTGCCCCC.CCC.TT.GGCCTTC.GGCCAATTCA.CCCCCCCGGAGAtttttggcca                                                     >  2:3957956/1‑94 (MQ=18)
                                               cgccgcGACATCGGCTGGCGCGCGCTCTGAGACGAAACGGGG.GGGGCTGGCCCC.CTC.TT.GGCCTTT.CGCCAATTCA.CCCCCCCGGACATTTCGGGCAAGaa                                                  >  1:3980246/1‑100 (MQ=255)
                                               cgccgcGACATCGGCTGGCGCGCGCTCTGAGACGAAACGGGG.GGCGCTGCCCCC.CTC.TT.GGCCTTC.GGCCAATTCA.CCCCCCGGGATATTTCGGGCAAGAt                                                  <  2:2060072/101‑1 (MQ=21)
                                               cgccgcGACATCGGCGGGCGCGCGCTCTGAGACGAAACGGGG.GGCGCTGCCCCC.CTC.TT.GGCCTTC.GGCCAATTCA.CCCCCCGGGATATTTCGGGCAAGAt                                                  <  2:523134/101‑1 (MQ=21)
                                                gccgcGACATCGGCTGGCGCGCGCTCTGAGACGAAACGGGG.GGCGCTGCCCCC.CTC.TT.GGCCTTC.GGCCAATTCA.CCCCCCGGGATATTTCGGGCAAGATa                                                 <  2:2066625/101‑1 (MQ=21)
                                                gccgcGACATCGGCTGGCGCGCGCTCTGAGACGAAACGGGG.GGCGCTGCCCCC.CTC.TT.GGCCTTC.GGCCAATTCA.CCCCCCGGGATATTTCGGGCAAGATa                                                 <  2:3108432/101‑1 (MQ=21)
                                                  cgcgACATCGGCTGGCGCGCGCTCTGAGACGAAACGGGG.GGCGCTGCCCCC.CTC.TT.GGCCTTC.GGCCCATTCA.CCCCCCCGGATATTTCGGGCAAGATaaa                                               >  1:3906941/1‑101 (MQ=21)
                                                   gcgACATCGGCTGGCGCGCGCTCTGAGACGAAACGGGG.GGCGCTGCCCCC.CTC.TT.GGCCTTC.GGCCAATTCA.CCCCCCGGGATATTTCGGGCAAGATAaag                                              <  1:430223/101‑1 (MQ=21)
                                                    cgACATCGGCTGGCGCGCGCTCTGAGACGAAACGGGG.GGCGCTGCCCCC.CTC.TT.GGCCTTC.GGCCAATTCA.CCCCCCGGGATATTTCGGGCAAGATAaaga                                             <  1:1811536/101‑1 (MQ=21)
                                                    cgACATCGGCTGGCGCGCGCTCTGAGACGAAACGGGG.GGCGCTGCCCCC.CTC.TT.GGCCTTC.GGCCAATTCA.CCCCCCGGGATATTTCGGGCAAGATAaaga                                             >  2:3176725/1‑101 (MQ=21)
                                                     gACATCGGCTGGCGCGCGCTCTGAGACGAAACGGGG.GGCGCTGCCCCC.CTC.TT.GGCCTTC.GGCCAATTCA.CCCCCCGGGATATTTCGGGCAAGATAaagaa                                            <  1:19493/101‑1 (MQ=21)
                                                      aCATCGGCTGGCGCTCGCTCTGAGACGAAACGGGG.GGCCCCGCCCCC.CTT.TTGG..CCTCGGGCCAATTCA.CCCCCCGCGCCGTTTCGGGGGCGATAaaaaag                                           >  1:4115095/1‑101 (MQ=9)
                                                        aTCGGCTGGCGCGCGCTCTGAGACGAAACGGGG.GGCGCTGCCCCC.CTC.TT.GGCCTTC.GGCCCATTCA.CCCCCCCGGAAATTTTGGGGAAGATAAAGAAggg                                         >  2:3039499/1‑100 (MQ=18)
                                                         tCGGCTGGCGCGCGCTCTGAGACGAAACGGGG.GGCGCTGCCCCC.CTC.TT.GGCCTTC.GGCCAATTCA.CCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTc                                        >  1:3848117/1‑101 (MQ=21)
                                                           ggCTGGGGCGCGCTCTGAGACGAAACGGGG.GGCGCTGCCCCC.CTC.TT.GGCCTTC.GGCCAATTCA.CCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCtt                                      <  1:3073594/101‑1 (MQ=21)
                                                           ggCTGGCGCGCCCTCTGAGACGAAACGGGG.GGCGCTGCCCCC.CTC.TT.GGCCTTC.GGCCAATTCA.CCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCtt                                      <  1:1943026/101‑1 (MQ=21)
                                                            gCTGGCGCGCGCTCTGAGACGAAACGGGG.GGCGCTGCCCCC.CTC.TT.GGCCTTC.GGCCAATTCA.CCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCttt                                     <  2:1078597/101‑1 (MQ=21)
                                                             cTGGCGCGCGCTCTGAGACGAAACGGGG.GGCGCTGCCCCC.CTC.TT.GGCCTTC.GGCCAATTCA.CCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCtttt                                    <  1:2169706/101‑1 (MQ=21)
                                                              tGGCGCGCGCTCTGAGACGAAACGGGG.GGCGCTGCCCCC.CTC.TT.GGCCTTC.GGCCAATTCA.CCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTc                                   >  2:1009588/1‑101 (MQ=21)
                                                              tGGCGCGCGCTCTGAGACGAAACGGGG.GGCGCTGCCCCC.CTC.TT.GGCCTTC.GGCCAATTCA.CCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTc                                   <  1:3248456/101‑1 (MQ=21)
                                                                gcgcgcgcTCTGAGACGAAACGGGG.GGCGCTGCCCCC.CTC.TT.GGCCTTC.GGCCAATTCA.CCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCt                                 <  1:3376059/101‑1 (MQ=21)
                                                                gcgcgcgcTCTGAGACGAAACGGGG.GGCGCTGCCCCC.CTC.TT.GGCCTTC.GGCCAATTCA.CCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCt                                 <  1:3727302/101‑1 (MQ=21)
                                                                 cgcgcgcTCTGAGACGAAACGGGG.GGCGCTGCCCCC.CTC.TT.GGCCTTC.GGCCAATTCA.CCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCtt                                <  1:2512408/101‑1 (MQ=21)
                                                                 cgcgcgcTCTGAGACGAAACGGGG.GGCGCTGCCCCC.CTC.TT.GGCCTTC.GGCCAATTCA.CCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCtt                                <  1:1533020/101‑1 (MQ=21)
                                                                 cgcgcgcTCTGAGACGAAACGGGG.GGCGCTGCCCCC.CTC.TT.GGCCTTC.GGCCAATTCA.CCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCtt                                <  2:3626458/101‑1 (MQ=21)
                                                                  gcgcgcTCTGAGACGAAACGGGG.GGCGCTGCCCCC.CTC.TT.GGCCTTC.GGCCAATTCC.CCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCttt                               <  2:823803/101‑1 (MQ=21)
                                                                  gcgcgcTCTGAGACGAAACGGGG.GGCGCTGCCCCC.CTC.TT.GGCCTTC.GGCCAATTCA.CCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCttt                               <  2:2428260/101‑1 (MQ=21)
                                                                  gcgcgcTCTGAGACGAAACGGGG.GGCGCTGCCCCC.CTC.TT.GGCCTTC.GGCCAATTCA.CCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCttt                               <  2:1830164/101‑1 (MQ=21)
                                                                    gcgcTCTGAGGCGAAACGGGG.GGCGCTGCCCCC.CTC.TT.GGCCTTC.GGCCAATTCA.CCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCa                             <  1:2096435/101‑1 (MQ=21)
                                                                     cgcTCTGAGACGAAACGGGG.GGCGCTGCCCCC.CTC.TT.GGCCTTC.GGCCAATTCA.CCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCAt                            <  1:4145111/101‑1 (MQ=21)
                                                                      gcTCTGAGACGAAACGGGG.GGCGCTGCCCCC.CTC.TT.GGCCTTC.GGCCAATTCA.CCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATc                           <  1:22833/101‑1 (MQ=21)
                                                                      gcTCTGAGACGAAACGGGG.GGCGCTGCCCCC.CTC.TT.GGCCTTC.GGCCAATTCA.CCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATc                           <  1:31132/101‑1 (MQ=21)
                                                                        tctGAGACGAAACGGGG.GGCGCTGCCCCC.CTC.TT.GGCCTTC.GGCCAATTCA.CCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCtt                         <  1:1173596/101‑1 (MQ=21)
                                                                         ctGAGACGAAACGGGG.GGCGCTGCCCCC.CTC.TT.GGCCTTC.GGCCAATTCA.CCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTg                        <  1:569431/101‑1 (MQ=21)
                                                                          tGAGACGAAACGGGG.GGCGCTGCCCCC.CTC.TT.GGCCTTC.GGCCAATTCA.CCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGc                       >  2:2272230/1‑101 (MQ=21)
                                                                           gagaCGAAACGGGG.GGCGCTGCCCCC.CTC.TT.GGCCTTC.GGCCAATTCA.CCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGcc                      >  1:1222553/1‑101 (MQ=21)
                                                                           gagaCGAAACGGGG.GGCGCTGCCCCC.CTC.TT.GGCCTTC.GGCCAATTCA.CCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGcc                      >  2:1826958/1‑101 (MQ=21)
                                                                           gagaCGAAACAGGG.GGCGCTACCCCC.CTC.TT.GGCCTTC.GGCCAATTCC.CCCCCCGGAATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGcc                      <  1:3307486/101‑1 (MQ=21)
                                                                               cGAAACGGGG.GGCGCTGCCCCC.CTC.TT.GGCCTTC.GGCCAATTCA.CCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGaa                  >  2:3077751/1‑101 (MQ=21)
                                                                               cGAAACGGGG.GGCGCTGCCCCC.CTC.TT.GGCCTTC.GGCCAATTCA.CCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGaa                  <  1:3404108/101‑1 (MQ=21)
                                                                                gAAACGGGG.GGCGCTGCTCCC.CTC.TT.GGCCTTC.GGCCAATTCA.CCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGaaa                 <  2:3806296/101‑1 (MQ=21)
                                                                                gAAACGGGG.GGCGCTGCCCCC.CTC.TT.GGCCTTC.GGCCAATTCA.CCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGaaa                 <  2:2702264/101‑1 (MQ=21)
                                                                                gAAACGGGG.GGCGCTGCCCCC.CTC.TT.GGCCTTC.GGCCAATTCA.CCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGaaa                 <  2:1593263/101‑1 (MQ=21)
                                                                                  aacgggg.ggCGCTGCCCCC.CTC.TT.GGCCTTC.GGCCAATTCACCCCCCCGGGAAATTT.GGGCCAGAAAAAGAAAGGCTTTTTCTTTCATCTTTcccccaaaa               >  1:2185736/5‑95 (MQ=2)
                                                                                  aacgggg.ggCGCTGCCCCC.CTC.TT.GGCCTTC.GGCCAATTCA.CCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAata               <  1:3581109/97‑1 (MQ=25)
                                                                                   acgggg.ggCGCTGCCCCC.CTC.TT.GGCCTTC.GGCCAATTCA.CCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAatat              <  2:4156408/98‑1 (MQ=25)
                                                                                    cgggg.ggCGCTGCCCCC.CTC.TT.GGCCTTC.GGCCAATTCA.CCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATc             <  2:1120311/99‑1 (MQ=25)
                                                                                    cgggg.ggCGCTGCCCCC.CTC.TT.GGCCTTC.GGCCAATTCA.CCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATc             <  2:4015211/99‑1 (MQ=25)
                                                                                     gggg.ggCGCTGCCCCC.CTC.TT.GGCCTTC.GGCCAATTCA.CCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATcc            <  2:1571825/100‑1 (MQ=25)
                                                                                     gggg.ggCGCTGCCCCC.CTC.TT.GGCCTTC.GGCCAATTCA.CCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATcc            <  1:949433/100‑1 (MQ=25)
                                                                                       gg.ggCGCTGCCCCCTCTT.GG.C..CTTT.GGCCAATTCA..CCCCTGGGATATTTCGGGCAAGAT...GAA.ATCATGCCCCTTCATCTTGCCGCAAATATCCCgggggga    <  1:3959563/101‑3 (MQ=2)
                                                                                        g.ggCGCTGCCCCCTCTT.GG.C..CTTT.GGCCAATTCA..CCCCTGGGATATTTCGGGCAAGAT...GAA.ATCATGCCCCTTCATCTTGCCGCAAATATCCCggggggag   >  1:601326/1‑98 (MQ=2)
                                                                                        g.ggCGCTG.CCCC.CTC.TT.GGCCTTT.GGCCAATTCA..CCCCTGGGATATTTCGGGCAAGAT...GAA.ATCATGCCCCTTCATCTTGCCGCAAATATCCCggggggag   >  1:1254032/1‑98 (MQ=2)
                                                                                          ggCGCTGCCCCCTCTT.GG.C..CTTT.GGCCAATTCA..CCCCTGGGATATTTCGGGCAAGAT...GAA.ATCATGCCCCTTCATCTTGCCGCAAATATCCCggggggagc  <  1:1048067/101‑5 (MQ=2)
                                                                                          ggCGCTGCCCCCTCTT.GG.C..CTTT.GGCCAATTCA..CCCCTGGGATATTTCGGGCAAGAT...GAA.ATCATGCCCCTTCATCTTGCCGCAAATATCCCggggggagc  <  1:2334609/101‑5 (MQ=2)
                                                                                            cgcTGCCCCC.CTC.TT.GGCCTTC.GGCCAATTCA.CCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCCAAATATCCCggggg      >  2:3022218/1‑101 (MQ=25)
                                                                                             gcTGCCCCC.CTC.TT.GGCCTTC.GGCCAATTCT.CCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCCGGGGGa     <  1:2866857/101‑1 (MQ=37)
                                                                                              ccgCCCCC.CTC.TT.GGCCTTC.GGCCAAGTCC.CCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCCGGGGGAg    <  2:3568447/99‑1 (MQ=37)
                                                                                                      |                                                                                                   
GCGCCTATGCGATGGTCGATGCGATCGACTGGCCCGAGGGGTTCTGCCGCCGCGACATCGGCTGGCGCGCGCTCTGAGACGAAAC.GGG.GGCGCTGCCCCC.CTC.CT.G..CTTC.GGCCAATTCA.CCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCCGGGGGAGCC  >  NC_014034/2566452‑2566644

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