Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_014034 2,566,599 Δ1 bp 5.4% intergenic (+71/+123) 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,5990A→. 5.4% 14.3 ‑0.4 ‑0.0129intergenic (+71/+123)purD/sdaAphosphoribosylamine‑‑glycine ligase/L‑serine ammonia‑lyase
Reads supporting (aligned to +/- strand):  new base (4/3):  ref base (55/67):  total (59/70)
Fisher's exact test strand distribution p-value = 7.02e-01
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 9.45e-01

Only 100 of 153 total aligned reads displayed.

CGCGACATCGGCTGGCGCGCGCTCTGAGACGAAAC.GGGGGCGCTGCCCCC.CTCCTG..CTTCGGCCAATTCACCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCCGGGGGAGCCCGCTTGCGGGCGGGGGCAGCGCCCCCGAAGCGCCCGACGCGCGGCGGGCTTGCCC  >  NC_014034/2566502‑2566699
                                                                                                     |                                                                                                    
cgcgACATCGGCTGGCGCGCGCTCTGAGACGAAACGGGGGGCGCTGCCCCC.CTCTTGGCCTTCGGCCCATTCACCCCCCCGGATATTTCGGGCAAGATaaa                                                                                                      >  1:3906941/1‑101 (MQ=21)
  cgACATCGGCTGGCGCGCGCTCTGAGACGAAACGGGGGGCGCTGCCCCC.CTCTTGGCCTTCGGCCAATTCACCCCCCGGGATATTTCGGGCAAGATAaaga                                                                                                    >  2:3176725/1‑101 (MQ=21)
  cgACATCGGCTGGCGCGCGCTCTGAGACGAAACGGGGGGCGCTGCCCCC.CTCTTGGCCTTCGGCCAATTCACCCCCCGGGATATTTCGGGCAAGATAaaga                                                                                                    <  1:1811536/101‑1 (MQ=21)
   gACATCGGCTGGCGCGCGCTCTGAGACGAAACGGGGGGCGCTGCCCCC.CTCTTGGCCTTCGGCCAATTCACCCCCCGGGATATTTCGGGCAAGATAaagaa                                                                                                   <  1:19493/101‑1 (MQ=21)
      aTCGGCTGGCGCGCGCTCTGAGACGAAACGGGGGGCGCTGCCCCC.CTCTTGGCCTTCGGCCCATTCACCCCCCCGGAAATTTTGGGGAAGATAAAGAAggg                                                                                                >  2:3039499/1‑100 (MQ=18)
       tCGGCTGGCGCGCGCTCTGAGACGAAACGGGGGGCGCTGCCCCC.CTCTTGGCCTTCGGCCAATTCACCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTc                                                                                               >  1:3848117/1‑101 (MQ=21)
         ggCTGGGGCGCGCTCTGAGACGAAACGGGGGGCGCTGCCCCC.CTCTTGGCCTTCGGCCAATTCACCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCtt                                                                                             <  1:3073594/101‑1 (MQ=21)
          gCTGGCGCGCGCTCTGAGACGAAACGGGGGGCGCTGCCCCC.CTCTTGGCCTTCGGCCAATTCACCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCttt                                                                                            >  1:2196690/1‑101 (MQ=21)
          gCTGGCGCGCGCTCTGAGACGAAACGGGGGGCGCTGCCCCC.CTCTTGGCCTTCGGCCAATTCACCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCttt                                                                                            <  2:1078597/101‑1 (MQ=21)
           cTGGCGCGCGCTCTGAGACGAAACGGGGGGCGCTGCCCCC.CTCTTGGCCTTCGGCCAATTCACCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCtttt                                                                                           <  1:2327547/101‑1 (MQ=21)
           cTGGCGCGCGCTCTGAGACGAAACGGGGGGCGCTGCCCCC.CTCTTGGCCTTCGGCCAATTCACCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCtttt                                                                                           <  1:2169706/101‑1 (MQ=21)
            tGGCGCGCGCTCTGAGACGAAACGGGGGGCGCTGCCCCC.CTCTTGGCCTTCGGCCAATTCACCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTc                                                                                          >  2:1009588/1‑101 (MQ=21)
              gcgcgcgcTCTGAGACGAAACGGGGGGCGCTGCCCCC.CTCTTGGCCTTCGGCCAATTCACCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCt                                                                                        <  1:3727302/101‑1 (MQ=21)
              gcgcgcgcTCTGAGACGAAACGGGGGGCGCTGCCCCC.CTCTTGGCCTTCGGCCAATTCACCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCt                                                                                        <  1:3376059/101‑1 (MQ=21)
               cgcgcgcTCTGAGACGAAACGGGGGGCGCTGCCCCC.CTCTTGGCCTTCGGCCAATTCACCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCtt                                                                                       <  1:1533020/101‑1 (MQ=21)
               cgcgcgcTCTGAGACGAAACGGGGGGCGCTGCCCCC.CTCTTGGCCTTCGGCCAATTCACCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCtt                                                                                       <  2:3626458/101‑1 (MQ=21)
               cgcgcgcTCTGAGACGAAACGGGGGGCGCTGCCCCC.CTCTTGGCCTTCGGCCAATTCACCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCtt                                                                                       <  1:2512408/101‑1 (MQ=21)
                gcgcgcTCTGAGACGAAACGGGGGGCGCTGCCCCC.CTCTTGGCCTTCGGCCAATTCCCCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCttt                                                                                      <  2:823803/101‑1 (MQ=21)
                gcgcgcTCTGAGACGAAACGGGGGGCGCTGCCCCC.CTCTTGGCCTTCGGCCAATTCACCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCttt                                                                                      <  2:1830164/101‑1 (MQ=21)
                gcgcgcTCTGAGACGAAACGGGGGGCGCTGCCCCC.CTCTTGGCCTTCGGCCAATTCACCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCttt                                                                                      <  2:2428260/101‑1 (MQ=21)
                gcgcgcTCTGAGACGAAACGGGGGGCGCTGCCCCC.CTCTTGGCCTTCGGCCAATTCACCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCttt                                                                                      >  2:3297536/1‑101 (MQ=21)
                 cgcgcTCTGAGACGAAACGGGGGGCGCTGCCCCC.CTCTTGGCCTTCGGCCAATTCACCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTc                                                                                     <  1:2076249/101‑1 (MQ=21)
                 cgcgcTCTGAGACGAAACGGGGGGCGCTGCCCCC.CTCTTGGCCTTCGGCCAATTCACCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTc                                                                                     <  1:509035/101‑1 (MQ=21)
                   cgcTCTGAGACGAAACGGGGGGCGCTGCCCCC.CTCTTGGCCTTCGGCCAATTCACCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCAt                                                                                   <  1:4145111/101‑1 (MQ=21)
                    gcTCTGAGACGAAACGGGGGGCGCTGCCCCC.CTCTTGGCCTTCGGCCAATTCACCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATc                                                                                  <  1:31132/101‑1 (MQ=21)
                      tctGAGACGAAACGGGGGGCGCTGCCCCC.CTCTTGGCCTTCGGCCAATTCACCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCtt                                                                                <  1:1173596/101‑1 (MQ=21)
                      tctGAGACGAAACGGGGGGCGCTGCCCCC.CTCTTGGCCTTCGGCCAATTCACCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCtt                                                                                <  1:298907/101‑1 (MQ=21)
                       ctGAGACGAAACGGGGGGCGCTGCCCCC.CTCTTGGCCTTCGGCCAATTCACCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTg                                                                               <  1:569431/101‑1 (MQ=21)
                        tGAGACGAAACGGGGGGCGCTGCCCCC.CTCTTGGCCTTCGGCCAATTCACCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTGTTTTTTTTTCTTCTTGc                                                                              >  2:853783/1‑101 (MQ=18)
                        tGAGACGAAACGGGGGGCGCTGCCCCC.CTCTTGGCCTTCGGCCAATTCACCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGc                                                                              >  2:2272230/1‑101 (MQ=21)
                         gagaCGAAACGGGGGGCGCTGCCCCC.CTCTTGGCCTTCGGCCAATTCACCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGcc                                                                             >  1:1222553/1‑101 (MQ=21)
                          agaCGAAACGGGGGGCGCTGCCCCC.CTCTTGGCCTTCGGCCAATTCACCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGccc                                                                            <  2:236725/101‑1 (MQ=21)
                            aCGAAACGGGGGGCGCTGCCCCC.CTCTTGGCCTTCGGCCAATTCACCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGa                                                                          >  1:1787489/1‑101 (MQ=21)
                              gAAACGGGGGGCGCTGCCCCC.CTCTTGGCCTTCGGCCAATTCACCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGaaa                                                                        <  2:2702264/101‑1 (MQ=21)
                              gAAACGGGGGGCGCTGCCCCC.CTCTTGGCCTTCGGCCAATTCACCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGaaa                                                                        <  2:1593263/101‑1 (MQ=21)
                               aaacggggggCGCTGCCCCC.CTCTTGGCCTTCGGCCAATTCACCCCCCGGGATATTTCGGGCAAGATAAAGAAGGGCTTTTCCTTTCATCTTGcccccaat                                                                       >  2:4223765/6‑96 (MQ=18)
                                aacggggggCGCTGCCCCC.CTCTTGGCCTTCGGCCAATTCACCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAata                                                                      <  1:3581109/97‑1 (MQ=25)
                                  cggggggCGCTGCCCCC.CTCTTGGCCTTCGGCCAATTCACCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATc                                                                    <  2:4015211/99‑1 (MQ=25)
                                   ggggggCGCTGCCCCC.CTCTTGGCCTTCGGCCAATTCACCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATcc                                                                   <  2:1571825/100‑1 (MQ=25)
                                   agggggCGCTG.CCCC.CTCTTGGCCTTTGGCCAATTCA.CCCCTGGGATATTTCGGGCAAGAT...GAA.ATCATGCCCCTTCATCTTGCCGCAAATATCCCggggg                                                             >  1:4076720/2‑101 (MQ=0)
                                     ggggCGCTGCCCCCTCTTGGC..CTTTGGCCAATTCA.CCCCTGGGATATTTCGGGCAAGAT...GAA.ATCATGCCCCTTCATCTTGCCGCAAATATCCCgggggga                                                           <  1:3959563/101‑3 (MQ=2)
                                      gggCGCTGCCCCCTCTTGGC..CTTTGGCCAATTCA.CCCCTGGGATATTTCGGGCAAGAT...GAA.ATCATGCCCCTTCATCTTGCCGCAAATATCCCggggggag                                                          >  1:601326/1‑98 (MQ=2)
                                       ggCGCTGCCCCCTCTTGGC..CTTTGGCCAATTCA.CCCCTGGGATATTTCGGGCAAGAT...GAA.ATCATGCCCCTTCATCTTGCCGCAAATATCCCggggggagc                                                         <  1:1048067/101‑5 (MQ=2)
                                       ggCGCTGCCCCCTCTTGGC..CTTTGGCCAATTCA.CCCCTGGGATATTTCGGGCAAGAT...GAA.ATCATGCCCCTTCATCTTGCCGCAAATATCCCggggggagc                                                         <  1:2334609/101‑5 (MQ=2)
                                         cgcTGCCCCC.CTCTTGGCCTTCGGCCAATTCACCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCCAAATATCCCggggg                                                             >  2:3022218/1‑101 (MQ=25)
                                          gcTGCCCCC.CTCTTGGCCTTCGGCCAATTCTCCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCCGGGGGa                                                            <  1:2866857/101‑1 (MQ=37)
                                           cTGCCCCC.CTCTTGGCCTTCGGCCAATTCACCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCCGGGGGAg                                                           <  2:1142532/101‑1 (MQ=37)
                                                  ccctcttggccTTCGGCCAATTCCCCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCCGGGGGAGCCCGCt                                                     <  1:1052427/91‑1 (MQ=38)
                                                  ccctcttggccTTCGGCCAATTCACCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCCGGGGGAGCCCGCt                                                     <  1:93659/91‑1 (MQ=255)
                                                    ctcttggccTTCGGCCAATTCCCCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCCGGGGGAGCCCGCTTg                                                   <  2:2117738/93‑1 (MQ=255)
                                                      cttggccTTCGGCCAATTCACCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCCGGGGGAGCCCGCTTgcg                                                 >  2:3746717/7‑101 (MQ=255)
                                                       ttggccTTCGGCCAATTCACCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCCGGGGGAGCCCGCTTgcgg                                                <  1:3070202/96‑1 (MQ=255)
                                                       ttggccTTCGGCCAATTCACCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCCGGGGGAGCCCGCTTgcgg                                                <  2:2553261/96‑1 (MQ=255)
                                                        tggccTTCGGCCAATTCACCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCCGGGGGAGACCGCTTgccgg                                               >  2:1792567/5‑101 (MQ=255)
                                                         ggccTTCGGCCAATTCACCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCCGGGGGAGCCCGCTTgcgggc                                              >  1:1210820/4‑101 (MQ=255)
                                                          gccTTCGGCCAATTCACCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCCGGGGGAGCCCGCTTgcgggag                                             >  2:3850891/3‑99 (MQ=255)
                                                           ccTTCGGCCAATTCACCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCCGGGGGAGCCCGCTTgcgggcag                                            >  2:1717474/2‑99 (MQ=255)
                                                               cGGCCAATTCACCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCCGGGGGAGCCCGCTTGCGGGCGGGGGc                                        <  2:2144179/101‑1 (MQ=255)
                                                                ggCCAATTCACCCCCCGGGATATTTTGGGCGAGATAAAGAAGGGGTTTTTCTTTTCTCTTTCCCGAAAAATCACGGGGGGGGCCGCCTCGTGGCggggggg                                       >  2:756706/1‑99 (MQ=16)
                                                                 gCCAATTCACCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCCGGGGGAGCCCGCTTGCGGGCGGGGGCAg                                      <  2:1958095/101‑1 (MQ=255)
                                                                 gCCAATTCACCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCCGGGGGAGCCCGCTTGCGGGCGGGGGCAg                                      <  1:3377084/101‑1 (MQ=255)
                                                                   cAAGTCACCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCCGGGGGAGCCCGCTTGCGGGCGGGGGCAgcg                                    <  1:1876199/101‑1 (MQ=255)
                                                                     aTTCACCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCCGGGGGAGCCCGCTTGCGGGCGGGGGCAGCGcc                                  >  2:1662111/1‑101 (MQ=255)
                                                                     aTTCACCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCCGGGGGAGCCCGCTTGCGGGCGGGGGCAGCGcc                                  >  1:3091680/1‑101 (MQ=255)
                                                                     aTTCACCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCCGGGGGAGCCCGCTTGCGGGCGGGGGCAGCGcc                                  >  1:2704649/1‑101 (MQ=255)
                                                                      ttCACCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCCGGGGGAGCCCGCTTGCGGGCGGGGGCAGCGccc                                 >  1:652542/1‑101 (MQ=255)
                                                                          ccccccGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCCGGGGGAGCCCGCTTGCGGGCGGGGGCCGCGCCCCCGa                             >  2:2440648/1‑101 (MQ=255)
                                                                           cccccGAGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCCGGGGGAGCCCGCTTGCGGGCGGGGGCAGCGCCCCCGaa                            <  2:2040485/101‑1 (MQ=255)
                                                                             cccGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCCGGGGGAGCCCGCTTGCGGGCGGGGGCAGCGCCCCCGAAgc                          <  1:88752/101‑1 (MQ=255)
                                                                              ccGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCCGGGGGAGCCCGCTTGCGGGCGGGGGCAGCGCCCCCGAAgcg                         >  2:1190061/1‑101 (MQ=255)
                                                                              ccGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCCGGGGGAGCCCGCTTGCGGGCGGGGGCAGCGCCCCCGAAgcg                         <  1:3962109/101‑1 (MQ=255)
                                                                               cGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCCGGGGGAGCCCGCTTGCGGGCGGGGGCAGCGCCCCCGAAgcgc                        >  1:1676698/1‑101 (MQ=255)
                                                                               cGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCCGGGGGAGCCCGCTTGCGGGCGGGGGCAGCGCCCCCGAAgcgc                        >  1:946245/1‑101 (MQ=255)
                                                                                gggATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCCGGGGGAGCCCGCTTGCGGGCGGGGGCAGCGCCCCCGAAGCGcc                       <  1:3779274/101‑1 (MQ=255)
                                                                                 ggATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCCGGGGGAGCCCGCTTGCGGGCGGGGGCAGCGCCCCCGAAGCGccc                      >  1:3557518/1‑101 (MQ=255)
                                                                                 ggATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCCGGGGAAGCCCGCTTGCGGGCGGGGGCAGCGCCCCCGAAGCGccc                      <  2:680417/101‑1 (MQ=255)
                                                                                  gATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCCGGGGGAGCCCGCTTGCGGGCGGGGGCAGCGCCCCCGAAGCGCCCg                     >  2:579481/1‑101 (MQ=255)
                                                                                   atatTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCCGGGGGAGCCCGCTTGCGGGCGGGGGCAGCGCCCCCGAAGCGCCCGa                    >  1:3851558/1‑101 (MQ=255)
                                                                                     atTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCCGGGGGAGCCCGCTTGCGGGCGGGGGCAGCGCCCCCGAAGCGCCCGAcg                  >  1:2516173/1‑101 (MQ=255)
                                                                                     atTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCCGGGGGAGCCCGCTTGCGGGCGGGGGCAGCGCCCCCGAAGCGCCCGAcg                  >  1:323970/1‑101 (MQ=255)
                                                                                     atTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCCGGGGGAGCCCGCTTGCGGGCGGGGGCAGCGCCCCCGAAGCGCCCGAcg                  >  1:1927029/1‑101 (MQ=255)
                                                                                      tttCGGACAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCCGGGGGAGCCCGCTTGCGGGCGGGGGCAGCGCCCCCGAAGCGCCCGAcgc                 <  2:436788/101‑1 (MQ=255)
                                                                                       ttCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCCGGGGGAGCCCGCTTGCGGGCGGGGGCAGCGCCCCCGAAGCGCCCGAcgcg                <  2:3428342/101‑1 (MQ=255)
                                                                                       ttCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCCGGGGGAGCCCGCTTGCGGGCGGGGGCAGCGCCCCCGAAGCGCCCGAcgcg                <  1:2683626/101‑1 (MQ=255)
                                                                                         cGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCCGGGGGAGCCCGCTTGCGGGCGGGGGCAGCGCCCCCGAAGCGCCCGAcgcgcg              >  1:3450870/1‑101 (MQ=255)
                                                                                          gggCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCCGGGGGAGCCCGCTTGCGGGCGGGGGCAGCGCCCCCGAAGCGCCCGACGCgcgg             >  2:3547478/1‑101 (MQ=255)
                                                                                           ggCAAGTTAAAGAAGGTTTTTTCCTTTCTTCTTGCCCGAAATATCCCGGGGGATCCCGCTTGCGGGCGGGGGCAGCGCCCCCGAAGCGCCCGACGCgcggc            <  1:3282303/101‑1 (MQ=39)
                                                                                           ggCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCCGGGGGAGCCCGCTTGCGGGCGGGGGCAGCGCCCCCGAAGCGCCCGACGCgcggc            >  2:185544/1‑101 (MQ=255)
                                                                                              aaGATAAAGAAGGTCTTTTCTTTTCATCTTGCCCGAAATATCCCGGGGGAGCCCGCTTGCGGGCGGGGGCAGCGCCCCCGAAGCGCCCGACGCGCGGCggg         <  1:4021231/101‑1 (MQ=255)
                                                                                              aaGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCCGGGGGAGCCCGCTTGCGGGCGGGGGCAGCGCCCCCGAAGCGCCCGACGCGCGGCggg         >  2:2700590/1‑101 (MQ=255)
                                                                                              aaGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCCGGGGGAGCCCGCTTGCGGGCGGGGGCAGCGCCCCCGAAGCGCCCGACGCGCGGCggg         >  1:117707/1‑101 (MQ=255)
                                                                                               aGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCCGGGGGAGCCCGCTTGCGGGCGGGGGCAGCGCCCCCGAAGCGCCCGACGCGCGGCGGGc        >  1:203708/1‑101 (MQ=255)
                                                                                                  tAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCCGGGGGAGCCCGCTTGCGGGCGGGGGCAGCGCCCCCGAAGCGCCCGACGCGCGGCGGGGttc     >  1:3181312/1‑100 (MQ=255)
                                                                                                   aaaGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCCGGGGGAGCCCGCTTGCGGGCGGGGGCAGCGCCCCCGAAGCGCCCGACGCGCGGCGGGCTTgc    >  2:3147613/1‑101 (MQ=255)
                                                                                                   aaaGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCCGGGGGAGCCCGCTTGCGGGCGGGGGCAGCGCCCCCGAAGCGCCCGACGCGCGGCGGGCTTgc    >  1:2248743/1‑101 (MQ=255)
                                                                                                   aaaGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCCGGGGGAGCCCGCTTGCGGGCGGGGGCAGCGCCCCCGAAGCGCCCGACGCGCGGCGGGCTTgc    >  1:794897/1‑101 (MQ=255)
                                                                                                   aaaGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCCGGGGGAGCCCGCTTGCGGGCGGGGGCAGCGCCCCCGAAGCGCCCGACGCGCGGCGGGCTTgc    >  1:1983891/1‑101 (MQ=255)
                                                                                                    aagaagGTCTTTTCCTTTCATCTTGCCCGAAATATCCCGGGGGAGCCCGCTTGCGGGCGGGGGCAGCGCCCCCGAAGCGCCCGACGCGCGGCGGGCTTgcc   >  1:3421357/1‑101 (MQ=255)
                                                                                                    aagaagGTCTTTTCCTTTCATCTTGCCCGAAATATCCCGGGGGAGCCCGCTTGCGGGCGGGGGCAGCGCCCCCGAAGCGCCCGACGCGCGGCGGGCTTccc   >  2:2561971/1‑101 (MQ=255)
                                                                                                    aagaagGTCTTTTCCTTTCATCTTGCCCGAAATATCCCGGGGGAGCCCGCTTGCGGGCGGGGGCAGCGCCCCCGAAGCGCCCGACGCGCGGCGGGCTTccc   >  1:3903167/1‑101 (MQ=255)
                                                                                                     |                                                                                                    
CGCGACATCGGCTGGCGCGCGCTCTGAGACGAAAC.GGGGGCGCTGCCCCC.CTCCTG..CTTCGGCCAATTCACCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCCGGGGGAGCCCGCTTGCGGGCGGGGGCAGCGCCCCCGAAGCGCCCGACGCGCGGCGGGCTTGCCC  >  NC_014034/2566502‑2566699

Base quality scores:  ATCG  < 3 ≤  ATCG  < 17 ≤  ATCG  < 26 ≤  ATCG  < 33 ≤  ATCG  < 41 ≤  ATCG