Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_014034 2,566,557 +C 96.2% intergenic (+29/+165) 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,5572.→C96.2% 5.3 ‑0.6 ‑1.2105intergenic (+29/+165)purD/sdaAphosphoribosylamine‑‑glycine ligase/L‑serine ammonia‑lyase
Reads supporting (aligned to +/- strand):  new base (18/83):  ref base (1/3):  total (19/86)
Fisher's exact test strand distribution p-value = 5.56e-01
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 6.77e-02

Only 100 of 129 total aligned reads displayed.

ATGCGATGGTCGATGCGATCGACTGGCCCGAGGGGTTCTGCCGCCGCGACATCGGCTGGCGCGCGCTCTGAGACGAAAC.GGG.GGCGCTGCCCCC.CTC.CT..G..CTTC.GGCC.AATTCACCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCCGGGGGAGCC  >  NC_014034/2566458‑2566644
                                                                                                          |                                                                                          
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)
    gATGGTCGATGCGATCGACTGGCCCGAGGGGTTCTGCCGCCGCGACATCGGCTGGCGCGCGCTCTGAGACGAAACGGGG.GGCGCTGCCCCC.CTC.TT..ggcct                                                                                         <  1:1113287/101‑5 (MQ=39)
        gTCGATGCGATCGACTGGCCCGAGGGGTTCTGCCGCCGCGACATCGGCTGGCGCGCGCTCTGAGACGAAACGGGG.GGCGCTGCCCCC.CTC.TT..ggccttcgg                                                                                     <  2:2429512/101‑9 (MQ=38)
           gATGCGCTCGGCTGGCCCGAGGGGTTCTGCCGCCGCGACATCGGCTGGCGCGCGCTCTGAGACGAAACGGGG.GGCGCTGCCCCC.CTC.TT..GGCCTTC.GGCC.a                                                                                <  1:705356/101‑1 (MQ=25)
             tGCGATCGACTGGCCCGAGGGGTTCTGCCGCCGCGACATCGGCTGGCGCGCGCTCTGAGACGAAACGGGG.GGCGCTGCCCCC.CTC.TT..GGCCTTC.GGCC.AAt                                                                              <  2:3129681/101‑1 (MQ=37)
              gCGATCGACTGGCCCGAGGGGTTCTGCCGCCGCGACATCGGCTGGCGCGCGCTCTGAGACGAAACGGGG.GGCGCTGCCCCC.CTC.TT..GGCCTTC.GGCC.AAtt                                                                             <  1:3672165/101‑1 (MQ=37)
              gCGATCGACTGGCCCGAGGGGTTCTGCCGCCGCGACATCGGCTGGCGCGCGCTCTGAGACGAAACGGGG.GGCGCTGCCCCC.CTC.TT..GGCCTTC.GGCC.AAtt                                                                             <  1:1545824/101‑1 (MQ=37)
              gCGATCGACTGGCCCGAGGGGTTCTGCCGCCGCGACATCGGCTGGCGCGCGCTCTGAGACGAAACGGGG.GGCGCTGCCCCC.CTC.TA..GGCCTTC.GGCC.AAtt                                                                             <  1:3566316/101‑1 (MQ=37)
               cGATCGACTGGCCCGAGGGGTTCTGCCGCCGCGACATCGGCTGGCGCGCGCTCTGAGACGAAACGGGG.GGCGCTGCCCCC.CTC.TT..GGCCTTC.GGCC.AATTc                                                                            <  1:3392828/101‑1 (MQ=37)
               cGATCGACTGGCCCGAGGGGTTCTGCCGCCGCGACATCGGCTGGCGCGCGCTCTGAGACGAAACGGGG.GGCGCTGCCCCC.CTC.TT..GGCCTTC.GGCC.AATTc                                                                            <  1:2286059/101‑1 (MQ=37)
                  tCGACTGGCCCGAGGGGTTCTGCCGCCGCGACATCGGCTGGCGCGCGCTCTGAGACGAAACGGGG.GGCGCTGCCCCC.CTC.TT..GGCCTTC.GGCC.AATTCAcc                                                                         <  1:486815/101‑1 (MQ=37)
                  tCGACTGGCCCGAGGGGTTCTGCCGCCGCGACATCGGCTGGCGCGCGCTCTGAGACGAAACGGGG.GGCGCTGCCCCC.CCC.TT..GGCCCTC.GGCC.AATTCCcc                                                                         >  1:2553452/1‑101 (MQ=25)
                   cGACTGGCCCGAGGGGTTCTGCCGCCGCGACATCGGCTGGCGCGCGCTCTGAGACGAAACGGGG.GGCGCTGCCCCC.CTC.TT..GGCCTTC.GGCC.AATTCAccc                                                                        <  2:659481/101‑1 (MQ=37)
                     ccTGGCCCGAGGGGTTCTGCCGCCGCGACATCGGCTGGCGCGCGCTCTGAGACGAAACGGGG.GGCGCTGCCCCC.CTC.TT..GGCCTTC.GGCC.AATTCAccccc                                                                      <  2:203294/100‑1 (MQ=25)
                     aCTGGCCCGAGGGGTTCTGCCGCCGCGACATCGGCTGGCGCGCGCTCTGAGACGAAACGGGG.GGCGCTGCCCCC.CTC.TT..GGCCTTC.GGCC.AATTCAccccc                                                                      <  1:1187826/101‑1 (MQ=25)
                      cTGGCCCGAGGGGTTCTGCCGCCGCGACATCGGCTGGCGCGCGCTCTGAGACGAAACGGGG.GGCGCTGCCCCC.CTC.TT..GGCCTTC.GGCC.AATTCAcccccc                                                                     <  2:702261/101‑1 (MQ=25)
                      cTGGCCCGAGGGGTTCTGCCGCCGCGACATCGGCTGGCGCGCGCTCTGAGACGAAACGGGG.GGCGCTGCCCCC.CTC.TT..GGCCTTC.GGCC.AATTCAcccccc                                                                     <  2:2895798/101‑1 (MQ=25)
                       tGGCCCGAGGGGTTCTGCCGCCGCGACATCGGCTGGCGCGCGCTCTGAGACGAAACGGGG.GGCGCTGCCCCC.CTC.TT..GGCCTTC.GGCC.AATTCACCCCCCg                                                                    <  1:2181231/101‑1 (MQ=25)
                        ggCCCGAGGGGTTCTGCCGCCGCGACATCGGCTGGCGCGCGCTCTGAGACGAAACGGGG.GGCGCTGCCCCC.CTC.TT..GGCCTTC.GGCC.AATTCACCCCCCgg                                                                   <  1:503621/101‑1 (MQ=25)
                        ggCCCGAGGGGTTCTGCCGCCGCGACATCGGCTGGCGCGCGCTCTGAGACGAAACGGGG.GGCGCTGCCCCC.CTC.TT..GGCCTTC.GGCC.AATTCACCCCCCgg                                                                   <  1:2944444/101‑1 (MQ=25)
                        ggCCCGAGGGGTTCTGCCGCCGCGACATCGGCTGGCGCGCGCTCTGAGACGAAACGGGG.GGCGCTGCCCCC.CTC.TT..GGCCTTC.GGCC.AATTCACCCCCCgg                                                                   <  1:1099121/101‑1 (MQ=25)
                         gCCCGAGGGGTTCTGCCGCCGCGACATCGGCTGGCGCGCGCTCTGAGACGAAAC.GGGTGGGGCTGCCCCC.CTC.TT..GGCCCTC.CGCC.AATTCACCCCCCcgg                                                                  >  2:3490303/1‑101 (MQ=25)
                          cccGAGGGGTTCTGCCGCCGCGACATCGGCTGGCGCGCGCTCTGAGACGAAACGGGG.GGCGCTGCCCCC.CTC.TT..GGCCTTC.GGCC.AATTCACCCCCCGGGa                                                                 <  1:1829286/101‑1 (MQ=25)
                           ccGAGGGGTTCTGCCGCCGCGACATCGGCTGGCGCGCGCTCTGAGACGAAACGGGG.GGGGCCGCCCCC.CCCTTT..GC.CTTC.GGCC.AAATCACCCCCCcggga                                                                >  1:1978096/1‑99 (MQ=255)
                           ccGAGGGGTTCTGCCGCCGCGACATCGGCTGGCGCGCGCTCTGAGACGAAACGGGG.GGCGCTGCCCCC.CTC.TT..GGCCTTC.GGCC.AATTCACCCCCCGGGat                                                                <  2:4180170/101‑1 (MQ=25)
                           ccGAGGGGTTCTGCCGCCGCGACATCGGCTGGCGCGCGCTCTGAGACGAAACGGGG.GGCGCTGCCCCC.CTC.TT..GGCCTTC.GGCC.AATTCACCCCCCGGGat                                                                <  2:2168907/101‑1 (MQ=25)
                           ccGAGGGGTTCTGCCGCCGCGACATCGGCTGGCGCGCGCTCTGAGACGAAACGGGG.GGCGCTGCCCCC.CTC.TT..GGCCTTC.GGCC.AATTCACCCCCCGGGat                                                                <  1:3194193/101‑1 (MQ=25)
                            cGAGGGGTTCTGCCGCCGCGACATCGGCTGGCGCGCGCTCTGAGACGAAACGGGG.GGGGCTGCCCCC.CCC.TT..GGCCTTC.CGCC.AAATCCCCCCCCCGGgta                                                               >  2:3168341/1‑101 (MQ=255)
                            cGAGGGGTTCTGCCGCCGCGACATCGGCTGGCGCGCGCTCTGAGACGAAACGGGG.GGCGCTGCCCCC.CTC.TT..GGCCTTC.GGCC.AATTCACCCCCCGGGata                                                               <  2:1018617/101‑1 (MQ=25)
                             gAGGGGTTCTGCCGCCGCGACATCGGCTGGCGCGCGCTCTGAGACGAAACGGGG.GGCGCTGCCCCC.CTC.TT..GGCCTTC.GGCC.AATTCACCCCCCGGGatat                                                              <  2:4131168/101‑1 (MQ=25)
                             gAGGGGTTCTGCCGCCGCGACATCGGCTGGCGCGCGCTCTGAGACGAAACGGGG.GGCGCTGCCCCC.CTC.TT..GGCCTTC.GGCC.AATTCACCCCCCCGGatat                                                              >  1:2119448/1‑101 (MQ=25)
                              aGGGGTTCTGCCGCCGCGACATCGGCTGGCGCGCGCTCTGAGACGAAACGGGG.GGCGCTGCCCCC.CTC.TT..GGCCTTC.GGCC.AATTCAACCCCCGGGATAtt                                                             >  2:1052333/1‑101 (MQ=21)
                               ggggTTCTGCCGCCGCGACATCGGCTGGCGCGCGCTCTGAGACGAAACGGGG.GGCGCTGCCCCC.CTC.TT..GGCCTTC.GGCC.AATTCACCCCCCGGGATAttt                                                            <  1:1697153/101‑1 (MQ=21)
                                gggTTCTGCCGCCGCGACATCGGCTGGCGCGCGCTCTGAGACGAAACGGGG.GGGGGTGCCCCC.CCC.TTGGG..CTTC.GGCC.AATTCCCCCCCCGGGATATTTc                                                           >  1:1249184/1‑101 (MQ=25)
                                    tCTGCCGCCGCGACATCGGCTGGCGCGCGCTCTGAGACGAAACGGGG.GGCGCTGCCCCC.CTC.TT..GGCCTTC.GGCC.CATTAAACCCCCCGGATATTTTgggg                                                       >  2:1592017/1‑100 (MQ=21)
                                    tCTGCCGCCGCGACATCGGCTGGCGCGCGCTCTGAGACGAAACGGGG.GGCGCTGCCCCC.CTC.TT..GGCCTTC.GGCC.AATTCACCCCCCGGGATATTTCGGGc                                                       <  2:229061/101‑1 (MQ=21)
                                    tCTGCCGCCGCGACATCGGCTGGCGCGCGCTCTGAGACGAAACGGGG.GGCGCTGCCCCC.CTC.TT..GGCCTTC.GGCC.AATTCACCCCCCGGGATATTTCGGGc                                                       <  1:1268244/101‑1 (MQ=21)
                                    tCTGCCGCCGCGACATCGGCTGGCGCGCGCTCTGAGACGAAACGGGG.GGCGCCGCCCCC.CTC.TTTGG..CCTC.GGCCAAATCCACCCCCCcggtacttttggg                                                        >  2:3685661/1‑91 (MQ=255)
                                     cTGCCGCCGCGACATCGGCTGGCGCGCGCTCTGAGACGAAACGGGG.GGCGCTGCCCCC.CTC.TT..GGCCTTC.GGCC.AATTCACCCCCCGGGATATTTCGGGCa                                                      <  2:2900729/101‑1 (MQ=21)
                                     cTGCCGCCGCGACATCGGCTGGCGCGCGCTCTGAGACGAAACGGGG.GGCGCTGCCCCC.CTC.TT..GGCCTTC.GGCC.AATTCACCCCCCCGGCTATTTCGGGCa                                                      >  2:861613/1‑101 (MQ=21)
                                      tGCCGCCGCGACATCGGCTGGCGCGCGCTCTGAGACGAAGCGGGG.GGCGCTGCCCCC.CTC.TT..GGCCTTC.GGCC.AATTCACCCCCCGGGATATTTCGGGCaa                                                     <  1:473092/101‑1 (MQ=21)
                                         cgccgcGACATCGGCTGGCGCGCGCTCTGAGACGAAACGGGG.GGCGCTGCCCCC.CTC.TT..GGCCTTC.GGCC.AATTCACCCCCCGGGATATTTCGGGCAAGAt                                                  <  2:1620144/101‑1 (MQ=21)
                                         cgccgcGACATCGGCTGGCGCGCGCTCTGAGACGAAACGGGG.GGCGCTGCCCCC.CTC.TT..GGCCTTC.GGCC.AATTCACCCCCCGGGATATTTCGGGCAAGAt                                                  <  2:2060072/101‑1 (MQ=21)
                                         cgccgcGACATCGGCGGGCGCGCGCTCTGAGACGAAACGGGG.GGCGCTGCCCCC.CTC.TT..GGCCTTC.GGCC.AATTCACCCCCCGGGATATTTCGGGCAAGAt                                                  <  2:523134/101‑1 (MQ=21)
                                          gccgcGACATCGGCTGGCGCGCGCTCTGAGACGAAACGGGG.GGCGCTGCCCCC.CTC.TT..GGCCTTC.GGCC.AATTCACCCCCCGGGATATTTCGGGCAAGATa                                                 <  2:3108432/101‑1 (MQ=21)
                                          gccgcGACATCGGCTGGCGCGCGCTCTGAGACGAAACGGGG.GGCGCTGCCCCC.CTC.TT..GGCCTTC.GGCC.AATTCACCCCCCGGGATATTTCGGGCAAGATa                                                 <  2:2066625/101‑1 (MQ=21)
                                            cgcgACATCGGCTGGCGCGCGCTCTGAGACGAAACGGGG.GGCGCTGCCCCC.CTC.TT..GGCCTTC.GGCC.CATTCACCCCCCCGGATATTTCGGGCAAGATaaa                                               >  1:3906941/1‑101 (MQ=21)
                                             gcgACATCGGCTGGCGCGCGCTCTGAGACGAAACGGGG.GGCGCTGCCCCC.CTC.TT..GGCCTTC.GGCC.AATTCACCCCCCGGGATATTTCGGGCAAGATAaag                                              <  1:430223/101‑1 (MQ=21)
                                              cgACATCGGCTGGCGCGCGCTCTGAGACGAAACGGGG.GGCGCTGCCCCC.CTC.TT..GGCCTTC.GGCC.AATTCACCCCCCGGGATATTTCGGGCAAGATAaaga                                             <  1:1811536/101‑1 (MQ=21)
                                              cgACATCGGCTGGCGCGCGCTCTGAGACGAAACGGGG.GGCGCTGCCCCC.CTC.TT..GGCCTTC.GGCC.AATTCACCCCCCGGGATATTTCGGGCAAGATAaaga                                             >  2:3176725/1‑101 (MQ=21)
                                               gACATCGGCTGGCGCGCGCTCTGAGACGAAACGGGG.GGCGCTGCCCCC.CTC.TT..GGCCTTC.GGCC.AATTCACCCCCCGGGATATTTCGGGCAAGATAaagaa                                            <  1:19493/101‑1 (MQ=21)
                                                aCATCGGCTGGCGCTCGCTCTGAGACGAAACGGGG.GGCCCCGCCCCC.CTT.TTG.G..CCTCGGGCC.AATTCACCCCCCGCGCCGTTTCGGGGGCGATAaaaaag                                           >  1:4115095/1‑101 (MQ=9)
                                                     ggCTGGGGCGCGCTCTGAGACGAAACGGGG.GGCGCTGCCCCC.CTC.TT..GGCCTTC.GGCC.AATTCACCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCtt                                      <  1:3073594/101‑1 (MQ=21)
                                                     ggCTGGCGCGCCCTCTGAGACGAAACGGGG.GGCGCTGCCCCC.CTC.TT..GGCCTTC.GGCC.AATTCACCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCtt                                      <  1:1943026/101‑1 (MQ=21)
                                                      gCTGGCGCGCGCTCTGAGACGAAACGGGG.GGCGCTGCCCCC.CTC.TT..GGCCTTC.GGCC.AATTCACCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCttt                                     >  1:2196690/1‑101 (MQ=21)
                                                       cTGGCGCGCGCTCTGAGACGAAACGGGG.GGCGCTGCCCCC.CTC.TT..GGCCTTC.GGCC.AATTCACCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCtttt                                    <  1:2169706/101‑1 (MQ=21)
                                                        tGGCGCGCGCTCTGAGACGAAACGGGG.GGCGCTGCCCCC.CTC.TT..GGCCTTC.GGCC.AATTCACCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTc                                   <  1:3248456/101‑1 (MQ=21)
                                                        tGGCGCGCGCTCTGAGACGAAACGGGG.GGCGCTGCCCCC.CTC.TT..GGCCTTC.GGCC.AATTCACCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTc                                   >  2:1009588/1‑101 (MQ=21)
                                                          gcgcgcgcTCTGAGACGAAACGGGG.GGCGCTGCCCCC.CTC.TT..GGCCTTC.GGCC.AATTCACCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCt                                 <  1:3727302/101‑1 (MQ=21)
                                                           cgcgcgcTCTGAGACGAAACGGGG.GGCGCTGCCCCC.CTC.TT..GGCCTTC.GGCC.AATTCACCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCtt                                <  1:2512408/101‑1 (MQ=21)
                                                           cgcgcgcTCTGAGACGAAACGGGG.GGCGCTGCCCCC.CTC.TT..GGCCTTC.GGCC.AATTCACCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCtt                                <  2:3626458/101‑1 (MQ=21)
                                                           cgcgcgcTCTGAGACGAAACGGGG.GGCGCTGCCCCC.CTC.TT..GGCCTTC.GGCC.AATTCACCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCtt                                <  1:1533020/101‑1 (MQ=21)
                                                            gcgcgcTCTGAGACGAAACGGGG.GGCGCTGCCCCC.CTC.TT..GGCCTTC.GGCC.AATTCCCCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCttt                               <  2:823803/101‑1 (MQ=21)
                                                            gcgcgcTCTGAGACGAAACGGGG.GGCGCTGCCCCC.CTC.TT..GGCCTTC.GGCC.AATTCACCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCttt                               >  2:3297536/1‑101 (MQ=21)
                                                            gcgcgcTCTGAGACGAAACGGGG.GGCGCTGCCCCC.CTC.TT..GGCCTTC.GGCC.AATTCACCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCttt                               <  2:1830164/101‑1 (MQ=21)
                                                            gcgcgcTCTGAGACGAAACGGGG.GGCGCTGCCCCC.CTC.TT..GGCCTTC.GGCC.AATTCACCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCttt                               <  2:2428260/101‑1 (MQ=21)
                                                             cgcgcTCTGAGACGAAACGGGG.GGCGCTGCCCCC.CTC.TT..GGCCTTC.GGCC.AATTCACCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTc                              <  1:2076249/101‑1 (MQ=21)
                                                               cgcTCTGAGACGAAACGGGG.GGCGCTGCCCCC.CTC.TT..GGCCTTC.GGCC.AATTCACCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCAt                            <  1:4145111/101‑1 (MQ=21)
                                                                gcTCTGAGACGAAACGGGG.GGCGCTGCCCCC.CTC.TT..GGCCTTC.GGCC.AATTCACCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATc                           <  1:902741/101‑1 (MQ=21)
                                                                gcTCTGAGACGAAACGGGG.GGCGCTGCCCCC.CTC.TT..GGCCTTC.GGCC.AATTCACCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATc                           <  1:22833/101‑1 (MQ=21)
                                                                  tctGAGACGAAACGGGG.GGCGCTGCCCCC.CTC.TT..GGCCTTC.GGCC.AATTCACCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCtt                         <  1:298907/101‑1 (MQ=21)
                                                                   ctGAGACGAAACGGGG.GGCGCTGCCCCC.CTC.TT..GGCCTTC.GGCC.AATTCACCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTg                        <  1:569431/101‑1 (MQ=21)
                                                                    tGAGACGAAACGGGG.GGCGCTGCCCCC.CTC.TT..GGCCTTC.GGCC.AATTCACCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTGTTTTTTTTTCTTCTTGc                       >  2:853783/1‑101 (MQ=18)
                                                                    tGAGACGAAACGGGG.GGCGCTGCCCCC.CTC.TT..GGCCTTC.GGCC.AATTCACCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGc                       >  2:2272230/1‑101 (MQ=21)
                                                                     gagaCGAAACGGGG.GGCGCTGCCCCC.CTC.TT..GGCCTTC.GGCC.AATTCACCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGcc                      >  2:1826958/1‑101 (MQ=21)
                                                                     gagaCGAAACGGGG.GGCGCTGCCCCC.CTC.TT..GGCCTTC.GGCC.AATTCACCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGcc                      >  1:1222553/1‑101 (MQ=21)
                                                                     gagaCGAAACAGGG.GGCGCTACCCCC.CTC.TT..GGCCTTC.GGCC.AATTCCCCCCCCGGAATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGcc                      <  1:3307486/101‑1 (MQ=21)
                                                                      agaCGAAACGGGG.GGCGCTGCCCCC.CTC.TT..GGCCTTC.GGCC.AATTCACCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGccc                     <  2:236725/101‑1 (MQ=21)
                                                                        aCGAAACGGGG.GGCGCTGCCCCC.CTC.TT..GGCCTTC.GGCC.AATTCACCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGa                   >  1:1787489/1‑101 (MQ=21)
                                                                         cGAAACGGGG.GGCGCTGCCCCC.CTC.TT..GGCCTTC.GGCC.AATTCACCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGaa                  >  2:3077751/1‑101 (MQ=21)
                                                                         cGAAACGGGG.GGCGCTGCCCCC.CTC.TT..GGCCTTC.GGCC.AATTCACCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGaa                  <  1:3404108/101‑1 (MQ=21)
                                                                          gAAACGGGG.GGCGCTGCTCCC.CTC.TT..GGCCTTC.GGCC.AATTCACCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGaaa                 <  2:3806296/101‑1 (MQ=21)
                                                                          gAAACGGGG.GGCGCTGCCCCC.CTC.TT..GGCCTTC.GGCC.AATTCACCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGaaa                 <  2:2702264/101‑1 (MQ=21)
                                                                            aacgggg.ggCGCTGCCCCC.CTC.TT..GGCCTTC.GGCC.AATTCACCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAata               <  1:3581109/97‑1 (MQ=25)
                                                                             acgggg.ggCGCTGCCCCC.CTC.TT..GGCCTTC.GGCC.AATTCACCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAatat              <  2:4156408/98‑1 (MQ=25)
                                                                              cgggg.ggCGCTGCCCCC.CTC.TT..GGCCTTC.GGCC.AATTCACCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATc             <  2:1120311/99‑1 (MQ=25)
                                                                              cgggg.ggCGCTGCCCCC.CTC.TT..GGCCTTC.GGCC.AATTCACCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATc             <  2:4015211/99‑1 (MQ=25)
                                                                              caggg.ggCGCTG.CCCC.CTC.TT..GGCCTTT.GGCC.AATTCA.CCCCTGGGATATTTCGGGCAAGAT...GAA.ATCATGCCCCTTCATCTTGCCGCAAATATCCCgggg       >  1:3951563/3‑101 (MQ=0)
                                                                               gggg.ggCGCTGCCCCC.CTC.TT..GGCCTTC.GGCC.AATTCACCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATcc            <  2:1571825/100‑1 (MQ=25)
                                                                               aggg.ggCGCTG.CCCC.CTC.TT..GGCCTTT.GGCC.AATTCA.CCCCTGGGATATTTCGGGCAAGAT...GAA.ATCATGCCCCTTCATCTTGCCGCAAATATCCCggggg      >  1:4076720/2‑101 (MQ=0)
                                                                                 gg.ggCGCTGCCCCCTCTT.GG..C..CTTT.GGCC.AATTCA.CCCCTGGGATATTTCGGGCAAGAT...GAA.ATCATGCCCCTTCATCTTGCCGCAAATATCCCgggggga    <  1:3959563/101‑3 (MQ=2)
                                                                                  g.ggCGCTGCCCCCTCTT.GG..C..CTTT.GGCC.AATTCA.CCCCTGGGATATTTCGGGCAAGAT...GAA.ATCATGCCCCTTCATCTTGCCGCAAATATCCCggggggag   >  1:601326/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.AATTCACCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCCAAATATCCCggggg      >  2:3022218/1‑101 (MQ=25)
                                                                                       gcTGCCCCC.CTC.TT..GGCCTTC.GGCC.AATTCTCCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCCGGGGGa     <  1:2866857/101‑1 (MQ=37)
                                                                                        cTGCCCCC.CTC.TT..GGCCTTC.GGCC.AATTCACCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCCGGGGGAg    <  2:1142532/101‑1 (MQ=37)
                                                                                                          |                                                                                          
ATGCGATGGTCGATGCGATCGACTGGCCCGAGGGGTTCTGCCGCCGCGACATCGGCTGGCGCGCGCTCTGAGACGAAAC.GGG.GGCGCTGCCCCC.CTC.CT..G..CTTC.GGCC.AATTCACCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCCGGGGGAGCC  >  NC_014034/2566458‑2566644

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