Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_014034 2,566,561 C→T 5.8% intergenic (+33/+161) 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,5610C→T 5.8% 11.2 ‑3.0 ‑0.8120intergenic (+33/+161)purD/sdaAphosphoribosylamine‑‑glycine ligase/L‑serine ammonia‑lyase
Reads supporting (aligned to +/- strand):  new base (4/3):  ref base (24/89):  total (28/92)
Fisher's exact test strand distribution p-value = 5.07e-02
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 1.43e-01

Only 100 of 138 total aligned reads displayed.

GATGCGATCGACTGGCCCGAGGGGTTCTGCCGCCGCGACATCGGCTGGCGCGCGCTCTGAGACGAAAC.GGG.GGCGCTGCCCCC.CTC.CT..G..CTTC.GGCC.AATTCACCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCCGGGGGAGCCCGCTTGCGGGCGGGGGC  >  NC_014034/2566469‑2566661
                                                                                                    |                                                                                                      
gATGCGCTCGGCTGGCCCGAGGGGTTCTGCCGCCGCGACATCGGCTGGCGCGCGCTCTGAGACGAAACGGGG.GGCGCTGCCCCC.CTC.TT..GGCCTTC.GGCC.a                                                                                                 <  1:705356/101‑1 (MQ=25)
 aTGCGATCGACTGGCCCGAGGGGTTCTGCCGCCGCGACATCGGCTGGCGCGCGCTCTGAGACGAAACGGGG.GGCGCTGCCCCC.CTC.TT..GGCCTTC.GGCC.aa                                                                                                >  2:2028863/1‑101 (MQ=37)
   gCGATCGACTGGCCCGAGGGGTTCTGCCGCCGCGACATCGGCTGGCGCGCGCTCTGAGACGAAACGGGG.GGCGCTGCCCCC.CTC.TT..GGCCTTC.GGCC.AAtt                                                                                              <  1:1545824/101‑1 (MQ=37)
   gCGATCGACTGGCCCGAGGGGTTCTGCCGCCGCGACATCGGCTGGCGCGCGCTCTGAGACGAAACGGGG.GGCGCTGCCCCC.CTC.TT..GGCCTTC.GGCC.AAtt                                                                                              <  1:3672165/101‑1 (MQ=37)
   gCGATCGACTGGCCCGAGGGGTTCTGCCGCCGCGACATCGGCTGGCGCGCGCTCTGAGACGAAACGGGG.GGCGCTGCCCCC.CTC.TA..GGCCTTC.GGCC.AAtt                                                                                              <  1:3566316/101‑1 (MQ=37)
    cGATCGACTGGCCCGAGGGGTTCTGCCGCCGCGACATCGGCTGGCGCGCGCTCTGAGACGAAACGGGG.GGCGCTGCCCCC.CTC.TT..GGCCTTC.GGCC.AATTc                                                                                             <  1:2286059/101‑1 (MQ=37)
    cGATCGACTGGCCCGAGGGGTTCTGCCGCCGCGACATCGGCTGGCGCGCGCTCTGAGACGAAACGGGG.GGCGCTGCCCCC.CTC.TT..GGCCTTC.GGCC.AATTc                                                                                             <  1:3392828/101‑1 (MQ=37)
        cGACTGGCCCGAGGGGTTCTGCCGCCGCGACATCGGCTGGCGCGCGCTCTGAGACGAAACGGGG.GGCGCTGCCCCC.CTC.TT..GGCCTTC.GGCC.AATTCAccc                                                                                         <  2:659481/101‑1 (MQ=37)
          gcgggCCCGAGGGGTTCTGCCGCCGCGACATCGGCTGGCGCGCGCTCTGAGACGAAACGGGG.GGCGCTGCCCCC.CTC.TT..GGCCTTC.GGCC.AATTCAccccc                                                                                       <  2:2391411/98‑1 (MQ=25)
          aCTGGCCCGAGGGGTTCTGCCGCCGCGACATCGGCTGGCGCGCGCTCTGAGACGAAACGGGG.GGCGCTGCCCCC.CTC.TT..GGCCTTC.GGCC.AATTCAccccc                                                                                       <  1:1187826/101‑1 (MQ=25)
           cTGGCCCGAGGGGTTCTGCCGCCGCGACATCGGCTGGCGCGCGCTCTGAGACGAAACGGGG.GGCGCTGCCCCC.CTC.TT..GGCCTTC.GGCC.AATTCAcccccc                                                                                      <  2:2895798/101‑1 (MQ=25)
           cTGGCCCGAGGGGTTCTGCCGCCGCGACATCGGCTGGCGCGCGCTCTGAGACGAAACGGGG.GGCGCTGCCCCC.CTC.TT..GGCCTTC.GGCC.AATTCAcccccc                                                                                      <  2:702261/101‑1 (MQ=25)
            tGGCCCGAGGGGTTCTGCCGCCGCGACATCGGCTGGCGCGCGCTCTGAGACGAAACGGGG.GGCGCTGCCCCC.CTC.TT..GGCCTTC.GGCC.AATTCACCCCCCg                                                                                     <  1:2181231/101‑1 (MQ=25)
             ggCCCGAGGGGTTCTGCCGCCGCGACATCGGCTGGCGCGCGCTCTGAGACGAAACGGGG.GGCGCTGCCCCC.CTC.TT..GGCCTTC.GGCC.AATTCACCCCCCgg                                                                                    <  1:2944444/101‑1 (MQ=25)
             ggCCCGAGGGGTTCTGCCGCCGCGACATCGGCTGGCGCGCGCTCTGAGACGAAACGGGG.GGCGCTGCCCCC.CTC.TT..GGCCTTC.GGCC.AATTCACCCCCCgg                                                                                    <  1:1099121/101‑1 (MQ=25)
             ggCCCGAGGGGTTCTGCCGCCGCGACATCGGCTGGCGCGCGCTCTGAGACGAAACGGGG.GGCGCTGCCCCC.CTC.TT..GGCCTTC.GGCC.AATTCACCCCCCgg                                                                                    <  1:503621/101‑1 (MQ=25)
              gCCCGAGGGGTTCTGCCGCCGCGACATCGGCTGGCGCGCGCTCTGAGACGAAAC.GGGTGGGGCTGCCCCC.CTC.TT..GGCCCTC.CGCC.AATTCACCCCCCcgg                                                                                   >  2:3490303/1‑101 (MQ=25)
               cccGAGGGGTTCTGCCGCCGCGACATCGGCTGGCGCGCGCTCTGAGACGAAACGGGG.GGCGCTGCCCCC.CTC.TT..GGCCTTC.GGCC.AATTCACCCCCCGGGa                                                                                  <  1:1829286/101‑1 (MQ=25)
                ccGAGGGGTTCTGCCGCCGCGACATCGGCTGGCGCGCGCTCTGAGACGAAACGGGG.GGGGCCGCCCCC.CCCTTT..GC.CTTC.GGCC.AAATCACCCCCCcggga                                                                                 >  1:1978096/1‑99 (MQ=255)
                ccGAGGGGTTCTGCCGCCGCGACATCGGCTGGCGCGCGCTCTGAGACGAAACGGGG.GGCGCTGCCCCC.CTC.TT..GGCCTTC.GGCC.AATTCACCCCCCGGGat                                                                                 <  2:2168907/101‑1 (MQ=25)
                ccGAGGGGTTCTGCCGCCGCGACATCGGCTGGCGCGCGCTCTGAGACGAAACGGGG.GGCGCTGCCCCC.CTC.TT..GGCCTTC.GGCC.AATTCACCCCCCGGGat                                                                                 <  1:3194193/101‑1 (MQ=25)
                ccGAGGGGTTCTGCCGCCGCGACATCGGCTGGCGCGCGCTCTGAGACGAAACGGGG.GGCGCTGCCCCC.CTC.TT..GGCCTTC.GGCC.AATTCACCCCCCGGGat                                                                                 <  2:4180170/101‑1 (MQ=25)
                 cGAGGGGTTCTGCCGCCGCGACATCGGCTGGCGCGCGCTCTGAGACGAAACGGGG.GGGGCTGCCCCC.CCC.TT..GGCCTTC.CGCC.AAATCCCCCCCCCGGgta                                                                                >  2:3168341/1‑101 (MQ=255)
                 cGAGGGGTTCTGCCGCCGCGACATCGGCTGGCGCGCGCTCTGAGACGAAACGGGG.GGCGCTGCCCCC.CTC.TT..GGCCTTC.GGCC.AATTCACCCCCCGGGata                                                                                <  1:1708987/101‑1 (MQ=25)
                  gAGGGGTTCTGCCGCCGCGACATCGGCTGGCGCGCGCTCTGAGACGAAACGGGG.GGCGCTGCCCCC.CTC.TT..GGCCTTC.GGCC.AATTCACCCCCCGGGatat                                                                               <  2:4131168/101‑1 (MQ=25)
                  gAGGGGTTCTGCCGCCGCGACATCGGCTGGCGCGCGCTCTGAGACGAAACGGGG.GGCGCTGCCCCC.CTC.TT..GGCCTTC.GGCC.AATTCACCCCCCCGGatat                                                                               >  1:2119448/1‑101 (MQ=25)
                   aGGGGTTCTGCCGCCGCGACATCGGCTGGCGCGCGCTCTGAGACGAAACGGGG.GGCGCTGCCCCC.CTC.TT..GGCCTTC.GGCC.AATTCAACCCCCGGGATAtt                                                                              >  2:1052333/1‑101 (MQ=21)
                    ggggTTCTGCCGCCGCGACATCGGCTGGCGCGCGCTCTGAGACGAAACGGGG.GGCGCTGCCCCC.CTC.TT..GGCCTTC.GGCC.AATTCACCCCCCGGGATAttt                                                                             <  1:1697153/101‑1 (MQ=21)
                    ggggTTCTGCCGCCGCGACATCGGCTGGCGCGCGCTCTGAGACGAAACGGGG.GGCGCTGCCCCC.CTC.TT..GGCCTTC.GGCC.AATTCACCCCCCGGGATAttt                                                                             <  2:77585/101‑1 (MQ=21)
                        ttCTGCCGCCGCGACATCGGCTGGCGCGCGCTCTGAGACGAAACGGGG.GGCGCTGCCCCC.CTC.TT..GGCCTTC.GGCC.AATTCACCCCCCGGGATATTTCggg                                                                         <  2:2145273/101‑1 (MQ=21)
                         tCTGCCGCCGCGACATCGGCTGGCGCGCGCTCTGAGACGAAACGGGG.GGCGCTGCCCCC.CTC.TT..GGCCTTC.GGCC.CATTAAACCCCCCGGATATTTTgggg                                                                        >  2:1592017/1‑100 (MQ=21)
                         tCTGCCGCCGCGACATCGGCTGGCGCGCGCTCTGAGACGAAACGGGG.GGCGCTGCCCCC.CTC.TT..GGCCTTC.GGCC.AATTCACCCCCCGGGATATTTCGGGc                                                                        <  1:1268244/101‑1 (MQ=21)
                         tCTGCCGCCGCGACATCGGCTGGCGCGCGCTCTGAGACGAAACGGGG.GGCGCCGCCCCC.CTC.TTTGG..CCTC.GGCCAAATCCACCCCCCcggtacttttggg                                                                         >  2:3685661/1‑91 (MQ=255)
                          cTGCCGCCGCGACATCGGCTGGCGCGCGCTCTGAGACGAAACGGGG.GGCGCTGCCCCC.CTC.TT..GGCCTTC.GGCC.AATTCACCCCCCGGGATATTTCGGGCa                                                                       <  2:2900729/101‑1 (MQ=21)
                          cTGCCGCCGCGACATCGGCTGGCGCGCGCTCTGAGACGAAACGGGG.GGCGCTGCCCCC.CTC.TT..GGCCTTC.GGCC.AATTCACCCCCCCGGCTATTTCGGGCa                                                                       >  2:861613/1‑101 (MQ=21)
                           tGCCGCCGCGACATCGGCTGGCGCGCGCTCTGAGACGAAGCGGGG.GGCGCTGCCCCC.CTC.TT..GGCCTTC.GGCC.AATTCACCCCCCGGGATATTTCGGGCaa                                                                      <  1:473092/101‑1 (MQ=21)
                           tGCCGCCGCGACATCGGCTGGCGCGCGCTCTGAGACGAAACGGGG.GGCGCTGCCCCC.CCC.TT..GGCCTTC.GGCC.AATTCACCCCCCCGGAGAtttttggcca                                                                      >  2:3957956/1‑94 (MQ=18)
                            gccgccGCGACATCGGCTGGCGCGCGCTCTGAGACGAAACGGGG.GGCGCTGCCCCC.CTC.TT..GGCCTTC.GGCC.AATTCACCCCCCGGGATATTTCGGGCAAg                                                                     <  1:45837/101‑1 (MQ=21)
                              cgccgcGACATCGGCTGGCGCGCGCTCTGAGACGAAACGGGG.GGCGCTGCCCCC.CTC.TT..GGCCTTC.GGCC.AATTCACCCCCCGGGATATTTCGGGCAAGAt                                                                   <  2:1620144/101‑1 (MQ=21)
                              cgccgcGACATCGGCTGGCGCGCGCTCTGAGACGAAACGGGG.GGCGCTGCCCCC.CTC.TT..GGCCTTC.GGCC.AATTCACCCCCCGGGATATTTCGGGCAAGAt                                                                   <  2:2060072/101‑1 (MQ=21)
                              cgccgcGACATCGGCTGGCGCGCGCTCTGAGACGAAACGGGG.GGCGCTGCCCCC.CTC.TT..GGCCTTC.GGCC.AATTCACCCCCCGGGATATTTCGGGCAAGAt                                                                   <  2:3164/101‑1 (MQ=21)
                               gccgcGACATCGGCTGGCGCGCGCTCTGAGACGAAACGGGG.GGCGCTGCCCCC.CTC.TT..GGCCTTC.GGCC.AATTCACCCCCCGGGATATTTCGGGCAAGATa                                                                  <  2:2066625/101‑1 (MQ=21)
                               gccgcGACATCGGCTGGCGCGCGCTCTGAGACGAAACGGGG.GGCGCTGCCCCC.CTC.TT..GGCCTTC.GGCC.AATTCACCCCCCGGGATATTTCGGGCAAGATa                                                                  <  2:3108432/101‑1 (MQ=21)
                                 cgcgACATCGGCTGGCGCGCGCTCTGAGACGAAACGGGG.GGCGCTGCCCCC.CTC.TT..GGCCTTC.GGCC.CATTCACCCCCCCGGATATTTCGGGCAAGATaaa                                                                >  1:3906941/1‑101 (MQ=21)
                                  gcgACATCGGCTGGCGCGCGCTCTGAGACGAAACGGGG.GGCGCTGCCCCC.CTC.TT..GGCCTTC.GGCC.AATTCACCCCCCGGGATATTTCGGGCAAGATAaag                                                               <  1:430223/101‑1 (MQ=21)
                                   cgACATCGGCTGGCGCGCGCTCTGAGACGAAACGGGG.GGCGCTGCCCCC.CTC.TT..GGCCTTC.GGCC.AATTCACCCCCCGGGATATTTCGGGCAAGATAaaga                                                              >  2:3176725/1‑101 (MQ=21)
                                   cgACATCGGCTGGCGCGCGCTCTGAGACGAAACGGGG.GGCGCTGCCCCC.CTC.TT..GGCCTTC.GGCC.AATTCACCCCCCGGGATATTTCGGGCAAGATAaaga                                                              <  1:1811536/101‑1 (MQ=21)
                                     aCATCGGCTGGCGCTCGCTCTGAGACGAAACGGGG.GGCCCCGCCCCC.CTT.TTG.G..CCTCGGGCC.AATTCACCCCCCGCGCCGTTTCGGGGGCGATAaaaaag                                                            >  1:4115095/1‑101 (MQ=9)
                                       aTCGGCTGGCGCGCGCTCTGAGACGAAACGGGG.GGCGCTGCCCCC.CTC.TT..GGCCTTC.GGCC.CATTCACCCCCCCGGAAATTTTGGGGAAGATAAAGAAggg                                                          >  2:3039499/1‑100 (MQ=18)
                                        tCGGCTGGCGCGCGCTCTGAGACGAAACGGGG.GGCGCTGCCCCC.CTC.TT..GGCCTTC.GGCC.AATTCACCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTc                                                         >  1:3848117/1‑101 (MQ=21)
                                        tCGGCTGGCGCGCGCTCTGAGACGAAACGGGG.GGCGCTGCCCCC.CTC.TT..GGCCTTC.GGCC.AATTCACCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTc                                                         <  1:3471783/101‑1 (MQ=21)
                                          ggCTGGCGCGCCCTCTGAGACGAAACGGGG.GGCGCTGCCCCC.CTC.TT..GGCCTTC.GGCC.AATTCACCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCtt                                                       <  1:1943026/101‑1 (MQ=21)
                                           gCTGGCGCGCGCTCTGAGACGAAACGGGG.GGCGCTGCCCCC.CTC.TT..GGCCTTC.GGCC.AATTCACCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCttt                                                      >  1:2196690/1‑101 (MQ=21)
                                           gCTGGCGCGCGCTCTGAGACGAAACGGGG.GGCGCTGCCCCC.CTC.TT..GGCCTTC.GGCC.AATTCACCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCttt                                                      <  2:1078597/101‑1 (MQ=21)
                                            cTGGCGCGCGCTCTGAGACGAAACGGGG.GGCGCTGCCCCC.CTC.TT..GGCCTTC.GGCC.AATTCACCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCtttt                                                     <  1:2169706/101‑1 (MQ=21)
                                             tGGCGCGCGCTCTGAGACGAAACGGGG.GGCGCTGCCCCC.CTC.TT..GGCCTTC.GGCC.AATTCACCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTc                                                    >  2:1009588/1‑101 (MQ=21)
                                               gcgcgcgcTCTGAGACGAAACGGGG.GGCGCTGCCCCC.CTC.TT..GGCCTTC.GGCC.AATTCACCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCt                                                  <  1:3376059/101‑1 (MQ=21)
                                               gcgcgcgcTCTGAGACGAAACGGGG.GGCGCTGCCCCC.CTC.TT..GGCCTTC.GGCC.AATTCACCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCt                                                  <  1:3727302/101‑1 (MQ=21)
                                                cgcgcgcTCTGAGACGAAACGGGG.GGCGCTGCCCCC.CTC.TT..GGCCTTC.GGCC.AATTCACCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCtt                                                 <  1:2512408/101‑1 (MQ=21)
                                                 gcgcgcTCTGAGACGAAACGGGG.GGCGCTGCCCCC.CTC.TT..GGCCTTC.GGCC.AATTCACCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCttt                                                <  2:2428260/101‑1 (MQ=21)
                                                 gcgcgcTCTGAGACGAAACGGGG.GGCGCTGCCCCC.CTC.TT..GGCCTTC.GGCC.AATTCACCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCttt                                                <  2:1830164/101‑1 (MQ=21)
                                                  cgcgcTCTGAGACGAAACGGGG.GGCGCTGCCCCC.CTC.TT..GGCCTTC.GGCC.AATTCACCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTc                                               <  1:2076249/101‑1 (MQ=21)
                                                  cgcgcTCTGAGACGAAACGGGG.GGCGCTGCCCCC.CTC.TT..GGCCTTC.GGCC.AATTCACCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTc                                               <  1:509035/101‑1 (MQ=21)
                                                     gcTCTGAGACGAAACGGGG.GGCGCTGCCCCC.CTC.TT..GGCCTTC.GGCC.AATTCACCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATc                                            <  1:22833/101‑1 (MQ=21)
                                                     gcTCTGAGACGAAACGGGG.GGCGCTGCCCCC.CTC.TT..GGCCTTC.GGCC.AATTCACCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATc                                            <  1:31132/101‑1 (MQ=21)
                                                      ctctGAGACGAAACGGGG.GGCGCTGCCCCC.CTC.TT..GGCCTTC.GGCC.AATTCACCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCt                                           <  1:2775027/101‑1 (MQ=21)
                                                         tGAGACGAAACGGGG.GGCGCTGCCCCC.CTC.TT..GGCCTTC.GGCC.AATTCACCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTGTTTTTTTTTCTTCTTGc                                        >  2:853783/1‑101 (MQ=18)
                                                         tGAGACGAAACGGGG.GGCGCTGCCCCC.CTC.TT..GGCCTTC.GGCC.AATTCACCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGc                                        >  2:2272230/1‑101 (MQ=21)
                                                          gagaCGAAACGGGG.GGCGCTGCCCCC.CTC.TT..GGCCTTC.GGCC.AATTCACCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGcc                                       >  1:1222553/1‑101 (MQ=21)
                                                          gagaCGAAACAGGG.GGCGCTACCCCC.CTC.TT..GGCCTTC.GGCC.AATTCCCCCCCCGGAATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGcc                                       <  1:3307486/101‑1 (MQ=21)
                                                           agaCGAAACGGGG.GGCGCTGCCCCC.CTC.TT..GGCCTTC.GGCC.AATTCACCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGccc                                      <  2:236725/101‑1 (MQ=21)
                                                              cGAAACGGGG.GGCGCTGCCCCC.CTC.TT..GGCCTTC.GGCC.AATTCACCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGaa                                   >  2:3077751/1‑101 (MQ=21)
                                                               gAAACGGGG.GGCGCTGCTCCC.CTC.TT..GGCCTTC.GGCC.AATTCACCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGaaa                                  <  2:3806296/101‑1 (MQ=21)
                                                               gAAACGGGG.GGCGCTGCCCCC.CTC.TT..GGCCTTC.GGCC.AATTCACCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGaaa                                  <  2:2702264/101‑1 (MQ=21)
                                                                aaacgggg.ggCGCTGCCCCC.CTC.TT..GGCCTTC.GGCC.AATTCACCCCCCGGGATATTTCGGGCAAGATAAAGAAGGGCTTTTCCTTTCATCTTGcccccaat                                 >  2:4223765/6‑96 (MQ=18)
                                                                   cgggg.ggCGCTGCCCCC.CTC.TT..GGCCTTC.GGCC.AATTCACCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATc                              <  2:4015211/99‑1 (MQ=25)
                                                                   caggg.ggCGCTG.CCCC.CTC.TT..GGCCTTT.GGCC.AATTCA.CCCCTGGGATATTTCGGGCAAGAT...GAA.ATCATGCCCCTTCATCTTGCCGCAAATATCCCgggg                        >  1:3951563/3‑101 (MQ=0)
                                                                    gggg.ggCGCTGCCCCC.CTC.TT..GGCCTTC.GGCC.AATTCACCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATcc                             <  2:1571825/100‑1 (MQ=25)
                                                                    gggg.ggCGCTGCCCCC.CTC.TT..GGCCTTC.GGCC.AATTCACCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATcc                             <  1:949433/100‑1 (MQ=25)
                                                                    aggg.ggCGCTG.CCCC.CTC.TT..GGCCTTT.GGCC.AATTCA.CCCCTGGGATATTTCGGGCAAGAT...GAA.ATCATGCCCCTTCATCTTGCCGCAAATATCCCggggg                       >  1:4076720/2‑101 (MQ=0)
                                                                      gg.ggCGCTGCCCCCTCTT.GG..C..CTTT.GGCC.AATTCA.CCCCTGGGATATTTCGGGCAAGAT...GAA.ATCATGCCCCTTCATCTTGCCGCAAATATCCCgggggga                     <  1:3959563/101‑3 (MQ=2)
                                                                       g.ggCGCTGCCCCCTCTT.GG..C..CTTT.GGCC.AATTCA.CCCCTGGGATATTTCGGGCAAGAT...GAA.ATCATGCCCCTTCATCTTGCCGCAAATATCCCggggggag                    >  1:601326/1‑98 (MQ=2)
                                                                       g.ggCGCTG.CCCC.CTC.TT..GGCCTTT.GGCC.AATTCA.CCCCTGGGATATTTCGGGCAAGAT...GAA.ATCATGCCCCTTCATCTTGCCGCAAATATCCCggggggag                    >  1:1254032/1‑98 (MQ=2)
                                                                         ggCGCTGCCCCCTCTT.GG..C..CTTT.GGCC.AATTCA.CCCCTGGGATATTTCGGGCAAGAT...GAA.ATCATGCCCCTTCATCTTGCCGCAAATATCCCggggggagc                   <  1:1048067/101‑5 (MQ=2)
                                                                         ggCGCTGCCCCCTCTT.GG..C..CTTT.GGCC.AATTCA.CCCCTGGGATATTTCGGGCAAGAT...GAA.ATCATGCCCCTTCATCTTGCCGCAAATATCCCggggggagc                   <  1:2334609/101‑5 (MQ=2)
                                                                            gcTGCCCCC.CTC.TT..GGCCTTC.GGCC.AATTCTCCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCCGGGGGa                      <  1:2866857/101‑1 (MQ=37)
                                                                             ccgCCCCC.CTC.TT..GGCCTTC.GGCC.AAGTCCCCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCCGGGGGAg                     <  2:3568447/99‑1 (MQ=37)
                                                                             cTGCCCCC.CTC.TT..GGCCTTC.GGCC.AATTCACCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCCGGGGGAg                     <  2:1142532/101‑1 (MQ=37)
                                                                                       ccctcttggccTTC.GGCC.AATTCCCCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCCGGGGGAGCCCGCt               <  1:1052427/91‑1 (MQ=38)
                                                                                       ccctcttggccTTC.GGCC.AATTCACCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCCGGGGGAGCCCGCt               <  1:93659/91‑1 (MQ=255)
                                                                                         ctcttggccTTC.GGCC.AATTCCCCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCCGGGGGAGCCCGCTTg             <  2:2117738/93‑1 (MQ=255)
                                                                                          tcttggccTTC.GGCC.AATTCACCCCCCGGGATATTTCGGGCAAGATAAAGAAGGGCTTTTCCTTTCATCTTGCCCGAAATATCCCGGGGGAGCCCGCTTgc            >  1:3011369/8‑101 (MQ=255)
                                                                                           cttggccTTC.GGCC.AATTCACCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCCGGGGGAGCCCGCTTgcg           >  2:3746717/7‑101 (MQ=255)
                                                                                            ttggccTTC.GGCC.AATTCACCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCCGGGGGAGCCCGCTTgcgg          <  1:3176964/96‑1 (MQ=255)
                                                                                            ttggccTTC.GGCC.AATTCACCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCCGGGGGAGCCCGCTTgcgg          <  1:3070202/96‑1 (MQ=255)
                                                                                            ttggccTTC.GGCC.AATTCACCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCCGGGGGAGCCCGCTTgcgg          <  2:2553261/96‑1 (MQ=255)
                                                                                             tggccTTC.GGCC.AATTCACCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCCGGGGGAGACCGCTTgccgg         >  2:1792567/5‑101 (MQ=255)
                                                                                              ggccTTC.GGCC.AATTCACCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCCGGGGGAGCCCGCTTgcgggc        >  1:721715/4‑101 (MQ=255)
                                                                                              ggccTTC.GGCC.AATTCACCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCCGGGGGAGCCCGCTTgcgggc        >  1:1210820/4‑101 (MQ=255)
                                                                                                ccTTC.GGCC.AATTCACCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCCGGGGGAGCCCGCTTgcgggcag      >  2:1717474/2‑99 (MQ=255)
                                                                                                    |                                                                                                      
GATGCGATCGACTGGCCCGAGGGGTTCTGCCGCCGCGACATCGGCTGGCGCGCGCTCTGAGACGAAAC.GGG.GGCGCTGCCCCC.CTC.CT..G..CTTC.GGCC.AATTCACCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCCGGGGGAGCCCGCTTGCGGGCGGGGGC  >  NC_014034/2566469‑2566661

Base quality scores:  ATCG  < 3 ≤  ATCG  < 17 ≤  ATCG  < 25 ≤  ATCG  < 33 ≤  ATCG  < 40 ≤  ATCG