Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_014034 2,566,597 Δ1 bp 2.5% intergenic (+69/+125) 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,5970A→. 2.5% 4.4 0.0 ‑0.2160intergenic (+69/+125)purD/sdaAphosphoribosylamine‑‑glycine ligase/L‑serine ammonia‑lyase
Reads supporting (aligned to +/- strand):  new base (2/2):  ref base (69/86):  total (72/88)
Fisher's exact test strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 5.64e-01

Only 100 of 187 total aligned reads displayed.

GCCGCGACATCGGCTGGCGCGCGCTCTGAGACG.AAAC.GGGGGCGCTGCCCCC.CTCCTG..CTTCGGCCAATTCA.CCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCCG.GGGGAGCCCGCTTGCGGGCGGGGGCAGCGCCCCCGAAGCGCCCGACGCGCGGCGGGCTTGC  >  NC_014034/2566500‑2566697
                                                                                                       |                                                                                                     
cccgcGACATCGGCTGGCGCGCCCTCTGAGACG.AAACGGGGGGCGCTGCCCCC.CTCTTGGCCTTCGGCCAATTCA.CCCCCCGGGATATTTCGGGCAAGATa                                                                                                       <  2:3598144/100‑1 (MQ=21)
 ccgcGACATCGGCTGGCGCGCGCTCTGAGACG.AAACGGGGGGCGCTGCCCCC.CTCTTGGCCTTCGGCCAATTCA.CCCCCCGGGATATTTCGGGCAAGATaa                                                                                                      <  2:2601768/101‑1 (MQ=21)
  cgcgACATCGGCTGGCGCGCGCTCTGAGACG.AAACGGGGGGCGCTGCCCCC.CTCTTGGCCTTCGGCCAATTCA.CCCCCCGGGATATTTCGGGCAAGATaaa                                                                                                     <  2:1980484/101‑1 (MQ=21)
  cgcgACATCGGCTGGCGCGCGCTCTGAGACG.AAACGGGGGGCGCTGCCCCC.CTCTTGGCCTTCGGCCAATTCA.CCCCCCGGGATATTTCGGGCAAGATaaa                                                                                                     <  2:3837105/101‑1 (MQ=21)
  cgcgACATCGGCTGGCGCGCGCTCTGAGACG.AAACGGGGGGCGCTGCCCCC.CTCTTGGCCTTCGGCCAATTCA.CCCCCCGGGATATTTCGGGCAAGATaaa                                                                                                     <  1:2582465/101‑1 (MQ=21)
  cgcgACATCGGCTGGCGCGCGCTCTGAGACG.AAACGGGGGGCGCTGCCCCC.CTCTTGGCCTTCGGCCAATTCA.CCCCCCGGGATATTTCGGGCAAAATaaa                                                                                                     <  1:295598/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)
     gACATCGGCTGGCGCGCGCTCTGAGACG.AAACGGGGGGCGCTGCCCCC.CTCTTGGCCTTCGGCCAATTCA.CCCCCCGGGATATTTCGGGCAAGATAaagaa                                                                                                  <  2:3920021/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)
           ggCTGGCGCGCGCTCTGAGACG.AAACGGGGGGCGCTGCCCCC.CTCTTGGCCTTCGGCCAATTCA.CCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCtt                                                                                            <  1:3569336/101‑1 (MQ=21)
           ggCTGGCGCGCGCTCTGAGACG.AAACGGGGGGCGCTGCCCCC.CTCTTGGCCTTCGGCCAATTCA.CCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCtt                                                                                            <  2:1394035/101‑1 (MQ=21)
           ggCTGGCGCGCGCTCTGAGACG.AAACGGGGGGCGCTGCCCCC.CTCTTGGCCTTCGGCCAATTCA.CCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCtt                                                                                            >  2:1299253/1‑101 (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                                                                                         >  1:3459804/1‑101 (MQ=21)
              tGGCGCGCGCTCTGAGACG.AAACGGGGGGCGCTGCCCCC.CTCTTGGCCTTCGGCCAATTCA.CCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTc                                                                                         <  1:2981142/101‑1 (MQ=21)
                gcgcgcgcTCTGAGACG.AAACGGGGGGCGCTGCCCCC.CTCTTGGCCTTCGGCCAATTCA.CCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCt                                                                                       <  1:4349531/101‑1 (MQ=21)
                  gcgcgcTCTGAGACG.AAACGGGGGGCGCTGCCCCC.CTCTTGGCCTTCGGCCAATTCA.CCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCttt                                                                                     <  1:3001361/101‑1 (MQ=21)
                  gcgcgcTCTGAGACG.AAACGGGGGGCGCTGCCCCC.CTCTTGGCCTTCGGCCAATTCA.CCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCttt                                                                                     <  1:2359893/101‑1 (MQ=21)
                   cgcgcTCTGAGACG.AAACGGGGGGCGCTGCCCCC.CTCTTGGCCTTCGGCCAATTCA.CCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTc                                                                                    <  1:1807017/101‑1 (MQ=21)
                    gcgcTCTGAGACG.AAACGGGGGGCGCTGCCCCC.CTCTTGGCCTTCGGCCAATTCA.CCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCa                                                                                   <  2:124700/101‑1 (MQ=21)
                     cgcTCTGAGACG.AAACGGGGGGCGCTGCCCCC.CTCTTGGCCTTCGGCCAATTCACCCCCCCGGGAAATTTCGGGCAAAAAAAAAAAGGTTTTTTCTttttt                                                                                   >  1:3164572/1‑99 (MQ=12)
                     cgcTCTGAGACG.AAACGGGGGGCGCTGCCCCC.CTCTTGGCCTTCGGCCAATTCA.CCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCAt                                                                                  <  1:3596700/101‑1 (MQ=21)
                      gcTCTGAGACG.AAACGGGGGGCGCTGCCCCC.CTCTTGGCCTTCGGCCAATTCC.CCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATc                                                                                 <  1:817442/101‑1 (MQ=21)
                         ctGAGACG.AAACGGGGGGCGCTGCCCCC.CTCTTGGCCTTCGGCCAATTCA.CCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTg                                                                              >  2:2689914/1‑101 (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)
                            agaCG.AAACGGGGGGCGCTGCCCCC.CTCTTGGCCTTCGGCCAATTCA.CCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGccc                                                                           <  2:4443182/101‑1 (MQ=21)
                             gaCG.AAACGGGGGGCGCTGCCCCC.CTCTTGGCCTTCGGCCAATTCA.CCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCg                                                                          <  2:2243119/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.CTCTTGGCCTTCGGCCAATTCA.CCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGaaa                                                                       <  2:4135017/101‑1 (MQ=21)
                                  aaacggggggCGCTGCCCCC.CTCTTGGCCTTCGGCCAATTCA.CCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAat                                                                      <  1:661025/96‑1 (MQ=21)
                                   aacggggggCGCTGCCCCC.CTCTTGGCCTTCGGCCAATTCA.CCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAata                                                                     >  1:2794101/5‑101 (MQ=25)
                                    acggggggCGCTGCCCCC.CTCTTGGCCTTCGGCCAATTCA.CCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAatat                                                                    <  1:1158219/98‑1 (MQ=25)
                                    acggggggCGCTGCCCCC.CTCTTGGCCTTCGGCCAATTCA.CCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAatat                                                                    >  2:4193579/4‑101 (MQ=25)
                                       gggggCGCTGCCCCCTCTTGGC..CTTTGGCCAATTCA..CCCCTGGGATATTTCGGGCAAGAT...GAA.ATCATGCCCCTTCATCTTGCCGCAAATATCCCg.ggggg                                                          <  1:3971574/101‑2 (MQ=0)
                                       gggggCGCTGCCCCC.CTCTTGGCCTTCGGCCAATTCA.CCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATccc                                                                 <  1:2729715/101‑1 (MQ=25)
                                       gggggCGCTGCCCCC.CTCTTGGCCTTCGGCCAATTCA.CCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATccc                                                                 <  1:3095488/101‑1 (MQ=25)
                                        ggggCGCTGCCCCCTCTTGGC..CTTTGGCCAATTCA..CCCCTGGGATATTTCGGGCAAGAT...GAA.ATCATGCCCCTTCATCTTGCCGCAAATATCCCg.ggggga                                                         <  2:2846502/101‑3 (MQ=2)
                                         gggCGCTGCCCCC.CTCTTGGCCTTCGGCCAATTCA.CCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCCg.g                                                              <  1:3108485/101‑1 (MQ=37)
                                           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:4128436/101‑1 (MQ=37)
                                              cTGCCCCC.CTCTTGGCCTTCGGCCAATTCA.CCCCCCGGGATACTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCCG.GGGGAg                                                         >  2:3937016/1‑101 (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)
                                                      cctcttggccTTCGGCCAATTCA.CCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGCAATATCCCG.GGGGAGCCCGCtt                                                  >  2:4177887/10‑101 (MQ=37)
                                                      cctcttggccTTCGGCCAATTCA.CCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCCG.GGGGAGCCCGCtt                                                  <  2:2849051/92‑1 (MQ=255)
                                                       ctcttggccTTCGGCCAATTCA.CCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCCG.GGGGAGCCCGCTTg                                                 <  1:281610/93‑1 (MQ=255)
                                                         cttggccTTCGGCCAATTCA.CCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCCG.GGGGAGCCCGCTTgcg                                               <  1:3958514/95‑1 (MQ=255)
                                                         cttggccTTCGGCCAATTCA.CCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCCG.GGGGAGCCCGCTTgcg                                               >  1:2922025/7‑101 (MQ=255)
                                                           tggccTTCGGCCAATTCA.CCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCCG.GGGGAGCCCGCTTgcggg                                             <  1:747769/97‑1 (MQ=255)
                                                            ggccTTCGGCCAATTCC.CCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCCG.GGGGAGCCCGCTTgcgggc                                            <  1:3289530/98‑1 (MQ=255)
                                                            ggccTTCGGCCAATTCA.CCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCCG.GGGGAGCCCGCCTgccggc                                            >  2:2330775/4‑101 (MQ=255)
                                                               cTTCGGGCAATTCA.CCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCCG.GGGGAGCCCGCTTGcgggcggg                                         <  1:725245/101‑1 (MQ=255)
                                                                      cAATTCA.CCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCCG.GGGGAGCCCGCTTGCGGGCGGGGGCAgcg                                  <  1:3361897/101‑1 (MQ=255)
                                                                       aaTTCA.CCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCCG.GGGGAGCCCGCTTGCGGGCGGGGGCAgcgc                                 >  1:1584280/1‑101 (MQ=255)
                                                                        aTTCA.CCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCCG.GGGGAGCCCGCTTGCGGGCGGGGGCAGCGcc                                >  2:1839637/1‑101 (MQ=255)
                                                                          tCA.CCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCCG.GGGGAGCCCGCTTGCGGGCGGGGGCAGCGcccc                              >  2:3359653/1‑101 (MQ=255)
                                                                           cA.CCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTTCCCGAAATATCCCG.GGGGACCCCGCCTGCGGGCGGGGGCCGGGccccc                             >  1:209796/1‑101 (MQ=37)
                                                                           cA.CCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCCG.GGGGAGCCCGCTTGCGGGCGGGGGCAGCGccccc                             <  2:2317730/101‑1 (MQ=255)
                                                                           cA.CCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCCG.GGGGAGCCCGCTTGCGGGCGGGGGCAGCGccccc                             >  2:1404180/1‑101 (MQ=255)
                                                                           cA.CCCCCCGGGATATTTCGGGCAAGATAAAGAAGGGCTTTTTCTTTCATCTTGCCCGAAATATCCCG.GGGGGGCCCGCTTGCGGGGGGGGGGCGCGccccc                             >  2:3072127/1‑101 (MQ=21)
                                                                            cgccccccGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCCG.GGGGAGCCCGCTTGCGGGCGGGGGCAGCGccccc                             <  1:1738269/99‑1 (MQ=255)
                                                                             ccccctcCGGATATTTCGGGCAGGATAAAGCAGGTCTTTTCCCTCCACGTTGTCGCAAACACCCCG.GGGGGGCCCGCCAGCGGGCGGGGGCAGCGCCCCCg                            <  2:1427563/100‑1 (MQ=255)
                                                                             cccccccGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCCG.GGGGAGCCCGCTTGCGGGCGGGGGCAGCGCCCCCg                            <  2:1374118/100‑1 (MQ=255)
                                                                               cccccGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCCG.GGGGAGCCCGCTTGCGGGCGGGGGCAGCGCCCCCGaa                          <  2:2318409/101‑1 (MQ=255)
                                                                                 cccGGGATATTTCGGGCTAGATAAAGATGGTCTTTTCCTTTCATCTTGCCCAAAATATCCCG.GGGGTTCCCGCTTGCGGGCGGGGGCAGCGCCCCCGAAgc                        <  2:3970547/101‑1 (MQ=38)
                                                                                    gggATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCCG.GGGGAGCCCGCTTGCGGGCGGGGGCAGCGCCCCCGAAGCGcc                     <  2:4105696/101‑1 (MQ=255)
                                                                                      gATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCCG.GGGGAGCCCGCTTGCGGGCGGGGGCAGCGCCCCCGAAGCGCCCg                   >  2:3995857/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)
                                                                                         atTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCCG.GGGGAGCCCGCTTGCGGGCGGGGGCAGCGCCCCCGAAGCGCCCGAcg                >  2:2809734/1‑101 (MQ=255)
                                                                                          tttCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCCG.GGGGAGCCCGCTTGCGGGCGGGGGCAGCGCCCCCGAAGCGCCCGAcgc               <  2:1164977/101‑1 (MQ=255)
                                                                                           ttCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCCG.GGGGAGCCCGCTTGCGGGCGGGGGCAGCGCCCCCGAAGCGCCCGAcgcg              >  1:2178200/1‑101 (MQ=255)
                                                                                             cGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCCG.GGGGAGCCCGCTTGCGGGCGGGGGCAGCGCCCCCGAAGCGCCCGAcgcgcg            >  1:1274712/1‑101 (MQ=255)
                                                                                               ggCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCCG.GGGGAGCCCGCTTGCGGGCGGGGGCAGCGCCCCCGAAGCGCCCGACGCgcggc          >  1:1614818/1‑101 (MQ=255)
                                                                                               ggCAAAATAAACAAGGTCTTTTCCTTTCATCTTTCCCGAAATATCCCG.GGGGAGCCCGCTTGCGGGCGGGGGCAGCGCCCCCGAAGCGCCCGACGCgcggc          <  1:3299003/101‑1 (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:3375116/1‑101 (MQ=255)
                                                                                                 cAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCCG.GGGGAGCCCGCTTGCGGGCGGGGGCAGCGCCCCCGAAGCGCCCGACGCGcggcgg        >  2:2515150/1‑101 (MQ=255)
                                                                                                 cAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCCG.GGGGAGCCCGCTTGCGGGCGGGGGCAGCGCCCCCGAAGCGCCCGACGCGcggcgg        >  1:2476696/1‑101 (MQ=255)
                                                                                                   aGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCCG.GGGGAGCCCGCTTGCGGGCGGGGGCAGCGCCCCCGAAGCGCCCGACGCGCGGCGGGc      >  1:2897534/1‑101 (MQ=255)
                                                                                                    gATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCCG.GGGGAGCCCGCTTGCGGGCGGGGGCAGCGCCCCCGAAGCGCCCGACGCGCGGCGGGCt     >  2:107087/1‑101 (MQ=255)
                                                                                                     aTAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCCG.GGGGAGCCCGCTTGCGGGCGGGGGCAGCGCCCCCGAAGCGCCCGACGCGCGGCGGGCtt    <  1:996566/101‑1 (MQ=255)
                                                                                                      tAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCCG.GGGGAGCCCGCTTGCGGGCGGGGGCAGCGCCCCCGAAGCGCCCGACGCGCGGCGGGCtt  >  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)
                                                                                                       |                                                                                                     
GCCGCGACATCGGCTGGCGCGCGCTCTGAGACG.AAAC.GGGGGCGCTGCCCCC.CTCCTG..CTTCGGCCAATTCA.CCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCCG.GGGGAGCCCGCTTGCGGGCGGGGGCAGCGCCCCCGAAGCGCCCGACGCGCGGCGGGCTTGC  >  NC_014034/2566500‑2566697

Base quality scores:  ATCG  < 3 ≤  ATCG  < 20 ≤  ATCG  < 28 ≤  ATCG  < 34 ≤  ATCG  < 41 ≤  ATCG