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

Only 100 of 132 total aligned reads displayed.

TATGCGATGGTCGATGCGATCGACTGGCCCGAGGGGTTCTGCCGCCGCGACATCGGCTGGCGCGCGCTCTGAGACGAAAC.GGG.GGCGCTGCCCCC.CTC.CT..G..CTTC.GGCC.AATTCACCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCCGGGGGAGCC  >  NC_014034/2566457‑2566644
                                                                                                       |                                                                                              
tATGCGATGGTCGATGCGATCGACTGGCCCGAGGGGTTCTGCCGCCGCGACATCGGCTGGCGCGCGCTCTGAGACGAAACGGGG.GGGGCTGCCCCc.ccc.ct                                                                                                >  1:1789161/1‑101 (MQ=255)
 aTGCGATGGTCGATGCGATCGACTGGCCCGAGGGGTTCTGCCGCCGCGACATCGGCTGGCGCGCGCTCTGAGACGAAACGGGG.GGCGCTGCCCCC.CTC.TT..g                                                                                             <  1:2279289/101‑1 (MQ=255)
  tGCGATGGTCGATGCGATCGACTGGCCCGAGGGGTTCTGCCGCCGCGACATCGGCTGGCGCGCGCTCTGAGACGAAACGGGG.GGCGCTGCCCCC.CTC.TT..gg                                                                                            <  1:1918000/101‑2 (MQ=255)
     gATGGTCGATGCGATCGACTGGCCCGAGGGGTTCTGCCGCCGCGACATCGGCTGGCGCGCGCTCTGAGACGAAACGGGG.GGCGCTGCCCCC.CTC.TT..ggcct                                                                                         <  1:1058341/101‑5 (MQ=39)
     gATGGTCGATGCGATCGACTGGCCCGAGGGGTTCTGCCGCCGCGACATCGGCTGGCGCGCGCTCTGAGACGAAACGGGG.GGCGCTGCCCCC.CTC.TT..ggcct                                                                                         <  1:1113287/101‑5 (MQ=39)
      aTGGTCGATGCGATCGACTGGCCCGAGGGGTTCTGCCGCCGCGACATCGGCTGGCGCGCGCTCTGAGACGAAACGGGG.GGCGCCGCCCCc.ctc.cttggcct                                                                                          >  1:2144324/1‑95 (MQ=255)
      aTGGTCGATGCGATCGACTGGCCCGAGGGGTTCTGCCGCCGCGACATCGGCTGGCGCGCGCTCTGAGACGAAAC.GGG.GGGGGCGGCCCc.ccc.ctttgcctt                                                                                         >  1:3714625/1‑94 (MQ=255)
            gATGCGCTCGGCTGGCCCGAGGGGTTCTGCCGCCGCGACATCGGCTGGCGCGCGCTCTGAGACGAAACGGGG.GGCGCTGCCCCC.CTC.TT..GGCCTTC.GGCC.a                                                                                <  1:705356/101‑1 (MQ=25)
            gATGCGATCGACTGGCCCGAGGGGTTCTGCCGCCGCGACATCGGCTGGCGCGCGCTCTGAGACGAAACGGGG.GGCGCTGCCCCC.CTC.TT..GGCCTTC.GGCC.a                                                                                <  2:2743559/101‑1 (MQ=37)
             aTGCGATCGACTGGCCCGAGGGGTTCTGCCGCCGCGACATCGGCTGGCGCGCGCTCTGAGACGAAACGGGG.GGCGCTGCCCCC.CTC.TT..GGCCTTC.GGCC.aa                                                                               >  2:2028863/1‑101 (MQ=37)
              tGCGATCGACTGGCCCGAGGGGTTCTGCCGCCGCGACATCGGCTGGCGCGCGCTCTGAGACGAAACGGGG.GGCGCTGCCCCC.CTC.TT..GGCCTTC.GGCC.AAt                                                                              <  2:3129681/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)
                   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)
                      gcgggCCCGAGGGGTTCTGCCGCCGCGACATCGGCTGGCGCGCGCTCTGAGACGAAACGGGG.GGCGCTGCCCCC.CTC.TT..GGCCTTC.GGCC.AATTCAccccc                                                                      <  2:2391411/98‑1 (MQ=25)
                      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:1099121/101‑1 (MQ=25)
                         ggCCCGAGGGGTTCTGCCGCCGCGACATCGGCTGGCGCGCGCTCTGAGACGAAACGGGG.GGCGCTGCCCCC.CTC.TT..GGCCTTC.GGCC.AATTCACCCCCCgg                                                                   <  1:2944444/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:2168907/101‑1 (MQ=25)
                            ccGAGGGGTTCTGCCGCCGCGACATCGGCTGGCGCGCGCTCTGAGACGAAACGGGG.GGCGCTGCCCCC.CTC.TT..GGCCTTC.GGCC.AATTCACCCCCCGGGat                                                                <  1:3194193/101‑1 (MQ=25)
                            ccGAGGGGTTCTGCCGCCGCGACATCGGCTGGCGCGCGCTCTGAGACGAAACGGGG.GGCGCTGCCCCC.CTC.TT..GGCCTTC.GGCC.AATTCACCCCCCGGGat                                                                <  2:4180170/101‑1 (MQ=25)
                             cGAGGGGTTCTGCCGCCGCGACATCGGCTGGCGCGCGCTCTGAGACGAAACGGGG.GGGGCTGCCCCC.CCC.TT..GGCCTTC.CGCC.AAATCCCCCCCCCGGgta                                                               >  2:3168341/1‑101 (MQ=255)
                             cGAGGGGTTCTGCCGCCGCGACATCGGCTGGCGCGCGCTCTGAGACGAAACGGGG.GGCGCTGCCCCC.CTC.TT..GGCCTTC.GGCC.AATTCACCCCCCGGGata                                                               <  1:1708987/101‑1 (MQ=25)
                             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)
                               aGGGGTTCTGCCGCCGCGACATCGGCTGGCGCGCGCTCTGAGACGAAACGGGG.GGCGCTGCCCCC.CTC.TT..GGCCTTC.GGCC.AATTCAACCCCCGGGATAtt                                                             >  2:1052333/1‑101 (MQ=21)
                                 gggTTCTGCCGCCGCGACATCGGCTGGCGCGCGCTCTGAGACGAAACGGGG.GGGGGTGCCCCC.CCC.TTGGG..CTTC.GGCC.AATTCCCCCCCCGGGATATTTc                                                           >  1:1249184/1‑101 (MQ=25)
                                    ttCTGCCGCCGCGACATCGGCTGGCGCGCGCTCTGAGACGAAACGGGG.GGCGCTGCCCCC.CTC.TT..GGCCTTC.GGCC.AATTCACCCCCCGGGATATTTCggg                                                        <  2:2145273/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)
                                       tGCCGCCGCGACATCGGCTGGCGCGCGCTCTGAGACGAAACGGGG.GGCGCTGCCCCC.CCC.TT..GGCCTTC.GGCC.AATTCACCCCCCCGGAGAtttttggcca                                                     >  2:3957956/1‑94 (MQ=18)
                                          cgccgcGACATCGGCTGGCGCGCGCTCTGAGACGAAACGGGG.GGGGCTGGCCCC.CTC.TT..GGCCTTT.CGCC.AATTCACCCCCCCGGACATTTCGGGCAAGaa                                                  >  1:3980246/1‑100 (MQ=255)
                                          cgccgcGACATCGGCTGGCGCGCGCTCTGAGACGAAACGGGG.GGCGCTGCCCCC.CTC.TT..GGCCTTC.GGCC.AATTCACCCCCCGGGATATTTCGGGCAAGAt                                                  <  2:2060072/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:3164/101‑1 (MQ=21)
                                          cgccgcGACATCGGCGGGCGCGCGCTCTGAGACGAAACGGGG.GGCGCTGCCCCC.CTC.TT..GGCCTTC.GGCC.AATTCACCCCCCGGGATATTTCGGGCAAGAt                                                  <  2:523134/101‑1 (MQ=21)
                                              gcgACATCGGCTGGCGCGCGCTCTGAGACGAAACGGGG.GGCGCTGCCCCC.CTC.TT..GGCCTTC.GGCC.AATTCACCCCCCGGGATATTTCGGGCAAGATAaag                                              <  1:430223/101‑1 (MQ=21)
                                                 aCATCGGCTGGCGCTCGCTCTGAGACGAAACGGGG.GGCCCCGCCCCC.CTT.TTG.G..CCTCGGGCC.AATTCACCCCCCGCGCCGTTTCGGGGGCGATAaaaaag                                           >  1:4115095/1‑101 (MQ=9)
                                                   aTCGGCTGGCGCGCGCTCTGAGACGAAACGGGG.GGCGCTGCCCCC.CTC.TT..GGCCTTC.GGCC.CATTCACCCCCCCGGAAATTTTGGGGAAGATAAAGAAggg                                         >  2:3039499/1‑100 (MQ=18)
                                                    tCGGCTGGCGCGCGCTCTGAGACGAAACGGGG.GGCGCTGCCCCC.CTC.TT..GGCCTTC.GGCC.AATTCACCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTc                                        >  1:3848117/1‑101 (MQ=21)
                                                    tCGGCTGGCGCGCGCTCTGAGACGAAACGGGG.GGCGCTGCCCCC.CTC.TT..GGCCTTC.GGCC.AATTCACCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTc                                        <  1:3471783/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:2327547/101‑1 (MQ=21)
                                                        cTGGCGCGCGCTCTGAGACGAAACGGGG.GGCGCTGCCCCC.CTC.TT..GGCCTTC.GGCC.AATTCACCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCtttt                                    <  1:2169706/101‑1 (MQ=21)
                                                         tGGCGCGCGCTCTGAGACGAAACGGGG.GGCGCTGCCCCC.CTC.TT..GGCCTTC.GGCC.AATTCACCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTc                                   >  2:1009588/1‑101 (MQ=21)
                                                         tGGCGCGCGCTCTGAGACGAAACGGGG.GGCGCTGCCCCC.CTC.TT..GGCCTTC.GGCC.AATTCACCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTc                                   <  1:3248456/101‑1 (MQ=21)
                                                           gcgcgcgcTCTGAGACGAAACGGGG.GGCGCTGCCCCC.CTC.TT..GGCCTTC.GGCC.AATTCACCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCt                                 <  1:3376059/101‑1 (MQ=21)
                                                           gcgcgcgcTCTGAGACGAAACGGGG.GGCGCTGCCCCC.CTC.TT..GGCCTTC.GGCC.AATTCACCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCt                                 <  1:3727302/101‑1 (MQ=21)
                                                            cgcgcgcTCTGAGACGAAACGGGG.GGCGCTGCCCCC.CTC.TT..GGCCTTC.GGCC.AATTCACCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCtt                                <  1:1533020/101‑1 (MQ=21)
                                                            cgcgcgcTCTGAGACGAAACGGGG.GGCGCTGCCCCC.CTC.TT..GGCCTTC.GGCC.AATTCACCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCtt                                <  1:2512408/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:2428260/101‑1 (MQ=21)
                                                              cgcgcTCTGAGACGAAACGGGG.GGCGCTGCCCCC.CTC.TT..GGCCTTC.GGCC.AATTCACCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTc                              <  1:2076249/101‑1 (MQ=21)
                                                              cgcgcTCTGAGACGAAACGGGG.GGCGCTGCCCCC.CTC.TT..GGCCTTC.GGCC.AATTCACCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTc                              <  1:509035/101‑1 (MQ=21)
                                                               gcgcTCTGAGGCGAAACGGGG.GGCGCTGCCCCC.CTC.TT..GGCCTTC.GGCC.AATTCACCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCa                             <  1:2096435/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:31132/101‑1 (MQ=21)
                                                                 gcTCTGAGACGAAACGGGG.GGCGCTGCCCCC.CTC.TT..GGCCTTC.GGCC.AATTCACCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATc                           <  1:22833/101‑1 (MQ=21)
                                                                 gcTCTGAGACGAAACGGGG.GGCGCTGCCCCC.CTC.TT..GGCCTTC.GGCC.AATTCACCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATc                           <  1:902741/101‑1 (MQ=21)
                                                                   tctGAGACGAAACGGGG.GGCGCTGCCCCC.CTC.TT..GGCCTTC.GGCC.AATTCACCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCtt                         <  1:298907/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)
                                                                       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)
                                                                           gAAACGGGG.GGCGCTGCCCCC.CTC.TT..GGCCTTC.GGCC.AATTCACCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGaaa                 <  2:1593263/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)
                                                                               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)
                                                                                  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:2334609/101‑5 (MQ=2)
                                                                                     ggCGCTGCCCCCTCTT.GG..C..CTTT.GGCC.AATTCA.CCCCTGGGATATTTCGGGCAAGAT...GAA.ATCATGCCCCTTCATCTTGCCGCAAATATCCCggggggagc  <  1:1048067/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)
                                                                                         ccgCCCCC.CTC.TT..GGCCTTC.GGCC.AAGTCCCCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCCGGGGGAg    <  2:3568447/99‑1 (MQ=37)
                                                                                         cTGCCCCC.CTC.TT..GGCCTTC.GGCC.AATTCACCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCCGGGGGAg    <  2:1142532/101‑1 (MQ=37)
                                                                                                       |                                                                                              
TATGCGATGGTCGATGCGATCGACTGGCCCGAGGGGTTCTGCCGCCGCGACATCGGCTGGCGCGCGCTCTGAGACGAAAC.GGG.GGCGCTGCCCCC.CTC.CT..G..CTTC.GGCC.AATTCACCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCCGGGGGAGCC  >  NC_014034/2566457‑2566644

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