Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_014034 2,566,561 C→T 2.0% intergenic (+33/+161) purD → / ← sdaA phosphoribosylamine‑‑glycine ligase/L‑serine ammonia‑lyase

Read alignment evidence...
  seq id position change freq score reads annotation genes product
*NC_0140342,566,5610C→T 2.0% 2.7 ‑1.5 ‑0.8203intergenic (+33/+161)purD/sdaAphosphoribosylamine‑‑glycine ligase/L‑serine ammonia‑lyase
Reads supporting (aligned to +/- strand):  new base (2/2):  ref base (43/156):  total (45/158)
Fisher's exact test strand distribution p-value = 2.14e-01
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 1.52e-01

Only 100 of 224 total aligned reads displayed.

CGATGCGATCGACTGGCCCGAGGGGTTCTGCCGCCGCGACATCGGCTGGCGCGCGCTCTGAGACGAAAC.GG.GGG.CGCTGCCCCC.CTCCT..G..CTTCGGCCAATTCACCCCCCGGG.ATATTTCGGGCAAGATAAAGAAG.GTCTTTTCCTTTCATC.TTGCCCGAAATATCCCGGGGGAGCCCGCTTGCGGGCGGGGGC  >  NC_014034/2566468‑2566661
                                                                                                     |                                                                                                       
cGATGCGATCGACTGGCCCGAGGGGTTCTGCCGCCGCGACATCGGCTGGCGCGCGCTCTGAGACGAAACGGG.GGG.CGCTGCCCCC.CTCTT..GGCCTTCGGcc                                                                                                     <  1:936316/101‑1 (MQ=37)
cGATGCGATCGACTGGCCCGAGGGGTTCTGCCGCCGCGACATCGGCTGGCGCGCGCTCTGAGACGAAACGGG.GGG.CGCTGCCCCC.CTCTT..GGCCTTCGGcc                                                                                                     <  1:3815222/101‑1 (MQ=37)
     cGATCGACTGGCCCGAGGGGTTCTGCCGCCGCGACATCGGCTGGCGCGCGCTCTGAGACGAAACGGG.GGG.CGCTGCCCCC.CTCTT..GGCCTTCGGCCAATTc                                                                                                <  1:1152670/101‑1 (MQ=37)
      ggtCGACTGGCCCGAGGGGTTCTGCCGCCGCGACATCGGCTGGCGCGCGCTCTGAGACGAAACGGG.GGG.CGCTGCCCCC.CTCTT..GGCCTTCGGCCAATTCa                                                                                               <  2:1842999/99‑1 (MQ=25)
      gATCGACTGGCCCGAGGGGTTCTGCCGCCGCGACATCGGCTGGCGCGCGCTCTGAGACGAAACGGG.GGG.CGCTGCCCCC.CTCTT..GGCCTTCGGCCAATTCa                                                                                               <  2:3126198/101‑1 (MQ=37)
        tCGACTGGCCCGAGGGGTTCTGCCGCCGCGACATCGGCTGGCGCGCGCTCTGAGACGAAACGGG.GGG.CGCTGCCCCC.CTCTT..GGCCTTCGGCCAATTCAcc                                                                                             <  2:3286119/101‑1 (MQ=37)
          gACTGTCCCGAGGGGTTCTGCCGCCGCGACATCGGCTGGCGCGCGCTCTGAGACGAAACGGG.GGG.CGCTGCCCCC.CTCTT..GGCCTTCGGCCAATTCAcccc                                                                                           <  2:1103212/101‑1 (MQ=25)
          gACTGGCCCGAGGGGTTCTGCCGCCGCGACATCGGCTGGCGCGCGCTCTGAGACGAAACGGG.GGG.CGCTGCCCCC.CTCTTGGG..CCTCGGGCACTTCAcccc                                                                                           >  2:546213/1‑101 (MQ=255)
           aCTGGCCCGAGGGGTTCTGCCGCCGCGACATCGGCTGGCGCGCGCTCTGAGACGAAACGGG.GGG.CGCTGCCCCC.CTCTT..GGCCTTCGGCCAATTCAccccc                                                                                          <  2:4887723/101‑1 (MQ=25)
             tGGCCCGAGGGGTTCTGCCGCCGCGACATCGGCTGGCGCGCGCTCTGAGACGAAACGGG.GGG.CGCTGCCCCC.CTCTT..GGCCTTCGGCCAATTCACCCCCCg                                                                                        <  2:2580723/101‑1 (MQ=25)
             tGGCCCGAGGGGTTCTGCCGCCGCGACATCGGCTGGCGCGCGCTCTGAGACGAAACGGG.GGG.CGCTGCCCCC.CTCTT..GGCCTTCGGCCAATTCACCCCCCg                                                                                        <  2:4553621/101‑1 (MQ=25)
             tGGCCCGAGGGGTTCTGCCGCCGCGACATCGGCTGGCGCGCGCTCTGAGACGAAACGGG.GGG.CGCTGCCCCC.CCCCT..GGCCCTCCGCCGATTCACCCCCCg                                                                                        >  1:3059487/1‑101 (MQ=255)
               gCCCGAGGGGTTCTTACGGCGCGGCAACGGCTGGCGCGCGCTCTGAGAAGAAAC.GGAGGG.CGCTGCGCCC.CTCTT..GGCCTTCGGCCAATTCACCCCCCggg                                                                                      <  2:3343907/101‑1 (MQ=16)
               gCCCGAGGGGTTCTGCCGCCGCGACATCGGCTGTCGCGCGCTCTGAGACGAAACGGG.GGG.CGCTGCCCCC.CTCTT..GGCCTTCGGCCAATTCACCCCCCggg                                                                                      <  1:4309561/101‑1 (MQ=25)
                  cGAGGGGTTCTGCCGCCGCGACATCGGCTGGCGCGCGCTCTGAGACGAAACGGG.GGG.CGCTGCCCCC.CTCTT..GGCCTTCGGCCAATTCACCCCCCGGG.ata                                                                                  <  1:2547606/101‑1 (MQ=25)
                    aGGGGTTCTGCCGCCGCGACATCGGCTGGCGCGCGCTCTGAGACGAAACGGG.GGG.CGCTGCCCCC.CTCTT..GGCCTTCGGCCAATTCACCCCCCGGG.ATAtt                                                                                <  1:15296/101‑1 (MQ=21)
                    aGGGGTTCTGCCGCCGCGACATCGGCTGGCGCGCGCTCTGAGACGAAACGGG.GGG.CGCTGCCCCC.CTCTT..GGCCTTCGGCCAATTCACCCCCCGGG.ATAtt                                                                                <  1:5046351/101‑1 (MQ=21)
                     ggggTTCTGCCGCCGCGACATCGGCTGGCGCGCGCTCTGAGACGAAACGGG.GGG.CGCTGCCCCC.CTCTT..GGCCTTCGGCCAATTCACCCCCCGGG.ATAttt                                                                               <  1:4823843/101‑1 (MQ=21)
                      gggTTCTGCCGCCGCGACATCGGCTGGCGCGCGCTCTGAGACGAAACGGG.GGG.CGCTGCCCCC.CTCTT..GGCCTTCGGCCAATTCACCCCCCGGG.ATATTTc                                                                              <  2:3172939/101‑1 (MQ=21)
                        gTTCTGCCGCCGCGACATCGGCTGGCGCGCGCTCTGAGACGAAACGGG.GGG.CGCTGCCCCC.CTCTT..GGCCTTCGGCCCAATCACCCCCCGGG.ATATTTTgg                                                                            >  2:4745212/1‑101 (MQ=21)
                         ttCTGCCGCCGCGACATCGGCTGGCGCGCGCTCTGAGACGAAACGGG.GGG.CGCTGCCCCC.CTCTT..GGCCTTCGGCCAATTCACCCCCCGGG.ATATTTCggg                                                                           >  2:2644884/1‑101 (MQ=21)
                             gcctccGCGACATCGGCTGGCGCGCGTTCTGAGACGAAACGGG.GGG.CGCTGCCCCC.CTCTT..GGCCTTCGGCCAATTCACCCCCCGGG.ATATTTCGGGCAAg                                                                       <  2:558382/101‑1 (MQ=18)
                              ccgccgCGACATCGGCTGGCGCGCGCTCTGAGACGAAACGGG.GGG.CGCTGCCCCC.CTCTT..GGCCTTCGGCCAATTCACCCCCCGGG.ATATTTCGGGCAAGa                                                                      <  1:3738925/101‑1 (MQ=21)
                               cgccgcGACATCGGCTGGCGCGCGCTCTGAGACGAAACGGG.GGG.CGCTGCCCCC.CTCTT..GGCCTTCGGCCAATTCACCCCCCGGG.ATATTTCGGGCAAGAt                                                                     <  1:3134645/101‑1 (MQ=21)
                                gccgcGACATCGGCTGGCGCGCGCTCTGAGACGAAACGGG.GGG.CGCTGCCCCC.CTCTT..GGCCTTCGGCCAATTCACCCCCCGGG.ATATTTCGGGCAAGATa                                                                    <  1:758907/101‑1 (MQ=21)
                                gccgcGACATCGGCTGGCGCGCGCTCTGAGACGAAACGGG.GGG.CGCTGCCCCC.CTCTT..GGCCTTCGGCCAATTCACCCCCCGGG.ATATTTCGGGCAAGATa                                                                    <  2:150726/101‑1 (MQ=21)
                                gccgcGACATCGGCTGGCGCGCGCTCTGAGACGAAACGGG.GGG.CGCTGCCCCC.CTCTT..GGCCTTCGGCCAATTCACCCCCCGGG.ATATTTCGGGCAAGATa                                                                    <  2:2373027/101‑1 (MQ=21)
                                 ccgcGACATCGGCTGGCGCGCGCTCTGAGACGAAACGGG.GGG.GGGCGCCCCC.CTCCT..GGCCTTTCGCCCATTCACCCCCCGGG.ATACTTCGGGCGAGACaa                                                                   >  2:812208/1‑101 (MQ=16)
                                 ccgcGACATCGGCTGGCGCGCGCTCTGAGACGAAACGGG.GGG.CGCTGCCCCC.CTCTT..GGCCTTCGGCCAATTCACCCCCCGGG.ATATTTCGGGCAAGATaa                                                                   <  2:1821522/101‑1 (MQ=21)
                                 ccgcGACATCGGCTGGCGCGCGCTCTGAGACGAAACGGG.GGG.CGCTGCCCCC.CTCTT..GGCCTTCGGCCAATTCACCCCCCGGG.ATATTTCGGGCAAGATaa                                                                   <  2:2102822/101‑1 (MQ=21)
                                 ccgcGACATCGGCTGGCGCGCGCTCTGAGACGAAACGGG.GGG.CGCTGCCCCC.CTCTT..GGCCTTCGGCCAATTCACCCCCCGGG.ATATTTCGGGCAAGATaa                                                                   <  2:2169155/101‑1 (MQ=21)
                                 ccgcGACATCGGCTGGCGCGCGATCTGAGACGAAACGGG.GGG.CGCTGCCCCC.CTCTT..GGCCTTCGGCCAATTCACCCCCCGGG.ATATTTCGGGCAAGATaa                                                                   <  1:2766840/101‑1 (MQ=21)
                                  cgcgACATCGGCTGGCGCGCGCTCTGAGACGAAACGGG.GGG.CGCTGCCCCC.CTCTT..GGCCTTCGGCCAATTCACCCCCCGGG.ATATTTCGGGCAAGATaaa                                                                  <  1:4025861/101‑1 (MQ=21)
                                   gtgACTTCGGCGGGCGCGCGCTCTGAGACAAAACGGG.GGG.CGCTGCCCCC.CTCTT..GGCCTTCGGCCAATTCCCCCCCCGGG.ATATTTCGGGCAAGATAaag                                                                 <  1:2206603/99‑1 (MQ=18)
                                   gcgACATCGGCTGGCGCGCGCTCTGAGACGAAACGGG.GGG.CGCTGCCCCC.CTCTT..GGCCTTCGGCCAATTCACCCCCCGGG.ATATTTCGGGCAAGATAaag                                                                 <  1:3001494/101‑1 (MQ=21)
                                    cgACATCGGCTGGCGCGCGCTCTGAGACGAAACGGG.GGG.CGCTGCCCCC.CTCTT..GGCCTTCGGCCAATTCACCCCCCGGG.ATATTTCGGGCAAGATAaaga                                                                <  2:2991831/101‑1 (MQ=21)
                                     gACATCGGCTGGCGCTCGCTCTGAGACGAAACGGG.GGG.CGCTGCCCCC.CTCTT..GGCCTTCGGGCCAATCACCCCCCCGG.GTAATTCGGGCCAGATAaaaaa                                                               >  1:2354133/1‑101 (MQ=12)
                                     gACATCGGCTGGCGCGCGCTCTGAGACGAAACGGG.GGG.CGCTGCCCCC.CTCTT..GGCCTTCGGCCAATTCAACCCCCGGG.ATATTTCGGGCAAGATAaagaa                                                               <  1:1375187/101‑1 (MQ=21)
                                      aCATCGGCTGGCGCGCGCTCTGAGACGAAACGGG.GGG.CGCTGCCCCC.CTCTT..GGCCTTCGGCCAATTCACCCCCCGGG.ATATTACGGGCAAGATAaagaag                                                              <  2:2525570/101‑1 (MQ=21)
                                         tCGGCTGGCGCGCGCTCTGAGACGAAACGGG.GGG.CGCTGCCCCC.CTCTT..GGCCTTCGGCCAATTCACCCCCCGGG.ATATTTCGGGCAAGATAAAGAAG.GTc                                                          <  2:2440729/101‑1 (MQ=21)
                                          cGGCTGGCGCGCGCTCTGAGACGAAACGGG.GGG.CGCTGCCCCC.CTCTT..GGCCTTCGGCCAATTCACCCCCCGGG.ATATTTCGGGCAAGATAAAGAAG.GTCt                                                         <  1:535576/101‑1 (MQ=21)
                                           ggCTGGCGCGCGCTCTGAGACGAAACGGG.GGG.CGCTGCCCCC.CTCTT..GGCCTTCGGCCAATTCACCCCCCGGG.ATATTTCGGGCAAGATAAAGAAG.GTCtt                                                        <  1:182174/101‑1 (MQ=21)
                                           ggCTGGCGCGCGCTCTGAGACGAAACGGG.GGG.CGCTGCCCCC.CTCTT..GGCCTTCGGCCAATTCACCCCCCGGG.ATATTTCGGGCAAGATAAAGAAG.GTCtt                                                        <  2:3402325/101‑1 (MQ=21)
                                           ggCTGGCGCGCCCTCTGAGACGAAACGGG.GGG.CGCTGCCCCC.CTCTT..GGCCTTCGGCCAATTCACCCCCCGGG.ATATTTCGGGCAAGATAAAGAAG.GTCtt                                                        <  2:4794865/101‑1 (MQ=21)
                                             cTGGCGCGCGCTCTGAGACGAAACGGG.GGG.CGCTGCCCCC.CTCTT..GGCCTTCGGCCAATTCACCCCCCGGG.ATATTTCGGGCAAGATAAAGAAG.GTCtttt                                                      <  2:174030/101‑1 (MQ=21)
                                             cTGGCGCGCGCTCTGAGACGAAACGGG.GGG.CGCTGCCCCC.CTCTT..GGCCTTCGGCCAATTCACCCCCCGGG.ATATTTCGGGCAAGATAAAGAAG.GTCtttt                                                      <  1:1041996/101‑1 (MQ=21)
                                              tGGCGCGCGCTCTGAGACGAAACGGG.GGG.CGCTGCCCCC.CTCTT..GGCCTTCGGCCAATTCACCCCCCGGG.ATATTTCGGGCAAGATAAAGAAG.GTCTTTTc                                                     <  2:3994153/101‑1 (MQ=21)
                                              tGGCGCGCGCTCTGAGACGAAACGGG.GGG.CGCTGCCCCC.CTCTT..GGCCTTCGGCCAATTCACCCCCCGGG.ATATTTCGGGCAAGATAAAGAAG.GTCTTTTc                                                     <  2:5020865/101‑1 (MQ=21)
                                               ggCGCGCGCTCTGAGACGAAACGGG.GGG.CGCTGCCCCC.CTCTT..GGCCTTCGGCCAATTCACCCCCCGGG.ATATTTCGGGCAAGATAAAGAAG.GTCTTTTcc                                                    >  1:4956338/1‑101 (MQ=21)
                                                gcgcgcgcTCTGAGACGAAACGGG.GGG.CGCTGCCCCC.CTCTT..GGCCTTCGGCCAATTCACCCCCCGGG.ATATTTCGGGCAAGATAAAGAAG.GTCTTTTCCt                                                   <  2:2861298/101‑1 (MQ=21)
                                                 cgcgcgcTCTGAGACGAAACGGG.GGG.CGCTGCCCCC.CTCTT..GGCCTTCGGCCAATTCACCCCCCGGG.ATATTTCGGGCAAGATAAAGAAG.GTCTTTTCCtt                                                  <  1:2093234/101‑1 (MQ=21)
                                                  gcgcgcTCTGAGACGAAACGGG.GGG.CGCTGCCCCC.CTCTT..GGCCTTCGGCCAATTCACCCCCCGGG.ATATTTCGGGCAAGATAAAGAAG.GTCTTTTCCttt                                                 <  1:4426888/101‑1 (MQ=21)
                                                  gcgcgcTCTGAGACGAAACGGG.GGG.CGCTGCCCCC.CTCTT..GGCCTTCGGCCAATTCACCCCCCGGG.ATATTTCGGGCAAGATAAAGAAG.GTCTTTTCCttt                                                 <  1:1700501/101‑1 (MQ=21)
                                                   cgcgcTCTGAGACGAAACGGG.GGG.CGCTGCCCCC.CTCTT..GGCCTTCGGCCAATTCACCCCCCGGG.ATATTTCGGGCAAGATAAAGAAG.GTCTTTTCCTTTc                                                <  2:1918172/101‑1 (MQ=21)
                                                   cgcgcTCTGAGACGAAACGGG.GGG.CGCTGCCCCC.CTCTT..GGCCTTCGGCCAATTCACCCCCCGGG.ATATTTCGGGCAAGATAAAGAAG.GTCTTTTCCTTTc                                                <  2:406436/101‑1 (MQ=21)
                                                    gcgcTCTGAGACGAAACGGG.GGG.CGCTGCCCCC.CTCTT..GGCCTTCGGCCAATTCACCCCCCGGG.ATATTTCGGGCAAGATAAAGAAG.GTCTTTTCCTTTCa                                               <  1:2952744/101‑1 (MQ=21)
                                                    gcgcTCTGAGACGAAACGGG.GGG.CGCTGCCCCC.CTCTT..GGCCTTCGGCCAATTCACCCCCCGGG.ATATTTCGGGCAAGATAAAGAAG.GTCTTTTCCTTTCa                                               <  2:426189/101‑1 (MQ=21)
                                                     cgcTCTGAGACGAAACGGG.GGG.CGCTGCCCCC.CTCTT..GGCCTTCGGCCAATTCACCCCCCGGG.ATATTTCGGGCAAGATAAAGAAG.GTCTTTTCCTTTCAt                                              <  1:3208980/101‑1 (MQ=21)
                                                     cgcTCTGAGACGAAAC.GG.GGGTCGCTTCCCCC.CTCTT..GGCCTTCGGCCAATTCACCCCCCGGG.ATATTTCGGGCAAGATAAAGAAG.GTCTTTTCCTTTCAt                                              <  1:1616734/101‑1 (MQ=25)
                                                      gcTCTGAGACGAAACGGG.GGG.CGCTGCCCCC.CTCTT..GGCCTTCGGCCAATTCACCCCCCGGG.ATATTTCGGGCAAGATAAAGAAG.GTCTTTTCCTTTCATc                                             <  1:2154159/101‑1 (MQ=21)
                                                      gcTCTGAGACGAAACGGG.GGG.CGCTGCCCCC.CTCTT..GGCCTTCGGCCAATTCACCCCCCGGG.ATATTTCGGGCAAGATAAAGAAG.GTCTTTTCCTTTCATc                                             <  1:2565977/101‑1 (MQ=21)
                                                      gcTCTGAGACGAAACGGG.GGG.CGCTGCCCCC.CTCTT..GGCCTTCGGCCAATTCACCCCCCGGG.ATATTTCGGGCAAGATAAAGAAG.GTCTTTTCCTTTCATc                                             <  1:3191357/101‑1 (MQ=21)
                                                       ctctGAGACGAAACGGG.GGG.CGCTGCCCCC.CTCTT..GGCCTTCGGCCAATTCACCCCCCGGG.ATATTTCGGGCAAGATAAAGAAG.GTCTTTTCCTTTCATC.t                                           >  2:2387314/1‑101 (MQ=21)
                                                       ctctGAGACGAAACGGG.GGG.CGCTGCCCCC.CTCTT..GGCCTTCGGCCAATTCACCCCCCGGG.ATATTTCGGGCAAGATAAAGAAG.GTCTTTTCCTTTCATC.t                                           >  2:4678419/1‑101 (MQ=21)
                                                       ctctGAGACGAAACGGG.GGG.CGCTGCCCCC.CTCTT..GGCCTTCGGCCAATTCACCCCCCGGG.ATATTTCGGGCAAGATAAAGAAG.GTCTTTTCCTTTCATC.t                                           <  2:3411039/101‑1 (MQ=21)
                                                       ctctGAGACGAAACGGG.GGG.CGCTGCCCCC.CTCTT..GGCCTTCGGCCAATTCACCCCCCGGG.ATATTTCAGGCAAGATAAAGAAG.GTCTTTTCCTTTCATC.t                                           <  2:1093383/101‑1 (MQ=21)
                                                        tctGAGACGAAACGGG.GGG.CGCTGCCCCC.CTCTT..GGCCTTCGGCCAATTCACCCCCCGGG.ATATTTCGGGCAAGATAAAGAAG.GTCTTTTCCTTTCATC.tt                                          >  2:612440/1‑101 (MQ=21)
                                                        tctGAGACGAAACGGG.GGG.CGCTGCCCCC.CTCTT..GGCCTTCGGCCAATTCACCCCCCGGG.ATATTTCGGGCAAGATAAAGAAG.GTCTTTTCCTTTCATC.tt                                          <  1:88760/101‑1 (MQ=21)
                                                         ctGAGACGAAACGGG.GGG.CGCTGCCCCC.CTCTT..GGCCTTCGGCCAATTCACCCCCCGGG.ATATTTCGGGCAAGATAAAGAAG.GTCTTTTCCTTTCATC.TTg                                         <  1:582768/101‑1 (MQ=21)
                                                         ctGAGACGAAACGGG.GGG.CGCTGCCCCC.CTCTT..GGCCTTCGGCCAATTCACCCCCCGGG.ATATTTCGGGCAAGATAAAGAAG.GTCTTATCCTTTCATC.TTg                                         <  2:1974418/101‑1 (MQ=21)
                                                          tGAGACGAAACGGG.GGG.CGCTGCCCCC.CTCTT..GGCCTTCGGCCAATTCACCCCCCGGG.ATATTTCGGGCAAGATAAAGAAG.GTCTTTTCCTTTCATC.TTGc                                        <  2:1073410/101‑1 (MQ=21)
                                                           gagaCGAAACGGG.GGG.CGCTGCCCCC.CTCTT..GGCCTTCGGCCAATTCACCCCCCGGG.ATATTTCGGGCAAGATAAAGAAG.GTCTTTTCCTTTCATC.TTGcc                                       >  1:4110983/1‑101 (MQ=21)
                                                           gagaCGAAACGGG.GGG.CGCTGCCCCC.CTCTT..GGCCTTCGGCCAATTCACCCCCCCGG.AGATTTTGGGGCAGATAAAAAAG.GGGTTTTTCTTTTATT.ttttc                                       >  1:3572030/1‑98 (MQ=2)
                                                            agaCGAAACGGG.GGG.CGCTGCCCCC.CTCTT..GGCCTTCGGCCAATTCACCCCCCGGG.ATATTTCGGGCAAGATAAAGAAG.GTCTTTTCCTTTCATC.TTGccc                                      <  1:1748694/101‑1 (MQ=21)
                                                            agaCGAAACGGG.GGG.CGCTGCCCCC.CTCTT..GGCCTTCGGCCAATTCACCCCCCGGG.ATATTTCGGGCAAGATAAAGAAG.GTCTTTTCCTTTCATC.TTGccc                                      <  1:191171/101‑1 (MQ=21)
                                                               cGAAACGGG.GGG.CGCTGCCCCC.CTCTT..GGCCTTCGGCCAATTCACCCCCCGGG.ATATTTCGGGCAAGATAAAGAAG.GTCTTTTCCTTTCATC.TTGCCCGaa                                   >  1:3357414/1‑101 (MQ=21)
                                                               cGAAACGGG.GGG.CGCTGCCCCC.CTCTT..GGCCTTCGGCCAATTCACCCCCCGGG.ACATTTCGGGCAAGATAAAGAAG.GTCTTTTCCTTTCATC.TTGCCCGaa                                   >  2:613597/1‑101 (MQ=21)
                                                                 aaacggg.ggg.CGCTGCCCCC.CTCTT..GGCCTTCGGCCAATTCACCCCCCGGG.ATATTTCGGGCAAGATAAAGAAG.GTCTTTTCCTTTCATC.TTGCCCGAAat                                 <  1:3998860/96‑1 (MQ=21)
                                                                 aaacggg.ggg.CGCTGCCCCC.CTCTT..GGCCTTCGGCCAATTCACCCCCCGGG.AAATTCCGGCCAAAAAAAGAAAG.GACTTTTCTTTTCTTC.TTCCCCGaaaa                                 >  2:4359411/6‑100 (MQ=2)
                                                                   acggg.ggg.CGCTGCCCCC.CTCTT..GGCCTTCGGCCAATTCTCCCCCCGGG.ATATTTCGGGCAAGATAAAGAAG.GTCTTTTCCTTTCATC.TTGCCCGAAatat                               <  2:3198269/98‑1 (MQ=25)
                                                                   acggg.ggg.CGCTGCCCCC.CTCTT..GGCCTTCGGCCAATTCACCCCCCGGG.ATATTTCGGGCAAGATAAAGAAG.GTCTTTTCCTTTCATC.TTGCCCGAAatat                               <  2:908820/98‑1 (MQ=25)
                                                                   acggg.ggg.CGCTGCCCCC.CTCTT..GGCCTTCGGCCAATTCACCCCCCGGG.ATATTTCGGGCAAGATAAAGAAG.GTCTTTTCCTTTCATC.TTGCCCGAAatat                               <  1:2129400/98‑1 (MQ=25)
                                                                    cggg.ggg.CGCTGCCCCC.CTCTT..GGCCTTCGGCCAATTCACCCCCCGGGTAAATTTGGGGC.AGATAAAAAAGAGCTTTTTCTTTTCTTT.TTGCCCCaaaatc                               >  1:226992/3‑98 (MQ=0)
                                                                      gg.ggg.CGCTGCCCCCTCTTGG..C..CTTTGGCCAATTCA.CCCCTGGG.ATATTTCGGGCAAGAT...GAA..ATCATGCCCCTTCATC.TTGCCGCAAATATCCCgggggg                      <  2:1478495/101‑2 (MQ=2)
                                                                       g.ggg.CGCTGCCCCC.CTCTT..GGCCTTCGGCCAATTCACCCCCCGGG.ATATTTCGGGCAAGATAAAGAAG.GTCTTTTCCTTTCATC.TTGCCCGAAATATCCCg                           <  2:2456950/101‑1 (MQ=25)
                                                                       g.ggg.CCCTGCCCCC.TTTTG..GC.CTTTGGCCAATTCA.CCCCTGGG.ATATTTCGGGCAAGAT...GAA..ATCATGCCCCTTCATC.TTGCCGCAAATATCCCgggggga                     <  2:359618/101‑3 (MQ=2)
                                                                                 gCCCCC.CTCTT..GGCCTTCGGCCAATTCTCCCCCCGGG.ATATTTCGGTCAAGATAAAGAAG.GTCTTTTCCTTTCATC.TTTCCCGAAATATCCCGGGGGAGcc                   <  1:709533/101‑1 (MQ=25)
                                                                                 gCCCCC.CTCTT..GGCCTTCGGCCAATTCACCCCCCGGG.ATATTTCGGGCAAGATAAAGAAG.GTCTTTTCCTTTCATC.TTGGCCGAAATATCCCGGGGGAGcc                   >  2:4973558/1‑101 (MQ=25)
                                                                                 gCCCCC.CTCTT..GGCCTTCGGCCAATTCACCCCCCGGG.ATATTTCGGGCAAGATAAAGAAG.GTCTTTTCCTTTCATC.TTGCCCGAAATATCCCGGGGGAGcc                   <  2:4036376/101‑1 (MQ=37)
                                                                                        ccctcttggccTTCGGCCAATTCACCCCCCGGG.ATATTTCGGGCAAGATAAAGAAG.GTCTTTTCCTTTCATC.TTGCCCGAAATATCCCGGGGGAGCCCGCt               <  1:2645069/91‑1 (MQ=255)
                                                                                         cctcttggccTTCGGCCAATTCACCCCCCGGG.ATATTTCGGGCAAGATAAAGAAG.GTCTTTTCCTTTCATC.TTGCCCGAAATATCCCGGGGGAGCCCGCtt              <  2:3833724/92‑1 (MQ=255)
                                                                                         cctcctggccTTCGCCCAATTCCCCCCCCGGG.ATATTTCGGGCAAGATAAAGAAG.GTCTTTTCCTTTCATC.TTGCCCGAAATATCCCGGGGGAGCCCGCtt              <  1:3930720/92‑1 (MQ=37)
                                                                                          ctcttggccTTCGGCCAATTCACCCCCCGGG.ATATTTCGGGCAAGATAAAGAAG.GTCTTTTCCTTTCATC.TTGCCCGAAATATCCCGGGGGAGCCCGCTTg             <  1:3269055/93‑1 (MQ=255)
                                                                                          ctcttggccTTCGGCCAATTCACCCCCCGGG.ATATTTCGGGCAAGATAAAGAAG.GTCTTTTCCTTTCATC.TTGCCCGAAATATCCCGGGGGAGCCCGCTTg             <  2:2205523/93‑1 (MQ=255)
                                                                                              tggccTTCGGCCAATTCCCCCCCCGGG.ATATTTCGGGCAAGATAAAGAAG.GTCTTTTCCTTTCATC.TTGCCCGAAATATCCCGGGGGAGCCCGCTTgcggg         <  2:4504973/97‑1 (MQ=255)
                                                                                               ggccTTCGGCCAATTCACCCCCCGGG.ATATTTCGGGCAAGATAAAGAAG.GTCTTTTCCTTTCATC.TTTCCCGAAATATCCCGGGGGAGCCCGCTTgcgggc        >  1:2724557/4‑101 (MQ=255)
                                                                                                 ccTTCGGCCAATTCACCCCCCGGG.ATAATTTGGGCAAGATAAAAAAA.GGGTTTTTCTTTTATCTTTTCCCCAAAATCCCCGGGGGGGCCCCCTGgcggggg       >  1:4235308/2‑99 (MQ=12)
                                                                                                  cTTCGGCCAATTCACCCCCCGGG.ATATTTCGGGCAAGATAAAGAAG.GTCTTTTCCTTTCATC.TTGCCCGAAATATCCCGGGGGAGCCCGCTTGcgggcggg     >  1:3400121/1‑101 (MQ=255)
                                                                                                  cTTCGGCCAATTCACCCCCCGGG.ATATTTCGGGCAAGATAAAGAAG.GTCTTTTCCTTTCATC.TTGCCCGAAATATCCCGGGGGAGCCCGCTTGcgggcggg     <  1:4678742/101‑1 (MQ=255)
                                                                                                    tCGGCCAATTCACCCCCCGGG.ATATTTCGGGCAAGATAAAGAAG.GTCTTTTCCTTTCATC.TTGCCCGAAATATCCCGGGGGAGCCCGCTTGCGGGCggggg   >  1:2669553/1‑101 (MQ=255)
                                                                                                     |                                                                                                       
CGATGCGATCGACTGGCCCGAGGGGTTCTGCCGCCGCGACATCGGCTGGCGCGCGCTCTGAGACGAAAC.GG.GGG.CGCTGCCCCC.CTCCT..G..CTTCGGCCAATTCACCCCCCGGG.ATATTTCGGGCAAGATAAAGAAG.GTCTTTTCCTTTCATC.TTGCCCGAAATATCCCGGGGGAGCCCGCTTGCGGGCGGGGGC  >  NC_014034/2566468‑2566661

Base quality scores:  ATCG  < 3 ≤  ATCG  < 15 ≤  ATCG  < 26 ≤  ATCG  < 33 ≤  ATCG  < 40 ≤  ATCG