Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_014034 2,566,599 Δ1 bp 3.2% 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→. 3.2% 7.5 ‑0.4 ‑0.0156intergenic (+71/+123)purD/sdaAphosphoribosylamine‑‑glycine ligase/L‑serine ammonia‑lyase
Reads supporting (aligned to +/- strand):  new base (3/2):  ref base (69/82):  total (72/84)
Fisher's exact test strand distribution p-value = 6.63e-01
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 9.94e-01

Only 100 of 187 total aligned reads displayed.

CGCGACATCGGCTGGCGCGCGCTCTGAGACG.AAAC.GGGGGCGCTGCCCCC.CTCCTG..CTTCGGCCAATTCA.CCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCCG.GGGGAGCCCGCTTGCGGGCGGGGGCAGCGCCCCCGAAGCGCCCGACGCGCGGCGGGCTTGCCC  >  NC_014034/2566502‑2566699
                                                                                                       |                                                                                                     
cgcgACATCGGCTGGCGCGCGCTCTGAGACG.AAACGGGGGGCGCTGCCCCC.CTCTTGGCCTTCGGCCAATTCA.CCCCCCGGGATATTTCGGGCAAGATaaa                                                                                                       <  2:1980484/101‑1 (MQ=21)
cgcgACATCGGCTGGCGCGCGCTCTGAGACG.AAACGGGGGGCGCTGCCCCC.CTCTTGGCCTTCGGCCAATTCA.CCCCCCGGGATATTTCGGGCAAGATaaa                                                                                                       <  1:2582465/101‑1 (MQ=21)
 tcgACATCGGCTGGCGCGCGCTCTGAGACG.AAACGGGGGGCGCTGCCCCC.CTCTTGGCCTTCGGCCAATTCA.CCCCCCGGGATATCTCGGGCAAGATAaag                                                                                                      <  1:3750265/100‑1 (MQ=21)
 gcgACATCGGCTGGCGCGCGCTCTGAGACGAAAAG.GGGGGCGCTGCCCCC.CTCTTGGCCTTCCGCCAATTCA.CCCCCCCGGGTATTTCGGGGAAGATaaaa                                                                                                      >  2:3461737/1‑100 (MQ=18)
 gcgACATCGGCTGGCGCGCGCTCTGAGACG.AAACGGGGGGCGCTGCCCCC.CTCTTGGCCTTCGGTCAATTCA.CCCCCCGGGATATTTCGGGCAAGATAaag                                                                                                      <  1:1358392/101‑1 (MQ=21)
    aCATCGGCTGGCGCGCGCTCTGAGACG.AAACGGGGGGCGCTGCCCCC.CTCTTGGCCTTCGGCCAATTCA.CCCCCCGGGATATTTCGGGCAAGATAaaaaag                                                                                                   >  1:1271813/1‑101 (MQ=21)
    aCATCGGCTGGCGCGCGCTCTGAGACG.AAACGGGGGGCGCTGCCCCC.CTCTTGGCCTTCGGCCAATTCA.CCCCCCCGGAAAATTCGGGCAAGATAaagaag                                                                                                   >  1:2317880/1‑101 (MQ=21)
     cATCGGCTGGCGCGCGCTCTGAGACG.AAACGGGGGGCGCTGCCCCC.CTCTTGGCCTTCGGCCAATTCA.CCCCCCGGGATATTTCGGGCAAGATAAAGAAgg                                                                                                  <  1:215566/101‑1 (MQ=21)
     cATCGGCTGGCGCGCGCTCTGAGACG.AAACGGGGGGCGCTGCCCCC.CTCTTGGCCTTCGGCCAATTCA.CCCCCCGGGATATTTCGGGCAAGATAAAGAAgg                                                                                                  <  1:2985799/101‑1 (MQ=21)
        cGGCTGGCGCGCGCTCTGAGACG.AAACGGGGGGCGCTGCCCCC.CTCTTGGCCTTCGGCCAATTCA.CCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCt                                                                                               <  2:1414892/101‑1 (MQ=21)
         ggCTGGCGCGCGCTCTGAGACG.AAACGGGGGGCGCTGCCCCC.CTCTTGGCCTTCGGCCAATTCA.CCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCtt                                                                                              <  2:1394035/101‑1 (MQ=21)
           cTGGCGCGCGCTCTGAGACG.AAACGGGGGGCGCTGCCCCC.CTCTTGGCCTTCGGCCAATTCA.CCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCtttt                                                                                            <  2:2134379/101‑1 (MQ=21)
           cTGGCGCGCGCTCTGAGACG.AAACGGGGGGCGCTGCCCCC.CTCTTGGCCTTCGGCCAATTCA.CCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCtttt                                                                                            <  2:327401/101‑1 (MQ=21)
           cTGGCGCGCGCTCTGAGACG.AAACGGGGGGCGCTGCCCCC.CTCTTGGCCTTCGGCCAATTCA.CCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCtttt                                                                                            <  1:4081669/101‑1 (MQ=21)
            tGGCGCGCGCTCTGAGACG.AAACGGGGGGCGCTGCCCCC.CTCTTGGCCTTCGGCCAATTCA.CCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTc                                                                                           >  2:3361660/1‑101 (MQ=21)
              gcgcgcgcTCTGAGACG.AAACGGGGGGCGCTGCCCCC.CTCTTGGCCTTCGGCCAATTCA.CCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCt                                                                                         <  1:4349531/101‑1 (MQ=21)
               cgcgcgcTCTGAGACG.AAACGGGGGGCGCTGCCCCC.CTCTTGGCCTTCGGCCAATTCA.CCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCtt                                                                                        <  2:3926459/101‑1 (MQ=21)
                gcgcgcTCTGAGACG.AAACGGGGGGCGCTGCCCCC.CTCTTGGCCTTCGGCCAATTCA.CCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCttt                                                                                       <  1:2359893/101‑1 (MQ=21)
                 cgcgcTCTGAGACG.AAACGGGGGGCGCTGCCCCC.CTCTTGGCCTTCGGCCAATTCA.CCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTc                                                                                      <  1:1807017/101‑1 (MQ=21)
                 cgcgcTCTGAGACG.AAACGGGGGGCGCTGCCCCC.CTCTTGGCCTTCGGCCAATTCA.CCCCCCGGGATATTTCGGGCAAGATAAAGAAGGGCTTTTCCTTTc                                                                                      >  2:413218/1‑101 (MQ=21)
                  gcgcTCTGAGACG.AAACGGGGGGCGCTGCCCCC.CTCTTGGCCTTCGGCCAATTCA.CCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCa                                                                                     <  2:124700/101‑1 (MQ=21)
                   cgcTCTGAGACG.AAACGGGGGGCGCTGCCCCC.CTCTTGGCCTTCGGCCAATTCACCCCCCCGGGAAATTTCGGGCAAAAAAAAAAAGGTTTTTTCTttttt                                                                                     >  1:3164572/1‑99 (MQ=12)
                    gcTCTGAGACG.AAACGGGGGGCGCTGCCCCC.CTCTTGGCCTTCGGCCAATTCC.CCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATc                                                                                   <  1:817442/101‑1 (MQ=21)
                       ctGAGACG.AAACGGGGGGCGCTGCCCCC.CTCTTGGCCTTCGGCAAATTCG.CCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTg                                                                                <  1:4265530/101‑1 (MQ=21)
                        tGAGACG.AAACGGGGGGCGCTGCCCCC.CTCTTGGCCTTCGGCCAATTCA.CCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGc                                                                               <  2:1574485/101‑1 (MQ=21)
                        tGAGACG.AAACGGGGGGCGCTGCCCCC.CTCTTGGCCTTCGGCCAATTCA.CCCCCCGGGAAATTTCCTGCAAGATAAAGAAGGTCTTTTTCTTTCATCTTGc                                                                               >  1:495608/1‑101 (MQ=18)
                         gagaCG.AAACGGGGGGCGCTGCCCCC.CTCTTGGCCTTCGGCCAATTCA.CCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGcc                                                                              <  2:271613/101‑1 (MQ=21)
                          agaCG.AAACGGGGGGCGCTGCCCCC.CTCTTGGCCTTCGGCCAATTCA.CCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGccc                                                                             <  2:4443182/101‑1 (MQ=21)
                           gaCG.AAACGGGGGGCGCTGCCCCC.CTCTTGGCCTTCGGCCAATTCA.CCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCg                                                                            >  1:1777787/1‑101 (MQ=21)
                             cG.AAACGGGGGGCGCTGCCCCC.CTCTTGGCCTTCGGCCAATTCA.CCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGaa                                                                          >  1:3643487/1‑101 (MQ=21)
                              g.AAACGGGGGGCGCTGCCCCC.CTCTTGGCCTTCGGCCCATTCA.CCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGaaa                                                                         <  1:3107135/101‑1 (MQ=21)
                              g.AAACGGGGGGCGCTGCCCCC.CTCTTGGCCTTCGGCCAATTCA.CCCCCCGGGATATTTCGGGAAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGaaa                                                                         >  2:1377123/1‑101 (MQ=21)
                                aaacggggggCGCTGCCCCC.CTCTTGGCCTTCGGCCAATTCA.CCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAat                                                                        <  1:661025/96‑1 (MQ=21)
                                 aacggggggCGCTGCCCCC.CTCTTGGCCTTCGGCCAATTCA.CCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAata                                                                       <  1:3403555/97‑1 (MQ=25)
                                 aacggggggCGCTGCCCCC.CTCTTGGCCTTCGGCCAATTCA.CCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAata                                                                       >  1:2794101/5‑101 (MQ=25)
                                  acggggggCGCTGCCCCC.CTCTTGGCCTTCGGCCAATTCA.CCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAatat                                                                      >  2:4193579/4‑101 (MQ=25)
                                  acggggggCGCTGCCCCC.CTCTTGGCCTTCGGCCAATTCA.CCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAatat                                                                      <  1:1158219/98‑1 (MQ=25)
                                    agggggCGCTG.CCCC.CTCTTGGCCTTTGGCCAATTCA..CCCCTGGGATATTTCGGGCAAGAT...GAA.ATCATGCCCCTTCATCTTGCCGCAAATATCCCg.gggg                                                             >  1:3465033/2‑101 (MQ=0)
                                     gggggCGCTGCCCCCTCTTGGC..CTTTGGCCAATTCA..CCCCTGGGATATTTCGGGCAAGAT...GAA.ATCATGCCCCTTCATCTTGCCGCAAATATCCCg.ggggg                                                            <  1:3971574/101‑2 (MQ=0)
                                     gggggCGCTGCCCCC.CTCTTGGCCTTCGGCCAATTCA.CCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATccc                                                                   <  1:2729715/101‑1 (MQ=25)
                                      ggggCGCTGCCCCCTCTTGGC..CTTTGGCCAATTCA..CCCCTGGGATATTTCGGGCAAGAT...GAA.ATCATGCCCCTTCATCTTGCCGCAAATATCCCg.ggggga                                                           <  2:2846502/101‑3 (MQ=2)
                                         gcgcTGCCCCC.CTCTTGGCCTTCGGCCAATTCA.CCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCCg.ggg                                                              <  1:550421/101‑1 (MQ=37)
                                          cgcTGCCCCC.CTCTTGGCCTTCGGCCAATTCA.CCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCCg.gggg                                                             <  2:4078109/101‑1 (MQ=37)
                                           gcTGCCCCC.CTCTTGGCCTTCGGCCAATTCA.CCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCGTTTCATCTTGCCCGAAATATCCCG.GGGGa                                                            <  2:1102050/101‑1 (MQ=37)
                                           gcTGCCCCC.CTCTTGGCCTTCGGCCAATTCA.CCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCCG.GGGGa                                                            <  1:2932405/101‑1 (MQ=37)
                                           gcTGCCCCC.CTCTTGGCCTTCGGCCAATTCA.CCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCCG.GGGGa                                                            <  1:4128436/101‑1 (MQ=37)
                                             tGCCCCC.CTCTTGGCCTTCGGCCAATTCA.CCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCCG.GGGGAGc                                                          >  2:1008137/1‑101 (MQ=37)
                                             tGCCCCC.CTCTTGGCCTTCGGCCAATTCA.CCCCCCGGGAAATTTCGGGGCAGATAAAGAGGGTTTTTTCCTTTCTTCTTTCCCGGAATATCCCC.CGGGacc                                                          >  1:192838/1‑99 (MQ=16)
                                               ccccc.cTCTGG..CTTC.GCC.ATTCA.CCCCCCGGGATATTTCGGGCACGATGA..AAGG.CATCCGGTTTCATCTTGCTGCAAATATCCCGAGGGGAGGCCGATTgc                                                  >  2:2946351/1‑101 (MQ=18)
                                                   ccctcttggccTTCGGCCAATTCC.CCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCCG.GGGGAGCCCGCt                                                     <  1:3591338/91‑1 (MQ=38)
                                                     ctcttggccTTCGGCCAATTCA.CCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCCG.GGGGAGCCCGCTTg                                                   <  2:1699336/93‑1 (MQ=255)
                                                     ctcttggccTTCGGCCAATTCA.CCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCCG.GGGGAGCCCGCTTg                                                   <  1:281610/93‑1 (MQ=255)
                                                       cttggccTTCGGCCAATTCA.CCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCCG.GGGGAGCCCGCTTgcg                                                 <  1:3958514/95‑1 (MQ=255)
                                                         tggccTTCGGCCAATTCA.CCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCCG.GGGGAGCCCGCTTgcggg                                               >  1:2972992/5‑101 (MQ=255)
                                                         tggccTTCGGCCAATTCA.CCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCCG.GGGGAGCCCGCTTgcggg                                               <  1:747769/97‑1 (MQ=255)
                                                          ggccTTCGGCCAATTCA.CCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCCG.GGGGAGCCCGCCTgccggc                                              >  2:2330775/4‑101 (MQ=255)
                                                             cTTCGGGCAATTCA.CCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCCG.GGGGAGCCCGCTTGcgggcggg                                           <  1:725245/101‑1 (MQ=255)
                                                                cGGCCAATTCA.CCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCCG.GGGGAGCCCGCTTGCGGGCGGGGGc                                        <  1:329854/101‑1 (MQ=255)
                                                                 ggCCAATTCG.CCCCCCGGGATATTTCGGGCGAGATAAAGAAGGTCTGTTCCTTTCATCTTGCCCGAAATATCCCG.GGGGAGCCCGCTTGCGGGCGGGGGCa                                       <  2:3689381/101‑1 (MQ=255)
                                                                     aaTTCA.CCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCCG.GGGGAGCCCGCTTGCGGGCGGGGGCAgcgc                                   >  1:1584280/1‑101 (MQ=255)
                                                                        tCA.CCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCCG.GGGGAGCCCGCTTGCGGGCGGGGGCAGCGcccc                                <  2:2266890/101‑1 (MQ=255)
                                                                        tCA.CCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCCG.GGGGAGCCCGCTTGCGGGCGGGGGCAGCGcccc                                >  2:3359653/1‑101 (MQ=255)
                                                                         cA.CCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTTCCCGAAATATCCCG.GGGGACCCCGCCTGCGGGCGGGGGCCGGGccccc                               >  1:209796/1‑101 (MQ=37)
                                                                         cA.CCCCCCGGGATATTTCGGGCAAGATAAAGAAGGGCTTTTTCTTTCATCTTGCCCGAAATATCCCG.GGGGGGCCCGCTTGCGGGGGGGGGGCGCGccccc                               >  2:3072127/1‑101 (MQ=21)
                                                                          cgccccccGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCCG.GGGGAGCCCGCTTGCGGGCGGGGGCAGCGccccc                               <  1:1738269/99‑1 (MQ=255)
                                                                             cccccGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCCG.GGGGAGCCCGCTTGCGGGCGGGGGCAGCGCCCCCGaa                            <  2:2318409/101‑1 (MQ=255)
                                                                             cccccGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCCG.GGGGAGCCCGCTTGCGGGCGGGGGCAGCGCCCCCGaa                            >  2:2669688/1‑101 (MQ=255)
                                                                                 cGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCCG.GGGGAGCCCGCTTGCGGGCGGGGGCAGCGCCCCCGAAgcgc                        >  2:3884936/1‑101 (MQ=255)
                                                                                 cGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCCG.GGGGAGCCCGCTTGCGGGCGGGGGCAGCGCCCCCGAAgcgc                        >  2:4149359/1‑101 (MQ=255)
                                                                                    gATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCCG.GGGGAGCCCGCTTGCGGGCGGGGGCAGCGCCCCCGAAGCGCCCg                     >  2:3288299/1‑101 (MQ=255)
                                                                                    gATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCCG.GGGGAGCCCGCTTGCGGGCGGGGGCAGCGCCCCCGAAGCGCCCg                     >  1:4395774/1‑101 (MQ=255)
                                                                                      tatTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCCG.GGGGAGCCCGCTTGCGGGCGGGGGCAGCGCCCCCGAAGCGCCCGAc                   >  2:1474047/1‑101 (MQ=255)
                                                                                      tatTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATACCCCG.GGGGAGCCCGCTTGCGGGCGGGGGCAGCGCCCCCGAAGCGCCCGAc                   <  2:3980756/101‑1 (MQ=255)
                                                                                       atTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCCG.GGGGAGCCCGCTTGCGGGCGGGGGCAGCGCCCCCGAAGCGCCCGAcg                  >  2:2809734/1‑101 (MQ=255)
                                                                                        tttCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCCG.GGGGAGCCCGCTTGCGGGCGGGGGCAGCGCCCCCGAAGCGCCCGAcgc                 <  2:1164977/101‑1 (MQ=255)
                                                                                         ttCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCCG.GGGGAGCCCGCTTGCGGGCGGGGGCAGCGCCCCCGAAGCGCCCGAcgcg                >  1:2178200/1‑101 (MQ=255)
                                                                                          tCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCCG.GGGGAGCCCGCTTGCGGGCGGGGGCAGCGCCCCCGAAGCGCCCGAcgcgc               >  1:3069158/1‑101 (MQ=255)
                                                                                           cGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCCG.GGGGAGCCCGCTTGCGGGCGGGGGCAGCGCCCCCGAAGCGCCCGAcgcgcg              <  1:3331345/101‑1 (MQ=255)
                                                                                           cGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCCG.GGGGAGCCCGCTTGCGGGCGGGGGCAGCGCCCCCGAAGCGCCCGAcgcgcg              >  1:1274712/1‑101 (MQ=255)
                                                                                            gggCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCCG.GGGGAGCCCGCTTGCGGGCGGGGGCAGCGCCCCCGAAGCGCCCGACGCgcgg             >  2:1340149/1‑101 (MQ=255)
                                                                                             ggCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCCG.GGGGAGCCCGCTTGCGGGCGGGGGCAGCGCCCCCGAAGCGCCCGACGCgcggc            >  1:1614818/1‑101 (MQ=255)
                                                                                             ggCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCCG.GGGGAGCCCGCTTGCGGGCGGGGGCAGCGCCCCCGAAGAGCCCGACGCgcggc            >  1:3119114/1‑101 (MQ=255)
                                                                                              gCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCCG.GGGGAGCCCGCTTGCGGGCGGGGGCAGCGCCCCCGAAGCGCCCGACGCgcggcg           <  2:3201781/101‑1 (MQ=255)
                                                                                              gCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCCG.GGGGAGCCCGCTTGCGGGCGGGGGCAGCGCCCCCGAAGCGCCCGACGCgcggcg           >  1:1149322/1‑101 (MQ=255)
                                                                                               cAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCCG.GGGGAGCCCGCTTGCGGGCGGGGGCAGCGCCCCCGAAGCGCCCGACGCGcggcgg          >  2:853590/1‑101 (MQ=255)
                                                                                               cAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCCG.GGGGAGCCCGCTTGCGGGCGGGGGCAGCGCCCCCGAAGCGCCCGACGCGcggcgg          >  2:2515150/1‑101 (MQ=255)
                                                                                               cAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCCG.GGGGAGCCCGCTTGCGGGCGGGGGCAGCGCCCCCGAAGCGCCCGACGCGcggcgg          >  1:2476696/1‑101 (MQ=255)
                                                                                                aaGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCCG.GGGGAGCCCGCTTGCGGGCGGGGGCAGCGCCCCCGAAGCGCCCGACGCGCGGCggg         <  2:3966646/101‑1 (MQ=255)
                                                                                                 aGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCCG.GGGGAGCCCGCTTGCGGGCGGGGGCAGCGCCCCCGAAGCGCCCGACGCGCGGCGGGc        >  1:2897534/1‑101 (MQ=255)
                                                                                                  gATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCCG.GGGGAGCCCGCTTGCGGGCGGGGGCAGCGCCCCCGAAGCTCCCGGCGCGCGGCGGGCt       <  2:401558/101‑1 (MQ=255)
                                                                                                  gATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCCG.GGGGAGCCCGCTTGCGGGCGGGGGCAGCGCCCCCGAAGCGCCCGACGCGCGGCGGGCt       >  2:107087/1‑101 (MQ=255)
                                                                                                   aTAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCCG.GGGGAGCCCGCTTGCGGGCGGGGGCAGCGCCCCCGAAGCGCCCGACGCGCGGCGGGCtt      <  1:996566/101‑1 (MQ=255)
                                                                                                    tAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCCG.GGGGAGCCCGCTTGCGGGCGGGGGCAGCGCCCCCGAAGCGCCCGACGCGCGGCGGGCttc     >  1:3744093/1‑100 (MQ=255)
                                                                                                    tAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCCG.GGGGAGCCCGCTTGCGGGCGGGGGCAGCGCCCCCGAAGCGCCCGACGCGCGGCGGGCtcg     >  1:1054335/1‑99 (MQ=255)
                                                                                                     aaaGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCCG.GGGGAGCCCGCTTGCGGGCGGGGGCAGCGCCCCCGAAGCGCCCGACGCGCGGCGGGGTTgc    >  1:160862/1‑101 (MQ=255)
                                                                                                     aaaGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCCG.GGGGAGCCCGCTTGCGGGCGGGGGCAGCGCCCCCGAAGCGCCCGACGCGCGGCGGGCTTgc    >  2:55519/1‑101 (MQ=255)
                                                                                                      aagaagGTCTTTTCCTTTCATCTTGCCCGAAATATCCCG.GGGGAGCCCGCTTGCGGGCGGGGGCAGCGCCCCCGAAGCGCCCGACGCGCGGCGGGGTTgcc   >  1:2823153/1‑101 (MQ=255)
                                                                                                      aagaagGTCTTTTCCTTTCATCTTGCCCGAAATATCCCG.GGGGAGCCCGCTTGCGGGCGGGGGCAGCGCCCCCGAAGCGCCCGACGCGCGGCGGGCTTgcc   >  2:1850171/1‑101 (MQ=255)
                                                                                                      aagaagGTCTTTTCCTTTCATCTTGCCCGAAATATCCCG.GGGGAGCCCGCTTGCGGGCGGGGGCAGCGCCCCCGAAGCGCCCGACGCGCGGCGGGCTTgcc   >  2:1433106/1‑101 (MQ=255)
                                                                                                       agaagGTCTTTTCCTTTCATCTTGCCCGAAATATCCCG.GGGGAGCCCGCTTGCGGGCGGGGGCAGCGCCCCCGAAGCGCCCGACGCGCGGCGGGCCTgccc  >  2:1264771/1‑101 (MQ=255)
                                                                                                       |                                                                                                     
CGCGACATCGGCTGGCGCGCGCTCTGAGACG.AAAC.GGGGGCGCTGCCCCC.CTCCTG..CTTCGGCCAATTCA.CCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCCG.GGGGAGCCCGCTTGCGGGCGGGGGCAGCGCCCCCGAAGCGCCCGACGCGCGGCGGGCTTGCCC  >  NC_014034/2566502‑2566699

Base quality scores:  ATCG  < 3 ≤  ATCG  < 18 ≤  ATCG  < 27 ≤  ATCG  < 33 ≤  ATCG  < 41 ≤  ATCG