Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_014034 2,566,577 C→T 2.5% intergenic (+49/+145) 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,5770C→T 2.5% 4.6 ‑0.5 ‑0.0202intergenic (+49/+145)purD/sdaAphosphoribosylamine‑‑glycine ligase/L‑serine ammonia‑lyase
Reads supporting (aligned to +/- strand):  new base (2/3):  ref base (55/142):  total (57/145)
Fisher's exact test strand distribution p-value = 6.22e-01
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 9.91e-01

Only 100 of 224 total aligned reads displayed.

CTGGCCCGAGGGGTTCTGCCGCCGCGACATCGGCTGGCGCGCGCTCTGAGACGAAAC.GGGGG.CGCTGCCCCC.CTCCTG..CTTCGGCCAATTCACCCCCCGGG.ATATTTCGGGCAAGATAAAGAAG.GTCTTTTCCTTTCATCTTGCCCG.AAATATCCCGGGGGAGCCCGCTTGCGGGCGGGGGCAGCGCCCCCGAAGCGC  >  NC_014034/2566480‑2566677
                                                                                                      |                                                                                                       
cTGGCCCGAGGGGTTCTGCCGCCGCGACATCGGCTGGCGCGCGCTCTGAGACGAAACGGGGGG.CGCTGCCCCC.CTCTTGGCCTTCGGCCAATTCAcccccc                                                                                                         <  1:3407201/101‑1 (MQ=25)
cTGGCCCGAGGGGTTCTGCCGCCGCGACATCGGCTGGCGCGCGCTCTGAGACGAAACGGGGGG.CGCTGCCCCC.CTCTTGGCCTTCGGCCAATTCAcccccc                                                                                                         <  1:2895076/101‑1 (MQ=25)
 tGGCCCGAGGGGTTCTGCCGCCGCGACATCGGCTGGCGCGCGCTCTGAGACGAAACGGGGGG.CGCTGCCCCC.CCCCTGGCCCTCCGCCGATTCACCCCCCg                                                                                                        >  1:3059487/1‑101 (MQ=255)
        aGGGGTTCTGCCGCCGCGACATCGGCTGGCGCGCGCTCTGAGACGAAACGGGGGG.CGCTGCCCCC.CTCTTGGCCTTCGGCCAATTCACCCCCCGGG.ATAtt                                                                                                <  1:15296/101‑1 (MQ=21)
         ggggTTCTGCCGCCGCGACATCGGCTGGCGCGCGCTCTGAGACGAAACGGGGGG.CGCTGCCCCC.CTCTTGGCCTTCGGCCAATTCACCCCCCGGG.ATAttt                                                                                               <  2:82969/101‑1 (MQ=21)
          gggTTCTGCCGCCGCGACATCGGCTGGCGCGCGCTCTGAGACGAAACGGGGGG.CGCTGCCCCC.CTCTTGGCCTTCGGCCAATTCACCCCCCGGG.ATATTTc                                                                                              <  2:3172939/101‑1 (MQ=21)
           ggTTCTGCCGCCGCGACATCGGCTGGCGCGCGCTCTGAGACGAAACGGGGGG.CGCTGCCCCC.CTCTTGGCCTTCGGCCAATTCACCTCCCGGG.ATATTTCg                                                                                             <  2:4988316/101‑1 (MQ=21)
            gTTCTGCCGCCGCGACATCGGCTGGCGCGCGCTCTGAGACGAAACGGGGGG.GGCTGCCCCC.CTCTTGGCCTTCGGCCCATTTCCCCCCCGGG.GGATTTCgg                                                                                            >  1:2398914/1‑101 (MQ=21)
              tCTGCCGCCGCGACATCGGCTGGCGCGCGCTCTGAGACGAAACGGGGGG.CGCTGCCCCC.CTCTTGGCCTTCGGCCAATTCACCCCCCGGG.ATATTTCGGGc                                                                                          <  2:2196847/101‑1 (MQ=21)
                 gcctccGCGACATCGGCTGGCGCGCGTTCTGAGACGAAACGGGGGG.CGCTGCCCCC.CTCTTGGCCTTCGGCCAATTCACCCCCCGGG.ATATTTCGGGCAAg                                                                                       <  2:558382/101‑1 (MQ=18)
                 gccgccGCGACATCGGCTGGCGCGCGCTCTGAGACGAAACGGGGGG.CGCTGCCCCC.CTCTTGGCCTTCGGCCAATTCACCCCCCGGG.ATATTTCGGGCAAg                                                                                       <  2:404170/101‑1 (MQ=21)
                 gccgccGCGACATCGGCTGGCGCGCGCTCTGAGACGAAACGGGGGG.CGCTGCCCCC.CTCTTGGCCTTCGGCCAATTCACCCCCCGGG.ATATTTCGGGCAAg                                                                                       <  1:1811994/101‑1 (MQ=21)
                  ccgccgCGACATCGGCTGGCGCGCGCTCTGAGACGAAACGGGGGG.CGCTGCCCCC.CTCTTGGCCTTCGGCCAATTCACCCCCCGGG.ATATTTCGGGCAAGa                                                                                      <  1:4252492/101‑1 (MQ=21)
                  ccgccgCGACATCGGCTGGCGCGCGCTCTGAGACGAAACGGGGGG.CGCTGCCCCC.CTCTTGGCCTTCGGCCAATTCACCCCCCGGG.ATATTTCGGGCAAGa                                                                                      <  1:3738925/101‑1 (MQ=21)
                   cgccgcGACATCGGCTGGCGCGCGCTCTGAGACGAAACGGGGGG.CGCTGCCCCC.CTCTTGGCCTTCGGCCAATTCACCCCCCGGG.ATATTTCGGGCAAGAt                                                                                     <  1:4733697/101‑1 (MQ=21)
                   cgccgcGACATCGGCTGGCGCGCGCTCTGAGACGAAACGGGGGG.CGCTGCCCCC.CTCTTGGCCTTCGGCCAATTCACCCCCCGGG.ATATTTCGGGCAAGAt                                                                                     <  2:475921/101‑1 (MQ=21)
                    gccgcGACATCGGCTGGCGCGCGCTCTGAGACGAAACGGGGGG.CGCTGCCCCC.CTCTTGGCCTTCGGCCAATTCACCCCCCGGG.ATATTTCGGGCAAGATa                                                                                    <  1:758907/101‑1 (MQ=21)
                     ccgcGACATCGGCTGGCGCGCGCTCTGAGACGAAACGGGGGG.CGCTGCCCCC.CTCTTGGCCTTCGGCCAATTCACCCCCCGGG.ATATTTCGGGCAAGATaa                                                                                   <  2:4311421/101‑1 (MQ=21)
                     ccgcGACATCGGCTGGCGCGCGCTCTGAGACGAAACGGGGGG.CGCTGCCCCC.CTCTTGGCCTTCGGCCAATTCACCCCCCGGG.ATATTTCGGGCAAGATaa                                                                                   <  2:2102822/101‑1 (MQ=21)
                      cgcgACATCGGCTGGCGCGCGCTCTGAGACGAAACGGGGGG.CGCTGCCCCC.CTCTTGGCCTTCGGCCAATTCACCCCCCGGG.ATATTTCGGGCAAGATaaa                                                                                  <  1:4025861/101‑1 (MQ=21)
                      cgcgACATCGGCTGGCGCGCGCTCTGAGACGAAACGGGGGG.CGCTGCCCCC.CTCTTGGCCTTCGGCCAATTCACCCCCCGGG.ATATTTCGGGCAAGATaaa                                                                                  <  1:431612/101‑1 (MQ=21)
                       gcgACATCGGCTGGCGCGCGCTCTGAGACGAAACGGGGGG.CGCTGCCCCC.CTCTTGGCCTTCGGCCAATTCACCCCCCGGG.ATATTTCGGGCAAGATAaag                                                                                 <  1:2556790/101‑1 (MQ=21)
                       gcgACATCGGCTGGCGCGCGCTCTGAGACGAAACGGGGGG.CGCTGCCCCC.CTCTTGGCCTTCGGCCAATTCACCCCCCGGG.ATATTTCGGGCAAGATAaag                                                                                 <  1:3001494/101‑1 (MQ=21)
                       gcgACATCGGCTGGCGCGCGCTCTGAGACGAAACGGGGGG.CGCTGCCCCC.CTCTTGGCCTTCGGCCAATTCACCCCCCGGG.ATATTTCGGGCAAGATAaag                                                                                 <  1:3206906/101‑1 (MQ=21)
                        cgACATCGGCTGGCGCGCGCTCTGAGACGAAACGGGGGG.CGCTGCCCCC.CTCTTGGCCTTCGGCCAATTCACCCCCCGGG.ATATTTCGGGCAAGATAaaga                                                                                <  2:2589035/101‑1 (MQ=21)
                         gACATCGGCTGGCGCGCGCTCTGGGACGAAACGGGGGG.CGCTGCCCCC.CTCTTGGCCTTCGGCCAATTCACCCCCCGGG.ATATTTCGGGCAAGATAaagaa                                                                               <  2:1840456/101‑1 (MQ=21)
                         gACATCGGCTGGCGCGCGCTCTGAGACGAAACGGGGGG.CGCTGCCCCC.CTCTTGGCCTTCGGCCAATTCAACCCCCGGG.ATATTTCGGGCAAGATAaagaa                                                                               <  1:1375187/101‑1 (MQ=21)
                          aCATCGGCTGGCGCGCGCTCTGAGACGAAACGGGGGG.CGCTGCCCCC.CTCTTGGCCTTCGGCCAATTCACCCCCCGGG.ATATTTTGGGCAAGATAaaaaag                                                                              >  1:1738075/1‑101 (MQ=18)
                          aCATCGGCTGGCGCGCGCTCTGAGACGAAACGGGGGG.CGCTGCCCCC.CTCTTGGCCTTCGGCCAATTCACCCCCCGGG.ATATTTCGGGCAAGATAaagaag                                                                              <  1:3942432/101‑1 (MQ=21)
                           cATCGGCTGGCGCGCGCTCTGAGACGAAACGGGGGG.CGCTGCCCCC.CTCTTGGCCTTCGGCCAATTCACCCCCCGGG.ATATTTCGGGCAAGATAAAGAAg.g                                                                            <  1:48548/101‑1 (MQ=21)
                             tCGGCTGGCGCGCGCTCTGAGACGAAACGGGGGG.CGCTGCCCCC.CTCTTGGCCTTCGGCCAATTCACCCCCCGGG.ATATTTCGGGCAAGATAAAGAAG.GTc                                                                          <  2:2440729/101‑1 (MQ=21)
                              cGGCTGGCGCGCGCTCTGAGACGAAACGGGGGG.CGCTGCCCCC.CTCTTGGCCTTCGGCCAATTCACCCCCCGGG.ATATTTCGGGCAAGATAAAGAAG.GGCt                                                                         >  1:598450/1‑101 (MQ=21)
                               ggCTGGCGCGCGCTCTGAGACGAAACGGGGGG.CGCTGCCCCC.CTCTTGGCCTTCGGCCAATTCACCCCCCGGG.ATATTTCGGGCAAGATAAAGAAG.GTCtt                                                                        <  2:337048/101‑1 (MQ=21)
                                  tGGCGCGCGCTCTGAGACGAAACGGGGGG.CGCTGCCCCC.CTCTTGGCCTTCGGCCAATTCACCCCCCGGG.ATATTTCGGGCAAGATAAAGAAG.GTCTTTTc                                                                     <  2:1768660/101‑1 (MQ=21)
                                  tGGCGCGCGCTCTGAGACGAAACGGGGGG.CGCTGCCCCC.CTCTTGGCCTTCGGCCAATTCACCCCCCGGG.ATATTTCGGGCAAGATAAAGAAG.GTCTTTTc                                                                     >  2:4917854/1‑101 (MQ=21)
                                   ggCGCGCGCTCTGAGACGAAACGGGGGG.CGCTGCCCCC.CTCTTGGCCTTCGGCCAATTCACCCCCCGGG.ATATTTCGGGCAAGATAAAGAAG.GTCTTTTcc                                                                    <  2:4673142/101‑1 (MQ=21)
                                   ggCGCGCGCTCTGAAACGAAACGGGGGG.CGCTGCCCCC.CTCTTGGCCTTCGGCCAATTCACCCCCCGGG.ATATTTCGGGCAAGATAAAGAAG.GTCTTTTcc                                                                    <  1:3651329/101‑1 (MQ=21)
                                    gcgcgcgcTCTGAGACGAAACGGGGGG.CGCTGCCCCC.CTCTTGGCCTTCGGCCAATTCACCCCCCGGG.ATATTTCGGGCAAGATAAAGAAG.GTCTTTTCCt                                                                   <  2:2861298/101‑1 (MQ=21)
                                    gcgcgcgcTCTGAGACGAAACGGGGGG.CGCTGCCCCC.CTCTTGGCCTTCGGCCAATTCACCCCCCGGG.ATATTTCGGGCAAGATAAAGAAG.GGCTTTTCCt                                                                   >  1:4505293/1‑101 (MQ=21)
                                      gcgcgcTCTGAGACGAAACGGGGGG.CGCTGCCCCC.CTCTTGGCCTTCGGCCAATTCACCCCCCGGG.ATATTTCGGGCAAGATAAAGAAG.GTCTTTTCCttt                                                                 <  1:3285597/101‑1 (MQ=21)
                                      gcgcgcTCTGAGACGAAACGGGGGG.CGCTGCCCCC.CTCTTGGCCTTCGGCCAATTCACCCCCCGGG.ATATTTCGGGCAAGATAAAGAAG.GTCTTTTCCttt                                                                 <  2:4016807/101‑1 (MQ=21)
                                      gcgcgcTCTGAGACGAAACGGGGGG.CGCTGCCCCC.CTCTTGGCCTTCGGCCAATTCACCCCCCGGG.ATATTTCGGGCAAGATAAAGAAG.GTCTTTTCCttt                                                                 <  1:3186209/101‑1 (MQ=21)
                                       cgcgcTCTGAGACGAAACGGGGGG.CGCTGCCCCC.CTCTTGGCCTTCGGCCAATTCACCCCCCGGG.ATATTTCGGGCAAGATAAAGAAG.GTCTTTTCCTTTc                                                                <  2:3302022/101‑1 (MQ=21)
                                         cgcTCTGAGACGAAAC.GGGGGTCGCTTCCCCC.CTCTTGGCCTTCGGCCAATTCACCCCCCGGG.ATATTTCGGGCAAGATAAAGAAG.GTCTTTTCCTTTCAt                                                              <  1:1616734/101‑1 (MQ=25)
                                          gcTCTGAGACGAACCGGGGGG.CGCTGCCCCC.TTCTTGGCCTTCGGCCAATTCACCCCCCGGG.ATATTTCGGGCAAGATAAAGAAG.GTCTTTTCCTTTCATc                                                             <  1:2652893/101‑1 (MQ=21)
                                          gcTCTGAGACGAAACGGGGGG.CGCTGCCCCC.CTCTTGGCCTTCGGCCAATTCACCCCCCGGG.ATATTTCGGGCAAGATAAAGAAG.GTCTTTTCCTTTCATc                                                             <  2:2891000/101‑1 (MQ=21)
                                           ctctGAGACGAAACGGGGGG.CGCTGCCCCC.CTCTTGGCCTTCGGCCAATTCACCCCCCGGG.ATATTTCGGGCAAGATAAAGAAG.GTCTTTTCCTTTTATCt                                                            >  1:2119601/1‑101 (MQ=21)
                                           ctctGAGACGAAACGGGGGG.CGCTGCCCCC.CTCTTGGCCTTCGGCCAATTCACCCCCCGGG.ATATTTCGGGCAAGATAAAGAAG.GTCTTTTCCTTTCATCt                                                            >  2:2387314/1‑101 (MQ=21)
                                            tctGAGACGAAACGGGGGG.CGCTGCCCCC.CTCTTGGCCTTCGGCCAATTCACCCCCCGGG.ATATTTCGGGCAAGATAAAGAAG.GTCTTTTCCTTTCATCtt                                                           >  1:3091542/1‑101 (MQ=21)
                                            tctGAGACGAAACGGGGGG.CGCTGCCCCC.CTCTTGGCCTTCGGCCAATTCACCCCCCGGG.ATATTTCGGGCAAGATAAAGAAG.GTCTTTTCCTTTCATCtt                                                           >  2:4397614/1‑101 (MQ=21)
                                            tctGAGACGAAACGGGGGG.CGCTGCCCCC.CTCTTGGCCTTCGGCCAATTCACCCCCCGGG.ATATTTCGGGCAAGATAAAGAAG.GTCTTTTCCTTTCATCtt                                                           <  1:88760/101‑1 (MQ=21)
                                            tctGAGACGAAACGGGGGG.CGCTGCCCCC.CTCTTGGCCTTCGGCCAATTCACCCCCCGGG.ATATTTCGGGCAAGATAAAGAAG.GTCTTTTCCTTTCATCtt                                                           >  2:612440/1‑101 (MQ=21)
                                            tctGAGACGAAACGGGGGG.CGCTGCCCCC.CTCTTGGCCTTCGGCCAATTCACCCCCCGGG.ATATTTCGGGCAAGATAAAGAAG.GTCTTCTCCTTTCATCtt                                                           <  1:2728105/101‑1 (MQ=21)
                                             ctGAGGCGAAACGGGGGG.CGCTGCCCCC.CTCTTGGCCTTCGGCCAATTCTCCCCCCGGG.ATATTTCGGGCAAGATAAAGAAG.GTCTTTTCCTTTCATCTTg                                                          <  1:2566418/101‑1 (MQ=21)
                                             ctGAGACGAAACGGGGGG.CGCTGCCCCC.CTCTTGGCCTTCGGCCAATTCACCCCCCGGG.ATATTTCGGGCAAGATAAAGAAG.GTCTTATCCTTTCATCTTg                                                          <  2:1974418/101‑1 (MQ=21)
                                                agaCGAAACGGGGGG.CGCTGCCCCC.CTCTTGGCCTTCGGCCAATTCACCCCCCGGG.ATATTTCGGGCAAGATAAAGAAG.GTCTTTTCCTTTCATCTTGccc                                                       <  1:191171/101‑1 (MQ=21)
                                                agaCGAAACGGGGGG.CGCTGCCCCC.CTCTTGGCCTTCGGCCAATTCACCCCCCGGG.ATATTTCGGGCAAGATAAAGAAG.GTCTTTTCCTTTCATCTTGccc                                                       <  2:1699276/101‑1 (MQ=21)
                                                 gaCGAAACGGGGGG.CGCTGCCCCC.CTCTTGGCCTTCGGCCAATTCACCCCCCGGG.ATATTTCGGGCAAGATAAAGAAG.GTCTTTTCCTTTCATCTTGCCCg                                                      <  2:3407080/101‑1 (MQ=21)
                                                 gaCGAAACGGGGGG.CGCTGCCCCC.CTCTTGGCCTTCGGCCAATTCACCCCCCGGG.ATATTTCGGGCAAGATAAAGAAG.GTCTTTTCCTTTCATCTTGCCCg                                                      <  1:514012/101‑1 (MQ=21)
                                                  aCGAAACGGGGGG.CGCTGCCCCC.CTCTTGGCCTTCGGCCAATTCACCCCCCGGG.ATATTTCGGGCAAGATAAAGAAG.GTCTTTTCCTTTCATCTTGCCCG.a                                                    <  2:539174/101‑1 (MQ=21)
                                                  aCGAAACGGGGGG.CGCTGCCCCC.CTCTTGGCCTTCGGCCAATTCACCCCCCGGG.ATATTTCGGGCAAGATAAAGAAG.GTCTTTTCCTTTCATCTTGCCCG.a                                                    <  2:4727329/101‑1 (MQ=21)
                                                   cGAAACGGGGGG.CGCTGCCCCC.CTCTTGGCCTTCGGCCAATTCACCCCCCGGG.ACATTTCGGGCAAGATAAAGAAG.GTCTTTTCCTTTCATCTTGCCCG.aa                                                   >  2:613597/1‑101 (MQ=21)
                                                     aaacgggggg.CGCTGCCCCC.CTCTTGGCCTTCGGCCAATTCACCCCCCGGG.ATATTTCGGGCAAGATAAAGAAG.GTCTTTTCCTTTCATCTTGCCCG.AAat                                                 <  2:107585/96‑1 (MQ=21)
                                                     aaacgggggg.CGCTGCCCCC.CTCTTGGCCTTCGGCCAATTCACCCCCCGGG.ATATTTCGGGCAAGATAAAGAAG.GTCTTTTCCTTTCATCTTGCCCG.AAat                                                 <  1:3998860/96‑1 (MQ=21)
                                                     aaacgggggg.CGCTGCCCCC.CTCTTGGCCTTCGGCCAATTCACCCCCCGGG.ATATTTCGGGCAAGATAAAGAAG.GTCTTTTCCTTTCATCTTGCCCG.AAat                                                 >  1:3827214/6‑101 (MQ=21)
                                                     aaacgggggg.CGCTGCCCCC.CTCTTGGCCTTCGGCCAATTCACCCCCCCGG.ATATTTCGGGCAAGATAAAGAAG.GGCTTTTCCTTTCATCTTGCCCC.AAat                                                 >  1:2033105/6‑101 (MQ=21)
                                                      ggcaggggg.CGCTG.CCCC.CTCTTGGCCTTTGGCCAATTCA.CCCCTGGG.ATATTTCGGGCAAGAT...GAA..ATCATGCCCCTTCATCTTGCCGC.AAATATCCCgg                                          >  1:4977195/5‑101 (MQ=0)
                                                        cgggggg.CGCTGCCCCC.CTCTTGGCCTTCGGCCAATTCACCCCCCGGGTAAATTTGGGGC.AGATAAAAAAGAGCTTTTTCTTTTCTTTTTGCCCC.aaaatc                                               >  1:226992/3‑98 (MQ=0)
                                                        cgggggg.CGCTGCCCCC.CTCTTGGCCTTCGGCCAATTCACCCCCCGGG.ATATTTCGGGCAAGATAAAGAAG.GTCTTTTCCTTTCATCTTGCCCG.AAATATc                                              >  2:3038455/3‑101 (MQ=25)
                                                        cgggggg.CGCTGCCCCC.CTCTTGGCCTTCGGCCAATTCACCCCCCGGG.ATATTTCGGGCAAGATAAAGAAG.GTCTTTTCCTTTCATCTTGCCCG.AAATATc                                              <  2:237318/99‑1 (MQ=25)
                                                        cgggggg.CGCTGCCCCC.CTCTTGGCCTTCGGCCAATTCACCCCCCGGG.ATATTTCGGGCAAGATAAAGAAG.GTCTTTTCCTTTCATCTTGCCCG.AAATATc                                              <  2:1023372/99‑1 (MQ=25)
                                                        cgggggg.CGCTGCCCCC.CTCTTGGCCTTCGGCCAATTCACCCCCCGGG.ATATTTCGGGCAAGATAAAGAAG.GTCTTTTACTTTCATCTTGCCCG.AAATATc                                              >  2:1100723/3‑101 (MQ=21)
                                                          ggggg.CGCTGCCCCCTCTTGGC..CTTTGGCCAATTCA.CCCCTGGG.ATATTTCGGGCAAGAT...GAA..ATCATGCCCCTTCATCTTGCCGC.AAATATCCCgggggg                                      <  2:4231871/101‑2 (MQ=0)
                                                          ggggg.CGCTGCCCCC.CTCTTGGCCTTCGGCCAATTCACCCCCCGGG.ATATTTCGGGCAAGATAAAGAAG.GTCTTTTCCTTTCATCTTGCCCG.AAATATccc                                            <  1:2137791/101‑1 (MQ=25)
                                                          ggggg.CGCTG.CCCC.CTCTTGGCCTTTGGCCAATTCA.CCCCTGGG.ATATTTCGGGCAAGAT...GAA..ATCATGCCCCTTCATCTTGCCGC.AAATATCCCgggggg                                      >  2:1480212/1‑100 (MQ=0)
                                                           gggg.CCCTGCCCCC.TTTTGGC.CTTTGGCCAATTCA.CCCCTGGG.ATATTTCGGGCAAGAT...GAA..ATCATGCCCCTTCATCTTGCCGC.AAATATCCCgggggga                                     <  2:359618/101‑3 (MQ=2)
                                                            ggg.CGCTGCCCCC.CTCTTGGCCTTCGGCCAATTCACCCCCCGGG.ATATTTCGGGCAAGATAAAGAAG.GTCTTTTCCTTTCATCTTGCCCG.AAATATCCCgg                                          <  2:2353969/101‑1 (MQ=37)
                                                              g.cgcTGCCCCC.CTCTTGGCCTTCGGCCAATTCACCCCCCGGG.ATATTTCGGGCAAGATAAAGAAG.GTCTTTTCCTTTCATCTTGCCCG.AAATATCCCgggg                                        >  1:1121439/1‑101 (MQ=37)
                                                                    gCCCCC.CTCTTGGCCTTCGGCCAATTCTCCCCCCGGG.ATATTTCGGTCAAGATAAAGAAG.GTCTTTTCCTTTCATCTTTCCCG.AAATATCCCGGGGGAGcc                                   <  1:709533/101‑1 (MQ=25)
                                                                          cctcttggccTTCGGCCAATTCACCCCCCGGG.ATATTTCGGGCAAGATAAAGAAG.GTCTTTTCCTTTCATCTTGCCCG.AAATATCCCGGGGGAGCCCGCtt                              <  2:3833724/92‑1 (MQ=255)
                                                                             cttggccTTCGGCCAATTCACCCCCCGGG.ATATTTCGGGCAAGATAAAGAAG.GTCTTTTCCTTTCATCTTGCCCG.AAATATCCCGGGGGAGCCCGCTTgcg                           >  2:2497216/7‑101 (MQ=255)
                                                                             cttggccTTCGGCCAATTCACCCCCCGGG.ATATTTCGGGCAAGATAAAGAAG.GTCTTTTCCTTTCATCTTGCCCG.AAATATCCCGGGGGAGCCCGCTTgcg                           >  1:3232992/7‑101 (MQ=255)
                                                                               tggccTTCGGCCAATTCACCCCCCGGG.ATAATTCGGGCAAGATAAAGAAG.GTCTTTTCCTTTCATCTTGCCCG.AGATATCCCGGGGGAGCCCGCTTgcggg                         >  2:1612481/5‑101 (MQ=255)
                                                                                ggccTTCGGCCAATTCACCCCCCGGG.ATATTTCGGGCAAGATAAAGAAG.GTCTTTTCCTTTCATCTTGCCCG.AAATATCCCGGGGGAGCCCGCTTgcgggc                        >  1:2694411/4‑101 (MQ=255)
                                                                                ggccTTCGGCCAATTCACCCCCCGGG.ATATTTCGGGCAAGATAAAGAAG.GTCTTTTCCTTTCATCTTGCCCG.AAATATCCCGGGGGAGCCCGCTTgcgggc                        >  2:2538902/4‑101 (MQ=255)
                                                                                   cTTCGGCCAATTCACCCCCCGGG.ATATTTCGGGCAAGATAAAGAAG.GTCTTTTCCTTTCATCTTGCCCG.AAATATCCCGGGGGAGCCCGCTTGcgggcggg                     >  2:4518154/1‑101 (MQ=255)
                                                                                      cGGCCAATTCACCCCCCGGG.ATATTTCGGGCAAGATAAAGAAG.GTCTTTTCCTTTCATCTTGCCCG.AAATATCCCGGGGGAGCCCGCTTGCGGGCGGGGGc                  <  1:4162774/101‑1 (MQ=255)
                                                                                      cGGCCAATTCACCCCCCGGG.ATATTTCGGGCAAGATAAAGAAG.GTCTTTTCCTTTCATCTTGCCCG.AAATATCCCGGGGGAGCCCGCTTGCGGGCGGGGGc                  >  1:2183399/1‑101 (MQ=255)
                                                                                       ggCCAATTCACCCCCCGGG.ATATTTCGGGCAAGATAAAGAAG.GTCTTTTCCTTTCATCTTGCCCG.AAATATCCCGGGGGAGCCCGCTTGCGGGCGGGGGCa                 >  2:1382929/1‑101 (MQ=255)
                                                                                        ggcAAATCCCCCCCCGGG.GTTTTTTCGGGAAAAAAAAAAAG.GGTTTTTTCTTTCCTTTTTCCCGAAAATTCCCCGGGGGAGCCCGCTTGCGGGCGGGGGCa                 <  1:1020567/99‑1 (MQ=255)
                                                                                        gCCAATTCCCCCCCCGGG.ATATTTCGGGCAAGATAAAGAAG.GTCTTTTCCTTTCATCTTGCCCG.AAATATCCCGGGGGAGCCCGCTTGCGGGCGGGGGCAg                <  2:4848740/101‑1 (MQ=255)
                                                                                            aTTCACCCCCCGGG.ATATTTCGGGCAAGATAAAGAAG.GTCTTTTCCTTTCATCTTGCCCG.AAATATCCCGGGGGAGCCCGCTTGCGGGCGGGGGGAGGGcc            >  1:2966624/1‑101 (MQ=255)
                                                                                            aTTCACCCCCCGGG.ATATTTCGGGCAAGATAAAGAAG.GTCTTTTCCTTTCATCTTGCCCG.AAATATCCCGGGGGAGCCCGCTTGCGGGCGGGGGCAGCGcc            >  1:1842656/1‑101 (MQ=255)
                                                                                            aTTCACCCCCCGGG.ATATTTCGGGCAAGATAAAGAAG.GGCTTTTTCTTTCATCTTGCCCG.AAATATCCCGGGGGAGCCCGCTTGCGGGCGGGGGCAGCGcc            >  2:2525011/1‑101 (MQ=255)
                                                                                              tCACCCCCCGGC.ATATTTCGGGCAAGATAAAGAAG.GTCTTTTCCTTTCAGCTTGCCCG.AAATATCCCGGGGGAGCCCGCTTGCGGGCGGGGGCAGCGcccc          <  2:3964271/101‑1 (MQ=255)
                                                                                               cACCCCCCGGG.ATATTTCGGGCAAGATAAAGAAG.GTCTTTTCCTTTCATCTTGCCCG.AAATATCCCGGGGGAGCCCGCTTGCGGGCGGGGGCAGCGccccc         >  1:3361587/1‑101 (MQ=255)
                                                                                               cACCCCCCGGG.ATATTTCGGGCAAGATAAAGAAG.GTCTTTTCCTTTCATCTTGCCCG.AAATATCCCGGGGGAGCCCGCTTGCGGGCGGGGGCAGCGccccc         >  2:2158328/1‑101 (MQ=255)
                                                                                                gccccccGGG.ATATTTCGGGCAAGATAAAGAAG.GTCTTTTCCTTTCATCTTGCCCG.AAATATCCCGGGGGAGCCCGCTTGCGGGCGGGGGCAGCGCCCCCg        <  1:554365/100‑1 (MQ=255)
                                                                                                aCCCCCCGGG.ATATTTCGGGCAAGATAAAGAAG.GTCTTTTCCTTTCATCTTGCCCG.AAATATCCCGGGGGAGCCCGCTTGCGGGCGGGGGCAGCGCCCCCg        <  2:3232741/101‑1 (MQ=255)
                                                                                                      cGGG.ATATTTCGGGCAAGATAAAGAAG.GTCTTTTCCTTTCATCTTGCCCG.AAATATCCCGGGGGAGCCCGCTTGCGGGCGGGGGCAGCGCCCCCGCAgcgc  >  2:753876/1‑101 (MQ=255)
                                                                                                      |                                                                                                       
CTGGCCCGAGGGGTTCTGCCGCCGCGACATCGGCTGGCGCGCGCTCTGAGACGAAAC.GGGGG.CGCTGCCCCC.CTCCTG..CTTCGGCCAATTCACCCCCCGGG.ATATTTCGGGCAAGATAAAGAAG.GTCTTTTCCTTTCATCTTGCCCG.AAATATCCCGGGGGAGCCCGCTTGCGGGCGGGGGCAGCGCCCCCGAAGCGC  >  NC_014034/2566480‑2566677

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