Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_014034 2,566,536 +G 97.9% intergenic (+8/+186) 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,5361.→G97.9% 3.9 ‑0.6 ‑1.5193intergenic (+8/+186)purD/sdaAphosphoribosylamine‑‑glycine ligase/L‑serine ammonia‑lyase
Reads supporting (aligned to +/- strand):  new base (34/155):  ref base (1/3):  total (35/158)
Fisher's exact test strand distribution p-value = 5.54e-01
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 3.48e-02

Only 100 of 231 total aligned reads displayed.

GCCGAGGCGCAACAGCGCGCCTATGCGATGGTCGATGCGATCGACTGGCCCGAGGGGTTCTGCCGCCGCGACATCGGCTGGCGCGCGCTCTGAGACG.AAAC.GGGGG.CGCTGCCCCCCTCCT..G..CTTCGGCCAATTCACCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAA  >  NC_014034/2566436‑2566629
                                                                                                      |                                                                                                  
gCCGAGGCGCTACAGCGCGCCTATGCGAGGGTCTATGCGATCGATTGGCCCGAGGGGTTCTTCCGCCGCGACATCGGCTGGCGCGCGCTCTGAGACG.AAAc                                                                                                     <  1:2417699/101‑1 (MQ=255)
gCCGAGGCGCAACAGCGCGCCTATGCGATGGTCGATGCGATCGACTGGCCCGAGGGGTTCTGCCGCCGCGACATCGGCTGGCGCGCGCTCTGAGACG.AAAc                                                                                                     <  2:2663805/101‑1 (MQ=255)
 ccGAGGCGCAACAGCGCGCCTATGCGATGGTCGATGCGATCGACTGGCCCGAGGGGTTCTGCCGCCGCGACATCGGCTGGCGCGCGCTCTGAGACG.AAAC.g                                                                                                   <  1:2901635/101‑1 (MQ=255)
 ccGAGGCGCAACAGCGCGCCTATGCGATGGTCGATGCGATCGACTGGCCCGAGGGGTTCTGCCGCCGCGACATCGGCTGGCGCGCGCTCTGAGACG.AAAC.g                                                                                                   <  2:1288962/101‑1 (MQ=255)
   gAGGCGCAACAGCGCGCCTATGCGATGGTCGATGCGATCGACTGGCCCGAGGGGTTCTGCCGCCGCGACATCGGCTGGGGCGCGCTCTGAGACG.AAAC.ggg                                                                                                 >  2:2742679/1‑101 (MQ=255)
   gAGGCGCAACAGCGCGCCTATGCGATGGTCGATGCGATCGACTGGCCCGAGGGGTTCTGCCGCCGCGACATCGGCTGGCGCGCGCTCTGAGACG.AAAC.ggg                                                                                                 >  1:4085221/1‑101 (MQ=255)
    aGGCGCAACAGCGCGCCTATGCGATGGTCGATGGGATCGACTGGCCCGAGGGGGTCTGCCGCCGCGACATCGGCTGGCGCGCGCTCCGAGAAG.AAAA.gggg                                                                                                >  2:3930428/1‑101 (MQ=255)
       cgcAACAGCGCGCCTATGCGATGGTCGATGCGATCGACTGGCCCGAGGGGTTCTGCCGCCGCGACATCGGCTGGCGCGCGCTCTGAGAAG.AAAC.ggggggg                                                                                             >  1:35328/1‑99 (MQ=255)
            cAGCGCGCCTATGCGATGGTCGATGCGATCGACTGGCCCGAGGGGTTCTGCCGCCGCGACATCGGCTGGCGCGCGCTCTGAGACG.AAAC.cgggcgcccgca                                                                                        >  2:2952516/1‑93 (MQ=255)
              gcgcgcCTATGCGATGGTCGATGCGATCGACTGGCCCGAGGGGTTCTGCCGCCGCGACATCGGCTGGCGCGCGCTCTGAGACG.AAACGGGGGG.CGCTGccc                                                                                      <  1:1643011/101‑1 (MQ=255)
                gcgcCTATGCGATGGTCGATGCGATCGACTGGCCCGAGGGGTTCTGCCGCCGCGACATCGGCTGGCGCGCGCTCTGAGACG.AAACGGGGGG.CGCTGccccc                                                                                    <  1:330057/101‑1 (MQ=255)
                  gcCTATGCGATGGTCGATGCGATCGACTGGCCCGAGGGGTTCTGCCGCCGCGACATCGGCTGGCGCGCGCTCTGAGACG.AAACGGGGGG.CGCTGCCCCcct                                                                                  <  1:361247/101‑1 (MQ=255)
                     tATGCGATGGTCGATGCGATCGACTGGCCCGAGGGGTTCTGCCGCCGCGACATCGGCTGGCGCGCGCTCTGAGACGAAAAG.GGGGG.CGGTGCCCCccccct                                                                               >  2:967782/1‑101 (MQ=255)
                     tATGCGATGGTCGATGCGATCGACTGGCCCGAGGGGTTCTGCCGCCGCGACATCGGCTGGCGCGCGCTCTGAGACG.AAACGGGGGG.GGCTGCCCCcctctt                                                                               >  2:1577610/1‑99 (MQ=255)
                      aTGTGATGGTCGATGCGATCGACTGGCGCGAGGGGTTCTGCCTCCGCGACATCGGCTGGCGCGCGCTCTGAGACG.AAACGGGGGG.CGCTGCCCCCCTCTT..g                                                                            <  2:4231738/101‑1 (MQ=255)
                      aTGCGATGGTCGATGCGATCGACTGGCCCGAGGGGTTCTGCCGCCGCGACATCGGCTGGCGCGCGCTCTGAGACG.AAACGGGGGG.CGCTGCCCCCCTCTT..g                                                                            <  2:1108243/101‑1 (MQ=255)
                      aTGCGATGGTCGATGCGATCGACTGGCCCGAGGGGTTCTGCCGCCGCGACATCGGCTGGCGCGCGCTCTGAGACG.AAACGGGGGG.CGCTGCCCCCCTCTT..g                                                                            <  2:3022714/101‑1 (MQ=255)
                      aTGCGATGGTCGATGCGATCGACTGGCCCGAGGGGTTCTGCCGACGCGACATCGGCTGGCGCGAGCGCTGAGACG.AAACGGGGGG.CGCCGCCCCCCCCCT..g                                                                            >  1:4319592/1‑101 (MQ=255)
                        gCGATGGTCGATGCGATCGACTGGCCCGAGGGGTTCTGCCGCCGCGACATCGGCTGGGGCGCGCCCTGAGACGAAAAG.GGGGG.CGGCGCCCCcctctttgc                                                                            >  2:591367/1‑96 (MQ=255)
                        gCGATGGTCGATGCGATCGACTGGCCCGAGGGGTTCTGCCGCCGCGACATCGGCTGGCGCGCGCTCTGAGACG.AAACGGGGGG.CGCTGCCCCCCTCTT..ggc                                                                          <  1:146573/101‑3 (MQ=255)
                         cGATGGTCGATGCGATCGACTGGCCCGAGGGGTTCTGCCGCCGCGACATCGGCTGGCGCGCGCTCTGAGACG.AAACGGGGGG.CGCTGCCCCCCTCTT..ggcc                                                                         <  1:2053214/101‑4 (MQ=39)
                          gATGGTCGATGCGATCGACTGGCCCGAGGGGTTCTGCCGCCGCGACATCGGCTGGGGCGCGCCCTGGGACGAAAAC.GGGGG.GGCTGCCCCCCCCCT..ggccc                                                                        >  2:2812003/1‑97 (MQ=255)
                           aTGGTCGATGCGATCGACTGGCCCGAGGGGTTCTGCCGCCGCGACATCGGCTGGCGCGCGCTCTGAGACG.AAACGGGGGG.CGGGGCCCCCCCCCT..gggcct                                                                       >  1:3057443/1‑96 (MQ=255)
                            tGGTCGATGCGATCGACTGGCCCGAGGGGTTCTGCCGCCGCGACATCGGCTGGCGCGCGCTCTGAGACG.AAACGGGGGG.CGCTGCCCCCCTCTT..ggccttc                                                                      <  1:4899935/101‑7 (MQ=38)
                             ggTCGATGCGATCGACTGGCCCGAGGGGTTCTGCCGCCGCGACATCGGCTGGCGCGCGCTCTGAGACG.AAACGGGGGG.CGCTGCCCCCCTCTT..ggccttcg                                                                     <  2:3651458/101‑8 (MQ=38)
                                cGATGCGATCGACTGGCCCGAGGGGTTCTGCCGCCGCGACATCGGCTGGCGCGCGCTCTGAGACG.AAACGGGGGG.CGCTGCCCCCCTCTT..GGCCTTCGGcc                                                                  <  1:936316/101‑1 (MQ=37)
                                cGATGCGATCGACTGGCCCGAGGGGTTCTGCCGCCGCGACATCGGCTGGCGCGCGCTCTGAGACG.AAACGGGGGG.CGCTGCCCCCCTCTT..GGCCTTCGGcc                                                                  <  1:3815222/101‑1 (MQ=37)
                                  aTGCGATCGACTGGCCCGAGGGGTTCTGCCGCCGCGACATCGGCTGGCGCGCGCTCTGAGACG.AAACGGGGGG.CGCTGCCCCCCTCTT..GGCCTTCGGCCaa                                                                <  1:2214065/101‑1 (MQ=37)
                                    gCGATCGACTGGCCCGAGGGGTTCTGCCGCCGCGACATCGGCTGGCGCGCGCTCTGAGACG.AAACGGGGGG.CGCTGCCCCCCTCTT..GGCCTTCGGCCAAtt                                                              <  2:3415495/101‑1 (MQ=37)
                                     cGATCGACTGGCCCGAGGGGTTCTGCCGCCGCGACATCGGCTGGCGCGCGCTCTGAGACG.AAACGGGGGG.CGCTGCCCCCCTCTT..GGCCTTCGGCCAATTc                                                             <  1:1152670/101‑1 (MQ=37)
                                      ggtCGACTGGCCCGAGGGGTTCTGCCGCCGCGACATCGGCTGGCGCGCGCTCTGAGACG.AAACGGGGGG.CGCTGCCCCCCTCTT..GGCCTTCGGCCAATTCa                                                            <  2:1842999/99‑1 (MQ=25)
                                      gATCGACTGGCCCGAGGGGTTCTGCCGCCGCGACATCGGCTGGCGCGCGCTCTGAGACG.AAACGGGGGG.CGCTGCCCCCCTCTT..GGCCTTCGGCCAATTCa                                                            <  2:3126198/101‑1 (MQ=37)
                                        tCGACTGGCCCGGGGGGTTCTGCCGCCGCGACATCGGCTGGCGCGCGCTCTGAGACG.AAACGGGGGG.CGCTGCCCCCCTCTT..GGCCTTCGGCCAATTCAcc                                                          <  1:2288985/101‑1 (MQ=25)
                                        tCGACTGGCCCGAGGGGTTCTGCCGCCGCGACATCGGCTGGCGCGCGCTCTGAGACG.AAACGGGGGG.CGCTGCCCCCCTCTT..GGCCTTCGGCCAATTCAcc                                                          <  2:3286119/101‑1 (MQ=37)
                                         cGACTGGCCCGAGGGGTTCTGCCGCCGCGACATCGGCTGGCGCGCGCTCTGAGACG.AAACGGGGGG.CGCTGCCCCCCTCTT..GGCCTTCGGCCAATTCAccc                                                         <  2:3178069/101‑1 (MQ=37)
                                          gACTGGCCCGAGGGGTTCTGCCGCCGCGACATCGGCTGGCGCGCGCTCTGAGACG.AAACGGGGGG.CGGTGCCCCCCTCCTTGG..CCTCGGCCCATTTAcccc                                                        >  2:195032/1‑101 (MQ=255)
                                           aCTGGCCCGAGGGGTTCTGCCGCCGCGACATCGGCTGGCGCGCGCTCTGAGACG.AAACGGGGGG.CGCTGCCCCCCTCTT..GGCCTTCGGCCAATTCAccccc                                                       <  2:4887723/101‑1 (MQ=25)
                                             tGGCCCGAGGGGTTCTGCCGCCGCGACATCGGCTGGCGCGCGCTCTGAGACG.AAACGGGGGG.CGCTGCCCCCCTCTT..GGCCTTCGGCCAATTCACCCCCCg                                                     <  2:2580723/101‑1 (MQ=25)
                                               gCCCGAGGGGTTCTGCCGCCGCGACATCGGCTGTCGCGCGCTCTGAGACG.AAACGGGGGG.CGCTGCCCCCCTCTT..GGCCTTCGGCCAATTCACCCCCCggg                                                   <  1:4309561/101‑1 (MQ=25)
                                                  cGAGGGGTTCTGCCGCCGCGACATCGGCTGGCGCGCGCTCTGAGACG.AAACGGGGGG.CGGTGCCCCCCTCTT..GGCCTTCGGCCCATTCACCCCCCGGGaca                                                >  2:554286/1‑99 (MQ=25)
                                                    aGGGGTTCTGCCGCCGCGACATCGGCTGGCGCGCGCTCTGAGACG.AAACGGGGGG.CGCTGCCCCCCTCTT..GGCCTTCGGCCAATTCACCCCCCGGGATAtt                                              <  1:5046351/101‑1 (MQ=21)
                                                     ggggTTCTGCCGCCGCGACATCGGCTGGCGCGCGCTCTGAGACG.AAACGGGGGG.CGCTGCCCCCCTCTT..GGCCTTCGGCCAATTCACCCCCCGGGATAttt                                             <  1:4823843/101‑1 (MQ=21)
                                                     ggggTTCTGCCGCCGCGACATCGGCTGGCGCGCGCTCTGAGACG.AAACGGGGGG.CGCTGCCCCCCTCTT..GGCCTTCGGCCAATTCACCCCCCGGGATAttt                                             <  1:3399827/101‑1 (MQ=21)
                                                       ggTTCTGCCGCCGCGACATCGGCTGGCGCGCGCTCTGAGACG.AAACGGGGGG.CGCTGCCCCCCCCTT..GGCCTTCGGCCCATTCACCCCCCCGGAGATTTct                                           >  2:201523/1‑100 (MQ=21)
                                                        gTTCTGCCGCCGCGACATCGGCTGGCGCGCGCTCTGAGACG.AAACGGGGGG.CGCTGCCCCCCTCTT..GGCCTTCGGCCCAATCACCCCCCGGGATATTTTgg                                          >  2:4745212/1‑101 (MQ=21)
                                                          tCTGCCGCCGCGACATCGGCTGGCGCGCGCTCTGAGACG.AAACGGGGGG.CGCTGCCCCCCTCTT..GGCCTTCGGCCAATTCACCCCCCGGGATATTTCGGGc                                        <  2:2196847/101‑1 (MQ=21)
                                                           cTGCCGCCGCGACTTCGGCTGGCGCGCGCTCTGAGACG.AAACGGGGGG.CGCTGCCCCCCTCTT..GGCCTTCGGCCAATTCACCCCCCGGGATATTTCGGGCa                                       <  2:213462/101‑1 (MQ=21)
                                                             gccgccGCGACATCGGCTGGCGCGCGCTCTGAGACG.AAACGGGGGG.CGCTGCCCCCCTCTT..GGCCTTCGGCCAATTCACCCCCCGGGATATTTCGGGCAAg                                     <  1:4004963/101‑1 (MQ=21)
                                                              ccgccgCGACATCGGCTGGCGCGCGCTCTGAGACG.AAACGGGGGG.CGCTGCCCCCCTCTT..GGCCTTCGGCCAATTCACCCCCCGGGATATTTCGGGCAAGa                                    <  1:3738925/101‑1 (MQ=21)
                                                               cgccgcGACATCGGCTGGCGCGCGCTCTGAGACG.AAACGGGGGG.CGCTGCCCCCCTCTT..GGCCTTCGGCCAATTCACCCCCCGGGATATTTCGGGCAAGAt                                   <  1:4733697/101‑1 (MQ=21)
                                                                gccgcGACATCGGCTGGCGCGCGCTCTGAGACG.AAACGGGGGG.CGCTGCCCCCCTCTT..GGCCTTCGGCCAATTCACCCCCCGGGATATTTCGGGCAAGATa                                  <  2:2373027/101‑1 (MQ=21)
                                                                gccgcGACATCGGCTGGCGCGCGCTCTGAGACG.AAACGGGGGG.CGCTGCCCCCCTCTT..GGCCTTCGGCCAATTCACCCCCCGGGATATTTCGGGCAAGATa                                  <  1:3158213/101‑1 (MQ=21)
                                                                 ccgcGACATCGGCTGGCGCGCGCTCTGAGACG.AAACGGGGGG.GGGCGCCCCCCTCCT..GGCCTTTCGCCCATTCACCCCCCGGGATACTTCGGGCGAGACaa                                 >  2:812208/1‑101 (MQ=16)
                                                                 ccgcGACATCGGCTGGCGCGCGCTCTGAGACG.AAACGGGGGG.CGCTGCCCCCCTCTT..GGCCTTCGGCCAATTCACCCCCCGGGATATTTCGGGCAAGATaa                                 <  2:1821522/101‑1 (MQ=21)
                                                                 ccgcGACATCGGCTGGCGCGCGCTCTGAGACG.AAACGGGGGG.CGCTGCCCCCCTCTT..GGCCTTCGGCCAATTCACCCCCCGGGATATTTCGGGCAAGATaa                                 <  2:4311421/101‑1 (MQ=21)
                                                                 ccgcGACATCGGCTGGCGCGCGCTCTGAGACG.AAACGGGGGG.CGCTGCCCCCCTCTT..GGCCTTCGGCCAATTCACCCCCCGGGATATTTCGGGCAAGATaa                                 <  1:747985/101‑1 (MQ=21)
                                                                 ccgcGACATCGGCTGGCGCGCGCTCTGAGACG.AAACGGGGGG.CGCTGCCCCCCTCTT..GGCCTTCGGCCAATTCACCCCCCGGGATATTTCGGGCAAGATaa                                 <  1:655313/101‑1 (MQ=21)
                                                                 ccgcGACATCGGCTGGCGCGCGATCTGAGACG.AAACGGGGGG.CGCTGCCCCCCTCTT..GGCCTTCGGCCAATTCACCCCCCGGGATATTTCGGGCAAGATaa                                 <  1:2766840/101‑1 (MQ=21)
                                                                   gcgACATCGGCTGGCGCGCGCTCTGAGACG.AAACGGGGGG.CGCTGCCCCCCTCTT..GGCCTTCGGCCAATTCACCCCCCGGGATATTTCGGGCAAGATAaag                               <  1:2556790/101‑1 (MQ=21)
                                                                   gcgACATCGGCTGGCGCGCGCTCTGAGACG.AAACGGGGGG.CGCTGCCCCCCTCTT..GGCCTTCGGCCAATTCACCCCCCGGGATATTTCGGGCAAGATAaag                               <  1:3206906/101‑1 (MQ=21)
                                                                   gcgACATCGGCTGGCGCGCGCTCTGAGACG.AAACGGGGGG.CGCTGCCCCCCTCTT..GGCCTTCGGCCAATTCACCCCCCGGGATATTTCGGGCAAGATAaag                               <  2:641627/101‑1 (MQ=21)
                                                                    cgACATCGGCTGGCGCGCGCTCTGAGACG.AAACGGGGGG.CGCTGCCCCCCTCTT..GGCCTTCGGCCAATTCACCCCCCGGGCTATATAGGGGAAGATaaaaa                              >  2:3690992/1‑99 (MQ=18)
                                                                    cgACATCGGCTGGCGCGCGCTCTGAGACG.AAACGGGGGG.CGCTGCCCCCCTCTT..GGCCTTCGGCCAATTCACCCCCCGGGATATTTCGGGCAAGATAaaga                              <  1:1381444/101‑1 (MQ=21)
                                                                     gACATCGGCTGGCGCTCGCTCTGAGACG.AAACGGGGGG.CGCTGCCCCCCTCTT..GGCCTTCGGGCCAATCACCCCCCCGGGTAATTCGGGCCAGATAaaaaa                             >  1:2354133/1‑101 (MQ=12)
                                                                      aCATCGGCTGGCGCGCGCTCTGAGACG.AAACGGGGGG.CGCTGCCCCCCTCTT..GGCCTTCGGCCAATTCACCCCCCGGGATATTACGGGCAAGATAaagaag                            <  2:2525570/101‑1 (MQ=21)
                                                                         tCGGCTGGCGCGCGCTCTGAGACG.AAACGGGGGG.CGCTGCCCCCCTCTT..GGCCTTCGGCCAATTCACCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTc                         <  2:2440729/101‑1 (MQ=21)
                                                                          cGGCTGGCGCGCGCTCTGAGACG.AAACGGGGGG.CGCTGCCCCCCTCTT..GGCCTTCGGCCAATTCACCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCt                        <  1:535576/101‑1 (MQ=21)
                                                                          cGGCTGGCGCGCGCTCTGAGACG.AAACGGGGGG.CGCTGCCCCCCTCTT..GGCCTTCGGCCAATTCACCCCCCGGGATATTTCGGGCAAGATAAAGAAGGGCt                        >  1:598450/1‑101 (MQ=21)
                                                                           ggCTGGCGCGCGCTCTGAGACG.AAACGGGGGG.CGCTGCCCCCCTCTT..GGCCTTCGGCCAATTCACCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCtt                       >  2:4666595/1‑101 (MQ=21)
                                                                           ggCTGGCGCGCGCTCTGAGACG.AAACGGGGGG.CGCTGCCCCCCTCTT..GGCCTTCGGCCAATTCACCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCtt                       <  2:337048/101‑1 (MQ=21)
                                                                           ggCTGGCGCGCGCTCTGAGACG.AAACGGGGGG.CGCTGCCCCCCTCTT..GGCCTTCGGCCAATTCACCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCtt                       <  1:182174/101‑1 (MQ=21)
                                                                             cTGGCGCGCGCTCTGAGACG.AAACGGGGGG.CGCTGCCCCCCTCTT..GGCCTTCGGCCAATTCACCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCtttt                     <  1:1041996/101‑1 (MQ=21)
                                                                              tGGCGCGCGCTCTGAGACG.AAACGGGGGG.CGCTGCCCCCCTCTT..GGCCTTCGGCCAATTCACCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTc                    >  2:4917854/1‑101 (MQ=21)
                                                                              tGGCGCGCGCTCTGAGACG.AAACGGGGGG.CGCTGCCCCCCTCTT..GGCCTTCGGCCAATTCACCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTc                    <  2:5020865/101‑1 (MQ=21)
                                                                              tGGCGCGCGCTCTGAGACG.AAACGGGGGG.CGCTGCCCCCCTCTT..GGCCTTCGGCCAATTCACCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTc                    <  2:1768660/101‑1 (MQ=21)
                                                                               ggCGCGCGCTCTGAGACG.AAACGGGGGG.CGCTGCCCCCCTCTT..GGCCTTCGGCCAATTCACCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTcc                   >  1:4956338/1‑101 (MQ=21)
                                                                                gcgcgcgcTCTGAGACG.AAACGGGGGG.CGCTGCCCCCCTCTT..GGCCTTCGGCCAATTCACCCCCCGGGATATTTCGGGCAAGATAAAGAAGGGCTTTTCCt                  >  1:4505293/1‑101 (MQ=21)
                                                                                  gcgcgcTCTGAGACG.AAACGGGGGG.CGCTGCCCCCCTCTT..GGCCTTCGGCCAATTCACCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCttt                <  2:4016807/101‑1 (MQ=21)
                                                                                  gcgcgcTCTGAGACG.AAACGGGGGG.CGCTGCCCCCCTCTT..GGCCTTCGGCCAATTCACCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCttt                <  1:1700501/101‑1 (MQ=21)
                                                                                  gcgcgcTCTGAGACG.AAACGGGGGG.CGCTGCCCCCCTCTT..GGCCTTCGGCCAATTCACCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCttt                <  1:3285597/101‑1 (MQ=21)
                                                                                   cgcgcTCTGAGACG.AAACGGGGGG.CGCTGCCCCCCTCTT..GGCCTTCGGCCAATTCACCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTc               <  2:192013/101‑1 (MQ=21)
                                                                                   cgcgcTCTGAGACG.AAACGGGGGG.CGCTGCCCCCCTCTT..GGCCTTCGGCCAATTCACCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTc               <  2:406436/101‑1 (MQ=21)
                                                                                   cgcgcTCTGAGACG.AAACGGGGGG.CGCTGCCCCCCTCTT..GGCCTTCGGCCAATTCACCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTc               <  2:2870611/101‑1 (MQ=21)
                                                                                    gcgcTCTGAGACG.AAACGGGGGG.CGCTGCCCCCCTCTT..GGCCTTCGGCCAATTCACCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCa              <  1:2952744/101‑1 (MQ=21)
                                                                                     cgcTCTGAGACG.AAAC.GGGGGTCGCTTCCCCCCTCTT..GGCCTTCGGCCAATTCACCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCAt             <  1:1616734/101‑1 (MQ=25)
                                                                                      gcTCTGAGACG.AACCGGGGGG.CGCTGCCCCCTTCTT..GGCCTTCGGCCAATTCACCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATc            <  1:2652893/101‑1 (MQ=21)
                                                                                       ctctGAGACG.AAACGGGGGG.CGCTGCCCCCCTCTT..GGCCTTCGGCCAATTCACCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTTATCt           >  1:2119601/1‑101 (MQ=21)
                                                                                       ctctGAGACG.AAACGGGGGG.CGCTGCCCCCCTCTT..GGCCTTCGGCCAATTCACCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCt           >  2:4678419/1‑101 (MQ=21)
                                                                                        tctGAGACG.AAACGGGGGG.CGCTGCCCCCCTCTT..GGCCTTCGGCCAATTCACCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCtt          >  2:612440/1‑101 (MQ=21)
                                                                                        tctGAGACG.AAACGGGGGG.CGCTGCCCCCCTCTT..GGCCTTCGGCCAATTCACCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCtt          >  1:3091542/1‑101 (MQ=21)
                                                                                        tctGAGACG.AAACGGGGGG.CGCTGCCCCCCTCTT..GGCCTTCGGCCAATTCACCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCtt          <  1:88760/101‑1 (MQ=21)
                                                                                         ctGAGGCG.AAACGGGGGG.CGCTGCCCCCCTCTT..GGCCTTCGGCCAATTCTCCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTg         <  1:2566418/101‑1 (MQ=21)
                                                                                         ctGAGACG.AAACGGGGGG.CGCTGCCCCCCTCTT..GGCCTTCGGCCAATTCACCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTg         <  1:1346547/101‑1 (MQ=21)
                                                                                          tGAGACG.AAACGGGGGG.CGCTGCCCCCCTCTT..GGCCTTCGGCCAATTCACCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGc        <  2:1073410/101‑1 (MQ=21)
                                                                                          tGAGACG.AAACGGGGGG.CGCTGCCCCCCTCTT..GGCCTTCGGCCAATTCACCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGc        <  1:4249700/101‑1 (MQ=21)
                                                                                          tGAGACG.AAACGGGGGG.CGCTGCCCCCCTCTT..GGCCTTCGGCCAATTCACCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGc        >  2:3223909/1‑101 (MQ=21)
                                                                                            agaCG.AAACGGGGGG.CGCTGCCCCCCTCTT..GGCCTTCGGCCAATTCACCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGccc      <  2:1699276/101‑1 (MQ=21)
                                                                                            agaCG.AAACGGGGGG.CGCTGCCCCCCTCTT..GGCCTTCGGCCAATTCACCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGccc      <  1:191171/101‑1 (MQ=21)
                                                                                            agaCG.AAACGGGGGG.CGCTGCCCCCCTCTT..GGCCTTCGGCCAATTCACCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGccc      <  1:1748694/101‑1 (MQ=21)
                                                                                             gaCG.AAACGGGGGG.CGCTGCCCCCCTCTT..GGCCTTCGGCCAATTCACCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCg     <  2:3407080/101‑1 (MQ=21)
                                                                                                      |                                                                                                  
GCCGAGGCGCAACAGCGCGCCTATGCGATGGTCGATGCGATCGACTGGCCCGAGGGGTTCTGCCGCCGCGACATCGGCTGGCGCGCGCTCTGAGACG.AAAC.GGGGG.CGCTGCCCCCCTCCT..G..CTTCGGCCAATTCACCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAA  >  NC_014034/2566436‑2566629

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