Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_014034 2,566,625 C→G 2.7% intergenic (+97/+97) 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,6250C→G 2.7% 4.1 ‑1.6 ‑0.0193intergenic (+97/+97)purD/sdaAphosphoribosylamine‑‑glycine ligase/L‑serine ammonia‑lyase
Reads supporting (aligned to +/- strand):  new base (2/4):  ref base (121/66):  total (123/70)
Fisher's exact test strand distribution p-value = 1.92e-01
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 9.92e-01

Only 100 of 214 total aligned reads displayed.

AGACGAAAC.GGGGGCGCTGCCCCCCTCCTG..CTTCGGCCAATTCACCCCCCGGG.ATATTTCGGGCAAGATAAAGAAG.GTCTTTTCCTTTCATC.TTGCCCG.AAATATCCCGGGGGAGCCCGCTTGCGGGCGGGGGCAGCGCCCCCGAAGCGCCCGACGCGCGGCGGGCTTGCCCGCCCAAGGAGGGCGCAACCGCCTCAG  >  NC_014034/2566528‑2566725
                                                                                                       |                                                                                                     
agaCGAAACGGGGGGCGCTGCCCCCCTCTTGGCCTTCGGCCAATTCACCCCCCGGG.ATATTTCGGGCAAGATAAAGAAG.GTCTTTTCCTTTCATC.TTGccc                                                                                                       <  2:1699276/101‑1 (MQ=21)
agaCGAAACGGGGGGCGCTGCCCCCCTCTTGGCCTTCGGCCAATTCACCCCCCGGG.ATATTTCGGGCAAGATAAAGAAG.GTCTTTTCCTTTCATC.TTGccc                                                                                                       <  1:191171/101‑1 (MQ=21)
 gaCGAAACGGGGGGCGCTGCCCCCCTCTTGGCCTTCGGCCAATTCACCCCCCGGG.ATATTTCGGGCAAGATAAAGAAG.GTCTTTTCCTTTCATC.TTGCCCg                                                                                                      <  1:514012/101‑1 (MQ=21)
 gaCGAAACGGGGGGCGCTGCCCCCCTCTTGGCCTTCGGCCAATTCACCCCCCGGG.ATATTTCGGGCAAGATAAAGAAG.GTCTTTTCCTTTCATC.TTGCCCg                                                                                                      <  2:3407080/101‑1 (MQ=21)
 gaCGAAACGGGGGGCGCTGCCCCCCTCTTGGCCTTCGGCCAATTCACCCCCCGGG.ATATTTCGGGCAAGATAAAGAAG.GTCTTTTCCTTTCATC.TTGCCCg                                                                                                      <  1:1894315/101‑1 (MQ=21)
  aCGAAACGGGGGGCGCTGCCCCCCTCTTGGCCTTCGGCCAATTCACCCCCCGGG.ATATTTCGGGCAAGATAAAGAAG.GTCTTTTCCTTTCATC.TTGCCCG.a                                                                                                    <  2:539174/101‑1 (MQ=21)
  aCGAAACGGGGGGCGCTGCCCCCCTCTTGGCCTTCGGCCAATTCACCCCCCGGG.ATATTTCGGGCAAGATAAAGAAG.GTCTTTTCCTTTCATC.TTGCCCG.a                                                                                                    <  2:4727329/101‑1 (MQ=21)
   cGAAACGGGGGGCGCTGCCCCCCTCTTGGCCTTCGGCCAATTCACCCCCCGGG.ATATTTCGGGCAAGATAAAGAAG.GTCTTTTCCTTTCATC.TTGCCCG.aa                                                                                                   >  1:3357414/1‑101 (MQ=21)
   cGAAACGGGGGGCGCTGCCCCCCTCTTGGCCTTCGGCCAATTCACCCCCCGGG.ACATTTCGGGCAAGATAAAGAAG.GTCTTTTCCTTTCATC.TTGCCCG.aa                                                                                                   >  2:613597/1‑101 (MQ=21)
    gaagcggggggCGCAGCCCCCCTCTTGGCCATCGGCCAAGTCAGCCACAGGG.ATATATCGGGCAAGATAAAGAAG.GTCTTTTCCTTTCATC.TTGCCCG.aaa                                                                                                  <  2:4183595/95‑1 (MQ=255)
    gAAACGGGGGGCGCTGCCCCCCTCTTGGCCTTCGGCCAATTCACCCCCCGGG.ATATTTCGGGCAAGATAAAGAAG.GTCTTTTCCTTTCATC.TTGCCCG.aaa                                                                                                  >  2:229319/1‑101 (MQ=21)
     aaacggggggCGCTGCCCCCCTCTTGGCCTTCGGCCAATTCACCCCCCGGG.ATATTTCGGGCAAGATAAAGAAG.GTCTTTTCCTTTCATC.TTGCCCG.AAat                                                                                                 <  1:3998860/96‑1 (MQ=21)
     aaacggggggCGCTGCCCCCCTCTTGGCCTTCGGCCAATTCACCCCCCGGG.ATATTTCGGGCAAGATAAAGAAG.GTCTTTTCCTTTCATC.TTGCCCG.AAat                                                                                                 <  2:107585/96‑1 (MQ=21)
     aaacggggggCGCTGCCCCCCTCTTGGCCTTCGGCCAATTCACCCCCCGGG.AAATTCCGGCCAAAAAAAGAAAG.GACTTTTCTTTTCTTC.TTCCCCG.aaaa                                                                                                 >  2:4359411/6‑100 (MQ=2)
     aaacggggggCGCTGCCCCCCTCTTGGCCTTCGGCCAATTCACCCCCCCGG.ATATTTCGGGCAAGATAAAGAAG.GGCTTTTCCTTTCATC.TTGCCCC.AAat                                                                                                 >  1:2033105/6‑101 (MQ=21)
      ggcagggggCGCTG.CCCCCTCTTGGCCTTTGGCCAATTCA.CCCCTGGG.ATATTTCGGGCAAGAT...GAA..ATCATGCCCCTTCATC.TTGCCGC.AAATATCCCgg                                                                                          >  1:4977195/5‑101 (MQ=0)
       acggggggCGCTGCCCCCCTCTTGGCCTTCGGCCAATTCTCCCCCCGGG.ATATTTCGGGCAAGATAAAGAAG.GTCTTTTCCTTTCATC.TTGCCCG.AAatat                                                                                               <  2:3198269/98‑1 (MQ=25)
       acggggggCGCTGCCCCCCTCTTGGCCTTCGGCCAATTCACCCCCCGGG.ATATTTCGGGCAAGATAAAGAAG.GTCTTTTCCTTTCATC.TTGCCCG.AAatat                                                                                               <  1:2975659/98‑1 (MQ=25)
        cggggggCGCTGCCCCCCTCTTGGCCTTCGGCCAATTCACCCCCCGGGTAAATTTGGGGC.AGATAAAAAAGAGCTTTTTCTTTTCTTT.TTGCCCC.aaaatc                                                                                               >  1:226992/3‑98 (MQ=0)
          gggggCGCTGCCCCCCTCTTGGCCTTCGGCCAATTCACCCCCCGGG.ATATTTCGGGCAAGATAAAGAAG.GTCTTTTCCTTTCATC.TTGCCCG.AAATATccc                                                                                            <  1:2137791/101‑1 (MQ=25)
          gggggCGCTG.CCCCCTCTTGGCCTTTGGCCAATTCA.CCCCTGGG.ATATTTCGGGCAAGAT...GAA..ATCATGCCCCTTCATC.TTGCCGC.AAATATCCCgggggg                                                                                      >  2:1480212/1‑100 (MQ=0)
           ggggCGCTGCCCCCCTCTTGGCCTTCGGCCAATTCACCCCCCGGG.ATATTTCGGGCAAGATAAAGAAG.GTCTTTTCCTTTCATC.TTGCCCG.AAATATCCCg                                                                                           <  2:2456950/101‑1 (MQ=25)
           ggggCCCTGCCCCCTTTTGGC.CTTTGGCCAATTCA.CCCCTGGG.ATATTTCGGGCAAGAT...GAA..ATCATGCCCCTTCATC.TTGCCGC.AAATATCCCgggggga                                                                                     <  2:359618/101‑3 (MQ=2)
            gggCGCTGCCCCCCTCTTGGCCTTCGGCCAATTCACCCCCCGGG.ATATTTCGGGCAAGATAAAGAAG.GTCTTTTCCTTTCATC.TTGCCCG.AAATATCCCgg                                                                                          <  2:2353969/101‑1 (MQ=37)
               cgcTGCCCCCCTCTTGGCCTTCGGCCAATTCACCCCCCGGG.ATATTTCGGGCAAGATAAAGAAG.GTCTTTTCCTTTCATC.TTGCCCG.AAATATCCCggggg                                                                                       <  2:4869814/101‑1 (MQ=37)
                   gCCCCCCTCTTGGCCTTCGGCCAATTCTCCCCCCGGG.ATATTTCGGTCAAGATAAAGAAG.GTCTTTTCCTTTCATC.TTTCCCG.AAATATCCCGGGGGAGcc                                                                                   <  1:709533/101‑1 (MQ=25)
                   gCCCCCCTCTTGGCCTTCGGCCAATTCACCCCCCGGG.ATATTTCGGGCAAGATAAAGAAG.GTCTTTTCCTTTCATC.TTGGCCG.AAATATCCCGGGGGAGcc                                                                                   >  2:4973558/1‑101 (MQ=25)
                   gCCCCCCTCTTGGCCTTCGGCCAATTCACCCCCCGGG.ATATTTCGGGCAAGATAAAGAAG.GTCTTTTCCTTTCATC.TTGCCCG.AAATATCCCGGGGGAGcc                                                                                   <  2:4036376/101‑1 (MQ=37)
                         ctcttggccTTCGGCCAATTCACCCCCCGGG.ATATTTCGGGCAAGATAAAGAAG.GTCTTTTCCTTTCATC.TTGCCCG.AAATATCCCGGGGGAGCCCGCTTg                                                                             <  1:3269055/93‑1 (MQ=255)
                         ctcttggccTTCGGCCAATTCACCCCCCGGG.ATATTTCGGGCAAGATAAAGAAG.GTCTTTTCCTTTCATC.TTGCCCG.AAATATCCCGGGGGAGCCCGCTTg                                                                             <  1:1577742/93‑1 (MQ=255)
                              ggccTTCGGCCAATTCACCCCCCGGG.ATATTTCGGGCAAGATAAAGAAG.GTCTTTTCCTTTCATC.TTGCCCG.AAATATCCCGGGGGAGCCCGCTTgcgggc                                                                        >  1:2694411/4‑101 (MQ=255)
                                ccTTCGGCCAATTCACCCCCCGGG.ATAATTTGGGCAAGATAAAAAAA.GGGTTTTTCTTTTATCTTTTCCCC.AAAATCCCCGGGGGGGCCCCCTGgcggggg                                                                       >  1:4235308/2‑99 (MQ=12)
                                 cTTCGGCCAATTCACCCCCCGGG.ATATTTCGGGCAAGATAAAGAAG.GTCTTTTCCTTTCATC.TTGCCCG.AAATATCCCGGGGGAGCCCGCTTGcgggcggg                                                                     <  1:4678742/101‑1 (MQ=255)
                                 cTTCGGCCAATTCACCCCCCGGG.ATATTTCGGGCAAGATAAAGAAG.GTCTTTTCCTTTCATC.TTGCCCG.AAATATCCCGGGGGAGCCCGCTTGcgggcggg                                                                     >  1:3400121/1‑101 (MQ=255)
                                 cTTCGGCCAATTCACCCCCCGGG.ATATTTCGGGCAAGATAAAGAAG.GTCTTTTCCTTTCATC.TTGCCCG.AAATATCCCGGGGGAGCCCGCTTGcgggcggg                                                                     >  2:4518154/1‑101 (MQ=255)
                                      ggcAAATCCCCCCCCGGG.GTTTTTTCGGGAAAAAAAAAAAG.GGTTTTTTCTTTCCTT.TTTCCCGAAAATTCCCCGGGGGAGCCCGCTTGCGGGCGGGGGCa                                                                 <  1:1020567/99‑1 (MQ=255)
                                          aTTCACCCCCCGGG.ATATTTCGGGCAAGATAAAGAAG.GTCTTTTCCTTTCATC.TTGCCCG.AAATATCCCGGGGGAGCCCGCTTGCGGGGGGGGGCAGCGcc                                                            >  2:2375588/1‑101 (MQ=255)
                                          aTTCACCCCCCGGG.ATATTTCGGGCAAGATAAAGAAG.GGCTTTTTCTTTCATC.TTGCCCG.AAATATCCCGGGGGAGCCCGCTTGCGGGCGGGGGCAGCGcc                                                            >  2:2525011/1‑101 (MQ=255)
                                           ttCACCCCCCGGG.ATATTTCGGGCAAGATAAAGAAG.GTCTTTTCCTTTCATC.TTGCCCG.AAATATCCCGGGGGAGCCCGCTTGCGGGCGGGGGCAGCGccc                                                           >  2:1581810/1‑101 (MQ=255)
                                            tCACCCCCCGGC.ATATTTCGGGCAAGATAAAGAAG.GTCTTTTCCTTTCAGC.TTGCCCG.AAATATCCCGGGGGAGCCCGCTTGCGGGCGGGGGCAGCGcccc                                                          <  2:3964271/101‑1 (MQ=255)
                                             cACCCCCCGGG.ATATTTCGGGCAAGATAAAGAAG.GTCTTTTCCTTTCATC.TTGCCCG.AAATATCCCGGGGGAGCCCGCTTGCGGGCGGGGGCAGCGccccc                                                         >  2:2158328/1‑101 (MQ=255)
                                              gccccccGGG.ATATTTCGGGCAAGATAAAGAAG.GTCTTTTCCTTTCATC.TTGCCCG.AAATATCCCGGGGGAGCCCGCTTGCGGGCGGGGGCAGCGCCCCCg                                                        <  1:554365/100‑1 (MQ=255)
                                                    cGGG.ATATTTCGGGCAAGATAAAGAAG.GTCTTTTCCTTTCATC.TTGCCCG.AAATATCCCGGGGGAGCCCGCTTGCGGGCGGGGGCAGCGCCCCCGCAgcgc                                                  >  2:753876/1‑101 (MQ=255)
                                                     ggg.ATATTTCGGTCAATATAAAGAAG.GTCTTTTCCTTTCTTC.TTGCCCA.AAATATCCCGGGGGAGCCCGCTTGCGGGCGGGGGCATCGCCCCCGAAGCGcc                                                 <  1:1987380/101‑1 (MQ=255)
                                                     ggg.ATATTTCGGGCAAGATAAAGAAG.GTCTTTTCCTTTCATT.TTGCCCG.AAATATCCCGGGGAAGCCCGCTTGCGGGCGGGGGCAGCGCCCCCGAAGCGcc                                                 >  1:857155/1‑101 (MQ=255)
                                                     ggg.ATATTTCGGGCAAGATAAAGAAG.GTCTTTTCCTTTCATC.TTGCCCG.AAATATCCCGGGGGAGCCCGCTTGCGGGCGGGGGCAGCGCCCCCGAAGCGcc                                                 >  1:3596420/1‑101 (MQ=255)
                                                     ggg.ATATTTCGGGCAAGATAAAGAAG.GTCTTTTCCTTTCATC.TTGCCCG.AAATATCCCGGGGGAGCCCGCTTGCGGGCGGGGGCAGCGCCCCCGAAGCGcc                                                 >  1:2764816/1‑101 (MQ=255)
                                                       g.ATATTTCGGGCAAGATAAAGAAG.GTCTTTTCCTTTCATC.TTGCCCG.AAATATCCCGGGGGAGCCCGCTTGCGGGCGGGGGCAGCGCCCCCGAAGCGCCCg                                               <  2:2749651/101‑1 (MQ=255)
                                                       g.ATATTTCGGGCAAGATAAAGAAG.GTCTTTTCCTTTCATC.TTGCCCG.AAATATCCCGGGGGAGCCCGCTTGCGGGCGGGGGCAGCGCCCCCGAAGCGCCCg                                               <  1:2403758/101‑1 (MQ=255)
                                                         atatTTCGGGCAAGATAAAGAAG.GTCTTTTCCTTTCATC.TTGCCCG.AAATATCCCGGGGGAGCCCGCTTGCGGGCGGGGGCAGCGCCCCCGAAGCGCCCgc                                              >  2:733843/1‑100 (MQ=255)
                                                         atatTTCGGGCAAGATAAAGAAG.GTCTTTTCCTTTCATC.TTGCCCG.AAATATCCCGGGGGAGCCCGCTTGCGGGCGGGGGCAGCGCCCCCGAAGCGCCCGa                                              >  1:2907800/1‑101 (MQ=255)
                                                         atatTTCGGGCAAGATAAAGAAG.GTCTTTTCCTTTCATC.TTGCCCG.AAATATCCCGGGGGAGCCCGCTTGCGGGCGGGGGCAGCGCCCCCGAAGCGCCCGa                                              >  2:3544099/1‑101 (MQ=255)
                                                          tatTTCGGGCAAGATAAAGAAG.GTCTTTTCCTTTCATC.TTGCCCG.AAATATCCCGGGGGAGCCCGCTTGCGGGCGGGGGCAGCGCCCCCGAAGCGCCCGAc                                             >  2:3366566/1‑101 (MQ=255)
                                                           atTTCGTGCAAGATAAAGAAG.GTCTTTTCCTTTCATC.TTGCCCG.AAATATCCCGGGGGAGCCCGCTTGCGGGCGGGGGCAGCGCCCCCGAAGCGCCCGAcg                                            <  1:812310/101‑1 (MQ=255)
                                                               cGGGCAAGATAAAGAAG.GTCTTTTCCTTTCATC.TTGCCCG.AAATATCCCGGGGGAGCCCGCTTGCGGGCGGGGGCAGCGCCCCCGAAGCGCCCGAcgcgcg                                        >  1:2528990/1‑101 (MQ=255)
                                                                 ggCAAGATAAAGAAG.GTCTTTTCCTTTCATC.TTGCCCG.AAATATCCCGGGGGAGCCCGCTTGCGGGCGGGGGCAGCGCCCCCGAAGCGCCCGACGCgcggc                                      <  1:3030695/101‑1 (MQ=255)
                                                                 ggCAAGATAAAGAAG.GTCTTTTCCTTTCATC.TTGCCCG.AAATATCCCGGGGGAGCCCGCTTGCGGGCGGGGGCAGCGCCCCCGAAGCGCCCGACGCgcggc                                      <  1:3670577/101‑1 (MQ=255)
                                                                    aaGATAAAGAAG.GTCTTTTCCTTTCATC.TTGCCCG.AAATATCCCGGGGGAGCCCGCTTGCGGGCGGGGGCAGCGCCCCCGAAGCGCCCGACGCGCGGCggg                                   >  2:1638185/1‑101 (MQ=255)
                                                                    aaGATAAAGAAG.GTCTTTTCCTTTCATC.TTGCCCG.AAATATCCCGGGGGAGCCCGCTTGCGGGCGGGGGCAGCGCCCCCGAAGCGCCCGACGCGCGGCggg                                   >  2:88232/1‑101 (MQ=255)
                                                                     aGATAAAGAAG.GTCTTTTCCTTTCATC.TTGCCCG.AAATATCCCGGGGGAGCCCGCTTGCGGGCGGGGGCAGCGCCCCCGAAGCGCCCGACGCGCGGCGGGc                                  >  2:2044745/1‑101 (MQ=255)
                                                                     aGATAAAGAAG.GTCTTTTCCTTTCATC.TTGCCCG.AAATATCCCGGGGGAGCCCGCTTGCGGGCGGGGGCAGCGCCCCCGAAGCGCCCGACGCGCGGCGGGc                                  >  1:650738/1‑101 (MQ=255)
                                                                     aGATAAAGAAG.GTCTTTTCCTTTCATC.TTGCCCG.AAATATCCCGGGGGAGCCCGCTTGCGGGCGGGGGCAGCGCCCCCGAAGCGCCCGACGCGCGGCGGGc                                  >  1:3104665/1‑101 (MQ=255)
                                                                       aTAAAGAAG.GTCTTTTCCTTTCATC.TTGCCCG.AAATATCCCGGGGGAGCCCGCTTGCGGGCGGGGGCAGCGCCCCCGAAGCGCCCGACGCGCGGCGGGCtt                                >  2:47594/1‑101 (MQ=255)
                                                                       aTAAAGAAG.GTCTTTTCCTTTCATC.TTGCCCG.AAATATCCCGGGGGAGCCCGCTTGCGGGCGGGGGCAGCGCCCCCGAAGCGCCCGACGCGCGGCGGGCtt                                >  2:3760075/1‑101 (MQ=255)
                                                                       aTAAAGAAG.GTCTTTTCCTTTCATC.TTGCCCG.AAATATCCCGGGGGAGCCCGCTTGCGGGCGGGGGCAGCGCCCCCCAAGCGCCCGACGAGCGGCGGGGtt                                >  2:3006850/1‑101 (MQ=255)
                                                                        tAAAGAAG.GTCTTTTCCTTTCATC.TTGCCCG.AAATATCCCGGGGGAGCCCGCTTGCGGGCGGGGGCAGCGCCCCCGAAGCGCCCGACGCGCGGCGGGCTTg                               >  1:3944030/1‑101 (MQ=255)
                                                                        tAAAGAAG.GTCTTTTCCTTTCATC.TTGCCCG.AAATACCCCGGGGGAGCCCGCTTGCGGGCGGGGGCAGCGCCCCCGAAGCGCCCGACGCGCGGCGGGCTTg                               <  2:3496729/101‑1 (MQ=255)
                                                                         aaaGAAG.GTCTTTTCCTTTCATC.TTGCCCG.AAATATCCCGGGGGAGCCCGCTTGCGGGCGGGGGCAGCGCCCCCGAAGCGCCCGACGCGCGGCGGGCTTgc                              >  2:3912354/1‑101 (MQ=255)
                                                                         aaaGAAG.GTCTTTTCCTTTCATC.TTGCCCG.AAATATCCCGGGGGAGCCCGCTTGCGGGCGGGGGCAGCGCCCCCGAAGCGCCCGACGCGCGGCGGGCTTgc                              >  2:4333584/1‑101 (MQ=255)
                                                                         aaaGAAG.GTCTTTTCCTTTCATC.TTGCCCG.AAATATCCCGGGGGAGCCCGCTTGCGGGCGGGGGCAGCGCCCCCGAAGCGCCCGACGCGCGGCGGGCTTgc                              >  2:2262552/1‑101 (MQ=255)
                                                                          aagaag.GTCTTTTCCTTTCATC.TTGCCCG.AAATATCCCGGGGGAGCCCGCTTGCGGGCGGGGGCAGCGCCCCCGAAGCGCCCGACGCGCGGCGGGCTTgcc                             >  2:4710155/1‑101 (MQ=255)
                                                                          aagaag.GTCTTTTCCTTTCATC.TTGCCCG.AAATATCCCGGGGGAGCCCGCTTGCGGGCGGGGGCAGCGCCCCCGAAGCGCCCGACGCGCGGCGGGCTTgcc                             >  1:2454976/1‑101 (MQ=255)
                                                                          aagaag.GTCTTTTCCTTTCATC.TTGCCCG.AAATATCCCGGGGGAGCCCGCTTGCGGGCGGGGGCAGCGCCCCCGAAGCGCCCGACGCGCGGCGGGCTTgcc                             >  2:2240128/1‑101 (MQ=255)
                                                                              ag.GTCTTTTCCTTTCATC.TTGCCCG.AAATATCCCGGGGGAGCCCGCTTGCGGGCGGGGGCAGCGCCCCCGAAGCGCCCGACGCGCGGCGGGCTTgcccgcc                         <  2:3057205/101‑1 (MQ=255)
                                                                                  tCTTTTCCTTTCATC.TTGCCCG.AAATATCCCGGGGGAGCCCGCTTGCGGGCGGGGGCAGCGCCCCCGAAGCGCCCGACGCGCGGCGGGCTTGCCCGCCCaa                      >  2:4479791/1‑101 (MQ=255)
                                                                                   cTTTTCCTTTCATC.TTGCCCG.AAATATCCCGGGGGAGCCCGCTTGCGGGCGGGGGGAGCGCCCCCGAAGCGCCCCACGCCCCGCGGGGTTTCCCCCCCAag                     >  2:4505076/1‑101 (MQ=37)
                                                                                   cTTTTCCTTTCATC.TTGCCCG.AAATATCCCGGGGGAGCCCGCTTGCGGGCGGGGGCAGCGCCCCCGAAGCGCCCGACGCGCGGCGGGCTTGCCCGCCCAag                     >  2:3847062/1‑101 (MQ=255)
                                                                                    ttttCCTTTCATC.TTGCCCG.AAATATCCCGGGGGAGCCCGCTTGCGGGCGGGGGCAGCGCCCCCGAAGCGCCCGACGCGCGGCGGGCTTGCCCGCCCAagg                    >  2:4493705/1‑101 (MQ=255)
                                                                                     tttCCTTTCATC.TTGCCCG.AAATATCCCGGGGGAGCCCGCTTGCGGGCGGGGGCAGCGCCCCCGAAGCGCCCGACGCGCGGCGGGCTTGCCCGCCCAagga                   >  1:4260663/1‑101 (MQ=255)
                                                                                      ttCCTTTCATC.TTGCCCG.AAATATCCCGGGGGAGCCCGCTTGCGGGCGGGGGCAGCGCCCCCGAAGCGCCCGACGCGCGGCGGGCTTGCCCGCCCAaggag                  >  2:1330689/1‑101 (MQ=255)
                                                                                      ttCCTTTCATC.TTGCCCG.AAATATCCCGGGGGAGCCCGCTTGCGGGCGGGGGCAGCGCCCCCGAAGCGCCCGACGCGCGGCGGGCTTGCCCGCCCAaggag                  >  2:1395306/1‑101 (MQ=255)
                                                                                      ttCCTTTCATC.TTGCCCG.AAATATCCCGGGGGAGCCCGCTTGCGGGCGGGGGCAGCGCCCCCGAAGCGCCCGACGCGCGGCGGGCTTGCCCGCCCAaggag                  >  1:1832277/1‑101 (MQ=255)
                                                                                       tCCTTTCATC.TTGCCCG.AAATATCCCGGGGGAGCCCGCTTGCGGGCGGGGGCAGCGCCCCCGAAGCGCCCGACGCGCGGCGGGCTTGCCCGCCCAaggagg                 <  1:4710478/101‑1 (MQ=255)
                                                                                       tCCTTTCATC.TTGCCCG.AAATATCCCGGGGGAGCCCGCTTGCGGGCGGGGGCAGCGCCCCCGAAGCGCCCGACGCGCGGCGGGCTTGCCCGCCCAaggagg                 <  1:4481228/101‑1 (MQ=255)
                                                                                        ccTTTCATC.TTGCCCG.AAATATCCCGGGGGAGCCCGCTTGCGGGCGGGGGCAGCGCCCCCGAAGCGCCCGACGCGCGGCGGGCTTGCCCGCCCAAGGAggg                >  2:2446397/1‑101 (MQ=255)
                                                                                        ccTTTCATC.TTGCCCG.AAATACCCCGGGGGAGCCCGCTTGCGGGCGGGGGCAGCGCCCCCGAAGCGCCCGACGCGCGGCGGGCTTGCCCGCCCAAGGAggg                <  1:993428/101‑1 (MQ=255)
                                                                                         cTTTCATC.TTGCCCG.AAATATCCCGGGGGAGCCCGCTTGCGGGCGGGGGCAGCGCCCCCGAAGCGCCCGACGCGCGGCGGGCTTGCCCGCCCAAGGAgggg               >  1:1699350/1‑100 (MQ=255)
                                                                                          tttCATC.TTGCCCG.AAATATCCCGGGGGAGCCCGCTTGCGGGCGGGGGCAGCGCCCCCGAAGCGCCCGACGCGCGGCGGGCTTGCCCGCCCAAGGAGGgcg              >  1:3808134/1‑101 (MQ=255)
                                                                                          tttCATC.TTGCCCG.AAATATCCCGGGGGAGCCCGCTTGCGGGCGGGGGCAGCGCCCCCGAAGCGCCCGACGCGCGGCGGGCTTGCCCCCCCAAGGAGGgcg              >  1:3831648/1‑101 (MQ=255)
                                                                                           ttCATC.TTGCCCG.AAATATCCCGGGGGAGCCCGCTTGCGGGCGGGGGCAGCGCCCCCGAAGCGCCCGACGCGCGGCGGGCTTGCCCGCCCAAGGAGGgcgc             <  2:2694226/101‑1 (MQ=255)
                                                                                            tCATC.TTGCCCG.AAATATCCCGGGGGAGCCCGCTTGCGGGCGGGGGCAGCGCCCCCGAAGCGCCCGACGCGCGGCGGGCTTGCCCGCCCAAGGAGGGCGCa            <  2:2907587/101‑1 (MQ=255)
                                                                                            tCATC.TTGCCCG.AAATATCCCGGGGGAGCCCGCTTGCGGGCGGGGGCAGCGCCCCCGAAGCGCCCGACCCGCGGCGGGCTTGCCCGCCCAAGGAGGGCGCa            >  1:4762914/1‑101 (MQ=255)
                                                                                             cATC.TTGCCCG.AAATATCCCGGGGGAGCCCGCTTGCGGGCGGGGGCAGCGCCCCCGAAGCGCCCGACGCGCGGCGGGCTTGCCCGCCCAAGGAGGGCGCaa           >  1:4614807/1‑101 (MQ=255)
                                                                                              aTC.TTGCCCG.AAATATCCCGGGGGAGCCCGCTTGCGGGCGGGGGCAGCGCCCCCGAAGCGCCCGACGCGCGGCGGGCTTGCCCGCCCAAGGAGGGCGCAAc          >  2:1542318/1‑101 (MQ=255)
                                                                                                c.TTGCCCG.AAATATCCCGGGGGAGCCCGCTTGCGGGCGGGGGCAGCGCCCCCGAAGCGCCCGACGCGCGGCGGGCTTGCCCCCCCAAGGAGGGCGCAAccc        >  1:3887249/1‑100 (MQ=255)
                                                                                                  ttGCCCG.AAATATCCCGGGGGAGCCCGCTTGCGGGCGGGGGCAGCGCCCCCGAAGCGCCCGACGCGCGGCGGGCTTGCCCGCCCAAGGAGGGCGCAACCGc       >  1:1194678/1‑101 (MQ=255)
                                                                                                   tGCCCG.AAATATCCCGGGGGAGCCCGCTTGCGGGCGGGGGGAGCGCCCCCGAAGCGCCCGACGCGCGGCGGGGCTGCCCGCCCCAGGCGCGCGCAAACgtc      >  2:2610259/1‑99 (MQ=37)
                                                                                                   tGCCCG.AAATATCCCGGGGGAGCCCGCTTGCGGGCGGGGGCAGCGCCCCCGAAGCGCCCGACGCGCGGCGGGCTTGCCCGCCCAAGGAGGGCGCAACCGcc      <  2:2183244/101‑1 (MQ=255)
                                                                                                      ccG.AAATATCCCGGGGGAGCCCGCTTGCGGGCGGGGGCAGCGCCCCCGAAGCGCCCGACGCGCGGCGGGCTTGCCCGCCCAAGGAGGGGGCAACCGCCTca   >  1:1128085/1‑101 (MQ=255)
                                                                                                       cG.AAATATCCCGGGGGAGCCCGCTTGCGGGCGGGGGCAGCGCCCCCGAAGCGCCCGACGCGCGGCGGGCTTGCCCGCCCAAGGAGGGCGCAACCGCCTcag  <  1:3967466/101‑1 (MQ=255)
                                                                                                       |                                                                                                     
AGACGAAAC.GGGGGCGCTGCCCCCCTCCTG..CTTCGGCCAATTCACCCCCCGGG.ATATTTCGGGCAAGATAAAGAAG.GTCTTTTCCTTTCATC.TTGCCCG.AAATATCCCGGGGGAGCCCGCTTGCGGGCGGGGGCAGCGCCCCCGAAGCGCCCGACGCGCGGCGGGCTTGCCCGCCCAAGGAGGGCGCAACCGCCTCAG  >  NC_014034/2566528‑2566725

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