Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_014034 2,566,551 +T 1.9% intergenic (+23/+171) purD → / ← sdaA phosphoribosylamine‑‑glycine ligase/L‑serine ammonia‑lyase

Read alignment evidence...
  seq id position change freq score reads annotation genes product
*NC_0140342,566,5511.→T 1.9% 4.7 ‑1.6 ‑0.7162intergenic (+23/+171)purD/sdaAphosphoribosylamine‑‑glycine ligase/L‑serine ammonia‑lyase
Reads supporting (aligned to +/- strand):  new base (2/1):  ref base (45/114):  total (47/115)
Fisher's exact test strand distribution p-value = 2.02e-01
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 2.17e-01

Only 100 of 175 total aligned reads displayed.

GCGCCTATGCGATGGTCGATGCGATCGACTGGCCCGAGGGGTTCTGCCGCCGCGACATCGGCTGGCGCGCGCTCTGAGACG.AAAC.GGGG.GCGCTGCCCCC.CTC.CTG..CTTC.GGCCAATTCACCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCCGGGGGAGCC  >  NC_014034/2566452‑2566644
                                                                                                       |                                                                                                 
gcgcCTATGCGATGGTCGATGCGATCGACTGGCCCGAGGGGTTCTGCCGCCGCGACATCGGCTGGCGCGCGCTCTGAGACG.AAACGGGGG.GCGCTGccccc                                                                                                    >  1:1763875/1‑101 (MQ=255)
     tATGCGGTGGTCGATGCGATCGACTGGCCCGAGGGGTTCTGCCGCCGCGACATCGGCTGGCGCGCGCTCTGAGACG.AAACGGGGG.GCGCTGCCCCc.ctctt                                                                                              <  1:3729515/101‑3 (MQ=255)
     tATGCGATGGTCGATGCGATCGACTGGCCCGAGGGGTTCTGCCGCCGCGACATCGGCTGGCGCGCGCTCTGAGACG.AAACGGGGG.GCGCTGCCCCc.ctctt                                                                                              <  2:3209079/101‑3 (MQ=255)
      aTGCGATGGTCGATGCGATCGACTGGCCCGAGGGGTTCTGCCGCCGCGACATCGGCTGGCGCGCGCTCTGAGACG.AAACGGGGG.GCGCTGCCCCC.CTC.TTg                                                                                            <  2:4014949/101‑1 (MQ=255)
      aTGCGATGGTCGATGCGATCGACTGGCCCGAGGGGTTCTGCCGCCGCGACATCGGCTGGCGCGCGCTCTGAGACG.AAACGGGGG.GCGCTGCCCCC.CTC.TTg                                                                                            <  1:628652/101‑1 (MQ=255)
        gCGATGGTCGATGCGATCGACTGGCCCGAGGGGTTCTGCCGCCGCGACATCGGCTGGCGCGCGCTCTGAGACG.AAACGGGGG.GCGCTGCCCCC.CTC.TTggc                                                                                          <  1:1394923/101‑3 (MQ=255)
          gATGGTCGATGCGATCGACTGGCCCGAGGGGTTCTGCCGCCGCGACATCGGCTGGCGCGCGCTCTGAGACG.AAACGGGGG.GCGCTGCCCCC.CTC.TTggcct                                                                                        <  1:3895620/101‑5 (MQ=39)
           aTGGTCGATGCGATCGACTGGCCCGAGGGGTTCTGCCGCCGCGACATCGGCTGGCGCGCGCTCTGAGACG.AAACGGGGG.GCGCTGCCCCC.CTC.TTggcctt                                                                                       <  2:1659909/101‑6 (MQ=38)
               tCGATGCGATCGACTGGCCCGAGGGGTTCTGCCGCCGCGACATCGGCTGGCGCGCGCTCTGAGACG.AAACGGGGG.GCGCTGCCCCC.CTC.TTggccttcggc                                                                                   >  2:1480413/1‑92 (MQ=37)
                cGATGCGATCGACTGGCCCGAGGGGTTCTGCCGCCGCGACATCGGCTGGCGCGCGCTCTGAGACG.AAACGGGGG.GCGCTGCCCCC.CTC.TTGGCCTTC.GGcc                                                                                 <  1:4189647/101‑1 (MQ=37)
                cGATGCGATCGACTGGCCCGAGGGGTTCTGCCGCCGCGACATCGGCTGGCGCGCGCTCTGAGACG.AAACGGGGG.GCGCTGCCCCC.CTC.TTGGCCTTC.GGcc                                                                                 <  2:2103579/101‑1 (MQ=37)
                 gATGCGATCGACTGGCCCGAGGGGTTCTGCCGCCGCGACATCGGCTGGCGCGCGCTCTGAGACG.AAACGGGGG.GCGCTGCCCCC.CTC.TTGGCCTTC.GGCCa                                                                                <  2:373762/101‑1 (MQ=37)
                 gATGCGATCGACTGGCCCGAGGGGTTCTGCCGCCGCGACATCGGCTGGCGCGCGCTCTGAGACG.AAACGGGGG.GCGCTGCCCCC.CTC.TTGGCCTTC.GGCCa                                                                                >  2:1041278/1‑101 (MQ=37)
                  aTGCGATCGACTGGCCCGAGGGGTTCTGCCGCCGCGACATCGGCTAGCGCGCGCTCTGAGACG.AAACGGGGG.GCGCTGCCCCC.CTC.TTGGCCTTC.GGCCaa                                                                               <  1:3924085/101‑1 (MQ=25)
                   tGCGATCGACTGGCCCGAGGGGTTCTGCCGCCGCGACATCGGCTGGCGCGCGCTCTGAGACG.AAACGGGGG.GCGCTGCCCCC.CTC.TTGGCCTTC.GGCCAAt                                                                              <  1:1609935/101‑1 (MQ=37)
                    gCGATCGACTGGCCCGAGGGGTTCTGCCGCCGCGACATCGGCTGGCGCGCGCTCTGAGACG.AAACGGGGG.GCGCTGCCCCC.CTC.TTGGCCTTC.GGCCAAtt                                                                             <  2:2292218/101‑1 (MQ=37)
                     cGATCGACTGGCCCGAGGGGTTCTGCCGCCGCGACATCGGCTGGCGCGCGCTCTGAGACG.AAACGGGGG.GCGCTGCCCCC.CTC.TTGGCCTTC.GGCCAATTc                                                                            <  1:3216613/101‑1 (MQ=37)
                      gATCGACTGGCCCGAGGGGTTCTGCCGCCGCGACATCGGCTGGCGCGCGCTCTGAGACG.AAACGGGGG.GCGCTGCCCCC.CTC.TTGGCCTTC.GGCCAATTCa                                                                           <  2:2878599/101‑1 (MQ=37)
                      gATCGACTGGCCCGAGGGGTTCTGCCGCCGCGACATCGGCTGGCGCGCGCTCTGAGACG.AAACGGGGG.GCGCTGCCCCC.CTC.TTGGCCTTC.GGCCAATTCa                                                                           <  2:2545304/101‑1 (MQ=37)
                      gATCGACTGGCCCGAGGGGTTCTGCCGCCGCGACATCGGCTGGCGCGCGCTCTGAGACG.AAACGGGGG.GCGCTGCCCCC.CTC.TTGGCCTTC.GGCCAATTCa                                                                           <  1:2629475/101‑1 (MQ=37)
                           aCTGGCCCGAGGGGTTCTGCCGCCGCGACATCGGCTGGCGCGCGCTCTGAGACG.AAACGGGGG.GCGCTGCCCCC.CTC.TTGGCCTTC.GGCCAATTCAccccc                                                                      <  2:898937/101‑1 (MQ=25)
                            cTGGCCCGAGGGGTTCTGCCGCCGCGACATCGGCTGGCGCGCGCTCTGAGACG.AAACGGGGG.GCGCTGCCCCC.CTC.TTGGCCTTC.GGCCAATTCAcccccc                                                                     <  1:1662172/101‑1 (MQ=25)
                             tGGCCCGAGGGGTTCTGCCGCCGCGACATCGGCTGGCGCGCGCTCTGAGACG.AAACGGGGG.GCGCTGCCCCC.CTC.TTGGCCTTC.GGCCAATTCACCCCCCg                                                                    <  1:229977/101‑1 (MQ=25)
                             tGGCCCGAGGGGTTCTGCCGCCGCGACATCGGCTGGCGCGCGCTCTGAGACG.AAACGGGGG.GCGCTGCCCCC.CTC.TTGGCCTTC.GGCCAATTCACCCCCCg                                                                    <  1:1528177/101‑1 (MQ=25)
                             tGGCCCGAGGGGTTCTGCCGCCGCGACATCGGCTGGCGCGCGCTCTGAGACG.AAACGGGGG.GCGCTGCCCCC.CTC.TTGGCCTTC.GGCCAATTCACCCCCCg                                                                    <  1:730590/101‑1 (MQ=25)
                              ggCCCGAGGGGTTCTGCCGCCGCGACATCGGCTGGCGCGCGCTCTGAGACG.AAACGGGGG.GCGCTGCCCCC.CTC.TTGGCCTTC.GGCCAATTCACCCCCCgg                                                                   <  1:52255/101‑1 (MQ=25)
                              ggCCCGAGGGGTTCTGCCGCCGCGACATCGGCTGGCGCGCGCTCTGAGACG.AAACGGGGG.GCGCTGCCCCC.CTC.TTGGCCTTC.GGCCAATTCACCCCCCgg                                                                   <  1:4585449/101‑1 (MQ=25)
                                gccGGGGGTTTCTGCCGCCGCGACATCGGCTGGCGCGCGCTCTGAGACG.AAACGGGGG.GCGCTGCCCCC.CTC.TTGGCCTTC.GGCCAATTCACCCCCCGGGa                                                                 <  2:2989507/100‑1 (MQ=21)
                                cccGAGGGGTTCTGCCGCCGCGACATCGGCTGGCGCGCGCTCTGAGACG.AAACGGGGG.GCGCTGCCCCC.CCCTTTGC.CTTC.CGCCCATTCACCCCCCgggg                                                                 >  1:2437129/1‑100 (MQ=25)
                                 ccGAGGGGTTCTGCCGCCGCGGCATCGGCTGGCGCGCGCTCTGAGACG.AAACGGGGG.GCGCTGCCCCC.CTC.TTGGCCTTC.GGCCAATTCACCCCCCGGGat                                                                <  2:399086/101‑1 (MQ=21)
                                   gAGGGGTTCTGCCGCCGCGACATCGGCTGGCGGGCGCTCTTGGACG.AAACGGGGG.GCGGGGCCCCC.CTCTTTGC.CTTC.AGCCACTTCAGCCCCCCGGagat                                                              >  2:2732121/1‑101 (MQ=255)
                                    aGGGGTTCTGCCGCCGCGACATCGGCTGGCGCGCGCTCTGAGACG.AAACGGGGG.GCGCTGCCCCC.CTC.TTGGCCTTC.GGCCAATTCACCCCCCGGGATAtt                                                             <  1:2930665/101‑1 (MQ=21)
                                     ggggTTCTGCCGCCGCGACATCGGCTGGCGCGCGCTCTGAGACG.AAACGGGGG.GCGCTGCCCCC.CTC.TTGGCCTTC.GGCCAATTCACCCCCCGGGATAttt                                                            <  1:634884/101‑1 (MQ=21)
                                      gggTTCTGCCGCCGCGACATCGGCTGGCGCGCGCTCTGAGACG.AAACGGGGG.GCGCTGCCCCC.CTC.TTGGCCTTC.GGCCAATTCACCCCCCGGGATATTTc                                                           <  1:4294553/101‑1 (MQ=21)
                                      gggTTCTGCCGCCGCGACATCGGCTGGCGCGCGCTCTGAGACG.AAACGGGGG.GCGCTGCCCCC.CTC.TTGGCCTTC.GGCCAATTCACCCCCCGGGATACTTc                                                           <  1:3591918/101‑1 (MQ=21)
                                        gTTCTGCCGCCGCGACATCGGCTGGCGCGCGCTCTGAGACG.AAACGGGGG.GCGCTGCCCCC.CTC.TTGGCCTTC.GGCCAATTCACCCCCCGGGATATTTCgg                                                         <  1:1572687/101‑1 (MQ=21)
                                        gTTCTGCCGCCGCGACATCGGCTGGCGCGCGCTCTGAGACG.AAACGGGGG.GCGCTGCCCCC.CTC.TTGGCCCTC.TGCCCATTCACCCCCCCGGAGATTTCgg                                                         >  2:501699/1‑101 (MQ=21)
                                         ttCTGCCGCCGCGACATCGGCTGGCGCGCGCTCTGAGACG.AAACGGGGG.GCGCTGCCCCC.CTC.TTGGCCTTC.GGCCAATTCACCCCCCGGGATATTTCggg                                                        >  1:2827380/1‑101 (MQ=21)
                                         ttCTGCCGCCGCGACATCGGCTGGCGCGCGCTCTGAGACG.AAACGGGGG.GCGCTGCCCCC.CTC.TTGGCCTTC.GGCCAATTCACCCCCCGGGATATTTCggg                                                        <  2:555777/101‑1 (MQ=21)
                                         ttCTGCCGCCGCGACATCGGCTGGCGCGCGCTCTGAGACG.AAACGGGGG.GCGCCGCCCCC.CTC.TTGGCCTTC.GGCCAATTCACCCCCCCGGAAATTTTggg                                                        >  2:4079914/1‑101 (MQ=21)
                                          tCTGCCGCCGCGACATCGGCTGGCGCGCGCTCTGAGACG.AAACGGGGG.GCGCTGCCCCC.CTC.TTGGCCTTC.GGCCAATTCACCCCCCGGGATATTTCGGGc                                                       <  1:1433966/101‑1 (MQ=21)
                                           cTGCCGCCGCGACATCGGCTGGCGCGCGCTCTGAGACG.AAACGGGGG.GCGCTGCCCCC.CTCTTTGC.CTTC.GGCCAATTCACCCCCCCGGATATTTCGGGCa                                                      >  2:2244777/1‑101 (MQ=21)
                                           cTGCCGCCGCGACATCGGCTGGCGCGCGCTCTGAGACG.AAACGGGGG.GCGCTGCCCCC.CTC.TTGGCCTTC.GGCCAATTCACCCCCCCGGATATTTCGGGCa                                                      >  1:4470570/1‑101 (MQ=21)
                                            tGCCGCCGCGACATCGGCTGGCGCGCGCTCTGAGACGAAAAG.GGGG.GGGCTGCCCCC.CTCTTTG..CCTCTGGCCAATTCACCCCCCCGGGAATTTTCGGCaa                                                     >  1:3995048/1‑101 (MQ=16)
                                             gccgccGCGACATCGGCTGGCGCGCGCTCTGAGACG.AAACGGGGG.GCGCTGCCCCC.CTC.TTGGCCTTC.GGCCAATTCACCCCCCGGGATATTTCGGGCAAg                                                    <  2:2106335/101‑1 (MQ=21)
                                             gccgccGCGACATCGGCTGGCGCGCGCTCTGAGACG.AAACGGGGG.GCGCTGCCCCC.CTC.TTGGCCTTC.GGCCAATTCACCCCCCGGGATATTTCGGGCAAg                                                    <  1:4538410/101‑1 (MQ=21)
                                             gccgccGCGACATCGGCTGGCGCGCGCTCTGAGACG.AAACGGGGG.GCGCTGCCCCC.CTC.TTGGCCTTC.GGCCAATTCACCCCCCGGGATATTTCGGGCAAg                                                    <  1:1027741/101‑1 (MQ=21)
                                             gccgccGCGACATCGGCTGGCGCGCGCTCTGAGACG.AAACGGGGG.GCGCTGCCCCC.CTC.TTGGCCTTC.GGCCAATTCACCCCCCCGGATATTTCGGGCaaa                                                    >  1:1875235/1‑100 (MQ=21)
                                               cgccgcGACATCGGCTGGCGCGCGCTCTGAGACG.AAACGGGGG.GCGCTGCCCCC.CTC.TTGGCCTTC.GGCCAATTCACCCCCCGGGATATTTCGGGCAAGAt                                                  <  2:4296838/101‑1 (MQ=21)
                                                gccgcGACATCGGCTGGCGCGCGCTCTGAGACG.AAACGGGGG.GCGCTGCCCCC.CTC.TTGGCCTTC.GGCCAATTCACCCCCCGGGATATTTCGGGCAAGATa                                                 <  1:4086259/101‑1 (MQ=21)
                                                gccgcGACATCGGCTGGCGCGCGCTCTGAGACG.AAACGGGGG.GCGCTGCCCCC.CTC.TTGGCCTTC.GGCCAATTCACCCCCCGGGATATTTCGGGCAAGATa                                                 <  2:4475018/101‑1 (MQ=21)
                                                 ccgcGACATCGGCTGGCGCGCGCTCTGAGACG.AAACGGGGG.GCGCTGCCCCC.CTC.TTGGCCTTC.GGCCAATTCACCCCCCGGGATATTTCGGGCAAGATaa                                                <  1:1149516/101‑1 (MQ=21)
                                                  cgcgACATCGGCTGGCGCGCGCTCTGAGACG.AAACGGGGG.GCGCCGCCCCC.CCC.TTGGCCCTC.TGCCCAATCACCCCCCAGCGAATTTCCAGCAAAATaaa                                               >  1:3552382/1‑101 (MQ=255)
                                                    cgACATCGGCTGGCGCGCGCTCTGAGACG.AAACGGGGG.GCGCTGCCCCC.CTC.TTGGCCTTC.GGCCAATTCACCCCCCGGGATATTTCGGGCAAGATAaaga                                             >  2:3378349/1‑101 (MQ=21)
                                                    cgACATCGGCTGGCGCGCGCTCTGAGACG.AAACGGGGG.GCGCTGCCCCC.CTC.TTGGCCTTC.GGCCAATTCACCCCCCGGGATATTTCGGGCAAGATAaaga                                             <  1:4128806/101‑1 (MQ=21)
                                                    cgACATCGGCTGGCGCGCGCTCTGAGACG.AAACGGGGG.GCGCTGCCCCC.CTC.TTGGCCTTC.GGCCAATTCACCCCCCGGGATATTTCGGGCAAGATAaaga                                             >  1:2809323/1‑101 (MQ=21)
                                                      aCATCGGCTGGCGCGCGCTCTGAGACG.AAACGGGGG.GCGCTGCCCCC.CTC.TTGGCCTTC.GGCCAATTCACCCCCCGGGATATTTCGGGCAAGATAaagaag                                           <  2:3109085/101‑1 (MQ=21)
                                                          cGGCTGGCGCGCGCTCTGAGACG.AAACGGGGG.GCGCTGCCCCC.CTC.TTGGCCTTC.GGCCAATTCACCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCt                                       <  1:2317065/101‑1 (MQ=21)
                                                           ggCTGGCGCGCGCTCTGAGACG.AAACGGGGG.GCGCTGCCCCC.CTC.TTGGCCTTC.GGCCAATTCACCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCtt                                      <  1:823683/101‑1 (MQ=21)
                                                           ggCTGGCGCGCGCTCTGAGACG.AAACGGGGG.GCGCTGCCCCC.CTC.TTGGCCTTC.GGCCAATTCACCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCtt                                      >  2:1809397/1‑101 (MQ=21)
                                                            gCTGGCGCGCGCTCTGAGACG.AAACGGGGG.GCGCTGCCCCC.CTC.TTGGCCTTC.GGCCAATTCCCCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCttt                                     <  1:3311393/101‑1 (MQ=21)
                                                            gCTGGCGCGCGCGCTGATACG.AAAC.GGGGTGCGCTCCCCCC.CTC.TTGGCCTTC.GGCCAATTCACCCCCCGGGTTATTTCGGGCAAGATAAAGAAGGTCttt                                     <  2:776877/101‑1 (MQ=21)
                                                             cTGGCGCGCGCTCTGAGACG.AAACGGGGG.GCGCTGCCCCC.CTC.TTGGCCTTC.GGCCAATTCACCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCtttt                                    <  2:269921/101‑1 (MQ=21)
                                                             cTGGCGCGCGCTCTGAGACG.AAACGGGGG.GCGCTGCCCCC.CTC.TTGGCCTTC.GGCCAATTCACCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCtttt                                    <  2:3594595/101‑1 (MQ=21)
                                                              tGGCGCGCGCTCTGAGACG.AAACGGGGG.GCGCTGCCCCC.CTC.TTGGCCTTC.GGCCAATTCACCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTc                                   <  2:4561765/101‑1 (MQ=21)
                                                                gcgcgcgcTCTGAGACG.AAACGGGGG.GCGCTGCCCCC.CTC.TTGGCCTTC.GGCCAATTCACCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCt                                 <  1:3271464/101‑1 (MQ=21)
                                                                gcgcgcgcTCTGAGACG.AAACGGGGG.GCGCTGCCCCC.CTC.TTGGCCTTC.GGCCAATTCACCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCt                                 <  1:4637538/101‑1 (MQ=21)
                                                                  gcgcgcTCTGAGACG.AAACGGGGG.GCGCTGCCCCC.CTC.TTGGCCTTC.GGCCAATTCACCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCttt                               >  1:4404563/1‑101 (MQ=21)
                                                                  gcgcgcTCTGAGACG.AAACGGGGG.GCGCTGCCCCC.CTC.TTGGCCTTC.GGCCAATTCACCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCttt                               <  1:2863293/101‑1 (MQ=21)
                                                                   cgcgcTCTGAGACG.AAACGGGGG.GCGCTGCCCCC.CTC.TTGGCCTTC.GGCCAATTCACCCCCCGGGATATTTCGGGCAAGGTAAAGAAGGTCTTTTCCTTTc                              <  1:1909872/101‑1 (MQ=21)
                                                                   cgcgcTCTGAGACG.AAACGGGGG.GCGCTGCCCCC.CTC.TTGGCCTTC.GGCCAATTCACCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTc                              <  2:927318/101‑1 (MQ=21)
                                                                    gcgcTCTGAGACG.AAACGGGGG.GCGCTGCCCCC.CTC.TTGGCCTTC.GGCCAATTCACCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCa                             <  1:2887161/101‑1 (MQ=21)
                                                                    gcgcTCTGAGACG.AAACGGGGG.GCGCTGCCCCC.CTC.TTGGCCTTC.GGCCAATTCACCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCa                             <  2:1107032/101‑1 (MQ=21)
                                                                    gatctctGAGACG.AAACGGGGG.GCGCTGCCCCC.CTC.TTGGCCTTC.GGCCAATTCACCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCa                             <  1:949411/98‑1 (MQ=21)
                                                                     cgcTCTGAGACG.AAACGGGGG.GCGCTGCCCCC.CTC.TTGGCCTTC.GGCCAATTCACCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCAt                            >  1:2541633/1‑101 (MQ=21)
                                                                      gcTCTGAGACG.AAACGGGGG.GCGCTGCCCCC.CTC.TTGGCCTTC.GGCCAATTCACCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATc                           <  2:358972/101‑1 (MQ=21)
                                                                      gcTCTGAGACG.AAACGGGGG.GCGCTGCCCCC.CTC.TTGGCCTTC.GGCCAATTCACCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATc                           <  2:3829338/101‑1 (MQ=21)
                                                                      gcTCTGAGACG.AAACGGGGG.GCGCTGCCCCC.CTC.TTGGCCTTC.GGCCAATTCACCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATc                           <  1:2371067/101‑1 (MQ=21)
                                                                      gcTCTGAGACG.AAACGGGGG.GCGCTGCCCCC.CTC.TTGGCCTTC.GGCCAATTCACCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATc                           <  2:4332674/101‑1 (MQ=21)
                                                                            agaCG.AAACGGGGG.GCGCTGCCCCC.CTC.TTGCCCTTC.GGCCAATTCACCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGccc                     <  1:2658160/101‑1 (MQ=21)
                                                                              aCG.AAACGGGGG.GCGCTGCCCCC.CTC.TTGGCCTTC.GGCCAATTCACCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGa                   >  1:3476249/1‑101 (MQ=21)
                                                                                g.AAACGGGGG.GCGCTGCCCCC.CTC.TTGGCCTTC.GGCCAATTCACCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGaaa                 <  1:433325/101‑1 (MQ=21)
                                                                                  aaacggggg.gCGCTGCCCCC.CTC.TTGGCCTTC.GGCCAATTCACCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAat                <  1:4226586/96‑1 (MQ=21)
                                                                                   aacggggg.gCGCTGCCCCC.CTC.TTGGCCTTC.GGGCAATTCACCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAata               <  1:3970909/97‑1 (MQ=25)
                                                                                   aacggggg.gCGCTGCCCCC.CTC.TTGGCCTTC.GGCCAATTCACCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAata               >  2:2512064/5‑101 (MQ=25)
                                                                                   aacggggg.gCGCTGCCCCC.CTC.TTGGCCTTC.GGCCAATTCACCCCCCGGGATATTTCGGGCAAGATAAAGAAGGGCTTTTCCTTTCATCTTTCCCGAAata               >  1:3728697/5‑101 (MQ=21)
                                                                                     cggggg.gCGCTGCCCCC.CTC.TTGGCCTTC.GGCCAATTCACCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATc             <  1:2141213/99‑1 (MQ=25)
                                                                                     cagggg.gCGCTG.CCCC.CTC.TTGGCCTTT.GGCCAATTCA.CCCCTGGGATATTTCGGGCAAGAT...GAA.ATCATGCCCCTTCATCTTGCCGCAAATATCCCgggg       >  1:1296958/3‑101 (MQ=0)
                                                                                      ggggg.gCGCTGCCCCC.CTC.TTGGCCTTC.GGCCAATTCACCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATcc            <  2:3092070/100‑1 (MQ=25)
                                                                                      agggg.gCGCTGCCCCCTCTT.GGC..CTTT.GGCCAATTCA.CCCCTGGGATATTTCGGGCAAGAT...GAA.ATCATGCCCCTTCATCTTGCCGCAAATATCCCggggg      >  1:3136529/2‑101 (MQ=0)
                                                                                      agggg.gCGCTG.CCCC.CTC.TTGGCCTTT.GGCCAATTCA.CCCCTGGGATATTTCGGGCAAGAT...GAA.ATCATGCCCCTTCATCTTGCCGCAAATATCCCggggg      >  1:3636398/2‑101 (MQ=0)
                                                                                       gggg.gCGCTG.CCCC.CTC.TTGGCCTTT.GGCCAATTCA.CCCCTGGGATATTTCGGGCAAGAT...GAA.ATCATGCCCCTTCATCTTGCCGCAAATATCCCgggggg     >  2:1500855/1‑100 (MQ=0)
                                                                                        ggg.gCGCTGCCCCCTCTT.GGC..CTTT.GGCCAATTCA.CCCCTGGGATATTTCGGGCAAGAT...GAA.ATCATGCCCCTTCATCTTGCCGCAAATATCCCgggggga    <  1:1280857/101‑3 (MQ=0)
                                                                                        ggg.gCGCTGCCCCC.CTC.TTGGCCTTC.GGCCAATTCACCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCCg          <  1:4310501/101‑1 (MQ=25)
                                                                                         gg.gCGCTGCCCCC.CTC.TTGGCCTTC.GGCCAATTCACCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCCgg         <  2:1846739/101‑1 (MQ=37)
                                                                                         gg.gCGCTG.CCCC.CTC.TTGGCCTTT.GGCCAATTCA.CCCCTGGGATATTTCGGGCAAGAT...GAA.ATCATGCCCCTTCATCTTGCCGCAAATATCCCggggggag   >  2:2471942/1‑98 (MQ=2)
                                                                                         gg.gCGCTG.CCCC.CTC.TTGGCCTTT.GGCCAATTCA.CCCCTGGGATATTTCGGGCAAGAT...GAA.ATCATGCCCCTTCATCTTGCCGCAAATATCCCggggggag   >  2:2606470/1‑98 (MQ=2)
                                                                                          g.gCGCTGCCCCC.CTC.TTGGCCTTC.GGCCAATTCACCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCCggg        <  1:1041361/101‑1 (MQ=37)
                                                                                              gcTGCCCCC.CTC.TTGGCCTTC.GGCCAATTCACCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCCGGGGGa     >  1:3659540/1‑101 (MQ=37)
                                                                                                 gCCCCC.CTC.TTGGCCTTC.GGCCAATTCACCCCCCGGGAGATTTCGGGCCAGATAAAAAAGGGCTTTTTCTTTTATCTTTTCCGGAAAATCCCCGGGGGGcc  >  2:1755715/1‑101 (MQ=12)
                                                                                                       |                                                                                                 
GCGCCTATGCGATGGTCGATGCGATCGACTGGCCCGAGGGGTTCTGCCGCCGCGACATCGGCTGGCGCGCGCTCTGAGACG.AAAC.GGGG.GCGCTGCCCCC.CTC.CTG..CTTC.GGCCAATTCACCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCCGGGGGAGCC  >  NC_014034/2566452‑2566644

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