Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_014034 2,566,603 Δ1 bp 5.2% intergenic (+75/+119) 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,6030G→. 5.2% 14.6 ‑0.3 ‑0.1135intergenic (+75/+119)purD/sdaAphosphoribosylamine‑‑glycine ligase/L‑serine ammonia‑lyase
Reads supporting (aligned to +/- strand):  new base (4/3):  ref base (61/67):  total (65/70)
Fisher's exact test strand distribution p-value = 7.11e-01
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 7.97e-01

Only 100 of 158 total aligned reads displayed.

ACATCGGCTGGCGCGCGCTCTGAGACGAAAC.GGGGGCGCTGCCCCC.CTCCT.G..CTTC.GGCCAATTCA.CCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCCGGGGGAGCCCGCTTGCGGGCGGGGGCAGCGCCCCCGAAGCGCCCGACGCGCGGCGGGCTTGCCCGCCC  >  NC_014034/2566506‑2566703
                                                                                                        |                                                                                                    
aCATCGGCTGGCGCTCGCTCTGAGACGAAACGGGGGGCCCCGCCCCC.CTTTTGG..CCTCGGGCCAATTCA.CCCCCCGCGCCGTTTCGGGGGCGATAaaaaag                                                                                                      >  1:4115095/1‑101 (MQ=9)
  aTCGGCTGGCGCGCGCTCTGAGACGAAACGGGGGGCGCTGCCCCC.CTCTT.GGCCTTC.GGCCCATTCA.CCCCCCCGGAAATTTTGGGGAAGATAAAGAAggg                                                                                                    >  2:3039499/1‑100 (MQ=18)
   tCGGCTGGCGCGCGCTCTGAGACGAAACGGGGGGCGCTGCCCCC.CTCTT.GGCCTTC.GGCCAATTCA.CCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTc                                                                                                   <  1:3471783/101‑1 (MQ=21)
   tCGGCTGGCGCGCGCTCTGAGACGAAACGGGGGGCGCTGCCCCC.CTCTT.GGCCTTC.GGCCAATTCA.CCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTc                                                                                                   >  1:3848117/1‑101 (MQ=21)
      gCTGGCGCGCGCTCTGAGACGAAACGGGGGGCGCTGCCCCC.CTCTT.GGCCTTC.GGCCAATTCA.CCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCttt                                                                                                <  2:1078597/101‑1 (MQ=21)
       cTGGCGCGCGCTCTGAGACGAAACGGGGGGCGCTGCCCCC.CTCTT.GGCCTTC.GGCCAATTCA.CCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCtttt                                                                                               <  1:2169706/101‑1 (MQ=21)
        tGGCGCGCGCTCTGAGACGAAACGGGGGGCGCTGCCCCC.CTCTT.GGCCTTC.GGCCAATTCA.CCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTc                                                                                              <  1:3248456/101‑1 (MQ=21)
          gcgcgcgcTCTGAGACGAAACGGGGGGCGCTGCCCCC.CTCTT.GGCCTTC.GGCCAATTCA.CCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCt                                                                                            <  1:3376059/101‑1 (MQ=21)
           cgcgcgcTCTGAGACGAAACGGGGGGCGCTGCCCCC.CTCTT.GGCCTTC.GGCCAATTCA.CCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCtt                                                                                           <  1:2512408/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)
             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)
                  tctGAGACGAAACGGGGGGCGCTGCCCCC.CTCTT.GGCCTTC.GGCCAATTCA.CCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCtt                                                                                    <  1:298907/101‑1 (MQ=21)
                  tctGAGACGAAACGGGGGGCGCTGCCCCC.CTCTT.GGCCTTC.GGCCAATTCA.CCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCtt                                                                                    <  1:1173596/101‑1 (MQ=21)
                    tGAGACGAAACGGGGGGCGCTGCCCCC.CTCTT.GGCCTTC.GGCCAATTCA.CCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTGTTTTTTTTTCTTCTTGc                                                                                  >  2:853783/1‑101 (MQ=18)
                    tGAGACGAAACGGGGGGCGCTGCCCCC.CTCTT.GGCCTTC.GGCCAATTCA.CCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGc                                                                                  >  2:2272230/1‑101 (MQ=21)
                     gagaCGAAACGGGGGGCGCTGCCCCC.CTCTT.GGCCTTC.GGCCAATTCA.CCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGcc                                                                                 >  1:1222553/1‑101 (MQ=21)
                     gagaCGAAACGGGGGGCGCTGCCCCC.CTCTT.GGCCTTC.GGCCAATTCA.CCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGcc                                                                                 >  2:1826958/1‑101 (MQ=21)
                      agaCGAAACGGGGGGCGCTGCCCCC.CTCTT.GGCCTTC.GGCCAATTCA.CCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGccc                                                                                <  2:236725/101‑1 (MQ=21)
                        aCGAAACGGGGGGCGCTGCCCCC.CTCTT.GGCCTTC.GGCCAATTCA.CCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGa                                                                              >  1:1787489/1‑101 (MQ=21)
                         cGAAACGGGGGGCGCTGCCCCC.CTCTT.GGCCTTC.GGCCAATTCA.CCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGaa                                                                             <  1:3404108/101‑1 (MQ=21)
                          gAAACGGGGGGCGCTGCTCCC.CTCTT.GGCCTTC.GGCCAATTCA.CCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGaaa                                                                            <  2:3806296/101‑1 (MQ=21)
                           aaacggggggCGCTGCCCCC.CTCTT.GGCCTTC.GGCCAATTCA.CCCCCCGGGATATTTCGGGCAAGATAAAGAAGGGCTTTTCCTTTCATCTTGcccccaat                                                                           >  2:4223765/6‑96 (MQ=18)
                            aacggggggCGCTGCCCCC.CTCTT.GGCCTTC.GGCCAATTCACCCCCCCGGGAAATTT.GGGCCAGAAAAAGAAAGGCTTTTTCTTTCATCTTTcccccaaaa                                                                          >  1:2185736/5‑95 (MQ=2)
                             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)
                               ggggggCGCTGCCCCC.CTCTT.GGCCTTC.GGCCAATTCA.CCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATcc                                                                       <  1:949433/100‑1 (MQ=25)
                               ggggggCGCTGCCCCC.CTCTT.GGCCTTC.GGCCAATTCA.CCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATcc                                                                       <  2:1571825/100‑1 (MQ=25)
                                 ggggCGCTGCCCCCTCTTGG.C..CTTT.GGCCAATTCA..CCCCTGGGATATTTCGGGCAAGAT...GAA.ATCATGCCCCTTCATCTTGCCGCAAATATCCCgggggga                                                               <  1:3959563/101‑3 (MQ=2)
                                  gggCGCTGCCCCCTCTTGG.C..CTTT.GGCCAATTCA..CCCCTGGGATATTTCGGGCAAGAT...GAA.ATCATGCCCCTTCATCTTGCCGCAAATATCCCggggggag                                                              >  1:601326/1‑98 (MQ=2)
                                   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)
                                      gcTGCCCCC.CTCTT.GGCCTTC.GGCCAATTCT.CCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCCGGGGGa                                                                <  1:2866857/101‑1 (MQ=37)
                                       ccgCCCCC.CTCTT.GGCCTTC.GGCCAAGTCC.CCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCCGGGGGAg                                                               <  2:3568447/99‑1 (MQ=37)
                                               ccctcttggccTTC.GGCCAATTCC.CCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCCGGGGGAGCCCGCt                                                         <  1:1052427/91‑1 (MQ=38)
                                               ccctcttggccTTC.GGCCAATTCA.CCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCCGGGGGAGCCCGCt                                                         <  1:93659/91‑1 (MQ=255)
                                                 ctcttggccTTC.GGCCAATTCC.CCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCCGGGGGAGCCCGCTTg                                                       <  2:2117738/93‑1 (MQ=255)
                                                   cttggccTTC.GGCCAATTCA.CCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCCGGGGGAGCCCGCTTgcg                                                     >  2:3746717/7‑101 (MQ=255)
                                                    ttggccTTC.GGCCAATTCA.CCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCCGGGGGAGCCCGCTTgcgg                                                    <  1:3070202/96‑1 (MQ=255)
                                                    ttggccTTC.GGCCAATTCA.CCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCCGGGGGAGCCCGCTTgcgg                                                    <  2:2553261/96‑1 (MQ=255)
                                                    ttggccTTC.GGCCAATTCA.CCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCCGGGGGAGCCCGCTTgcgg                                                    <  1:3176964/96‑1 (MQ=255)
                                                     tggccTTC.GGCCAATTCA.CCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCCGGGGGAGACCGCTTgccgg                                                   >  2:1792567/5‑101 (MQ=255)
                                                      ggccTTC.GGCCAATTCA.CCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCCGGGGGAGCCCGCTTgcgggc                                                  >  1:721715/4‑101 (MQ=255)
                                                      ggccTTC.GGCCAATTCA.CCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCCGGGGGAGCCCGCTTgcgggc                                                  >  1:1210820/4‑101 (MQ=255)
                                                      ggccTTC.GGCCAATTCA.CCCCCCGGGATATTTCGGGCAAGATAAAGAAGGGCTTTTCCTTTCATCTTGCCCGAAAGATCCCGGGGGACCCCGCCTgggggc                                                  >  1:3885267/4‑101 (MQ=25)
                                                        ccTTC.GGCCAATTCA.CCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCCGGGGGAGCCCGCTTgcgggcag                                                >  2:1717474/2‑99 (MQ=255)
                                                            c.GGCCAATTCA.CCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCCGGGGGAGCCCGCTTGCGGGCGGGGGc                                            <  2:2144179/101‑1 (MQ=255)
                                                               gCCAATTCA.CCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCCGGGGGAGCCCGCTTGCGGGCGGGGGCAg                                          <  2:1958095/101‑1 (MQ=255)
                                                                 cAAGTCA.CCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCCGGGGGAGCCCGCTTGCGGGCGGGGGCAgcg                                        <  1:1876199/101‑1 (MQ=255)
                                                                   aTTCA.CCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCCGGGGGAGCCCGCTTGCGGGCGGGGGCAGCGcc                                      >  2:722881/1‑101 (MQ=255)
                                                                   aTTCA.CCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCCGGGGGAGCCCGCTTGCGGGCGGGGGCAGCGcc                                      >  1:2704649/1‑101 (MQ=255)
                                                                   aTTCA.CCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCCGGGGGAGCCCGCTTGCGGGCGGGGGCAGCGcc                                      >  1:3091680/1‑101 (MQ=255)
                                                                     tCA.CCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCCAGGGGAGCCCGCTTGCGGGGGGGGGGAGCGcccc                                    >  2:1633744/1‑101 (MQ=255)
                                                                      cA.CCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCCGGGGGAGCCCGCTTGCGGGCGGGGGCAGCGccccc                                   <  1:2029001/101‑1 (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)
                                                                            cccGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCCGGGGGAGCCCGCTTGCGGGCGGGGGCAGCGCCCCCGAAgc                              <  1:1742602/101‑1 (MQ=255)
                                                                             ccGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCCGGGGGAGCCCGCTTGCGGGCGGGGGCAGCGCCCCCGAAgcg                             >  2:1190061/1‑101 (MQ=255)
                                                                               gggATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCCGGGGGAGCCCGCTTGCGGGCGGGGGCAGCGCCCCCGAAGCGcc                           >  2:301716/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)
                                                                                    atTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCCGGGGGAGCCCGCTTGCGGGCGGGGGCAGCGCCCCCGAAGCGCCCGAcg                      >  1:2516173/1‑101 (MQ=255)
                                                                                    atTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCCGGGGGAGCCCGCTTGCGGGCGGGGGCAGCGCCCCCGAAGCGCCCGAcg                      >  1:1927029/1‑101 (MQ=255)
                                                                                      ttCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCCGGGGGAGCCCGCTTGCGGGCGGGGGCAGCGCCCCCGAAGCGCCCGAcgcg                    <  2:3428342/101‑1 (MQ=255)
                                                                                      ttCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCCGGGGGAGCCCGCTTGCGGGCGGGGGCAGCGCCCCCGAAGCGCCCGAcgcg                    <  2:2617717/101‑1 (MQ=255)
                                                                                      ttCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCCGGGGGAGCCCGCTTGCGGGCGGGGGCAGCGCCCCCGAAGCGCCCGAcgcg                    <  1:2683626/101‑1 (MQ=255)
                                                                                          ggCAAGTTAAAGAAGGTTTTTTCCTTTCTTCTTGCCCGAAATATCCCGGGGGATCCCGCTTGCGGGCGGGGGCAGCGCCCCCGAAGCGCCCGACGCgcggc                <  1:3282303/101‑1 (MQ=39)
                                                                                          ggCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCCGGGGGAGCCCGCTTGCGGGCGGGGGCAGCGCCCCCGAAGCGCCCGACGCgcggc                >  2:185544/1‑101 (MQ=255)
                                                                                             aaGATAAAGAAGGTCTTTTCTTTTCATCTTGCCCGAAATATCCCGGGGGAGCCCGCTTGCGGGCGGGGGCAGCGCCCCCGAAGCGCCCGACGCGCGGCggg             <  1:4021231/101‑1 (MQ=255)
                                                                                             aaGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCCGGGGGAGCCCGCTTGCGGGCGGGGGCAGCGCCCCCGAAGCGCCCGACGCGCGGCggg             >  2:2700590/1‑101 (MQ=255)
                                                                                             aaGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCCGGGGGAGCCCGCTTGCGGGCGGGGGCAGCGCCCCCGAAGCGCCCGACGCGCGGCggg             >  1:117707/1‑101 (MQ=255)
                                                                                                 tAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCCGGGGGAGCCCGCTTGCGGGCGGGGGCAGCGCCCCCGAAGCGCCCGACGCGCGGCGGGGttc         >  1:3181312/1‑100 (MQ=255)
                                                                                                 tAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCCGGGGGAGCCCGCTTGCGGGCGGGGGCAGCGCCCCCGAAGCGCCCGACGCGCGGCGGGCTTg         >  2:527904/1‑101 (MQ=255)
                                                                                                  aaaGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCCGGGGGAGCCCGCTTGCGGGCGGGGGCAGCGCCCCCGAAGCGCCCGACGCGCGGCGGGCTTgc        >  1:2248743/1‑101 (MQ=255)
                                                                                                  aaaGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCCGGGGGAGCCCGCTTGCGGGCGGGGGCAGCGCCCCCGAAGCGCCCGACGCGCGGCGGGCTTgc        >  2:3148243/1‑101 (MQ=255)
                                                                                                  aaaGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCCGGGGGAGCCCGCTTGCGGGCGGGGGCAGCGCCCCCGAAGCGCCCGACGCGCGGCGGGCTTgc        >  2:3147613/1‑101 (MQ=255)
                                                                                                  aaaGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCCGGGGGAGCCCGCTTGCGGGCGGGGGCAGCGCCCCCGAAGCGCCCGACGCGCGGCGGGCTTgc        >  1:1983891/1‑101 (MQ=255)
                                                                                                  aaaGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCCGGGGGAGCCCGCTTGCGGGCGGGGGCAGCGCCCCCGAAGCGCCCGACGCGCGGCGGGCTTgc        >  1:794897/1‑101 (MQ=255)
                                                                                                  aaaGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCCGGGGGAGCCCGCTTGCGGGCGGGGGCAGAGCCCCCGAAGCGCCCGACGCGCGGCGGGGTTgc        >  1:3609110/1‑101 (MQ=255)
                                                                                                   aagaagGTCTTTTCCTTTCATCTTGCCCGAAATATCCCGGGGGAGCCCGCTTGCGGGCGGGGGCAGCGCCCCCGAAGCGCCCGACGCGCGGCGGGCTTgcc       >  1:3421357/1‑101 (MQ=255)
                                                                                                   aagaagGTCTTTTCCTTTCATCTTGCCCGAAATATCCCGGGGGAGCCCGCTTGCGGGCGGGGGCAGCGCCCCCGAAGCGCCCGACGCGCGGCGGGCTTccc       >  1:3903167/1‑101 (MQ=255)
                                                                                                    agaagGTCTTTTCCTTTCATCTTGCCCGAAATATCCCGGGGGAGCCCGCTTGCGGGCGGGGGCAGCGCCCCCGAAGCGCCCGACGCGCGGCGGGCTTgccc      >  2:3338628/1‑101 (MQ=255)
                                                                                                    agaagGTCTTTTCCTTTCATCTTGCCCGAAATATCCCGGGGGAGCCCGCTTGCGGGCGGGGGCAGCGCCCCCGAAGCGCCCGACGCGCGGCGGGCTTgccc      >  1:3845393/1‑101 (MQ=255)
                                                                                                     gaagGTCTTTTCCTTTCATCTTGCCCGAAATATCCCGGGGGAGCCCGCTTGCGGGCGGGGGCAGCGCCCCCGAAGCGCCCGACGCGCGGCGGGCTTgcccg     <  2:1222457/101‑1 (MQ=255)
                                                                                                     gaagGTCTTTTCCTTTCATCTTGCCCGAAATATCCCGGGGGAGCCCGCTTGCGGGCGGGGGCAGCGCCCCCGAAGCGCCCGACGCGCGGCGGGCTTgcccg     >  2:2955844/1‑101 (MQ=255)
                                                                                                      aagGTCTTTTCCTTTCATCTTGCCCGAAATATCCCGGGGGAGCCCGCTTGCGGGCGGGGGCAGCGCCCCCGAAGCGCCCGACGCGCGGCGGGCTTgcccgc    >  2:1578716/1‑101 (MQ=255)
                                                                                                      aagGTCTTTTCCTTTCATCTTGCCCGAAATATCCCGGGGGAGCCCGCTTGCGGGCGGGGGCAGCGCCCCCGAAGCGCCCGACGCGCGGCGGGCTTgcccgc    >  2:3960384/1‑101 (MQ=255)
                                                                                                      aagGTCTTTTCCTTTCATCTTGCCCGAAATATCCCGGGGGAGCCCGCTTGCGGGCGGGGGCAGCGCCCCCGAAGCGCCCGACGCGCGGCGGGCTTgcccgc    >  2:1710587/1‑101 (MQ=255)
                                                                                                      aagGTCTTTTCCTTTCATCTTGCCCGAAATATCCCGGGGGAGCCCGCTTGCGGGCGGGGGCAGCGCCCCCGAAGCGCCCGACGCGCGGCGGGCTTgcccgc    >  1:3795794/1‑101 (MQ=255)
                                                                                                       agGTCTTTTCCTTTCATCTTGCCCGAAATATCCCGGGGGAGCCCGCTTGCGGGCGGGGGCAGCGCCCCCGAAGCGCCCGACGCGCGGCGGGCTTgcccgcc   >  2:800901/1‑101 (MQ=255)
                                                                                                        |                                                                                                    
ACATCGGCTGGCGCGCGCTCTGAGACGAAAC.GGGGGCGCTGCCCCC.CTCCT.G..CTTC.GGCCAATTCA.CCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCCGGGGGAGCCCGCTTGCGGGCGGGGGCAGCGCCCCCGAAGCGCCCGACGCGCGGCGGGCTTGCCCGCCC  >  NC_014034/2566506‑2566703

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