Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_014034 2,566,577 C→T 5.8% intergenic (+49/+145) 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,5770C→T 5.8% 11.7 ‑1.6 ‑0.1121intergenic (+49/+145)purD/sdaAphosphoribosylamine‑‑glycine ligase/L‑serine ammonia‑lyase
Reads supporting (aligned to +/- strand):  new base (4/3):  ref base (32/82):  total (36/85)
Fisher's exact test strand distribution p-value = 1.94e-01
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 8.20e-01

Only 100 of 146 total aligned reads displayed.

CTGGCCCGAGGGGTTCTGCCGCCGCGACATCGGCTGGCGCGCGCTCTGAGACGAAAC.GGG.GGCGCTGCCCCC.CTC.CT..G..CTTCGGCC.AATTCA.CCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCCGGGGGAGCCCGCTTGCGGGCGGGGGCAGCGCCCCCGAAGCGC  >  NC_014034/2566480‑2566677
                                                                                                           |                                                                                                    
cTGGCCCGAGGGGTTCTGCCGCCGCGACATCGGCTGGCGCGCGCTCTGAGACGAAACGGGG.GGCGCTGCCCCC.CTC.TT..GGCCTTCGGCC.AATTCA.cccccc                                                                                                      <  2:702261/101‑1 (MQ=25)
 tGGCCCGAGGGGTTCTGCCGCCGCGACATCGGCTGGCGCGCGCTCTGAGACGAAACGGGG.GGCGCTGCCCCC.CTC.TT..GGCCTTCGGCC.AATTCA.CCCCCCg                                                                                                     <  1:2181231/101‑1 (MQ=25)
  ggCCCGAGGGGTTCTGCCGCCGCGACATCGGCTGGCGCGCGCTCTGAGACGAAACGGGG.GGCGCTGCCCCC.CTC.TT..GGCCTTCGGCC.AATTCA.CCCCCCgg                                                                                                    <  1:503621/101‑1 (MQ=25)
  ggCCCGAGGGGTTCTGCCGCCGCGACATCGGCTGGCGCGCGCTCTGAGACGAAACGGGG.GGCGCTGCCCCC.CTC.TT..GGCCTTCGGCC.AATTCA.CCCCCCgg                                                                                                    <  1:1099121/101‑1 (MQ=25)
  ggCCCGAGGGGTTCTGCCGCCGCGACATCGGCTGGCGCGCGCTCTGAGACGAAACGGGG.GGCGCTGCCCCC.CTC.TT..GGCCTTCGGCC.AATTCA.CCCCCCgg                                                                                                    <  1:2944444/101‑1 (MQ=25)
   gCCCGAGGGGTTCTGCCGCCGCGACATCGGCTGGCGCGCGCTCTGAGACGAAAC.GGGTGGGGCTGCCCCC.CTC.TT..GGCCCTCCGCC.AATTCA.CCCCCCcgg                                                                                                   >  2:3490303/1‑101 (MQ=25)
    cccGAGGGGTTCTGCCGCCGCGACATCGGCTGGCGCGCGCTCTGAGACGAAACGGGG.GGCGCTGCCCCC.CTC.TT..GGCCTTCGGCC.AATTCA.CCCCCCGGGa                                                                                                  <  1:1829286/101‑1 (MQ=25)
     ccGAGGGGTTCTGCCGCCGCGACATCGGCTGGCGCGCGCTCTGAGACGAAACGGGG.GGGGCCGCCCCC.CCCTTT..GC.CTTCGGCC.AAATCA.CCCCCCcggga                                                                                                 >  1:1978096/1‑99 (MQ=255)
     ccGAGGGGTTCTGCCGCCGCGACATCGGCTGGCGCGCGCTCTGAGACGAAACGGGG.GGCGCTGCCCCC.CTC.TT..GGCCTTCGGCC.AATTCA.CCCCCCGGGat                                                                                                 <  2:2168907/101‑1 (MQ=25)
     ccGAGGGGTTCTGCCGCCGCGACATCGGCTGGCGCGCGCTCTGAGACGAAACGGGG.GGCGCTGCCCCC.CTC.TT..GGCCTTCGGCC.AATTCA.CCCCCCGGGat                                                                                                 <  2:4180170/101‑1 (MQ=25)
     ccGAGGGGTTCTGCCGCCGCGACATCGGCTGGCGCGCGCTCTGAGACGAAACGGGG.GGCGCTGCCCCC.CTC.TT..GGCCTTCGGCC.AATTCA.CCCCCCGGGat                                                                                                 <  1:3194193/101‑1 (MQ=25)
      cGAGGGGTTCTGCCGCCGCGACATCGGCTGGCGCGCGCTCTGAGACGAAACGGGG.GGGGCTGCCCCC.CCC.TT..GGCCTTCCGCC.AAATCC.CCCCCCCGGgta                                                                                                >  2:3168341/1‑101 (MQ=255)
      cGAGGGGTTCTGCCGCCGCGACATCGGCTGGCGCGCGCTCTGAGACGAAACGGGG.GGCGCTGCCCCC.CTC.TT..GGCCTTCGGCC.AATTCA.CCCCCCGGGata                                                                                                <  2:1018617/101‑1 (MQ=25)
       gAGGGGTTCTGCCGCCGCGACATCGGCTGGCGCGCGCTCTGAGACGAAACGGGG.GGCGCTGCCCCC.CTC.TT..GGCCTTCGGCC.AATTCA.CCCCCCCGGatat                                                                                               >  1:2119448/1‑101 (MQ=25)
         ggggTTCTGCCGCCGCGACATCGGCTGGCGCGCGCTCTGAGACGAAACGGGG.GGCGCTGCCCCC.CTC.TT..GGCCTTCGGCC.AATTCA.CCCCCCGGGATAttt                                                                                             <  1:1697153/101‑1 (MQ=21)
         ggggTTCTGCCGCCGCGACATCGGCTGGCGCGCGCTCTGAGACGAAACGGGG.GGCGCTGCCCCC.CTC.TT..GGCCTTCGGCC.AATTCA.CCCCCCGGGATAttt                                                                                             <  2:77585/101‑1 (MQ=21)
              tCTGCCGCCGCGACATCGGCTGGCGCGCGCTCTGAGACGAAACGGGG.GGCGCTGCCCCC.CTC.TT..GGCCTTCGGCC.AATTCA.CCCCCCGGGATATTTCGGGc                                                                                        <  1:1268244/101‑1 (MQ=21)
              tCTGCCGCCGCGACATCGGCTGGCGCGCGCTCTGAGACGAAACGGGG.GGCGCCGCCCCC.CTC.TTTGG..CCTCGGCCAAATCCA.CCCCCCcggtacttttggg                                                                                         >  2:3685661/1‑91 (MQ=255)
               cTGCCGCCGCGACATCGGCTGGCGCGCGCTCTGAGACGAAACGGGG.GGCGCTGCCCCC.CTC.TT..GGCCTTCGGCC.AATTCA.CCCCCCGGGATATTTCGGGCa                                                                                       <  2:2900729/101‑1 (MQ=21)
               cTGCCGCCGCGACATCGGCTGGCGCGCGCTCTGAGACGAAACGGGG.GGCGCTGCCCCC.CTC.TT..GGCCTTCGGCC.AATTCA.CCCCCCCGGCTATTTCGGGCa                                                                                       >  2:861613/1‑101 (MQ=21)
                tGCCGCCGCGACATCGGCTGGCGCGCGCTCTGAGACGAAGCGGGG.GGCGCTGCCCCC.CTC.TT..GGCCTTCGGCC.AATTCA.CCCCCCGGGATATTTCGGGCaa                                                                                      <  1:473092/101‑1 (MQ=21)
                tGCCGCCGCGACATCGGCTGGCGCGCGCTCTGAGACGAAACGGGG.GGCGCTGCCCCC.CCC.TT..GGCCTTCGGCC.AATTCA.CCCCCCCGGAGAtttttggcca                                                                                      >  2:3957956/1‑94 (MQ=18)
                   cgccgcGACATCGGCTGGCGCGCGCTCTGAGACGAAACGGGG.GGGGCTGGCCCC.CTC.TT..GGCCTTTCGCC.AATTCA.CCCCCCCGGACATTTCGGGCAAGaa                                                                                   >  1:3980246/1‑100 (MQ=255)
                   cgccgcGACATCGGCTGGCGCGCGCTCTGAGACGAAACGGGG.GGCGCTGCCCCC.CTC.TT..GGCCTTCGGCC.AATTCA.CCCCCCGGGATATTTCGGGCAAGAt                                                                                   <  2:1620144/101‑1 (MQ=21)
                   cgccgcGACATCGGCTGGCGCGCGCTCTGAGACGAAACGGGG.GGCGCTGCCCCC.CTC.TT..GGCCTTCGGCC.AATTCA.CCCCCCGGGATATTTCGGGCAAGAt                                                                                   <  2:2060072/101‑1 (MQ=21)
                   cgccgcGACATCGGCGGGCGCGCGCTCTGAGACGAAACGGGG.GGCGCTGCCCCC.CTC.TT..GGCCTTCGGCC.AATTCA.CCCCCCGGGATATTTCGGGCAAGAt                                                                                   <  2:523134/101‑1 (MQ=21)
                    gccgcGACATCGGCTGGCGCGCGCTCTGAGACGAAACGGGG.GGCGCTGCCCCC.CTC.TT..GGCCTTCGGCC.AATTCA.CCCCCCGGGATATTTCGGGCAAGATa                                                                                  <  2:3108432/101‑1 (MQ=21)
                      cgcgACATCGGCTGGCGCGCGCTCTGAGACGAAACGGGG.GGCGCTGCCCCC.CTC.TT..GGCCTTCGGCC.CATTCA.CCCCCCCGGATATTTCGGGCAAGATaaa                                                                                >  1:3906941/1‑101 (MQ=21)
                       gcgACATCGGCTGGCGCGCGCTCTGAGACGAAACGGGG.GGCGCTGCCCCC.CTC.TT..GGCCTTCGGCC.AATTCA.CCCCCCGGGATATTTCGGGCAAGATAaag                                                                               <  1:430223/101‑1 (MQ=21)
                        cgACATCGGCTGGCGCGCGCTCTGAGACGAAACGGGG.GGCGCTGCCCCC.CTC.TT..GGCCTTCGGCC.AATTCA.CCCCCCGGGATATTTCGGGCAAGATAaaga                                                                              >  2:3176725/1‑101 (MQ=21)
                             tCGGCTGGCGCGCGCTCTGAGACGAAACGGGG.GGCGCTGCCCCC.CTC.TT..GGCCTTCGGCC.AATTCA.CCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTc                                                                         >  1:3848117/1‑101 (MQ=21)
                             tCGGCTGGCGCGCGCTCTGAGACGAAACGGGG.GGCGCTGCCCCC.CTC.TT..GGCCTTCGGCC.AATTCA.CCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTc                                                                         <  1:3471783/101‑1 (MQ=21)
                               ggCTGGGGCGCGCTCTGAGACGAAACGGGG.GGCGCTGCCCCC.CTC.TT..GGCCTTCGGCC.AATTCA.CCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCtt                                                                       <  1:3073594/101‑1 (MQ=21)
                               ggCTGGCGCGCCCTCTGAGACGAAACGGGG.GGCGCTGCCCCC.CTC.TT..GGCCTTCGGCC.AATTCA.CCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCtt                                                                       <  1:1943026/101‑1 (MQ=21)
                                gCTGGCGCGCGCTCTGAGACGAAACGGGG.GGCGCTGCCCCC.CTC.TT..GGCCTTCGGCC.AATTCA.CCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCttt                                                                      <  2:1078597/101‑1 (MQ=21)
                                 cTGGCGCGCGCTCTGAGACGAAACGGGG.GGCGCTGCCCCC.CTC.TT..GGCCTTCGGCC.AATTCA.CCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCtttt                                                                     <  1:2327547/101‑1 (MQ=21)
                                  tGGCGCGCGCTCTGAGACGAAACGGGG.GGCGCTGCCCCC.CTC.TT..GGCCTTCGGCC.AATTCA.CCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTc                                                                    <  1:3248456/101‑1 (MQ=21)
                                  tGGCGCGCGCTCTGAGACGAAACGGGG.GGCGCTGCCCCC.CTC.TT..GGCCTTCGGCC.AATTCA.CCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTc                                                                    >  2:1009588/1‑101 (MQ=21)
                                    gcgcgcgcTCTGAGACGAAACGGGG.GGCGCTGCCCCC.CTC.TT..GGCCTTCGGCC.AATTCA.CCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCt                                                                  <  1:3376059/101‑1 (MQ=21)
                                     cgcgcgcTCTGAGACGAAACGGGG.GGCGCTGCCCCC.CTC.TT..GGCCTTCGGCC.AATTCA.CCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCtt                                                                 <  1:2512408/101‑1 (MQ=21)
                                      gcgcgcTCTGAGACGAAACGGGG.GGCGCTGCCCCC.CTC.TT..GGCCTTCGGCC.AATTCA.CCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCttt                                                                <  2:1830164/101‑1 (MQ=21)
                                      gcgcgcTCTGAGACGAAACGGGG.GGCGCTGCCCCC.CTC.TT..GGCCTTCGGCC.AATTCA.CCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCttt                                                                <  2:2428260/101‑1 (MQ=21)
                                        gcgcTCTGAGGCGAAACGGGG.GGCGCTGCCCCC.CTC.TT..GGCCTTCGGCC.AATTCA.CCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCa                                                              <  1:2096435/101‑1 (MQ=21)
                                         cgcTCTGAGACGAAACGGGG.GGCGCTGCCCCC.CTC.TT..GGCCTTCGGCC.AATTCA.CCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCAt                                                             <  1:4145111/101‑1 (MQ=21)
                                          gcTCTGAGACGAAACGGGG.GGCGCTGCCCCC.CTC.TT..GGCCTTCGGCC.AATTCA.CCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATc                                                            <  1:31132/101‑1 (MQ=21)
                                          gcTCTGAGACGAAACGGGG.GGCGCTGCCCCC.CTC.TT..GGCCTTCGGCC.AATTCA.CCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATc                                                            <  1:22833/101‑1 (MQ=21)
                                            tctGAGACGAAACGGGG.GGCGCTGCCCCC.CTC.TT..GGCCTTCGGCC.AATTCA.CCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCtt                                                          <  1:298907/101‑1 (MQ=21)
                                             ctGAGACGAAACGGGG.GGCGCTGCCCCC.CTC.TT..GGCCTTCGGCC.AATTCA.CCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTg                                                         <  1:569431/101‑1 (MQ=21)
                                              tGAGACGAAACGGGG.GGCGCTGCCCCC.CTC.TT..GGCCTTCGGCC.AATTCA.CCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTGTTTTTTTTTCTTCTTGc                                                        >  2:853783/1‑101 (MQ=18)
                                              tGAGACGAAACGGGG.GGCGCTGCCCCC.CTC.TT..GGCCTTCGGCC.AATTCA.CCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGc                                                        >  2:2272230/1‑101 (MQ=21)
                                               gagaCGAAACGGGG.GGCGCTGCCCCC.CTC.TT..GGCCTTCGGCC.AATTCA.CCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGcc                                                       >  1:1222553/1‑101 (MQ=21)
                                               gagaCGAAACGGGG.GGCGCTGCCCCC.CTC.TT..GGCCTTCGGCC.AATTCA.CCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGcc                                                       >  2:1826958/1‑101 (MQ=21)
                                                agaCGAAACGGGG.GGCGCTGCCCCC.CTC.TT..GGCCTTCGGCC.AATTCA.CCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGccc                                                      <  2:236725/101‑1 (MQ=21)
                                                  aCGAAACGGGG.GGCGCTGCCCCC.CTC.TT..GGCCTTCGGCC.AATTCA.CCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGa                                                    >  1:1787489/1‑101 (MQ=21)
                                                   cGAAACGGGG.GGCGCTGCCCCC.CTC.TT..GGCCTTCGGCC.AATTCA.CCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGaa                                                   >  2:3077751/1‑101 (MQ=21)
                                                    gAAACGGGG.GGCGCTGCTCCC.CTC.TT..GGCCTTCGGCC.AATTCA.CCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGaaa                                                  <  2:3806296/101‑1 (MQ=21)
                                                    gAAACGGGG.GGCGCTGCCCCC.CTC.TT..GGCCTTCGGCC.AATTCA.CCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGaaa                                                  <  2:1593263/101‑1 (MQ=21)
                                                    gAAACGGGG.GGCGCTGCCCCC.CTC.TT..GGCCTTCGGCC.AATTCA.CCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGaaa                                                  <  2:2702264/101‑1 (MQ=21)
                                                      aacgggg.ggCGCTGCCCCC.CTC.TT..GGCCTTCGGCC.AATTCACCCCCCCGGGAAATTT.GGGCCAGAAAAAGAAAGGCTTTTTCTTTCATCTTTcccccaaaa                                                >  1:2185736/5‑95 (MQ=2)
                                                      aacgggg.ggCGCTGCCCCC.CTC.TT..GGCCTTCGGCC.AATTCA.CCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAata                                                <  1:3581109/97‑1 (MQ=25)
                                                       acgggg.ggCGCTGCCCCC.CTC.TT..GGCCTTCGGCC.AATTCA.CCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAatat                                               <  2:4156408/98‑1 (MQ=25)
                                                        cgggg.ggCGCTGCCCCC.CTC.TT..GGCCTTCGGCC.AATTCA.CCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATc                                              <  2:1120311/99‑1 (MQ=25)
                                                        cgggg.ggCGCTGCCCCC.CTC.TT..GGCCTTCGGCC.AATTCA.CCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATc                                              <  2:4015211/99‑1 (MQ=25)
                                                         gggg.ggCGCTGCCCCC.CTC.TT..GGCCTTCGGCC.AATTCA.CCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATcc                                             <  2:1571825/100‑1 (MQ=25)
                                                         aggg.ggCGCTG.CCCC.CTC.TT..GGCCTTTGGCC.AATTCA..CCCCTGGGATATTTCGGGCAAGAT...GAA.ATCATGCCCCTTCATCTTGCCGCAAATATCCCggggg                                       >  1:4076720/2‑101 (MQ=0)
                                                           gg.ggCGCTGCCCCCTCTT.GG..C..CTTTGGCC.AATTCA..CCCCTGGGATATTTCGGGCAAGAT...GAA.ATCATGCCCCTTCATCTTGCCGCAAATATCCCgggggga                                     <  1:3959563/101‑3 (MQ=2)
                                                            g.ggCGCTG.CCCC.CTC.TT..GGCCTTTGGCC.AATTCA..CCCCTGGGATATTTCGGGCAAGAT...GAA.ATCATGCCCCTTCATCTTGCCGCAAATATCCCggggggag                                    >  1:1254032/1‑98 (MQ=2)
                                                              ggCGCTGCCCCCTCTT.GG..C..CTTTGGCC.AATTCA..CCCCTGGGATATTTCGGGCAAGAT...GAA.ATCATGCCCCTTCATCTTGCCGCAAATATCCCggggggagc                                   <  1:1048067/101‑5 (MQ=2)
                                                                cgcTGCCCCC.CTC.TT..GGCCTTCGGCC.AATTCA.CCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCCAAATATCCCggggg                                       >  2:3022218/1‑101 (MQ=25)
                                                                 gcTGCCCCC.CTC.TT..GGCCTTCGGCC.AATTCT.CCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCCGGGGGa                                      <  1:2866857/101‑1 (MQ=37)
                                                                  ccgCCCCC.CTC.TT..GGCCTTCGGCC.AAGTCC.CCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCCGGGGGAg                                     <  2:3568447/99‑1 (MQ=37)
                                                                            ccctcttggccTTCGGCC.AATTCC.CCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCCGGGGGAGCCCGCt                               <  1:1052427/91‑1 (MQ=38)
                                                                            ccctcttggccTTCGGCC.AATTCA.CCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCCGGGGGAGCCCGCt                               <  1:93659/91‑1 (MQ=255)
                                                                              ctcttggccTTCGGCC.AATTCC.CCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCCGGGGGAGCCCGCTTg                             <  2:2117738/93‑1 (MQ=255)
                                                                                cttggccTTCGGCC.AATTCA.CCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCCGGGGGAGCCCGCTTgcg                           >  2:3746717/7‑101 (MQ=255)
                                                                                 ttggccTTCGGCC.AATTCA.CCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCCGGGGGAGCCCGCTTgcgg                          <  2:2553261/96‑1 (MQ=255)
                                                                                 ttggccTTCGGCC.AATTCA.CCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCCGGGGGAGCCCGCTTgcgg                          <  1:3070202/96‑1 (MQ=255)
                                                                                   ggccTTCGGCC.AATTCA.CCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCCGGGGGAGCCCGCTTgcgggc                        >  1:1210820/4‑101 (MQ=255)
                                                                                   ggccTTCGGCC.AATTCA.CCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCCGGGGGAGCCCGCTTgcgggc                        >  1:721715/4‑101 (MQ=255)
                                                                                   ggccTTCGGCC.AATTCA.CCCCCCGGGATATTTCGGGCAAGATAAAGAAGGGCTTTTCCTTTCATCTTGCCCGAAAGATCCCGGGGGACCCCGCCTgggggc                        >  1:3885267/4‑101 (MQ=25)
                                                                                     ccTTCGGCC.AATTCA.CCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCCGGGGGAGCCCGCTTgcgggcag                      >  2:1717474/2‑99 (MQ=255)
                                                                                      cTTCGGCC.AATTCA.CCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCCGGGGGAGCCCGCTTGcgggcggg                     >  2:2407805/1‑101 (MQ=255)
                                                                                         cGGCC.AATTCA.CCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCCGGGGGAGCCCGCTTGCGGGCGGGGGc                  <  2:2144179/101‑1 (MQ=255)
                                                                                          ggCC.AATTCA.CCCCCCGGGATATTTTGGGCGAGATAAAGAAGGGGTTTTTCTTTTCTCTTTCCCGAAAAATCACGGGGGGGGCCGCCTCGTGGCggggggg                 >  2:756706/1‑99 (MQ=16)
                                                                                           gCC.AATTCA.CCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCCGGGGGAGCCCGCTTGCGGGCGGGGGCAg                <  1:3377084/101‑1 (MQ=255)
                                                                                             c.AAGTCA.CCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCCGGGGGAGCCCGCTTGCGGGCGGGGGCAgcg              <  1:1876199/101‑1 (MQ=255)
                                                                                               aaTTCA.CCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCCGGGGGAGCCCGCTTGCGGGCGGGGGCAgcgc             >  2:1742483/1‑101 (MQ=255)
                                                                                                aTTCA.CCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCCGGGGGAGCCCGCTTGCGGGCGGGGGCAGCGcc            >  1:3091680/1‑101 (MQ=255)
                                                                                                aTTCA.CCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCCGGGGGAGCCCGCTTGCGGGCGGGGGCAGCGcc            >  1:2704649/1‑101 (MQ=255)
                                                                                                 ttCA.CCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCCGGGGGAGCCCGCTTGCGGGCGGGGGCAGCGccc           >  1:652542/1‑101 (MQ=255)
                                                                                                  tCA.CCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCCAGGGGAGCCCGCTTGCGGGGGGGGGGAGCGcccc          >  2:1633744/1‑101 (MQ=255)
                                                                                                   cA.CCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCCGGGGGAGCCCGCTTGCGGGCGGGGGCAGCGccccc         <  1:2029001/101‑1 (MQ=255)
                                                                                                    a.CCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCCGGGGGAGCCCGCTTGCGGGCGGGGGCAGCGCCCCCg        >  1:584338/1‑101 (MQ=255)
                                                                                                       cccccGAGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCCGGGGGAGCCCGCTTGCGGGCGGGGGCAGCGCCCCCGaa      <  2:2040485/101‑1 (MQ=255)
                                                                                                         cccGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCCGGGGGAGCCCGCTTGCGGGCGGGGGCAGCGCCCCCGAAgc    <  1:88752/101‑1 (MQ=255)
                                                                                                         cccGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCCGGGGGAGCCCGCTTGCGGGCGGGGGCAGCGCCCCCGAAgc    <  1:1742602/101‑1 (MQ=255)
                                                                                                          ccGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCCGGGGGAGCCCGCTTGCGGGCGGGGGCAGCGCCCCCGAAgcg   <  1:3962109/101‑1 (MQ=255)
                                                                                                          ccGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCCGGGGGAGCCCGCTTGCGGGCGGGGGCAGCGCCCCCGAAgcg   >  2:1190061/1‑101 (MQ=255)
                                                                                                           cGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCCGGGGGAGCCCGCTTGCGGGCGGGGGCAGCGCCCCCGAAgcgc  >  1:946245/1‑101 (MQ=255)
                                                                                                           cGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCCGGGGGAGCCCGCTTGCGGGCGGGGGCAGCGCCCCCGAAgcgc  >  1:1676698/1‑101 (MQ=255)
                                                                                                           |                                                                                                    
CTGGCCCGAGGGGTTCTGCCGCCGCGACATCGGCTGGCGCGCGCTCTGAGACGAAAC.GGG.GGCGCTGCCCCC.CTC.CT..G..CTTCGGCC.AATTCA.CCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCCGGGGGAGCCCGCTTGCGGGCGGGGGCAGCGCCCCCGAAGCGC  >  NC_014034/2566480‑2566677

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