Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_014034 2,566,597 Δ1 bp 5.3% intergenic (+69/+125) 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,5970A→. 5.3% 14.2 ‑0.4 ‑0.0133intergenic (+69/+125)purD/sdaAphosphoribosylamine‑‑glycine ligase/L‑serine ammonia‑lyase
Reads supporting (aligned to +/- strand):  new base (4/3):  ref base (56/70):  total (60/73)
Fisher's exact test strand distribution p-value = 7.01e-01
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 9.79e-01

Only 100 of 151 total aligned reads displayed.

GCCGCGACATCGGCTGGCGCGCGCTCTGAGACGAAAC.GGGGGCGCTGCCCCC.CTCCT.G..CTTC.GGCCAATTCA.CCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCCGGGGGAGCCCGCTTGCGGGCGGGGGCAGCGCCCCCGAAGCGCCCGACGCGCGGCGGGCTTGC  >  NC_014034/2566500‑2566697
                                                                                                        |                                                                                                    
gccgcGACATCGGCTGGCGCGCGCTCTGAGACGAAACGGGGGGCGCTGCCCCC.CTCTT.GGCCTTC.GGCCAATTCA.CCCCCCGGGATATTTCGGGCAAGATa                                                                                                      <  2:906698/101‑1 (MQ=21)
gccgcGACATCGGCTGGCGCGCGCTCTGAGACGAAACGGGGGGCGCTGCCCCC.CTCTT.GGCCTTC.GGCCAATTCA.CCCCCCGGGATATTTCGGGCAAGATa                                                                                                      <  2:3108432/101‑1 (MQ=21)
gccgcGACATCGGCTGGCGCGCGCTCTGAGACGAAACGGGGGGCGCTGCCCCC.CTCTT.GGCCTTC.GGCCAATTCA.CCCCCCGGGATATTTCGGGCAAGATa                                                                                                      <  2:2066625/101‑1 (MQ=21)
    cgACATCGGCTGGCGCGCGCTCTGAGACGAAACGGGGGGCGCTGCCCCC.CTCTT.GGCCTTC.GGCCAATTCA.CCCCCCGGGATATTTCGGGCAAGATAaaga                                                                                                  >  2:3176725/1‑101 (MQ=21)
    cgACATCGGCTGGCGCGCGCTCTGAGACGAAACGGGGGGCGCTGCCCCC.CTCTT.GGCCTTC.GGCCAATTCA.CCCCCCGGGATATTTCGGGCAAGATAaaga                                                                                                  <  1:1811536/101‑1 (MQ=21)
     gACATCGGCTGGCGCGCGCTCTGAGACGAAACGGGGGGCGCTGCCCCC.CTCTT.GGCCTTC.GGCCAATTCA.CCCCCCGGGATATTTCGGGCAAGATAaagaa                                                                                                 <  1:19493/101‑1 (MQ=21)
      aCATCGGCTGGCGCTCGCTCTGAGACGAAACGGGGGGCCCCGCCCCC.CTTTTGG..CCTCGGGCCAATTCA.CCCCCCGCGCCGTTTCGGGGGCGATAaaaaag                                                                                                >  1:4115095/1‑101 (MQ=9)
        aTCGGCTGGCGCGCGCTCTGAGACGAAACGGGGGGCGCTGCCCCC.CTCTT.GGCCTTC.GGCCCATTCA.CCCCCCCGGAAATTTTGGGGAAGATAAAGAAggg                                                                                              >  2:3039499/1‑100 (MQ=18)
           ggCTGGGGCGCGCTCTGAGACGAAACGGGGGGCGCTGCCCCC.CTCTT.GGCCTTC.GGCCAATTCA.CCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCtt                                                                                           <  1:3073594/101‑1 (MQ=21)
             cTGGCGCGCGCTCTGAGACGAAACGGGGGGCGCTGCCCCC.CTCTT.GGCCTTC.GGCCAATTCA.CCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCtttt                                                                                         <  1:2169706/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)
                gcgcgcgcTCTGAGACGAAACGGGGGGCGCTGCCCCC.CTCTT.GGCCTTC.GGCCAATTCA.CCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCt                                                                                      <  1:3376059/101‑1 (MQ=21)
                 cgcgcgcTCTGAGACGAAACGGGGGGCGCTGCCCCC.CTCTT.GGCCTTC.GGCCAATTCA.CCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCtt                                                                                     <  1:1533020/101‑1 (MQ=21)
                  gcgcgcTCTGAGACGAAACGGGGGGCGCTGCCCCC.CTCTT.GGCCTTC.GGCCAATTCC.CCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCttt                                                                                    <  2:823803/101‑1 (MQ=21)
                  gcgcgcTCTGAGACGAAACGGGGGGCGCTGCCCCC.CTCTT.GGCCTTC.GGCCAATTCA.CCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCttt                                                                                    <  2:2428260/101‑1 (MQ=21)
                  gcgcgcTCTGAGACGAAACGGGGGGCGCTGCCCCC.CTCTT.GGCCTTC.GGCCAATTCA.CCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCttt                                                                                    <  2:1830164/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:902741/101‑1 (MQ=21)
                      gcTCTGAGACGAAACGGGGGGCGCTGCCCCC.CTCTT.GGCCTTC.GGCCAATTCA.CCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATc                                                                                <  1:22833/101‑1 (MQ=21)
                      gcTCTGAGACGAAACGGGGGGCGCTGCCCCC.CTCTT.GGCCTTC.GGCCAATTCA.CCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATc                                                                                <  1:31132/101‑1 (MQ=21)
                       ctctGAGACGAAACGGGGGGCGCTGCCCCC.CTCTT.GGCCTTC.GGCCAATTCA.CCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCt                                                                               <  1:2775027/101‑1 (MQ=21)
                        tctGAGACGAAACGGGGGGCGCTGCCCCC.CTCTT.GGCCTTC.GGCCAATTCA.CCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCtt                                                                              <  1:1173596/101‑1 (MQ=21)
                        tctGAGACGAAACGGGGGGCGCTGCCCCC.CTCTT.GGCCTTC.GGCCAATTCA.CCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCtt                                                                              <  1:298907/101‑1 (MQ=21)
                          tGAGACGAAACGGGGGGCGCTGCCCCC.CTCTT.GGCCTTC.GGCCAATTCA.CCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGc                                                                            >  2:2272230/1‑101 (MQ=21)
                           gagaCGAAACGGGGGGCGCTGCCCCC.CTCTT.GGCCTTC.GGCCAATTCA.CCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGcc                                                                           >  1:1222553/1‑101 (MQ=21)
                           gagaCGAAACAGGGGGCGCTACCCCC.CTCTT.GGCCTTC.GGCCAATTCC.CCCCCCGGAATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGcc                                                                           <  1:3307486/101‑1 (MQ=21)
                              aCGAAACGGGGGGCGCTGCCCCC.CTCTT.GGCCTTC.GGCCAATTCA.CCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGa                                                                        >  1:1787489/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)
                                  aacggggggCGCTGCCCCC.CTCTT.GGCCTTC.GGCCAATTCACCCCCCCGGGAAATTT.GGGCCAGAAAAAGAAAGGCTTTTTCTTTCATCTTTcccccaaaa                                                                    >  1:2185736/5‑95 (MQ=2)
                                  aacggggggCGCTGCCCCC.CTCTT.GGCCTTC.GGCCAATTCA.CCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAata                                                                    <  1:3581109/97‑1 (MQ=25)
                                   acggggggCGCTGCCCCC.CTCTT.GGCCTTC.GGCCAATTCA.CCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAatat                                                                   <  2:4156408/98‑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)
                                    cagggggCGCTG.CCCC.CTCTT.GGCCTTT.GGCCAATTCA..CCCCTGGGATATTTCGGGCAAGAT...GAA.ATCATGCCCCTTCATCTTGCCGCAAATATCCCgggg                                                            >  1:3951563/3‑101 (MQ=0)
                                     ggggggCGCTGCCCCC.CTCTT.GGCCTTC.GGCCAATTCA.CCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATcc                                                                 <  2:1571825/100‑1 (MQ=25)
                                     ggggggCGCTGCCCCC.CTCTT.GGCCTTC.GGCCAATTCA.CCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATcc                                                                 <  1:949433/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)
                                         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)
                                             cTGCCCCC.CTCTT.GGCCTTC.GGCCAATTCA.CCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCCGGGGGAg                                                         <  2:1142532/101‑1 (MQ=37)
                                                     ccctcttggccTTC.GGCCAATTCC.CCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCCGGGGGAGCCCGCt                                                   <  1:1052427/91‑1 (MQ=38)
                                                     ccctcttggccTTC.GGCCAATTCA.CCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCCGGGGGAGCCCGCt                                                   <  1:93659/91‑1 (MQ=255)
                                                       ctcttggccTTC.GGCCAATTCC.CCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCCGGGGGAGCCCGCTTg                                                 <  2:2117738/93‑1 (MQ=255)
                                                        tcttggccTTC.GGCCAATTCA.CCCCCCGGGATATTTCGGGCAAGATAAAGAAGGGCTTTTCCTTTCATCTTGCCCGAAATATCCCGGGGGAGCCCGCTTgc                                                >  1:3011369/8‑101 (MQ=255)
                                                         cttggccTTC.GGCCAATTCA.CCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCCGGGGGAGCCCGCTTgcg                                               >  2:3746717/7‑101 (MQ=255)
                                                          ttggccTTC.GGCCAATTCA.CCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCCGGGGGAGCCCGCTTgcgg                                              <  1:3176964/96‑1 (MQ=255)
                                                          ttggccTTC.GGCCAATTCA.CCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCCGGGGGAGCCCGCTTgcgg                                              <  2:2553261/96‑1 (MQ=255)
                                                           tggccTTC.GGCCAATTCA.CCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCCGGGGGAGACCGCTTgccgg                                             >  2:1792567/5‑101 (MQ=255)
                                                            ggccTTC.GGCCAATTCA.CCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCCGGGGGAGCCCGCTTgcgggc                                            >  1:721715/4‑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)
                                                                     gCCAATTCA.CCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCCGGGGGAGCCCGCTTGCGGGCGGGGGCAg                                    <  2:1958095/101‑1 (MQ=255)
                                                                       cAAGTCA.CCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCCGGGGGAGCCCGCTTGCGGGCGGGGGCAgcg                                  <  1:1876199/101‑1 (MQ=255)
                                                                        aaTTCA.CCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCCGGGGGAGCCCGCTTGCGGGCGGGGGCAgcgc                                 >  2:1742483/1‑101 (MQ=255)
                                                                         aTTCA.CCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCCGGGGGAGCCCGCTTGCGGGCGGGGGCAGCGcc                                >  2:1662111/1‑101 (MQ=255)
                                                                         aTTCA.CCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCCGGGGGAGCCCGCTTGCGGGCGGGGGCAGCGcc                                >  1:2704649/1‑101 (MQ=255)
                                                                         aTTCA.CCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCCGGGGGAGCCCGCTTGCGGGCGGGGGCAGCGcc                                >  2:722881/1‑101 (MQ=255)
                                                                           tCA.CCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCCAGGGGAGCCCGCTTGCGGGGGGGGGGAGCGcccc                              >  2:1633744/1‑101 (MQ=255)
                                                                            cA.CCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCCGGGGGAGCCCGCTTGCGGGCGGGGGCAGCGccccc                             <  1:3077981/101‑1 (MQ=255)
                                                                             a.CCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCCGGGGGAGCCCGCTTGCGGGCGGGGGCAGCGCCCCCg                            >  1:584338/1‑101 (MQ=255)
                                                                               ccccccGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCCGGGGGAGCCCGCTTGCGGGCGGGGGCAGCGCCCCCGa                           <  1:3964337/101‑1 (MQ=255)
                                                                                cccccGAGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCCGGGGGAGCCCGCTTGCGGGCGGGGGCAGCGCCCCCGaa                          <  2:2040485/101‑1 (MQ=255)
                                                                                    cGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCCGGGGGAGCCCGCTTGCGGGCGGGGGCAGCGCCCCCGAAgcgc                      >  1:1676698/1‑101 (MQ=255)
                                                                                     gggATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCCGGGGGAGCCCGCTTGCGGGCGGGGGCAGCGCCCCCGAAGCGcc                     <  1:3779274/101‑1 (MQ=255)
                                                                                      ggATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCCGGGGGAGCCCGCTTGCGGGCGGGGGCCGCGCCCCCGAAGCGccc                    >  2:2508330/1‑101 (MQ=255)
                                                                                      ggATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCCGGGGGAGCCCGCTTGCGGGCGGGGGCAGCGCCCCCGAAGCGccc                    >  1:3557518/1‑101 (MQ=255)
                                                                                      ggATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCCGGGGAAGCCCGCTTGCGGGCGGGGGCAGCGCCCCCGAAGCGccc                    <  2:680417/101‑1 (MQ=255)
                                                                                       gATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCCGGGGGAGCCCGCTTGCGGGCGGGGGCAGCGCCCCCGAAGCGCCCg                   <  1:2636550/101‑1 (MQ=255)
                                                                                        atatTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCCGGGGGAGCCCGCTTGCGGGCGGGGGCAGCGCCCCCGAAGCGCCCGa                  >  2:2199887/1‑101 (MQ=255)
                                                                                        atatTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCCGGGGGAGCCCGCTTGCGGGCGGGGGCAGCGCCCCCGAAGCGCCCGa                  >  1:3851558/1‑101 (MQ=255)
                                                                                          atTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCCGGGGGAGCCCGCTTGCGGGCGGGGGCAGCGCCCCCGAAGCGCCCGAcg                >  1:2516173/1‑101 (MQ=255)
                                                                                          atTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCCGGGGGAGCCCGCTTGCGGGCGGGGGCAGCGCCCCCGAAGCGCCCGAcg                >  1:323970/1‑101 (MQ=255)
                                                                                          atTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCCGGGGGAGCCCGCTTGCGGGCGGGGGCAGCGCCCCCGAAGCGCCCGAcg                >  1:1927029/1‑101 (MQ=255)
                                                                                           tttCGGACAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCCGGGGGAGCCCGCTTGCGGGCGGGGGCAGCGCCCCCGAAGCGCCCGAcgc               <  2:436788/101‑1 (MQ=255)
                                                                                              cGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCCGGGGGAGCCCGCTTGCGGGCGGGGGCAGCGCCCCCGAAGCGCCCGAcgcgcg            >  1:3450870/1‑101 (MQ=255)
                                                                                               gggCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCCGGGGGAGCCCGCTTGCGGGCGGGGGCAGCGCCCCCGAAGCGCCCGACGCgcgg           >  2:3547478/1‑101 (MQ=255)
                                                                                               gggCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCCGGGGGAGCCCGCTTGCGGGCGGGGGCAGCGCCCCCGAAGCGCCCGACGCgcgg           >  2:848799/1‑101 (MQ=255)
                                                                                                ggCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCCGGGGGAGCCCGCTTGCGGGCGGGGGCAGCGCCCCCGAAGCGCCCGACGCgcggc          >  2:185544/1‑101 (MQ=255)
                                                                                                   aaGATAAAGAAGGTCTTTTCTTTTCATCTTGCCCGAAATATCCCGGGGGAGCCCGCTTGCGGGCGGGGGCAGCGCCCCCGAAGCGCCCGACGCGCGGCggg       <  1:4021231/101‑1 (MQ=255)
                                                                                                   aaGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCCGGGGGAGCCCGCTTGCGGGCGGGGGCAGCGCCCCCGAAGCGCCCGACGCGCGGCggg       >  1:117707/1‑101 (MQ=255)
                                                                                                   aaGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCCGGGGGAGCCCGCTTGCGGGCGGGGGCAGCGCCCCCGAAGCGCCCGACGCGCGGCggg       >  2:2700590/1‑101 (MQ=255)
                                                                                                    aGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCCGGGGGAGCCCGCTTGCGGGCGGGGGCAGCGCCCCCGAAGCGCCCGACGCGCGGCGGGc      >  1:203708/1‑101 (MQ=255)
                                                                                                      aTAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCCGGGGGAGCCCGCTTGCGGGCGGGGGCAGCGCCCCCGAAGCGCCCGACGCGCGGCGGGCtt    <  2:1457617/101‑1 (MQ=255)
                                                                                                       tAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCCGGGGGAGCCCGCTTGCGGGCGGGGGCAGCGCCCCCGAAGCGCCCGACGCGCGGCGGGCTTg   >  1:247466/1‑101 (MQ=255)
                                                                                                       tAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCCGGGGGAGCCCGCTTGCGGGCGGGGGCAGCGCCCCCGAAGCGCCCGACGCGCGGCGGGCTTg   >  2:527904/1‑101 (MQ=255)
                                                                                                        aaaGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCCGGGGGAGCCCGCTTGCGGGCGGGGGCAGCGCCCCCGAAGCGCCCGACGCGCGGCGGGCTTgc  >  1:794897/1‑101 (MQ=255)
                                                                                                        aaaGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCCGGGGGAGCCCGCTTGCGGGCGGGGGCAGCGCCCCCGAAGCGCCCGACGCGCGGCGGGCTTgc  >  1:1983891/1‑101 (MQ=255)
                                                                                                        aaaGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCCGGGGGAGCCCGCTTGCGGGCGGGGGCAGCGCCCCCGAAGCGCCCGACGCGCGGCGGGCTTgc  >  1:2248743/1‑101 (MQ=255)
                                                                                                        aaaGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCCGGGGGAGCCCGCTTGCGGGCGGGGGCAGCGCCCCCGAAGCGCCCGACGCGCGGCGGGCTTgc  >  2:3147613/1‑101 (MQ=255)
                                                                                                        |                                                                                                    
GCCGCGACATCGGCTGGCGCGCGCTCTGAGACGAAAC.GGGGGCGCTGCCCCC.CTCCT.G..CTTC.GGCCAATTCA.CCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCCGGGGGAGCCCGCTTGCGGGCGGGGGCAGCGCCCCCGAAGCGCCCGACGCGCGGCGGGCTTGC  >  NC_014034/2566500‑2566697

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