Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_014034 2,566,554 C→T 3.6% intergenic (+26/+168) 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,5540C→T 3.6% 3.6 ‑0.6 ‑0.2111intergenic (+26/+168)purD/sdaAphosphoribosylamine‑‑glycine ligase/L‑serine ammonia‑lyase
Reads supporting (aligned to +/- strand):  new base (1/3):  ref base (18/89):  total (19/92)
Fisher's exact test strand distribution p-value = 5.33e-01
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 7.00e-01

Only 100 of 134 total aligned reads displayed.

CCTATGCGATGGTCGATGCGATCGACTGGCCCGAGGGGTTCTGCCGCCGCGACATCGGCTGGCGCGCGCTCTGAGACGAAAC.GGGGGCGCTGCCCCC.CTC.CT..G..CTTC.GGCC.AATTCA.CCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCCGGGGGAGCC  >  NC_014034/2566455‑2566644
                                                                                                     |                                                                                                  
ccTATGCGATGGTCGATGCGATCGACTGGCCCGAGGGGTTCTGCCGCCGCGACATCGGCTGGCGCGCGCTCTGAGACGAAACGGGGGGCGCTGCCCCc.ctc                                                                                                    <  2:878730/101‑1 (MQ=255)
   aTGCGATGGTCGATGCGATCGACTGGCCCGAGGGGTTCTGCCGCCGCGACATCGGCTGGCGCGCGCTCTGAGACGAAACGGGGGGCGCTGCCCCC.CTC.TT..g                                                                                              <  1:2279289/101‑1 (MQ=255)
    tGCGATGGTCGATGCGATCGACTGGCCCGAGGGGTTCTGCCGCCGCGACATCGGCTGGCGCGCGCTCTGAGACGAAACGGGGGGCGCTGCCCCC.CTC.TT..gg                                                                                             <  1:1918000/101‑2 (MQ=255)
     gCGATGGTCGATGCGATCGACTGGCCCGAGGGGTTCTGCCGCCGCGACATCGGCTGGCGCGCGCTCTGAGACGAAACGGGGGGCGCTGCCCCC.CTC.TT..ggc                                                                                            <  1:3437533/101‑3 (MQ=255)
      cGATGGTCGATGCCATCGACTGGCCCGAGGGGTTCTGCCGCCGCGACATCGGCTGGCGCGCGCTCTGAGACGAAACGGGGGGCGCTGCCCCC.CTC.TT..ggcc                                                                                           <  2:2937542/101‑4 (MQ=38)
       gATGGTCGATGCGATCGACTGGCCCGAGGGGTTCTGCCGCCGCGACATCGGCTGGCGCGCGCTCTGAGACGAAACGGGGGGCGCTGCCCCC.CTC.TT..ggcct                                                                                          <  1:1058341/101‑5 (MQ=39)
        aTGGTCGATGCGATCGACTGGCCCGAGGGGTTCTGCCGCCGCGACATCGGCTGGCGCGCGCTCTGAGACGAAACGGGGGGCGCCGCCCCc.ctc.cttggcct                                                                                           >  1:2144324/1‑95 (MQ=255)
        aTGGTCGATGCGATCGACTGGCCCGAGGGGTTCTGCCGCCGCGACATCGGCTGGCGCGCGCTCTGAGACGAAAC.GGGGGGGGCGGCCCc.ccc.ctttgcctt                                                                                          >  1:3714625/1‑94 (MQ=255)
              gATGCGCTCGGCTGGCCCGAGGGGTTCTGCCGCCGCGACATCGGCTGGCGCGCGCTCTGAGACGAAACGGGGGGCGCTGCCCCC.CTC.TT..GGCCTTC.GGCC.a                                                                                 <  1:705356/101‑1 (MQ=25)
              gATGCGATCGACTGGCCCGAGGGGTTCTGCCGCCGCGACATCGGCTGGCGCGCGCTCTGAGACGAAACGGGGGGCGCTGCCCCC.CTC.TT..GGCCTTC.GGCC.a                                                                                 <  2:2743559/101‑1 (MQ=37)
               aTGCGATCGACTGGCCCGAGGGGTTCTGCCGCCGCGACATCGGCTGGCGCGCGCTCTGAGACGAAACGGGGGGCGCTGCCCCC.CTC.TT..GGCCTTC.GGCC.aa                                                                                >  2:2028863/1‑101 (MQ=37)
                tGCGATCGACTGGCCCGAGGGGTTCTGCCGCCGCGACATCGGCTGGCGCGCGCTCTGAGACGAAACGGGGGGCGCTGCCCCC.CTC.TT..GGCCTTC.GGCC.AAt                                                                               <  2:3129681/101‑1 (MQ=37)
                 gCGATCGACTGGCCCGAGGGGTTCTGCCGCCGCGACATCGGCTGGCGCGCGCTCTGAGACGAAACGGGGGGCGCTGCCCCC.CTC.TT..GGCCTTC.GGCC.AAtt                                                                              <  1:3672165/101‑1 (MQ=37)
                  cGATCGACTGGCCCGAGGGGTTCTGCCGCCGCGACATCGGCTGGCGCGCGCTCTGAGACGAAACGGGGGGCGCTGCCCCC.CTC.TT..GGCCTTC.GGCC.AATTc                                                                             <  1:2286059/101‑1 (MQ=37)
                     tCGACTGGCCCGAGGGGTTCTGCCGCCGCGACATCGGCTGGCGCGCGCTCTGAGACGAAACGGGGGGCGCTGCCCCC.CTC.TT..GGCCTTC.GGCC.AATTCA.cc                                                                         <  1:486815/101‑1 (MQ=37)
                     tCGACTGGCCCGAGGGGTTCTGCCGCCGCGACATCGGCTGGCGCGCGCTCTGAGACGAAACGGGGGGCGCTGCCCCC.CCC.TT..GGCCCTC.GGCC.AATTCC.cc                                                                         >  1:2553452/1‑101 (MQ=25)
                        gcgggCCCGAGGGGTTCTGCCGCCGCGACATCGGCTGGCGCGCGCTCTGAGACGAAACGGGGGGCGCTGCCCCC.CTC.TT..GGCCTTC.GGCC.AATTCA.ccccc                                                                      <  2:2391411/98‑1 (MQ=25)
                        ccTGGCCCGAGGGGTTCTGCCGCCGCGACATCGGCTGGCGCGCGCTCTGAGACGAAACGGGGGGCGCTGCCCCC.CTC.TT..GGCCTTC.GGCC.AATTCA.ccccc                                                                      <  2:203294/100‑1 (MQ=25)
                        aCTGGCCCGAGGGGTTCTGCCGCCGCGACATCGGCTGGCGCGCGCTCTGAGACGAAACGGGGGGCGCTGCCCCC.CTC.TT..GGCCTTC.GGCC.AATTCA.ccccc                                                                      <  1:1187826/101‑1 (MQ=25)
                         cTGGCCCGAGGGGTTCTGCCGCCGCGACATCGGCTGGCGCGCGCTCTGAGACGAAACGGGGGGCGCTGCCCCC.CTC.TT..GGCCTTC.GGCC.AATTCA.cccccc                                                                     <  2:702261/101‑1 (MQ=25)
                         cTGGCCCGAGGGGTTCTGCCGCCGCGACATCGGCTGGCGCGCGCTCTGAGACGAAACGGGGGGCGCTGCCCCC.CTC.TT..GGCCTTC.GGCC.AATTCA.cccccc                                                                     <  2:2895798/101‑1 (MQ=25)
                           ggCCCGAGGGGTTCTGCCGCCGCGACATCGGCTGGCGCGCGCTCTGAGACGAAACGGGGGGCGCTGCCCCC.CTC.TT..GGCCTTC.GGCC.AATTCA.CCCCCCgg                                                                   <  1:503621/101‑1 (MQ=25)
                           ggCCCGAGGGGTTCTGCCGCCGCGACATCGGCTGGCGCGCGCTCTGAGACGAAACGGGGGGCGCTGCCCCC.CTC.TT..GGCCTTC.GGCC.AATTCA.CCCCCCgg                                                                   <  1:1099121/101‑1 (MQ=25)
                             cccGAGGGGTTCTGCCGCCGCGACATCGGCTGGCGCGCGCTCTGAGACGAAACGGGGGGCGCTGCCCCC.CTC.TT..GGCCTTC.GGCC.AATTCA.CCCCCCGGGa                                                                 <  1:1829286/101‑1 (MQ=25)
                              ccGAGGGGTTCTGCCGCCGCGACATCGGCTGGCGCGCGCTCTGAGACGAAACGGGGGGGGCCGCCCCC.CCCTTT..GC.CTTC.GGCC.AAATCA.CCCCCCcggga                                                                >  1:1978096/1‑99 (MQ=255)
                              ccGAGGGGTTCTGCCGCCGCGACATCGGCTGGCGCGCGCTCTGAGACGAAACGGGGGGCGCTGCCCCC.CTC.TT..GGCCTTC.GGCC.AATTCA.CCCCCCGGGat                                                                <  2:2168907/101‑1 (MQ=25)
                              ccGAGGGGTTCTGCCGCCGCGACATCGGCTGGCGCGCGCTCTGAGACGAAACGGGGGGCGCTGCCCCC.CTC.TT..GGCCTTC.GGCC.AATTCA.CCCCCCGGGat                                                                <  2:4180170/101‑1 (MQ=25)
                               cGAGGGGTTCTGCCGCCGCGACATCGGCTGGCGCGCGCTCTGAGACGAAACGGGGGGGGCTGCCCCC.CCC.TT..GGCCTTC.CGCC.AAATCC.CCCCCCCGGgta                                                               >  2:3168341/1‑101 (MQ=255)
                               cGAGGGGTTCTGCCGCCGCGACATCGGCTGGCGCGCGCTCTGAGACGAAACGGGGGGCGCTGCCCCC.CTC.TT..GGCCTTC.GGCC.AATTCA.CCCCCCGGGata                                                               <  2:1018617/101‑1 (MQ=25)
                               cGAGGGGTTCTGCCGCCGCGACATCGGCTGGCGCGCGCTCTGAGACGAAACGGGGGGCGCTGCCCCC.CTC.TT..GGCCTTC.GGCC.AATTCA.CCCCCCGGGata                                                               <  1:1708987/101‑1 (MQ=25)
                                gAGGGGTTCTGCCGCCGCGACATCGGCTGGCGCGCGCTCTGAGACGAAACGGGGGGCGCTGCCCCC.CTC.TT..GGCCTTC.GGCC.AATTCA.CCCCCCGGGatat                                                              <  2:4131168/101‑1 (MQ=25)
                                gAGGGGTTCTGCCGCCGCGACATCGGCTGGCGCGCGCTCTGAGACGAAACGGGGGGCGCTGCCCCC.CTC.TT..GGCCTTC.GGCC.AATTCA.CCCCCCCGGatat                                                              >  1:2119448/1‑101 (MQ=25)
                                 aGGGGTTCTGCCGCCGCGACATCGGCTGGCGCGCGCTCTGAGACGAAACGGGGGGCGCTGCCCCC.CTC.TT..GGCCTTC.GGCC.AATTCA.ACCCCCGGGATAtt                                                             >  2:1052333/1‑101 (MQ=21)
                                  ggggTTCTGCCGCCGCGACATCGGCTGGCGCGCGCTCTGAGACGAAACGGGGGGCGCTGCCCCC.CTC.TT..GGCCTTC.GGCC.AATTCA.CCCCCCGGGATAttt                                                            <  1:1697153/101‑1 (MQ=21)
                                  ggggTTCTGCCGCCGCGACATCGGCTGGCGCGCGCTCTGAGACGAAACGGGGGGCGCTGCCCCC.CTC.TT..GGCCTTC.GGCC.AATTCA.CCCCCCGGGATAttt                                                            <  2:77585/101‑1 (MQ=21)
                                      ttCTGCCGCCGCGACATCGGCTGGCGCGCGCTCTGAGACGAAACGGGGGGCGCTGCCCCC.CTC.TT..GGCCTTC.GGCC.AATTCA.CCCCCCGGGATATTTCggg                                                        <  2:2145273/101‑1 (MQ=21)
                                       tCTGCCGCCGCGACATCGGCTGGCGCGCGCTCTGAGACGAAACGGGGGGCGCTGCCCCC.CTC.TT..GGCCTTC.GGCC.AATTCA.CCCCCCGGGATATTTCGGGc                                                       <  1:1268244/101‑1 (MQ=21)
                                       tCTGCCGCCGCGACATCGGCTGGCGCGCGCTCTGAGACGAAACGGGGGGCGCCGCCCCC.CTC.TTTGG..CCTC.GGCCAAATCCA.CCCCCCcggtacttttggg                                                        >  2:3685661/1‑91 (MQ=255)
                                        cTGCCGCCGCGACATCGGCTGGCGCGCGCTCTGAGACGAAACGGGGGGCGCTGCCCCC.CTC.TT..GGCCTTC.GGCC.AATTCA.CCCCCCCGGCTATTTCGGGCa                                                      >  2:861613/1‑101 (MQ=21)
                                         tGCCGCCGCGACATCGGCTGGCGCGCGCTCTGAGACGAAACGGGGGGCGCTGCCCCC.CCC.TT..GGCCTTC.GGCC.AATTCA.CCCCCCCGGAGAtttttggcca                                                     >  2:3957956/1‑94 (MQ=18)
                                          gccgccGCGACATCGGCTGGCGCGCGCTCTGAGACGAAACGGGGGGCGCTGCCCCC.CTC.TT..GGCCTTC.GGCC.AATTCA.CCCCCCGGGATATTTCGGGCAAg                                                    <  1:45837/101‑1 (MQ=21)
                                            cgccgcGACATCGGCTGGCGCGCGCTCTGAGACGAAACGGGGGGGGCTGGCCCC.CTC.TT..GGCCTTT.CGCC.AATTCA.CCCCCCCGGACATTTCGGGCAAGaa                                                  >  1:3980246/1‑100 (MQ=255)
                                            cgccgcGACATCGGCTGGCGCGCGCTCTGAGACGAAACGGGGGGCGCTGCCCCC.CTC.TT..GGCCTTC.GGCC.AATTCA.CCCCCCGGGATATTTCGGGCAAGAt                                                  <  2:2060072/101‑1 (MQ=21)
                                            cgccgcGACATCGGCTGGCGCGCGCTCTGAGACGAAACGGGGGGCGCTGCCCCC.CTC.TT..GGCCTTC.GGCC.AATTCA.CCCCCCGGGATATTTCGGGCAAGAt                                                  <  2:1620144/101‑1 (MQ=21)
                                            cgccgcGACATCGGCGGGCGCGCGCTCTGAGACGAAACGGGGGGCGCTGCCCCC.CTC.TT..GGCCTTC.GGCC.AATTCA.CCCCCCGGGATATTTCGGGCAAGAt                                                  <  2:523134/101‑1 (MQ=21)
                                             gccgcGACATCGGCTGGCGCGCGCTCTGAGACGAAACGGGGGGCGCTGCCCCC.CTC.TT..GGCCTTC.GGCC.AATTCA.CCCCCCGGGATATTTCGGGCAAGATa                                                 <  2:3108432/101‑1 (MQ=21)
                                               cgcgACATCGGCTGGCGCGCGCTCTGAGACGAAACGGGGGGCGCTGCCCCC.CTC.TT..GGCCTTC.GGCC.CATTCA.CCCCCCCGGATATTTCGGGCAAGATaaa                                               >  1:3906941/1‑101 (MQ=21)
                                                gcgACATCGGCTGGCGCGCGCTCTGAGACGAAACGGGGGGCGCTGCCCCC.CTC.TT..GGCCTTC.GGCC.AATTCA.CCCCCCGGGATATTTCGGGCAAGATAaag                                              <  1:430223/101‑1 (MQ=21)
                                                 cgACATCGGCTGGCGCGCGCTCTGAGACGAAACGGGGGGCGCTGCCCCC.CTC.TT..GGCCTTC.GGCC.AATTCA.CCCCCCGGGATATTTCGGGCAAGATAaaga                                             >  2:3176725/1‑101 (MQ=21)
                                                 cgACATCGGCTGGCGCGCGCTCTGAGACGAAACGGGGGGCGCTGCCCCC.CTC.TT..GGCCTTC.GGCC.AATTCA.CCCCCCGGGATATTTCGGGCAAGATAaaga                                             <  1:1811536/101‑1 (MQ=21)
                                                  gACATCGGCTGGCGCGCGCTCTGAGACGAAACGGGGGGCGCTGCCCCC.CTC.TT..GGCCTTC.GGCC.AATTCA.CCCCCCGGGATATTTCGGGCAAGATAaagaa                                            <  1:19493/101‑1 (MQ=21)
                                                   aCATCGGCTGGCGCTCGCTCTGAGACGAAACGGGGGGCCCCGCCCCC.CTT.TTG.G..CCTCGGGCC.AATTCA.CCCCCCGCGCCGTTTCGGGGGCGATAaaaaag                                           >  1:4115095/1‑101 (MQ=9)
                                                     aTCGGCTGGCGCGCGCTCTGAGACGAAACGGGGGGCGCTGCCCCC.CTC.TT..GGCCTTC.GGCC.CATTCA.CCCCCCCGGAAATTTTGGGGAAGATAAAGAAggg                                         >  2:3039499/1‑100 (MQ=18)
                                                      tCGGCTGGCGCGCGCTCTGAGACGAAACGGGGGGCGCTGCCCCC.CTC.TT..GGCCTTC.GGCC.AATTCA.CCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTc                                        >  1:3848117/1‑101 (MQ=21)
                                                      tCGGCTGGCGCGCGCTCTGAGACGAAACGGGGGGCGCTGCCCCC.CTC.TT..GGCCTTC.GGCC.AATTCA.CCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTc                                        <  1:3471783/101‑1 (MQ=21)
                                                        ggCTGGGGCGCGCTCTGAGACGAAACGGGGGGCGCTGCCCCC.CTC.TT..GGCCTTC.GGCC.AATTCA.CCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCtt                                      <  1:3073594/101‑1 (MQ=21)
                                                         gCTGGCGCGCGCTCTGAGACGAAACGGGGGGCGCTGCCCCC.CTC.TT..GGCCTTC.GGCC.AATTCA.CCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCttt                                     >  1:2196690/1‑101 (MQ=21)
                                                          cTGGCGCGCGCTCTGAGACGAAACGGGGGGCGCTGCCCCC.CTC.TT..GGCCTTC.GGCC.AATTCA.CCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCtttt                                    <  1:2169706/101‑1 (MQ=21)
                                                          cTGGCGCGCGCTCTGAGACGAAACGGGGGGCGCTGCCCCC.CTC.TT..GGCCTTC.GGCC.AATTCA.CCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCtttt                                    <  1:2327547/101‑1 (MQ=21)
                                                           tGGCGCGCGCTCTGAGACGAAACGGGGGGCGCTGCCCCC.CTC.TT..GGCCTTC.GGCC.AATTCA.CCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTc                                   <  1:3248456/101‑1 (MQ=21)
                                                           tGGCGCGCGCTCTGAGACGAAACGGGGGGCGCTGCCCCC.CTC.TT..GGCCTTC.GGCC.AATTCA.CCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTc                                   >  2:1009588/1‑101 (MQ=21)
                                                             gcgcgcgcTCTGAGACGAAACGGGGGGCGCTGCCCCC.CTC.TT..GGCCTTC.GGCC.AATTCA.CCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCt                                 <  1:3376059/101‑1 (MQ=21)
                                                              cgcgcgcTCTGAGACGAAACGGGGGGCGCTGCCCCC.CTC.TT..GGCCTTC.GGCC.AATTCA.CCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCtt                                <  1:1533020/101‑1 (MQ=21)
                                                              cgcgcgcTCTGAGACGAAACGGGGGGCGCTGCCCCC.CTC.TT..GGCCTTC.GGCC.AATTCA.CCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCtt                                <  2:3626458/101‑1 (MQ=21)
                                                               gcgcgcTCTGAGACGAAACGGGGGGCGCTGCCCCC.CTC.TT..GGCCTTC.GGCC.AATTCC.CCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCttt                               <  2:823803/101‑1 (MQ=21)
                                                               gcgcgcTCTGAGACGAAACGGGGGGCGCTGCCCCC.CTC.TT..GGCCTTC.GGCC.AATTCA.CCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCttt                               >  2:3297536/1‑101 (MQ=21)
                                                               gcgcgcTCTGAGACGAAACGGGGGGCGCTGCCCCC.CTC.TT..GGCCTTC.GGCC.AATTCA.CCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCttt                               <  2:2428260/101‑1 (MQ=21)
                                                                cgcgcTCTGAGACGAAACGGGGGGCGCTGCCCCC.CTC.TT..GGCCTTC.GGCC.AATTCA.CCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTc                              <  1:509035/101‑1 (MQ=21)
                                                                cgcgcTCTGAGACGAAACGGGGGGCGCTGCCCCC.CTC.TT..GGCCTTC.GGCC.AATTCA.CCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTc                              <  1:2076249/101‑1 (MQ=21)
                                                                   gcTCTGAGACGAAACGGGGGGCGCTGCCCCC.CTC.TT..GGCCTTC.GGCC.AATTCA.CCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATc                           <  1:902741/101‑1 (MQ=21)
                                                                   gcTCTGAGACGAAACGGGGGGCGCTGCCCCC.CTC.TT..GGCCTTC.GGCC.AATTCA.CCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATc                           <  1:31132/101‑1 (MQ=21)
                                                                    ctctGAGACGAAACGGGGGGCGCTGCCCCC.CTC.TT..GGCCTTC.GGCC.AATTCA.CCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCt                          <  1:2775027/101‑1 (MQ=21)
                                                                     tctGAGACGAAACGGGGGGCGCTGCCCCC.CTC.TT..GGCCTTC.GGCC.AATTCA.CCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCtt                         <  1:1173596/101‑1 (MQ=21)
                                                                     tctGAGACGAAACGGGGGGCGCTGCCCCC.CTC.TT..GGCCTTC.GGCC.AATTCA.CCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCtt                         <  1:298907/101‑1 (MQ=21)
                                                                       tGAGACGAAACGGGGGGCGCTGCCCCC.CTC.TT..GGCCTTC.GGCC.AATTCA.CCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGc                       >  2:2272230/1‑101 (MQ=21)
                                                                        gagaCGAAACGGGGGGCGCTGCCCCC.CTC.TT..GGCCTTC.GGCC.AATTCA.CCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGcc                      >  1:1222553/1‑101 (MQ=21)
                                                                        gagaCGAAACGGGGGGCGCTGCCCCC.CTC.TT..GGCCTTC.GGCC.AATTCA.CCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGcc                      >  2:1826958/1‑101 (MQ=21)
                                                                        gagaCGAAACAGGGGGCGCTACCCCC.CTC.TT..GGCCTTC.GGCC.AATTCC.CCCCCCGGAATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGcc                      <  1:3307486/101‑1 (MQ=21)
                                                                         agaCGAAACGGGGGGCGCTGCCCCC.CTC.TT..GGCCTTC.GGCC.AATTCA.CCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGccc                     <  2:236725/101‑1 (MQ=21)
                                                                           aCGAAACGGGGGGCGCTGCCCCC.CTC.TT..GGCCTTC.GGCC.AATTCA.CCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGa                   >  1:1787489/1‑101 (MQ=21)
                                                                            cGAAACGGGGGGCGCTGCCCCC.CTC.TT..GGCCTTC.GGCC.AATTCA.CCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGaa                  >  2:3077751/1‑101 (MQ=21)
                                                                            cGAAACGGGGGGCGCTGCCCCC.CTC.TT..GGCCTTC.GGCC.AATTCA.CCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGaa                  <  1:3404108/101‑1 (MQ=21)
                                                                             gAAACGGGGGGCGCTGCCCCC.CTC.TT..GGCCTTC.GGCC.AATTCA.CCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGaaa                 <  2:2702264/101‑1 (MQ=21)
                                                                             gAAACGGGGGGCGCTGCCCCC.CTC.TT..GGCCTTC.GGCC.AATTCA.CCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGaaa                 <  2:1593263/101‑1 (MQ=21)
                                                                              aaacggggggCGCTGCCCCC.CTC.TT..GGCCTTC.GGCC.AATTCA.CCCCCCGGGATATTTCGGGCAAGATAAAGAAGGGCTTTTCCTTTCATCTTGcccccaat                >  2:4223765/6‑96 (MQ=18)
                                                                               aacggggggCGCTGCCCCC.CTC.TT..GGCCTTC.GGCC.AATTCACCCCCCCGGGAAATTT.GGGCCAGAAAAAGAAAGGCTTTTTCTTTCATCTTTcccccaaaa               >  1:2185736/5‑95 (MQ=2)
                                                                               aacggggggCGCTGCCCCC.CTC.TT..GGCCTTC.GGCC.AATTCA.CCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAata               <  1:3581109/97‑1 (MQ=25)
                                                                                 cggggggCGCTGCCCCC.CTC.TT..GGCCTTC.GGCC.AATTCA.CCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATc             <  2:1120311/99‑1 (MQ=25)
                                                                                 cggggggCGCTGCCCCC.CTC.TT..GGCCTTC.GGCC.AATTCA.CCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATc             <  2:4015211/99‑1 (MQ=25)
                                                                                 cagggggCGCTG.CCCC.CTC.TT..GGCCTTT.GGCC.AATTCA..CCCCTGGGATATTTCGGGCAAGAT...GAA.ATCATGCCCCTTCATCTTGCCGCAAATATCCCgggg       >  1:3951563/3‑101 (MQ=0)
                                                                                  ggggggCGCTGCCCCC.CTC.TT..GGCCTTC.GGCC.AATTCA.CCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATcc            <  2:1571825/100‑1 (MQ=25)
                                                                                  agggggCGCTG.CCCC.CTC.TT..GGCCTTT.GGCC.AATTCA..CCCCTGGGATATTTCGGGCAAGAT...GAA.ATCATGCCCCTTCATCTTGCCGCAAATATCCCggggg      >  1:4076720/2‑101 (MQ=0)
                                                                                     gggCGCTGCCCCCTCTT.GG..C..CTTT.GGCC.AATTCA..CCCCTGGGATATTTCGGGCAAGAT...GAA.ATCATGCCCCTTCATCTTGCCGCAAATATCCCggggggag   >  1:601326/1‑98 (MQ=2)
                                                                                     gggCGCTG.CCCC.CTC.TT..GGCCTTT.GGCC.AATTCA..CCCCTGGGATATTTCGGGCAAGAT...GAA.ATCATGCCCCTTCATCTTGCCGCAAATATCCCggggggag   >  1:1254032/1‑98 (MQ=2)
                                                                                      ggCGCTGCCCCCTCTT.GG..C..CTTT.GGCC.AATTCA..CCCCTGGGATATTTCGGGCAAGAT...GAA.ATCATGCCCCTTCATCTTGCCGCAAATATCCCggggggagc  <  1:1048067/101‑5 (MQ=2)
                                                                                      ggCGCTGCCCCCTCTT.GG..C..CTTT.GGCC.AATTCA..CCCCTGGGATATTTCGGGCAAGAT...GAA.ATCATGCCCCTTCATCTTGCCGCAAATATCCCggggggagc  <  1:2334609/101‑5 (MQ=2)
                                                                                        cgcTGCCCCC.CTC.TT..GGCCTTC.GGCC.AATTCA.CCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCCAAATATCCCggggg      >  2:3022218/1‑101 (MQ=25)
                                                                                         gcTGCCCCC.CTC.TT..GGCCTTC.GGCC.AATTCT.CCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCCGGGGGa     <  1:2866857/101‑1 (MQ=37)
                                                                                          ccgCCCCC.CTC.TT..GGCCTTC.GGCC.AAGTCC.CCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCCGGGGGAg    <  2:3568447/99‑1 (MQ=37)
                                                                                          cTGCCCCC.CTC.TT..GGCCTTC.GGCC.AATTCA.CCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCCGGGGGAg    <  2:1142532/101‑1 (MQ=37)
                                                                                                     |                                                                                                  
CCTATGCGATGGTCGATGCGATCGACTGGCCCGAGGGGTTCTGCCGCCGCGACATCGGCTGGCGCGCGCTCTGAGACGAAAC.GGGGGCGCTGCCCCC.CTC.CT..G..CTTC.GGCC.AATTCA.CCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCCGGGGGAGCC  >  NC_014034/2566455‑2566644

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