Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_014034 2,566,572 Δ1 bp 5.8% intergenic (+44/+150) 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,5720C→. 5.8% 15.6 ‑1.6 ‑0.6121intergenic (+44/+150)purD/sdaAphosphoribosylamine‑‑glycine ligase/L‑serine ammonia‑lyase
Reads supporting (aligned to +/- strand):  new base (4/3):  ref base (32/81):  total (37/84)
Fisher's exact test strand distribution p-value = 1.95e-01
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 2.74e-01

Only 100 of 145 total aligned reads displayed.

TCGACTGGCCCGAGGGGTTCTGCCGCCGCGACATCGGCTGGCGCGCGCTCTGAGACGAAAC.GGGGGCGCTGCCCCC.CTCCT..G..CTTC.GGCCAATTCA.CCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCCGGGGGAGCCCGCTTGCGGGCGGGGGCAGCGCCCCCGA  >  NC_014034/2566476‑2566672
                                                                                                        |                                                                                                    
 cGACTGGCCCGAGGGGTTCTGCCGCCGCGACATCGGCTGGCGCGCGCTCTGAGACGAAACGGGGGGCGCTGCCCCC.CTCTT..GGCCTTC.GGCCAATTCA.ccc                                                                                                    <  2:659481/101‑1 (MQ=37)
   aCTGGCCCGAGGGGTTCTGCCGCCGCGACATCGGCTGGCGCGCGCTCTGAGACGAAACGGGGGGCGCTGCCCCC.CTCTT..GGCCTTC.GGCCAATTCA.ccccc                                                                                                  <  1:1187826/101‑1 (MQ=25)
    cTGGCCCGAGGGGTTCTGCCGCCGCGACATCGGCTGGCGCGCGCTCTGAGACGAAACGGGGGGCGCTGCCCCC.CTCTT..GGCCTTC.GGCCAATTCA.cccccc                                                                                                 <  2:702261/101‑1 (MQ=25)
      ggCCCGAGGGGTTCTGCCGCCGCGACATCGGCTGGCGCGCGCTCTGAGACGAAACGGGGGGCGCTGCCCCC.CTCTT..GGCCTTC.GGCCAATTCA.CCCCCCgg                                                                                               <  1:1099121/101‑1 (MQ=25)
      ggCCCGAGGGGTTCTGCCGCCGCGACATCGGCTGGCGCGCGCTCTGAGACGAAACGGGGGGCGCTGCCCCC.CTCTT..GGCCTTC.GGCCAATTCA.CCCCCCgg                                                                                               <  1:2944444/101‑1 (MQ=25)
      ggCCCGAGGGGTTCTGCCGCCGCGACATCGGCTGGCGCGCGCTCTGAGACGAAACGGGGGGCGCTGCCCCC.CTCTT..GGCCTTC.GGCCAATTCA.CCCCCCgg                                                                                               <  1:503621/101‑1 (MQ=25)
        cccGAGGGGTTCTGCCGCCGCGACATCGGCTGGCGCGCGCTCTGAGACGAAACGGGGGGCGCTGCCCCC.CTCTT..GGCCTTC.GGCCAATTCA.CCCCCCGGGa                                                                                             <  1:1829286/101‑1 (MQ=25)
         ccGAGGGGTTCTGCCGCCGCGACATCGGCTGGCGCGCGCTCTGAGACGAAACGGGGGGCGCTGCCCCC.CTCTT..GGCCTTC.GGCCAATTCA.CCCCCCGGGat                                                                                            <  2:4180170/101‑1 (MQ=25)
         ccGAGGGGTTCTGCCGCCGCGACATCGGCTGGCGCGCGCTCTGAGACGAAACGGGGGGCGCTGCCCCC.CTCTT..GGCCTTC.GGCCAATTCA.CCCCCCGGGat                                                                                            <  2:2168907/101‑1 (MQ=25)
          cGAGGGGTTCTGCCGCCGCGACATCGGCTGGCGCGCGCTCTGAGACGAAACGGGGGGCGCTGCCCCC.CTCTT..GGCCTTC.GGCCAATTCA.CCCCCCGGGata                                                                                           <  2:1018617/101‑1 (MQ=25)
           gAGGGGTTCTGCCGCCGCGACATCGGCTGGCGCGCGCTCTGAGACGAAACGGGGGGCGCTGCCCCC.CTCTT..GGCCTTC.GGCCAATTCA.CCCCCCGGGatat                                                                                          <  2:4131168/101‑1 (MQ=25)
           gAGGGGTTCTGCCGCCGCGACATCGGCTGGCGCGCGCTCTGAGACGAAACGGGGGGCGCTGCCCCC.CTCTT..GGCCTTC.GGCCAATTCA.CCCCCCCGGatat                                                                                          >  1:2119448/1‑101 (MQ=25)
             ggggTTCTGCCGCCGCGACATCGGCTGGCGCGCGCTCTGAGACGAAACGGGGGGCGCTGCCCCC.CTCTT..GGCCTTC.GGCCAATTCA.CCCCCCGGGATAttt                                                                                        <  2:77585/101‑1 (MQ=21)
             ggggTTCTGCCGCCGCGACATCGGCTGGCGCGCGCTCTGAGACGAAACGGGGGGCGCTGCCCCC.CTCTT..GGCCTTC.GGCCAATTCA.CCCCCCGGGATAttt                                                                                        <  1:1697153/101‑1 (MQ=21)
              gggTTCTGCCGCCGCGACATCGGCTGGCGCGCGCTCTGAGACGAAACGGGGGGGGGTGCCCCC.CCCTTGGG..CTTC.GGCCAATTCC.CCCCCCGGGATATTTc                                                                                       >  1:1249184/1‑101 (MQ=25)
                 ttCTGCCGCCGCGACATCGGCTGGCGCGCGCTCTGAGACGAAACGGGGGGCGCTGCCCCC.CTCTT..GGCCTTC.GGCCAATTCA.CCCCCCGGGATATTTCggg                                                                                    <  2:2145273/101‑1 (MQ=21)
                  tCTGCCGCCGCGACATCGGCTGGCGCGCGCTCTGAGACGAAACGGGGGGCGCTGCCCCC.CTCTT..GGCCTTC.GGCCCATTAA.ACCCCCCGGATATTTTgggg                                                                                   >  2:1592017/1‑100 (MQ=21)
                  tCTGCCGCCGCGACATCGGCTGGCGCGCGCTCTGAGACGAAACGGGGGGCGCTGCCCCC.CTCTT..GGCCTTC.GGCCAATTCA.CCCCCCGGGATATTTCGGGc                                                                                   <  2:229061/101‑1 (MQ=21)
                   cTGCCGCCGCGACATCGGCTGGCGCGCGCTCTGAGACGAAACGGGGGGCGCTGCCCCC.CTCTT..GGCCTTC.GGCCAATTCA.CCCCCCGGGATATTTCGGGCa                                                                                  <  2:2900729/101‑1 (MQ=21)
                   cTGCCGCCGCGACATCGGCTGGCGCGCGCTCTGAGACGAAACGGGGGGCGCTGCCCCC.CTCTT..GGCCTTC.GGCCAATTCA.CCCCCCCGGCTATTTCGGGCa                                                                                  >  2:861613/1‑101 (MQ=21)
                    tGCCGCCGCGACATCGGCTGGCGCGCGCTCTGAGACGAAGCGGGGGGCGCTGCCCCC.CTCTT..GGCCTTC.GGCCAATTCA.CCCCCCGGGATATTTCGGGCaa                                                                                 <  1:473092/101‑1 (MQ=21)
                    tGCCGCCGCGACATCGGCTGGCGCGCGCTCTGAGACGAAACGGGGGGCGCTGCCCCC.CCCTT..GGCCTTC.GGCCAATTCA.CCCCCCCGGAGAtttttggcca                                                                                 >  2:3957956/1‑94 (MQ=18)
                     gccgccGCGACATCGGCTGGCGCGCGCTCTGAGACGAAACGGGGGGCGCTGCCCCC.CTCTT..GGCCTTC.GGCCAATTCA.CCCCCCGGGATATTTCGGGCAAg                                                                                <  1:45837/101‑1 (MQ=21)
                       cgccgcGACATCGGCTGGCGCGCGCTCTGAGACGAAACGGGGGGGGCTGGCCCC.CTCTT..GGCCTTT.CGCCAATTCA.CCCCCCCGGACATTTCGGGCAAGaa                                                                              >  1:3980246/1‑100 (MQ=255)
                       cgccgcGACATCGGCTGGCGCGCGCTCTGAGACGAAACGGGGGGCGCTGCCCCC.CTCTT..GGCCTTC.GGCCAATTCA.CCCCCCGGGATATTTCGGGCAAGAt                                                                              <  2:2060072/101‑1 (MQ=21)
                       cgccgcGACATCGGCTGGCGCGCGCTCTGAGACGAAACGGGGGGCGCTGCCCCC.CTCTT..GGCCTTC.GGCCAATTCA.CCCCCCGGGATATTTCGGGCAAGAt                                                                              <  2:1620144/101‑1 (MQ=21)
                       cgccgcGACATCGGCTGGCGCGCGCTCTGAGACGAAACGGGGGGCGCTGCCCCC.CTCTT..GGCCTTC.GGCCAATTCA.CCCCCCGGGATATTTCGGGCAAGAt                                                                              <  2:3164/101‑1 (MQ=21)
                        gccgcGACATCGGCTGGCGCGCGCTCTGAGACGAAACGGGGGGCGCTGCCCCC.CTCTT..GGCCTTC.GGCCAATTCA.CCCCCCGGGATATTTCGGGCAAGATa                                                                             <  2:2066625/101‑1 (MQ=21)
                        gccgcGACATCGGCTGGCGCGCGCTCTGAGACGAAACGGGGGGCGCTGCCCCC.CTCTT..GGCCTTC.GGCCAATTCA.CCCCCCGGGATATTTCGGGCAAGATa                                                                             <  2:3108432/101‑1 (MQ=21)
                          cgcgACATCGGCTGGCGCGCGCTCTGAGACGAAACGGGGGGCGCTGCCCCC.CTCTT..GGCCTTC.GGCCCATTCA.CCCCCCCGGATATTTCGGGCAAGATaaa                                                                           >  1:3906941/1‑101 (MQ=21)
                            cgACATCGGCTGGCGCGCGCTCTGAGACGAAACGGGGGGCGCTGCCCCC.CTCTT..GGCCTTC.GGCCAATTCA.CCCCCCGGGATATTTCGGGCAAGATAaaga                                                                         <  1:1811536/101‑1 (MQ=21)
                            cgACATCGGCTGGCGCGCGCTCTGAGACGAAACGGGGGGCGCTGCCCCC.CTCTT..GGCCTTC.GGCCAATTCA.CCCCCCGGGATATTTCGGGCAAGATAaaga                                                                         >  2:3176725/1‑101 (MQ=21)
                             gACATCGGCTGGCGCGCGCTCTGAGACGAAACGGGGGGCGCTGCCCCC.CTCTT..GGCCTTC.GGCCAATTCA.CCCCCCGGGATATTTCGGGCAAGATAaagaa                                                                        <  1:19493/101‑1 (MQ=21)
                              aCATCGGCTGGCGCTCGCTCTGAGACGAAACGGGGGGCCCCGCCCCC.CTTTTG.G..CCTCGGGCCAATTCA.CCCCCCGCGCCGTTTCGGGGGCGATAaaaaag                                                                       >  1:4115095/1‑101 (MQ=9)
                                   ggCTGGGGCGCGCTCTGAGACGAAACGGGGGGCGCTGCCCCC.CTCTT..GGCCTTC.GGCCAATTCA.CCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCtt                                                                  <  1:3073594/101‑1 (MQ=21)
                                   ggCTGGCGCGCCCTCTGAGACGAAACGGGGGGCGCTGCCCCC.CTCTT..GGCCTTC.GGCCAATTCA.CCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCtt                                                                  <  1:1943026/101‑1 (MQ=21)
                                    gCTGGCGCGCGCTCTGAGACGAAACGGGGGGCGCTGCCCCC.CTCTT..GGCCTTC.GGCCAATTCA.CCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCttt                                                                 >  1:2196690/1‑101 (MQ=21)
                                    gCTGGCGCGCGCTCTGAGACGAAACGGGGGGCGCTGCCCCC.CTCTT..GGCCTTC.GGCCAATTCA.CCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCttt                                                                 <  2:1078597/101‑1 (MQ=21)
                                     cTGGCGCGCGCTCTGAGACGAAACGGGGGGCGCTGCCCCC.CTCTT..GGCCTTC.GGCCAATTCA.CCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCtttt                                                                <  1:2169706/101‑1 (MQ=21)
                                     cTGGCGCGCGCTCTGAGACGAAACGGGGGGCGCTGCCCCC.CTCTT..GGCCTTC.GGCCAATTCA.CCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCtttt                                                                <  1:2327547/101‑1 (MQ=21)
                                      tGGCGCGCGCTCTGAGACGAAACGGGGGGCGCTGCCCCC.CTCTT..GGCCTTC.GGCCAATTCA.CCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTc                                                               >  2:1009588/1‑101 (MQ=21)
                                      tGGCGCGCGCTCTGAGACGAAACGGGGGGCGCTGCCCCC.CTCTT..GGCCTTC.GGCCAATTCA.CCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTc                                                               <  1:3248456/101‑1 (MQ=21)
                                        gcgcgcgcTCTGAGACGAAACGGGGGGCGCTGCCCCC.CTCTT..GGCCTTC.GGCCAATTCA.CCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCt                                                             <  1:3727302/101‑1 (MQ=21)
                                         cgcgcgcTCTGAGACGAAACGGGGGGCGCTGCCCCC.CTCTT..GGCCTTC.GGCCAATTCA.CCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCtt                                                            <  1:1533020/101‑1 (MQ=21)
                                         cgcgcgcTCTGAGACGAAACGGGGGGCGCTGCCCCC.CTCTT..GGCCTTC.GGCCAATTCA.CCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCtt                                                            <  1:2512408/101‑1 (MQ=21)
                                         cgcgcgcTCTGAGACGAAACGGGGGGCGCTGCCCCC.CTCTT..GGCCTTC.GGCCAATTCA.CCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCtt                                                            <  2:3626458/101‑1 (MQ=21)
                                          gcgcgcTCTGAGACGAAACGGGGGGCGCTGCCCCC.CTCTT..GGCCTTC.GGCCAATTCA.CCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCttt                                                           <  2:1830164/101‑1 (MQ=21)
                                          gcgcgcTCTGAGACGAAACGGGGGGCGCTGCCCCC.CTCTT..GGCCTTC.GGCCAATTCA.CCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCttt                                                           <  2:2428260/101‑1 (MQ=21)
                                          gcgcgcTCTGAGACGAAACGGGGGGCGCTGCCCCC.CTCTT..GGCCTTC.GGCCAATTCA.CCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCttt                                                           >  2:3297536/1‑101 (MQ=21)
                                           cgcgcTCTGAGACGAAACGGGGGGCGCTGCCCCC.CTCTT..GGCCTTC.GGCCAATTCA.CCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTc                                                          <  1:509035/101‑1 (MQ=21)
                                           cgcgcTCTGAGACGAAACGGGGGGCGCTGCCCCC.CTCTT..GGCCTTC.GGCCAATTCA.CCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTc                                                          <  1:2076249/101‑1 (MQ=21)
                                             cgcTCTGAGACGAAACGGGGGGCGCTGCCCCC.CTCTT..GGCCTTC.GGCCAATTCA.CCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCAt                                                        <  1:4145111/101‑1 (MQ=21)
                                              gcTCTGAGACGAAACGGGGGGCGCTGCCCCC.CTCTT..GGCCTTC.GGCCAATTCA.CCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATc                                                       <  1:31132/101‑1 (MQ=21)
                                              gcTCTGAGACGAAACGGGGGGCGCTGCCCCC.CTCTT..GGCCTTC.GGCCAATTCA.CCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATc                                                       <  1:902741/101‑1 (MQ=21)
                                               ctctGAGACGAAACGGGGGGCGCTGCCCCC.CTCTT..GGCCTTC.GGCCAATTCA.CCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCt                                                      <  1:2775027/101‑1 (MQ=21)
                                                tctGAGACGAAACGGGGGGCGCTGCCCCC.CTCTT..GGCCTTC.GGCCAATTCA.CCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCtt                                                     <  1:298907/101‑1 (MQ=21)
                                                  tGAGACGAAACGGGGGGCGCTGCCCCC.CTCTT..GGCCTTC.GGCCAATTCA.CCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTGTTTTTTTTTCTTCTTGc                                                   >  2:853783/1‑101 (MQ=18)
                                                   gagaCGAAACAGGGGGCGCTACCCCC.CTCTT..GGCCTTC.GGCCAATTCC.CCCCCCGGAATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGcc                                                  <  1:3307486/101‑1 (MQ=21)
                                                    agaCGAAACGGGGGGCGCTGCCCCC.CTCTT..GGCCTTC.GGCCAATTCA.CCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGccc                                                 <  2:236725/101‑1 (MQ=21)
                                                      aCGAAACGGGGGGCGCTGCCCCC.CTCTT..GGCCTTC.GGCCAATTCA.CCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGa                                               >  1:1787489/1‑101 (MQ=21)
                                                       cGAAACGGGGGGCGCTGCCCCC.CTCTT..GGCCTTC.GGCCAATTCA.CCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGaa                                              >  2:3077751/1‑101 (MQ=21)
                                                       cGAAACGGGGGGCGCTGCCCCC.CTCTT..GGCCTTC.GGCCAATTCA.CCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGaa                                              <  1:3404108/101‑1 (MQ=21)
                                                        gAAACGGGGGGCGCTGCTCCC.CTCTT..GGCCTTC.GGCCAATTCA.CCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGaaa                                             <  2:3806296/101‑1 (MQ=21)
                                                        gAAACGGGGGGCGCTGCCCCC.CTCTT..GGCCTTC.GGCCAATTCA.CCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGaaa                                             <  2:2702264/101‑1 (MQ=21)
                                                         aaacggggggCGCTGCCCCC.CTCTT..GGCCTTC.GGCCAATTCA.CCCCCCGGGATATTTCGGGCAAGATAAAGAAGGGCTTTTCCTTTCATCTTGcccccaat                                            >  2:4223765/6‑96 (MQ=18)
                                                          aacggggggCGCTGCCCCC.CTCTT..GGCCTTC.GGCCAATTCACCCCCCCGGGAAATTT.GGGCCAGAAAAAGAAAGGCTTTTTCTTTCATCTTTcccccaaaa                                           >  1:2185736/5‑95 (MQ=2)
                                                          aacggggggCGCTGCCCCC.CTCTT..GGCCTTC.GGCCAATTCA.CCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAata                                           <  1:3581109/97‑1 (MQ=25)
                                                            cggggggCGCTGCCCCC.CTCTT..GGCCTTC.GGCCAATTCA.CCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATc                                         <  2:1120311/99‑1 (MQ=25)
                                                            cggggggCGCTGCCCCC.CTCTT..GGCCTTC.GGCCAATTCA.CCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATc                                         <  2:4015211/99‑1 (MQ=25)
                                                             ggggggCGCTGCCCCC.CTCTT..GGCCTTC.GGCCAATTCA.CCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATcc                                        <  2:1571825/100‑1 (MQ=25)
                                                               ggggCGCTGCCCCCTCTTGG..C..CTTT.GGCCAATTCA..CCCCTGGGATATTTCGGGCAAGAT...GAA.ATCATGCCCCTTCATCTTGCCGCAAATATCCCgggggga                                <  1:3959563/101‑3 (MQ=2)
                                                                gggCGCTGCCCCCTCTTGG..C..CTTT.GGCCAATTCA..CCCCTGGGATATTTCGGGCAAGAT...GAA.ATCATGCCCCTTCATCTTGCCGCAAATATCCCggggggag                               >  1:601326/1‑98 (MQ=2)
                                                                gggCGCTG.CCCC.CTCTT..GGCCTTT.GGCCAATTCA..CCCCTGGGATATTTCGGGCAAGAT...GAA.ATCATGCCCCTTCATCTTGCCGCAAATATCCCggggggag                               >  1:1254032/1‑98 (MQ=2)
                                                                 ggCGCTGCCCCCTCTTGG..C..CTTT.GGCCAATTCA..CCCCTGGGATATTTCGGGCAAGAT...GAA.ATCATGCCCCTTCATCTTGCCGCAAATATCCCggggggagc                              <  1:1048067/101‑5 (MQ=2)
                                                                 ggCGCTGCCCCCTCTTGG..C..CTTT.GGCCAATTCA..CCCCTGGGATATTTCGGGCAAGAT...GAA.ATCATGCCCCTTCATCTTGCCGCAAATATCCCggggggagc                              <  1:2334609/101‑5 (MQ=2)
                                                                   cgcTGCCCCC.CTCTT..GGCCTTC.GGCCAATTCA.CCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCCAAATATCCCggggg                                  >  2:3022218/1‑101 (MQ=25)
                                                                     ccgCCCCC.CTCTT..GGCCTTC.GGCCAAGTCC.CCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCCGGGGGAg                                <  2:3568447/99‑1 (MQ=37)
                                                                              ccctcttggccTTC.GGCCAATTCC.CCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCCGGGGGAGCCCGCt                          <  1:1052427/91‑1 (MQ=38)
                                                                              ccctcttggccTTC.GGCCAATTCA.CCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCCGGGGGAGCCCGCt                          <  1:93659/91‑1 (MQ=255)
                                                                                  cttggccTTC.GGCCAATTCA.CCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCCGGGGGAGCCCGCTTgcg                      >  2:3746717/7‑101 (MQ=255)
                                                                                   ttggccTTC.GGCCAATTCA.CCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCCGGGGGAGCCCGCTTgcgg                     <  2:2553261/96‑1 (MQ=255)
                                                                                   ttggccTTC.GGCCAATTCA.CCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCCGGGGGAGCCCGCTTgcgg                     <  1:3070202/96‑1 (MQ=255)
                                                                                   ttggccTTC.GGCCAATTCA.CCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCCGGGGGAGCCCGCTTgcgg                     <  1:3176964/96‑1 (MQ=255)
                                                                                    tggccTTC.GGCCAATTCA.CCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCCGGGGGAGACCGCTTgccgg                    >  2:1792567/5‑101 (MQ=255)
                                                                                     ggccTTC.GGCCAATTCA.CCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCCGGGGGAGCCCGCTTgcgggc                   >  1:1210820/4‑101 (MQ=255)
                                                                                     ggccTTC.GGCCAATTCA.CCCCCCGGGATATTTCGGGCAAGATAAAGAAGGGCTTTTCCTTTCATCTTGCCCGAAAGATCCCGGGGGACCCCGCCTgggggc                   >  1:3885267/4‑101 (MQ=25)
                                                                                        cTTC.GGCCAATTCA.CCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCCGGGGGAGCCCGCTTGcgggcggg                >  2:2407805/1‑101 (MQ=255)
                                                                                           c.GGCCAATTCA.CCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCCGGGGGAGCCCGCTTGCGGGCGGGGGc             <  2:2144179/101‑1 (MQ=255)
                                                                                             ggCCAATTCA.CCCCCCGGGATATTTTGGGCGAGATAAAGAAGGGGTTTTTCTTTTCTCTTTCCCGAAAAATCACGGGGGGGGCCGCCTCGTGGCggggggg            >  2:756706/1‑99 (MQ=16)
                                                                                              gCCAATTCA.CCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCCGGGGGAGCCCGCTTGCGGGCGGGGGCAg           <  2:1958095/101‑1 (MQ=255)
                                                                                              gCCAATTCA.CCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCCGGGGGAGCCCGCTTGCGGGCGGGGGCAg           <  1:3377084/101‑1 (MQ=255)
                                                                                                 aaTTCA.CCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCCGGGGGAGCCCGCTTGCGGGCGGGGGCAgcgc        >  2:1742483/1‑101 (MQ=255)
                                                                                                  aTTCA.CCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCCGGGGGAGCCCGCTTGCGGGCGGGGGCAGCGcc       >  1:3091680/1‑101 (MQ=255)
                                                                                                  aTTCA.CCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCCGGGGGAGCCCGCTTGCGGGCGGGGGCAGCGcc       >  1:2704649/1‑101 (MQ=255)
                                                                                                  aTTCA.CCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCCGGGGGAGCCCGCTTGCGGGCGGGGGCAGCGcc       >  2:1662111/1‑101 (MQ=255)
                                                                                                   ttCA.CCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCCGGGGGAGCCCGCTTGCGGGCGGGGGCAGCGccc      >  1:652542/1‑101 (MQ=255)
                                                                                                    tCA.CCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCCAGGGGAGCCCGCTTGCGGGGGGGGGGAGCGcccc     >  2:1633744/1‑101 (MQ=255)
                                                                                                     cA.CCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCCGGGGGAGCCCGCTTGCGGGCGGGGGCAGCGccccc    <  1:2029001/101‑1 (MQ=255)
                                                                                                      a.CCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCCGGGGGAGCCCGCTTGCGGGCGGGGGCAGCGCCCCCg   >  1:584338/1‑101 (MQ=255)
                                                                                                        ccccccGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCCGGGGGAGCCCGCTTGCGGGCGGGGGCAGCGCCCCCGa  <  1:3964337/101‑1 (MQ=255)
                                                                                                        |                                                                                                    
TCGACTGGCCCGAGGGGTTCTGCCGCCGCGACATCGGCTGGCGCGCGCTCTGAGACGAAAC.GGGGGCGCTGCCCCC.CTCCT..G..CTTC.GGCCAATTCA.CCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCCGGGGGAGCCCGCTTGCGGGCGGGGGCAGCGCCCCCGA  >  NC_014034/2566476‑2566672

Base quality scores:  ATCG  < 3 ≤  ATCG  < 18 ≤  ATCG  < 26 ≤  ATCG  < 33 ≤  ATCG  < 40 ≤  ATCG