Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_014034 2,566,609 T→G 2.4% intergenic (+81/+113) 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,6090T→G 2.4% 2.3 0.0 ‑0.0167intergenic (+81/+113)purD/sdaAphosphoribosylamine‑‑glycine ligase/L‑serine ammonia‑lyase
Reads supporting (aligned to +/- strand):  new base (2/2):  ref base (86/76):  total (89/78)
Fisher's exact test strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 9.86e-01

Only 100 of 191 total aligned reads displayed.

GCTGGCGCGCGCTCTGAGACGAAAC.GGGGGCGCTGCCCCC.CTCCTG..CTTCGGCCAATTCACCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCCG.GGGGAGCCCGCTTGCGGGCGGGGGCAGCGCCCCCGAAGCGCCCGACGCGCGGCGGGCTTGCCCGCCCAAGGAG  >  NC_014034/2566512‑2566709
                                                                                                     |                                                                                                     
gCTGGCGCGCGCTCTGAGACGAAACGGGGGGCGCTGCCCCC.CTCTTGGCCTTCGGCCAATTCACCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCttt                                                                                                       >  1:4154432/1‑101 (MQ=21)
 cTGGCGCGCGCTCTGAGACGAAACGGGGGGCGCTGCCCCC.CTCTTGGCCTTCGGCCAATTCACCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCtttt                                                                                                      <  1:4081669/101‑1 (MQ=21)
 cTGGCGCGCGCTCTGAGACGAAACGGGGGGCGCTGCCCCC.CTCTTGGCCTTCGGCCAATTCACCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCtttt                                                                                                      <  2:2134379/101‑1 (MQ=21)
 cTGGCGCGCGCTCTGAGACGAAACGGGGGGCGCTGCCCCC.CTCTTGGCCTTCGGCCAATTCACCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCtttt                                                                                                      <  1:885005/101‑1 (MQ=21)
  tGGCGCGCGCTCTGAGACGAAACGGGGGGCGCTGCCCCC.CTCTTGGCCTTCGGCCAATTCACCGCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTc                                                                                                     <  1:2647695/101‑1 (MQ=21)
  tGGCGCGCGCTCTGAGACGAAACGGGGGGCGCTGCCCCC.CTCTTGGCCTTCGGCCAATTCACCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTc                                                                                                     >  2:3361660/1‑101 (MQ=21)
  tGGCGCGCGCTCTGAGACGAAACGGGGGGCGCTGCCCCC.CTCTTGGCCTTCGGCCAATTCACCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTc                                                                                                     <  1:2981142/101‑1 (MQ=21)
    gcgcgcgcTCTGAGACGAAACGGGGGGCGCTGCCCCC.CTCTTGGCCTTCGGCCAATTCACCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCt                                                                                                   <  2:53519/101‑1 (MQ=21)
     cgcgcgcTCTGAGACGAAACGGGGGGCGCTGCCCCC.CTCTTGGCCTTCGGCCAATTCCCCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCtt                                                                                                  <  1:2822833/101‑1 (MQ=21)
     cgcgcgcTCTGAGACGAAACGGGGGGCGCTGCCCCC.CTCTTGGCCTTCGGCCAATTCACCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCtt                                                                                                  <  2:1094175/101‑1 (MQ=21)
      gcgcgcTCTGAGACGAAACGGGGGGCGCTGCCCCC.CTCTTGGCCTTCGGCCAATTCACCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCttt                                                                                                 <  1:2359893/101‑1 (MQ=21)
       cgcgcTCTGAGACGAAACGGGGGGCGCTGCCCCC.CTCTTGGCCTTCGGCCAATTCACCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTc                                                                                                <  1:1807017/101‑1 (MQ=21)
       cgcgcTCTGAGACGAAACGGGGGGCGCTGCCCCC.CTCTTGGCCTTCGGCCAATTCACCCCCCGGGATATTTCGGGCAAGATAAAGAAGGGCTTTTCCTTTc                                                                                                >  2:413218/1‑101 (MQ=21)
        gcgcTCTGAGACGAAACGGGGGGCGCTGCCCCC.CTCTTGGCCTTCGGCCAATTCACCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCa                                                                                               <  2:124700/101‑1 (MQ=21)
             ctGAGACGAAACGGGGGGCGCTGCCCCC.CTCTTGGCCTTCGGCCAATTCACCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTg                                                                                          >  1:1274440/1‑101 (MQ=21)
              tGAGACGAAACGGGGGGCGCTGCCCCC.CTCTTGGCCTTCGGCCAATTCACCCCCCGGGAAATTTCCTGCAAGATAAAGAAGGTCTTTTTCTTTCATCTTGc                                                                                         >  1:495608/1‑101 (MQ=18)
               gagaCGAAACGGGGGGCGCTGCCCCC.CTCTTGGCCTTCGGCCAATTCACCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGcc                                                                                        <  2:271613/101‑1 (MQ=21)
                agaCGAAACGGGGGGCGCTGCCCCC.CTCTTGGCCTTCGGCCAATTCACCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGccc                                                                                       <  2:4443182/101‑1 (MQ=21)
                 gaCGAAACGGGGGGCGCTGCCCCC.CTCTTGGCCTTCGGCCAATTCACCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCg                                                                                      <  2:2243119/101‑1 (MQ=21)
                 gaCGAAACGGGGGGCGCTGCCCCC.CTCTTGGCCTTCGGCCAATTCACCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCg                                                                                      <  1:1785610/101‑1 (MQ=21)
                   cGAAACGGGGGGCGCTGCCCCC.CTCTTGGCCTTCGGCCAATTCACCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGaa                                                                                    >  1:3643487/1‑101 (MQ=21)
                   cGAAACGGGGGGCGCTGCCCCC.CTCTTGGCCTTCGGCCAATTCACCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGaa                                                                                    <  1:740473/101‑1 (MQ=21)
                    gAAACGGGGGGCGCTGCCCCC.CTCTTGGCCTTCGGCCCATTCACCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGaaa                                                                                   <  1:3107135/101‑1 (MQ=21)
                      aacggggggCGCTGCCCCC.CTCTTGGCCTTCGGCCAATTCACCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAata                                                                                 <  1:2176547/97‑1 (MQ=25)
                      aacggggggCGCTGCCCCC.CTCTTGGCCTTCGGCCAATTCACCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAata                                                                                 <  1:3403555/97‑1 (MQ=25)
                       acggggggCGCTGCCCCC.CTCTTGGCCTTCGGCCAATTCACCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAatat                                                                                <  1:1158219/98‑1 (MQ=25)
                       acggggggCGCTGCCCCC.CTCTTGGCCTTCGGCCAATTCACCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAatat                                                                                >  2:4193579/4‑101 (MQ=25)
                         agggggCGCTG.CCCC.CTCTTGGCCTTTGGCCAATTCA.CCCCTGGGATATTTCGGGCAAGAT...GAA.ATCATGCCCCTTCATCTTGCCGCAAATATCCCg.gggg                                                                       >  1:2760966/2‑101 (MQ=0)
                         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.CTCTTGGCCTTCGGCCAATTCACCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATccc                                                                             <  1:2802997/101‑1 (MQ=25)
                          gggggCGCTGCCCCC.CTCTTGGCCTTCGGCCAATGCACCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATccc                                                                             >  1:4105956/1‑101 (MQ=25)
                           ggggCGCTGCCCCCTCTTGGC..CTTTGGCCAATTCA.CCCCTGGGATATTTCGGGCAAGAT...GAA.ATCATGCCCCTTCATCTTGCCGCAAATATCCCg.ggggga                                                                     <  2:2846502/101‑3 (MQ=2)
                              gcgcTGCCCCC.CTCTTGGCCTTCGGCCAATTCACCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCCg.ggg                                                                        <  1:550421/101‑1 (MQ=37)
                               cgcTGCCCCC.CTCTTGGCCTTCGGCCAATTCACCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCCg.gggg                                                                       <  2:4078109/101‑1 (MQ=37)
                               cgcTGCCCCC.CTCTTGGCCTTCGGCCAATTCACCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCCAAATATCCCg.gggg                                                                       >  2:3331108/1‑101 (MQ=25)
                                gcTGCCCCC.CTCTTGGCCTTCGGCCAATTCACCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCCG.GGGGa                                                                      <  1:2932405/101‑1 (MQ=37)
                                 cTGCCCCC.CTCTTGGCCTTCGGCCAATTCACCCCCCGGGATACTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCCG.GGGGAg                                                                     >  2:3937016/1‑101 (MQ=37)
                                  tGCCCCC.CTCTTGGCCTTCGGCCAATTCACCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCCG.GGGGAGc                                                                    >  2:1008137/1‑101 (MQ=37)
                                  tGCCCCC.CTCTTGGCCTTCGGCCAATTCACCCCCCGGGAAATTTCGGGGCAGATAAAGAGGGTTTTTTCCTTTCTTCTTTCCCGGAATATCCCC.CGGGacc                                                                    >  1:192838/1‑99 (MQ=16)
                                   gCCCCC.CTCTTGGCCTTCGGCCAATTCCCCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCCG.GGGGAGcc                                                                   <  1:693658/101‑1 (MQ=38)
                                    ccccc.cTCTGG..CTTC.GCC.ATTCACCCCCCGGGATATTTCGGGCACGATGA..AAGG.CATCCGGTTTCATCTTGCTGCAAATATCCCGAGGGGAGGCCGATTgc                                                            >  2:2946351/1‑101 (MQ=18)
                                        ccctcttggccTTCGGCCAATTCCCCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCCG.GGGGAGCCCGCt                                                               <  1:3591338/91‑1 (MQ=38)
                                         cctcttggccTTCGTCCAATTCACCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCCG.GGGGAGCCCGCtt                                                              <  2:3267489/92‑1 (MQ=37)
                                               ggccTTCGGCCAATTCCCCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCCG.GGGGAGCCCGCTTgcgggc                                                        <  1:3289530/98‑1 (MQ=255)
                                               ggccTTCGGCCAATTCACCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCCG.GGGGAGCCCGCCTgccggc                                                        >  2:2330775/4‑101 (MQ=255)
                                                   ctCGGCCAATTCACCCTCCCGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCCG.GGGGAGCCCGCTTGCGGGCgggg                                                    <  2:563537/100‑1 (MQ=255)
                                                     cGGCCAATTCACCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCCG.GGGGAGCCCGCTTGCGGGCGGGGGc                                                  >  2:1158868/1‑101 (MQ=255)
                                                      ggCCAATTCGCCCCCCGGGATATTTCGGGCGAGATAAAGAAGGTCTGTTCCTTTCATCTTGCCCGAAATATCCCG.GGGGAGCCCGCTTGCGGGCGGGGGCa                                                 <  2:3689381/101‑1 (MQ=255)
                                                      ggCCAATTCACCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCCG.GGGGAGCCCGCTTGCGGGCGGGGGCa                                                 <  1:1299350/101‑1 (MQ=255)
                                                         cAATTCACCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCCG.GGGGAGCCCGCTTGCGGGCGGGGGCAgcg                                              <  1:3361897/101‑1 (MQ=255)
                                                          aaTTCACCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCCG.GGGGAGCCCGCTTGCGGGCGGGGGCAgcgc                                             >  1:1584280/1‑101 (MQ=255)
                                                           aTTCACCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCCG.GGGGAGCCCGCTTGCGGGCGGGGGCAGCGcc                                            >  2:1839637/1‑101 (MQ=255)
                                                             tCACCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCCG.GGGGAGCCCGCTTGCGGGCGGGGGCAGCGcccc                                          >  1:1776744/1‑101 (MQ=255)
                                                              cACCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTTCCCGAAATATCCCG.GGGGACCCCGCCTGCGGGCGGGGGCCGGGccccc                                         >  1:209796/1‑101 (MQ=37)
                                                              cACCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCCG.GGGGAGCCCGCTTGCGGGCGGGGGCAGCGccccc                                         >  2:1404180/1‑101 (MQ=255)
                                                               cccccccGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCCG.GGGGAGCCCGCTTGCGGGCGGGGGCAGCGCCCCCg                                        <  2:1374118/100‑1 (MQ=255)
                                                                 cccccGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCCG.GGGGAGCCCGCTTGCGGGCGGGGGCAGCGCCCCCGaa                                      <  2:2318409/101‑1 (MQ=255)
                                                                 cccccGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCCG.GGGGAGCCCGCTTGCGGGCGGGGGCAGCGCCCCCGaa                                      >  2:2669688/1‑101 (MQ=255)
                                                                   cccGGGATATTTCGGGCTAGATAAAGATGGTCTTTTCCTTTCATCTTGCCCAAAATATCCCG.GGGGTTCCCGCTTGCGGGCGGGGGCAGCGCCCCCGAAgc                                    <  2:3970547/101‑1 (MQ=38)
                                                                     cGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCCG.GGGGAGCCCGCTTGCGGGCGGGGGCAGCGCCCCCGAAgcgc                                  >  2:4149359/1‑101 (MQ=255)
                                                                     cGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCCG.GGGGAGCCCGCTTGCGGGCGGGGGCAGCGCCCCCGAAgcgc                                  >  2:3884936/1‑101 (MQ=255)
                                                                        gATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCCG.GGGGAGCCCGCTTGCGGGCGGGGGCAGCGCCCCCGAAGCGCCCg                               >  1:4395774/1‑101 (MQ=255)
                                                                        gATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCCG.GGGGAGCCCGCTTGCGGGCGGGGGCAGCGCCCCCGAAGCGCCCg                               >  2:3995857/1‑101 (MQ=255)
                                                                          tatTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCCG.GGGGAGCCCGCTTGCGGGCGGGGGCAGCGCCCCCGAAGCGCCCGAc                             >  2:2838482/1‑101 (MQ=255)
                                                                          tatTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCCG.GGGGAGCCCGCTTGCGGGCGGGGGCAGCGCCCCCGAAGCGCCCGAc                             >  2:1474047/1‑101 (MQ=255)
                                                                           atTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCCG.GGGGAGCCCGCTTGCGGGCGGGGGCAGCGCCCCCGAAGCGCCCGAcg                            >  2:3277120/1‑101 (MQ=255)
                                                                            tttCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCCG.GGGGAGCCCGCTTGCGGGCGGGGGCAGCGCCCCCGAAGCGCCCGAcgc                           <  2:1164977/101‑1 (MQ=255)
                                                                             ttCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCCG.GGGGAGCCCGCTTGCGGGCGGGGGCAGCGCCCCCGAAGCGCCCGAcgcg                          >  1:428516/1‑101 (MQ=255)
                                                                              tCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCCG.GGGGAGCCCGCTTGCGGGCGGGGGCAGCGCCCCCGAAGCGCCCGAcgcgc                         >  1:3069158/1‑101 (MQ=255)
                                                                               cGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCCG.GGGGAGCCCGCTTGCGGGCGGGGGCAGCGCCCCCGAAGCGCCCGAcgcgcg                        >  1:78253/1‑101 (MQ=255)
                                                                                 ggCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCCG.GGGGAGCCCGCTTGCGGGCGGGGGCAGCGCCCCCGAAGCGCCCGACGCgcggc                      <  2:3161282/101‑1 (MQ=255)
                                                                                 ggCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCCG.GGGGAGCCCGCTTGCGGGCGGGGGCAGCGCCCCCGAAGAGCCCGACGCgcggc                      >  1:3119114/1‑101 (MQ=255)
                                                                                 ggCAAAATAAACAAGGTCTTTTCCTTTCATCTTTCCCGAAATATCCCG.GGGGAGCCCGCTTGCGGGCGGGGGCAGCGCCCCCGAAGCGCCCGACGCgcggc                      <  1:3299003/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)
                                                                                    naGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCCG.GGGGAGCCCGCTTGCGGGCGGGGGCAGCGCCCCCGAAGCGCCCGACGCgnnnngnn                   >  1:612353/2‑99 (MQ=255)
                                                                                     aGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCCG.GGGGAGCCCGCTTGCGGGCGGGGGCAGCGCCCCCGAAGCGCCCGACGCGCGGCGGGc                  >  1:2897534/1‑101 (MQ=255)
                                                                                       aTAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCCG.GGGGAGCCCGCTTGCGGGCGGGGGCAGCGCCCCCGAAGCGCCCGACGCGCGGCGGGCtt                <  2:1852133/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              >  1:1732047/1‑101 (MQ=255)
                                                                                          aagaagGTCTTTTCCTTTCATCTTGCCCGAAATATCCCG.GGGGAGCCCGCTTGCGGGCGGGGGCAGCGCCCCCGAAGCGCCCGACGCGCGGCGGGCTTgcc             >  2:1433106/1‑101 (MQ=255)
                                                                                           agaagGTCTTTTCCTTTCATCTTGCCCGAAATATCCCG.GGGGAGCCCGCTTGCGGGCGGGGGCAGCGCCCCCGAAGCGCCCGACGCGCGGCGGGCTTgccc            <  2:1693902/101‑1 (MQ=255)
                                                                                           agaagGTCTTTTCCTTTCATCTTGCCCGAAATATCCCG.GGGGAGCCCGCTTGCGGGCGGGGGCAGCGCCCCCGAAGCGCCCGACGCGCGGCGGGCCTgccc            >  2:1264771/1‑101 (MQ=255)
                                                                                            gaagGTCTTTTCCTTTCATCTTGCCCGAAATATCCCG.GGGGAGCCCGCTTGCGGGCGGGGGCAGCGCCCCCGAAGCGCCCGACGCGCGGCGGGCTTgcccg           >  2:195873/1‑101 (MQ=255)
                                                                                               ggTCTTTTCCTTTCATCTTGCCCGAAATATCCCG.GGGGAGCCCGCTTGCGGGCGGGGGCAGCGCCCCCGAAGCGCCCGACGCGCGGCGGGCTTgcccgccc        <  1:1417379/101‑1 (MQ=255)
                                                                                                gTCTTTTCCTTTCATCTTGCCCGAAATATCCCG.GGGGAGCCCGCTTGCGGGCGGGGGCAGCGCCCCCGAAGCGCCCGACGCGCGGCGGGCTTGCCCGCCCa       >  1:823562/1‑101 (MQ=255)
                                                                                                 tCTTTTCCTTTCATCTTGCCCGAAATATCCCG.GGGGAGCCCGCTTGCGGGCGGGGGCAGCGCCCCCGAAGCGCCCGACGCGCGGCGGGCTTGCCCGCCCaa      >  1:2518637/1‑101 (MQ=255)
                                                                                                 tCTTTTCCTTTCATCTTGCCCGAAATATCCCG.GGGGAGCCCGCTTGCGGGCGGGGGCAGCGCCCCCGAAGCGCCCGACGCGCGGCGGGCTTGCCCGCCCaa      >  2:2313715/1‑101 (MQ=255)
                                                                                                  cTTTTCCTTTCATCTTGCCCGAAATATCCCG.GGGGAGCCCGCTTGCGGGCGGGGGCAGCGCCCCCGAAGCGCCCGACGCGCGGCGGGCTTGCCCGCCCAag     >  1:2858775/1‑101 (MQ=255)
                                                                                                  cTTTTCCTTTCATCTTGCCCGAAATATCCCG.GGGGAGCCCGCTTGCGGGCGGGGGCAGCGCCCCCGAAGCGCCCGACGCGCGGCGGGCTTGCCCGCCCAag     >  2:1510093/1‑101 (MQ=255)
                                                                                                   ttttCCTTTCATCTTGCCCGAAATATCCCG.GGGGAGCCCGCTTGCGGGCGGGGGCAGCGCCCCCGAAGCGCCCGACGCGCGGCGGGCTTGCCCGCCCCagg    >  2:3960001/1‑101 (MQ=255)
                                                                                                    tttCCTTTCATCTTGCCCGAAATATCCCG.GGGGAGCCCGCTTGCGGGCGGGGGCAGCGCCCCCGAAGCGCCCGACGCGCGGCGGGCTTGCCCGCCCAagga   >  1:1082312/1‑101 (MQ=255)
                                                                                                    tttCCTTTCATATGNCCCGAAANATCNGN.GGGGAGCCCGCTTGCGGGCGGGGGCAGCGCCCCCGAAGCGCCCGACGCGCGGCGGGCTTGCCCGCCCAagga   <  1:3885190/101‑1 (MQ=37)
                                                                                                     ttCCTTTCATCTTGCCCGAAATATCCCG.GGGGAGCCCGCTTGCGGGCGGGGGCAGCGCCCCCGAAGCGCCCGACGCGCGGCGGGCTTGCCCGCCCAaggag  >  1:4359464/1‑101 (MQ=255)
                                                                                                     ttCCTTTCATCTTGCCCGAAATATCCCG.GGGGAGCCCGCTTGCGGGCGGGGGCAGCGCCCCCGAAGCGCCCGACGCGCGGCGGGCTTGCCCGCCCAaggag  >  2:3572646/1‑101 (MQ=255)
                                                                                                     ttCCTTTCATCTTGCCCGAAATATCCCG.GGGGAGCCCGCTTGCGGGCGGGGGCAGCGCCCCCGAAGCGCCCGACGCGCGGCGGGCTTGCCCGCCCAaggag  >  2:2409043/1‑101 (MQ=255)
                                                                                                     |                                                                                                     
GCTGGCGCGCGCTCTGAGACGAAAC.GGGGGCGCTGCCCCC.CTCCTG..CTTCGGCCAATTCACCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCCG.GGGGAGCCCGCTTGCGGGCGGGGGCAGCGCCCCCGAAGCGCCCGACGCGCGGCGGGCTTGCCCGCCCAAGGAG  >  NC_014034/2566512‑2566709

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