Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_014034 2,566,572 Δ1 bp 2.5% 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→. 2.5% 6.9 ‑0.5 ‑0.3201intergenic (+44/+150)purD/sdaAphosphoribosylamine‑‑glycine ligase/L‑serine ammonia‑lyase
Reads supporting (aligned to +/- strand):  new base (2/3):  ref base (55/140):  total (57/144)
Fisher's exact test strand distribution p-value = 6.25e-01
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 4.47e-01

Only 100 of 231 total aligned reads displayed.

ATCGACTGGCCCGAGGGGTTCTGCCGCCGCGACATCGGCTGGCGCGCGCTCTGAGACGAAAC.GGGGGCGCTGCCCCCCTCCT..G..CTTCGGCCAATTCACCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCG.AAATATCCCGGGGGAGCCCGCTTGCGGGCGGGGGCAGCGCCCCCGA  >  NC_014034/2566475‑2566672
                                                                                                      |                                                                                                     
 tCGACTGGCCCGGGGGGTTCTGCCGCCGCGACATCGGCTGGCGCGCGCTCTGAGACGAAACGGGGGGCGCTGCCCCCCTCTT..GGCCTTCGGCCAATTCAcc                                                                                                      <  1:2288985/101‑1 (MQ=25)
   gACTGGCCCGAGGGGTTCTGCCGCCGCGACATCGGCTGGCGCGCGCTCTGAGACGAAACGGGGGGCGCTGCCCCCCTCTTGGG..CCTCGGGCACTTCAcccc                                                                                                    >  2:546213/1‑101 (MQ=255)
    aCTGGCCCGAGGGGTTCTGCCGCCGCGACATCGGCTGGCGCGCGCTCTGAGACGAAACGGGGGGCGCTGCCCCCCTCTT..GGCCTTCGGCCAATTCAccccc                                                                                                   <  2:4887723/101‑1 (MQ=25)
      tGGCCCGAGGGGTTCTGCCGCCGCGACATCGGCTGGCGCGCGCTCTGAGACGAAACGGGGGGCGCTGCCCCCCTCTT..GGCCTTCGGCCAATTCACCCCCCg                                                                                                 <  2:4553621/101‑1 (MQ=25)
        gCCCGAGGGGTTCTGCCGCCGCGACATCGGCTGTCGCGCGCTCTGAGACGAAACGGGGGGCGCTGCCCCCCTCTT..GGCCTTCGGCCAATTCACCCCCCggg                                                                                               <  1:4309561/101‑1 (MQ=25)
            gAGGGGTTCTGCCGCCGCGACATCGGCTGGCGCGCGCTCTGAGACGAAACGGGGGGCGCTGCCCCCCTCTT..GGCCTTCGGCCAATTCACCCCCCGGGatat                                                                                           <  1:3757308/101‑1 (MQ=25)
             aGGGGTTCTGCCGCCGCGACATCGGCTGGCGCGCGCTCTGAGACGAAACGGGGGGCGCTGCCCCCCTCTT..GGCCTTCGGCCAATTCACCCCCCGGGATAtt                                                                                          <  1:15296/101‑1 (MQ=21)
              ggggTTCTGCCGCCGCGACATCGGCTGGCGCGCGCTCTGAGACGAAACGGGGGGCGCTGCCCCCCTCTT..GGCCTTCGGCCAATTCACCCCCCGGGATAttt                                                                                         <  1:3399827/101‑1 (MQ=21)
               gggTTCTGCCGCCGCGACATCGGCTGGCGCGCGCTCTGAGACGAAACGGGGGGCGCTGCCCCCCTCTT..GGCCTTCGGCCAATTCACCCCCCGGGATATTTc                                                                                        <  2:3172939/101‑1 (MQ=21)
                ggTTCTGCCGCCGCGACATCGGCTGGCGCGCGCTCTGAGACGAAACGGGGGGCGCTGCCCCCCTCTT..GGCCTTCGGCCAATTCACCTCCCGGGATATTTCg                                                                                       <  2:4988316/101‑1 (MQ=21)
                ggTTCTGCCGCCGCGACATCGGCTGGCGCGCGCTCTGAGACGAAACGGGGGGCGCTGCCCCCCCCTT..GGCCTTCGGCCCATTCACCCCCCCGGAGATTTct                                                                                       >  2:201523/1‑100 (MQ=21)
                 gTTCTGCCGCCGCGACATCGGCTGGCGCGCGCTCTGAGACGAAACGGGGGGGGCTGCCCCCCTCTT..GGCCTTCGGCCCATTTCCCCCCCGGGGGATTTCgg                                                                                      >  1:2398914/1‑101 (MQ=21)
                 gTTCTGCCGCCGCGACATCGGCTGGCGCGCGCTCTGAGACGAAACGGGGGGCGCTGCCCCCCTCTT..GGCCTTCGGCCCAATCACCCCCCGGGATATTTTgg                                                                                      >  2:4745212/1‑101 (MQ=21)
                    cTGCCGCCGCGACTTCGGCTGGCGCGCGCTCTGAGACGAAACGGGGGGCGCTGCCCCCCTCTT..GGCCTTCGGCCAATTCACCCCCCGGGATATTTCGGGCa                                                                                   <  2:213462/101‑1 (MQ=21)
                    cTGCCGCCGCGACATCGGCTGGCGCGCGCTCTGAGACGAAACGGGGGGCGCTGCCCCCCTCTT..GGCCTTCGGCCAATTCACCCCCCGGGATATTTCGGGCa                                                                                   <  1:1928170/101‑1 (MQ=21)
                      gcctccGCGACATCGGCTGGCGCGCGTTCTGAGACGAAACGGGGGGCGCTGCCCCCCTCTT..GGCCTTCGGCCAATTCACCCCCCGGGATATTTCGGGCAAg                                                                                 <  2:558382/101‑1 (MQ=18)
                      gccgccGCGACATCGGCTGGCGCGCGCTCTGAGACGAAACGGGGGGCGCTGCCCCCCTCTT..GGCCTTCGGCCAATTCACCCCCCGGGATATTTCGGGCAAg                                                                                 <  1:4004963/101‑1 (MQ=21)
                      gccgccGCGACATCGGCTGGCGCGCGCTCTGAGACGAAACGGGGGGCGCTGCCCCCCTCTT..GGCCTTCGGCCAATTCACCCCCCGGGATATTTCGGGCAAg                                                                                 <  1:1811994/101‑1 (MQ=21)
                        cgccgcGACATCGGCTGGCGCGCGCTCTGAGACGAAACGGGGGGCGCTGCCCCCCTCTT..GGCCTTCGGCCAATTCACCCCCCGGGATATTTCGGGCAAGAt                                                                               <  1:4733697/101‑1 (MQ=21)
                        cgccgcGACATCGGCTGGCGCGCGCTCTGAGACGAAACGGGGGGCGCTGCCCCCCTCTT..GGCCTTCGGCCAATTCACCCCCCGGGATATTTCGGGCAAGAt                                                                               <  1:3134645/101‑1 (MQ=21)
                         gccgcGACATCGGCTGGCGCGCGCTCTGAGACGAAACGGGGGGCGCTGCCCCCCTCTT..GGCCTTCGGCCAATTCACCCCCCGGGATATTTCGGGCAAGATa                                                                              <  1:758907/101‑1 (MQ=21)
                         gccgcGACATCGGCTGGCGCGCGCTCTGAGACGAAACGGGGGGCGCTGCCCCCCTCTT..GGCCTTCGGCCAATTCACCCCCCGGGATATTTCGGGCAAGATa                                                                              <  2:2373027/101‑1 (MQ=21)
                          ccgcGACATCGGCTGGCGCGCGCTCTGAGACGAAACGGGGGGGGGCGCCCCCCTCCT..GGCCTTTCGCCCATTCACCCCCCGGGATACTTCGGGCGAGACaa                                                                             >  2:812208/1‑101 (MQ=16)
                          ccgcGACATCGGCTGGCGCGCGCTCTGAGACGAAACGGGGGGCGCTGCCCCCCTCTT..GGCCTTCGGCCAATTCACCCCCCGGGATATTTCGGGCAAGATaa                                                                             <  2:2169155/101‑1 (MQ=21)
                            gtgACTTCGGCGGGCGCGCGCTCTGAGACAAAACGGGGGGCGCTGCCCCCCTCTT..GGCCTTCGGCCAATTCCCCCCCCGGGATATTTCGGGCAAGATAaag                                                                           <  1:2206603/99‑1 (MQ=18)
                            gcgACATCGGCTGGCGCGCGCTCTGAGACGAAACGGGGGGCGCTGCCCCCCTCTT..GGCCTTCGGCCAATTCACCCCCCGGGATATTTCGGGCAAGATAaag                                                                           <  1:3001494/101‑1 (MQ=21)
                            gcgACATCGGCTGGCGCGCGCTCTGAGACGAAACGGGGGGCGCTGCCCCCCTCTT..GGCCTTCGGCCAATTCACCCCCCGGGATATTTCGGGCAAGATAaag                                                                           <  1:2556790/101‑1 (MQ=21)
                             cgACATCGGCTGGCGCGCGCTCTGAGACGAAACGGGGGGCGCTGCCCCCCTCTT..GGCCTTCGGCCAATTCACCCCCCGGGCTATATAGGGGAAGATaaaaa                                                                          >  2:3690992/1‑99 (MQ=18)
                             cgACATCGGCTGGCGCGCGCTCTGAGACGAAACGGGGGGCGCTGCCCCCCTCTT..GGCCTTCGGCCAATTCACCCCCCGGGATATTTCGGGCAAGATAaaga                                                                          <  2:2589035/101‑1 (MQ=21)
                              gACATCGGCTGGCGCGCGCTCTGGGACGAAACGGGGGGCGCTGCCCCCCTCTT..GGCCTTCGGCCAATTCACCCCCCGGGATATTTCGGGCAAGATAaagaa                                                                         <  2:1840456/101‑1 (MQ=21)
                              gACATCGGCTGGCGCGCGCTCTGAGACGAAACGGGGGGCGCTGCCCCCCTCTT..GGCCTTCGGCCAATTCAACCCCCGGGATATTTCGGGCAAGATAaagaa                                                                         <  1:1375187/101‑1 (MQ=21)
                               aCATCGGCTGGCGCGCGCTCTGAGACGAAACGGGGGGCGCTGCCCCCCTCTT..GGCCTTCGGCCAATTCACCCCCCGGGATATTTTGGGCAAGATAaaaaag                                                                        >  1:1738075/1‑101 (MQ=18)
                                cATCGGCTGGCGCGCGCTCTGAGACGAAACGGGGGGCGCTGCCCCCCTCTT..GGCCTTCGGCCAATTCACCCCCCGGGATATTTCGGGCAAGATAAAGAAgg                                                                       <  2:3278843/101‑1 (MQ=21)
                                cATCGGCTGGCGCGCGCTCTGAGACGAAACGGGGGGCGCTGCCCCCCTCTT..GGCCTTCGGCCAATTCACCCCCCGGGATATTTCGGGCAAGATAAAGAAgg                                                                       <  2:610399/101‑1 (MQ=21)
                                  tCGGCTGGCGCGCGCTCTGAGACGAAACGGGGGGCGCTGCCCCCCTCTT..GGCCTTCGGCCAATTCACCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTc                                                                     <  2:2440729/101‑1 (MQ=21)
                                    ggCTGGCGCGCGCTCTGAGACGAAACGGGGGGCGCTGCCCCCCTCTT..GGCCTTCGGCCAATTCACCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCtt                                                                   <  2:337048/101‑1 (MQ=21)
                                    ggCTGGCGCGCGCTCTGAGACGAAACGGGGGGCGCTGCCCCCCTCTT..GGCCTTCGGCCAATTCACCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCtt                                                                   <  1:2939631/101‑1 (MQ=21)
                                    ggCTGGCGCGCCCTCTGAGACGAAACGGGGGGCGCTGCCCCCCTCTT..GGCCTTCGGCCAATTCACCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCtt                                                                   <  2:4794865/101‑1 (MQ=21)
                                       tGGCGCGCGCTCTGAGACGAAACGGGGGGCGCTGCCCCCCTCTT..GGCCTTCGGCCAATTCACCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTc                                                                <  2:1768660/101‑1 (MQ=21)
                                       tGGCGCGCGCTCTGAGACGAAACGGGGGGCGCTGCCCCCCTCTT..GGCCTTCGGCCAATTCACCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTc                                                                <  2:3994153/101‑1 (MQ=21)
                                       tGGCGCGCGCTCTGAGACGAAACGGGGGGCGCTGCCCCCCTCTT..GGCCTTCGGCCAATTCACCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTc                                                                >  2:4917854/1‑101 (MQ=21)
                                         gcgcgcgcTCTGAGACGAAACGGGGGGCGCTGCCCCCCTCTT..GGCCTTCGGCCAATTCACCCCCCGGGATATTTCGGGCAAGATAAAGAAGGGCTTTTCCt                                                              >  1:4505293/1‑101 (MQ=21)
                                           gcgcgcTCTGAGACGAAACGGGGGGCGCTGCCCCCCTCTT..GGCCTTCGGCCAATTCACCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCttt                                                            <  1:3186209/101‑1 (MQ=21)
                                            cgcgcTCTGAGACGAAACGGGGGGCGCTGCCCCCCTCTT..GGCCTTCGGCCAATTCACCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTc                                                           <  2:406436/101‑1 (MQ=21)
                                            cgcgcTCTGAGACGAAACGGGGGGCGCTGCCCCCCTCTT..GGCCTTCGGCCAATTCACCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTc                                                           <  2:192013/101‑1 (MQ=21)
                                            cgcgcTCTGAGACGAAACGGGGGGCGCTGCCCCCCTCTT..GGCCTTCGGCCAATTCACCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTc                                                           <  2:3302022/101‑1 (MQ=21)
                                             gcgcTCTGAGACGAAACGGGGGGCGCTGCCCCCCTCTT..GGCCTTCGGCCAATTCACCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCa                                                          <  2:426189/101‑1 (MQ=21)
                                              cgcTCTGAGACGAAACGGGGGGCGCTGCCCCCCTCTT..GGCCTTCGGCCAATTCACCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCAt                                                         <  1:3208980/101‑1 (MQ=21)
                                               gcTCTGAGACGAACCGGGGGGCGCTGCCCCCTTCTT..GGCCTTCGGCCAATTCACCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATc                                                        <  1:2652893/101‑1 (MQ=21)
                                               gcTCTGAGACGAAACGGGGGGCGCTGCCCCCCTCTT..GGCCTTCGGCCAATTCACCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATc                                                        <  1:2154159/101‑1 (MQ=21)
                                               gcTCTGAGACGAAACGGGGGGCGCTGCCCCCCTCTT..GGCCTTCGGCCAATTCACCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATc                                                        <  2:2891000/101‑1 (MQ=21)
                                                ctctGAGACGAAACGGGGGGCGCTGCCCCCCTCTT..GGCCTTCGGCCAATTCACCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCt                                                       >  2:4678419/1‑101 (MQ=21)
                                                ctctGAGACGAAACGGGGGGCGCTGCCCCCCTCTT..GGCCTTCGGCCAATTCACCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCt                                                       <  2:3411039/101‑1 (MQ=21)
                                                 tctGAGACGAAACGGGGGGCGCTGCCCCCCTCTT..GGCCTTCGGCCAATTCACCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCtt                                                      <  1:88760/101‑1 (MQ=21)
                                                 tctGAGACGAAACGGGGGGCGCTGCCCCCCTCTT..GGCCTTCGGCCAATTCACCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTCTCCTTTCATCtt                                                      <  1:2728105/101‑1 (MQ=21)
                                                  ctGAGGCGAAACGGGGGGCGCTGCCCCCCTCTT..GGCCTTCGGCCAATTCTCCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTg                                                     <  1:2566418/101‑1 (MQ=21)
                                                  ctGAGACGAAACGGGGGGCGCTGCCCCCCTCTT..GGCCTTCGGCCAATTCACCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTg                                                     <  1:582768/101‑1 (MQ=21)
                                                   tGAGACGAAACGGGGGGCGCTGCCCCCCTCTT..GGCCTTCGGCCAATTCACCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGc                                                    <  1:4249700/101‑1 (MQ=21)
                                                   tGAGACGAAACGGGGGGCGCTGCCCCCCTCTT..GGCCTTCGGCCAATTCACCCCCCGGGAAATCTCGGGCAAGATAAAGAAGGTCGTTTCCTTTCATCTTGc                                                    >  1:4849067/1‑101 (MQ=21)
                                                    gagaCGAAACGGGGGGCGCTGCCCCCCTCTT..GGCCTTCGGCCAATTCACCCCCCCGGAGATTTTGGGGCAGATAAAAAAGGGGTTTTTCTTTTATTttttc                                                   >  1:3572030/1‑98 (MQ=2)
                                                     agaCGAAACGGGGGGCGCTGCCCCCCTCTT..GGCCTTCGGCCAATTCACCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGccc                                                  <  2:1699276/101‑1 (MQ=21)
                                                      gaCGAAACGGGGGGCGCTGCCCCCCTCTT..GGCCTTCGGCCAATTCACCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCg                                                 <  1:1894315/101‑1 (MQ=21)
                                                      gaCGAAACGGGGGGCGCTGCCCCCCTCTT..GGCCTTCGGCCAATTCACCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCg                                                 <  1:514012/101‑1 (MQ=21)
                                                        cGAAACGGGGGGCGCTGCCCCCCTCTT..GGCCTTCGGCCAATTCACCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCG.aa                                              <  1:1617993/101‑1 (MQ=21)
                                                        cGAAACGGGGGGCGCTGCCCCCCTCTT..GGCCTTCGGCCAATTCACCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCG.aa                                              >  1:3357414/1‑101 (MQ=21)
                                                        cGAAACGGGGGGCGCTGCCCCCCTCTT..GGCCTTCGGCCAATTCACCCCCCGGGACATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCG.aa                                              >  2:613597/1‑101 (MQ=21)
                                                         gAAACGGGGGGCGCTGCCCCCCTCTT..GGCCTTCGGCCAATTCACCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCG.aaa                                             >  2:229319/1‑101 (MQ=21)
                                                          aaacggggggCGCTGCCCCCCTCTT..GGCCTTCGGCCAATTCACCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCG.AAat                                            <  1:3944356/96‑1 (MQ=21)
                                                          aaacggggggCGCTGCCCCCCTCTT..GGCCTTCGGCCAATTCACCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCG.AAat                                            <  2:107585/96‑1 (MQ=21)
                                                          aaacggggggCGCTGCCCCCCTCTT..GGCCTTCGGCCAATTCACCCCCCGGGAAATTCCGGCCAAAAAAAGAAAGGACTTTTCTTTTCTTCTTCCCCG.aaaa                                            >  2:4359411/6‑100 (MQ=2)
                                                          aaacggggggCGCTGCCCCCCTCTT..GGCCTTCGGCCAATTCACCCCCCCGGATATTTCGGGCAAGATAAAGAAGGGCTTTTCCTTTCATCTTGCCCC.AAat                                            >  1:2033105/6‑101 (MQ=21)
                                                           ggcagggggCGCTG.CCCCCTCTT..GGCCTTTGGCCAATTCA.CCCCTGGGATATTTCGGGCAAGAT...GAA.ATCATGCCCCTTCATCTTGCCGC.AAATATCCCgg                                     >  1:4977195/5‑101 (MQ=0)
                                                            acggggggCGCTGCCCCCCTCTT..GGCCTTCGGCCAATTCTCCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCG.AAatat                                          <  2:3198269/98‑1 (MQ=25)
                                                            acggggggCGCTGCCCCCCTCTT..GGCCTTCGGCCAATTCACCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCG.AAatat                                          <  2:908820/98‑1 (MQ=25)
                                                            acggggggCGCTGCCCCCCTCTT..GGCCTTCGGCCAATTCACCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCG.AAatat                                          <  1:2129400/98‑1 (MQ=25)
                                                             cggggggCGCTGCCCCCCTCTT..GGCCTTCGGCCAATTCACCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCG.AAATATc                                         <  2:1023372/99‑1 (MQ=25)
                                                             cggggggCGCTGCCCCCCTCTT..GGCCTTCGGCCAATTCACCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCG.AAATATc                                         >  2:3038455/3‑101 (MQ=25)
                                                             cggggggCGCTGCCCCCCTCTT..GGCCTTCGGCCAATTCACCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTACTTTCATCTTGCCCG.AAATATc                                         >  2:1100723/3‑101 (MQ=21)
                                                               gggggCGCTG.CCCCCTCTT..GGCCTTTGGCCAATTCA.CCCCTGGGATATTTCGGGCAAGAT...GAA.ATCATGCCCCTTCATCTTGCCGC.AAATATCCCgggggg                                 >  2:1480212/1‑100 (MQ=0)
                                                                ggggCGCTGCCCCCCTCTT..GGCCTTCGGCCAATTCACCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCG.AAATATCCCg                                      <  2:2456950/101‑1 (MQ=25)
                                                                  ggCGCTGCCCCCCTCTT..GGCCTTCGGCCAATTCACCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCG.AAATATCCCggg                                    >  2:2674843/1‑101 (MQ=25)
                                                                   gcgcTGCCCCCCTCTT..GGCCTTCGGCCAATTCACCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCG.AAATATCCCgggg                                   >  1:1121439/1‑101 (MQ=37)
                                                                        gCCCCCCTCTT..GGCCTTCGGCCAATTCTCCCCCCGGGATATTTCGGTCAAGATAAAGAAGGTCTTTTCCTTTCATCTTTCCCG.AAATATCCCGGGGGAGcc                              <  1:709533/101‑1 (MQ=25)
                                                                        gCCCCCCTCTT..GGCCTTCGGCCAATTCACCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCG.AAATATCCCGGGGGAGcc                              <  2:4036376/101‑1 (MQ=37)
                                                                               cctcttggccTTCGGCCAATTCACCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCG.AAATATCCCGGGGGAGCCCGCtt                         <  2:3833724/92‑1 (MQ=255)
                                                                               cctcctggccTTCGCCCAATTCCCCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCG.AAATATCCCGGGGGAGCCCGCtt                         <  1:3930720/92‑1 (MQ=37)
                                                                                ctcttggccTTCGGCCAATTCACCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCG.AAATATCCCGGGGGAGCCCGCTTg                        <  1:3269055/93‑1 (MQ=255)
                                                                                  cttggccTTCGGCCAATTCACCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCG.AAATATCCCGGGGCAGCCGGCTTccg                      >  2:2502701/7‑101 (MQ=37)
                                                                                    tggccTTCGGCCAATTCCCCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCG.AAATATCCCGGGGGAGCCCGCTTgcggg                    <  2:4504973/97‑1 (MQ=255)
                                                                                     ggccTTCGGCCAATTCACCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCG.AAATATCCCGGGGGAGCCCGCTTgcgggc                   >  1:2694411/4‑101 (MQ=255)
                                                                                     ggccTTCGGCCAATTCACCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCG.AAATATCCCGGGGGAGCCCGCTTgcgggc                   >  2:2538902/4‑101 (MQ=255)
                                                                                           cGGCCAATTCACCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCG.AAATATCCCGGGGGAGCCCGCTTGCGGGCGGGGGc             <  1:4162774/101‑1 (MQ=255)
                                                                                            ggCCAATTCACCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCG.AAATATCCCGGGGGAGCCCGCTTGCGGGCGGGGGCa            >  2:1382929/1‑101 (MQ=255)
                                                                                             ggcAAATCCCCCCCCGGGGTTTTTTCGGGAAAAAAAAAAAGGGTTTTTTCTTTCCTTTTTCCCGAAAATTCCCCGGGGGAGCCCGCTTGCGGGCGGGGGCa            <  1:1020567/99‑1 (MQ=255)
                                                                                             gCCAATTCACCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCATTTCATCTTGCCCG.AAATATCCCGGGGGGGCCCGCTTGGGGGAggggggga           >  1:3076679/1‑98 (MQ=37)
                                                                                                 aTTCACCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCG.AAATATCCCGGGGGAGCCCGCTTGCGGGGGGGGGCAGCGcc       >  2:2375588/1‑101 (MQ=255)
                                                                                                 aTTCACCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCG.AAATATCCCGGGGGAGCCCGCTTGCGGGCGGGGGGAGGGcc       >  1:2966624/1‑101 (MQ=255)
                                                                                                  ttCACCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCG.AAATATCCCGGGGGAGCCCGCTTGCGGGCGGGGGCAGCGccc      >  2:1581810/1‑101 (MQ=255)
                                                                                                     gccccccGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCG.AAATATCCCGGGGGAGCCCGCTTGCGGGCGGGGGCAGCGCCCCCg   <  1:554365/100‑1 (MQ=255)
                                                                                                      ccccccGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCG.AAATATCCCGGGGGAGCCCGCTTGCGGTCGGGGGCAGCGCCCCCGa  <  1:3691270/101‑1 (MQ=255)
                                                                                                      |                                                                                                     
ATCGACTGGCCCGAGGGGTTCTGCCGCCGCGACATCGGCTGGCGCGCGCTCTGAGACGAAAC.GGGGGCGCTGCCCCCCTCCT..G..CTTCGGCCAATTCACCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCG.AAATATCCCGGGGGAGCCCGCTTGCGGGCGGGGGCAGCGCCCCCGA  >  NC_014034/2566475‑2566672

Base quality scores:  ATCG  < 3 ≤  ATCG  < 15 ≤  ATCG  < 25 ≤  ATCG  < 33 ≤  ATCG  < 41 ≤  ATCG