Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_014034 3,254,687 Δ1 bp 100% intergenic (+12/‑171) glpX → / → recJ fructose‑1,6‑bisphosphatase/single‑stranded‑DNA‑specific exonuclease RecJ

Read alignment evidence...
  seq id position change freq score reads annotation genes product
*NC_0140343,254,6870G→.100.0% 75.2 167intergenic (+12/‑171)glpX/recJfructose‑1,6‑bisphosphatase/single‑stranded‑DNA‑specific exonuclease RecJ
Reads supporting (aligned to +/- strand):  new base (3/139):  ref base (24/0):  total (28/139)
Fisher's exact test strand distribution p-value = 5.43e-26
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 1.12e-12

Only 100 of 224 total aligned reads displayed.

ACCCGAAATTCATCGAGACCGAAACGATCCTGATGCGCTCGAAAACCGGCTCGGTGCGCCGGGTGATCTATCGCAACCCGATCCGCTGAGGCAACGCT..CCGGGGGGTATGCCGTCCCCCGCACCCCCTGCGCGTATTTTGCCAAGGAAAACGGGGCCGTTTCACCTTGGCCGAAATATCCCGGGGGGGAATTTG  >  NC_014034/3254587‑3254780
                                                                                                      |                                                                                             
aCCCGAAATTCATCGAGACCGAAACGATCCTGATGCGCTCGAAAACCGGCTCGGTGCGCCGGGTGATCTATCGCAACCCGATCCGCTGAGGCAACGCT..ACg                                                                                               <  1:2366893/101‑1 (MQ=255)
    gAAATTCATCGAGACCGAAACGATCCTGATGCGCTCGAAAACCGGCTCGGTGCGCCGGGTGATCTATCGCAACCCGATCCGCTGAGGCAACGCT..ACggggg                                                                                           <  2:2334154/101‑1 (MQ=255)
      aaTTCATCGAGACCGAAACGATCCTGATGCGCTCGAAAACCGGCTCGGTGCGCCGGGTGATCTATCGCAACCCGATCCGCTGAGGCAACGCT..ACtgggggg                                                                                         >  2:4962161/1‑100 (MQ=255)
      aaTTCATCGAGACCGAAACGATCCTGATGCGCTCGAAAACCGGCTCGGTGCGCCGGGTGATCTATCGCAACCCGATCCGCTGAGGCAACGCT..ACgggggac                                                                                         <  1:1130860/101‑3 (MQ=255)
        ttCATCGAGACCGAAACGATCCTGGTGCGCTCGAAAACCGGCTCGGTGCGCCGGGTGATCTATCGCAACCCGATCCGCTGAGGCAACGCT..ACcgggggcgg                                                                                       >  2:1281054/1‑98 (MQ=255)
        ttCATCGAGACCGAAACGATCCTGATGCGCTCGAAAACCGGCTCGGTGCGCCGGGTGATCTATCGCAACCCGATCCGCTGAGGCAACGCT..ACgggggactg                                                                                       <  1:3687331/101‑5 (MQ=255)
         tCATCGAGACCGAAACGATCCTGATGCGCTCGAAAACCGGCTCGGTGCGCCGGGTGATCTATCGCAACCCGATCCGCTGAGGCAACGCT..ACggggggcggc                                                                                      >  2:4206397/1‑97 (MQ=255)
         tCATCGAGACCGAAACGATCCTGATGCGCTCGAAAACCGGCTCGGTGCGCCGGGTGATCTATCGCAACCCGATCCGCTGAGGCAACGCT..ACgggggactgc                                                                                      <  1:3559490/101‑6 (MQ=255)
         tCATCGAGACCGAAACGATCCTGATGCGCTCGAAAACCGGCTCGGTGCGCCGGGTGATCTATCGCAACCCGATCCGCTGAGGCAACGCT..ACgggggactgc                                                                                      <  2:1637883/101‑6 (MQ=255)
          cATCGAGACCGAAACGATCCTGATGCGCTCGAAAACCGGCTCGGTGCGCCGGGTGATCTATCGCAACCCGATCCGCTGAGGCAACGCT..ACggggggggccg                                                                                     >  1:5068418/1‑96 (MQ=255)
          cATCGAGACCGAAACGATCCTGATGCGCTCGAAAACCGGCTCGGTGCGCCGGGTGATCTATCGCAACCCGATCCGCTGAGGCAACGCT..ACgggggactgcc                                                                                     <  1:1150263/101‑7 (MQ=255)
          cATCGAGACCGAAACGATCCTGATGCGCTCGAAAACCGGCTCGGTGCGCCGGGTGATCTATCGCAACCCGATCCGCTGAGGCAACGCT..ACgggggactgcc                                                                                     <  1:5497523/101‑7 (MQ=255)
          cATCGAGACCGAAACGATCCTGATGCGCTCGAAAACCGGCTCGGTGCGCCGGGGGATCTATCGCAACCCGATCCGCTGAGGCAACGCTA.CCggggggcggc                                                                                      >  2:5247062/1‑97 (MQ=255)
           aTCGAGACCGAAACGATCCTGATGCGCTCGAAAACCGGCTCGGTGCGCCGGGTGATCTATCGCAACCCGATCCGCTGAGGCAACGCTACCCgggggggggg                                                                                      >  1:596882/1‑97 (MQ=255)
           aTCGAGACCGAAACGATCCTGATGCGCTCGAAAACCGGCTCGGTGCGCCGGGTGATCTATCGCAACCCGATCCGCTGAGGCAACGCTA.CCgggggggcggg                                                                                     >  2:405275/1‑96 (MQ=255)
           aTCGAGACCGAAACGATCCTGATGCGCTCGAAAACCGGCTCGGTGCGCCGGGTGATCTATCGCAACCCGATCCGCTGAGGCAACGCT..ACgggggggtgccg                                                                                    >  2:1598607/1‑95 (MQ=255)
           aTCGAGACCGAAACGATCCTGATGCGCTCGAAAACCGGCTCGGTGCGCCGGGTGATCTATCGCAACCCGATCCGCTGAGGCAACGCT..ACgggggactgccg                                                                                    <  1:387228/101‑8 (MQ=255)
            tCGAGACCGAAACGATCCTGATGCGCTCGAAAACCGGCTCGGTGCGCCGGGTGATCTATCGCAACCCGATCCGCTGAGGCAACGCTA.CCgggggggggcgc                                                                                    >  2:4076929/1‑95 (MQ=255)
            tCGAGACCGAAACGATCCTGATGCGCTCGAAAACCGGCTCGGTGCGCCGGGTGATCTATCGCAACCCGATCCGCTGAGGCAACGCTA.CCgggggggggccc                                                                                    >  1:2137635/1‑95 (MQ=255)
            tCGAGACCGAAACGATCCTGATGCGCTCGAAAACCGGCTCGGTGCGCCGGGTGATCTATCGCAACCCGATCCGCTGAGGCAACGCT..ACgggggactgccgt                                                                                   <  2:65727/101‑9 (MQ=255)
             cGAGACCGAAACGATCCTGATGCGCTCGAAAACCGGCTCGGTGCGCCGGGTGATCTATCGCAACCCGATCCGCTGAGGCAACGCT..ACgggggggggccggc                                                                                  >  2:460757/1‑93 (MQ=255)
              gagaCCGAAACGATCCTGATGCGCTCGAAAACCGGCTCGGTGCGCCGGGTGATCTATCGCAACCCGATCCGCTGAGGCAACGCTACCCgggggggggcggc                                                                                   >  1:1109250/1‑94 (MQ=255)
               agaCCGAAACGATCCTGATGCGCTCGAAAACCGGCTCGGTGCGCCGGGTGATCTATCGCAACCCGATCCGCTGAGGCAACGCT..AC.GGGGGACTGCCGTccc                                                                               <  2:2041644/101‑1 (MQ=255)
                gaCCGAAACGATCCTGATGCGCTCGAAAACCGGCTCGGTGCGCCGGGTGATCTATCGCAACCCGATCCGCTGAGGCAACGCT..AC.GGGGGACTGCCGTcccc                                                                              <  1:1714052/101‑1 (MQ=255)
                 aCCGAAACGATCCTGATGCGCTCGAAAACCGGCTCGGTGCGCCGGGTGATCTATCGCAACCCGATCCGCTGAGGCAACGCT..AC.GGGGGACTGCCGTccccc                                                                             <  1:4765207/101‑1 (MQ=255)
                 aCCGAAACGATCCTGATGCGCTCGAAAACCGGCTCGGTGCGCCGGGTGATCTATCGCAACCCGATCCGCTGAGGCAACGCT..AC.GGGGGACTGCCGTccccc                                                                             <  2:860192/101‑1 (MQ=255)
                 aCCGAAACGATCCTGATGCGCTCGAAAACCGGCTCGGTGCGCCGGGTGATCTATCGCAACCCGATCCGCTGAGGCAACGCT..AC.GGGGGACTGCCGTccccc                                                                             <  2:2549772/101‑1 (MQ=255)
                    gAAACGATCCTGATGCGCTCGAAAACCGGCTCGGTGCGCCGGGTGATCTATCGCAACCCGATCCGCTGAGGCAACGCT..ACGGGGGG.CCGCCGTCCCCCGcc                                                                          >  2:5787717/1‑100 (MQ=255)
                    gAAACGATCCTGATGCGCTCGAAAACCGGCTCGGTGCGCCGGGTGATCTATCGCAACCCGATCCGCTGAGGCAACGCT..AC.GGGGGACTGCCGTcccagatc                                                                          >  1:2227933/1‑96 (MQ=255)
                     aaaCGATCCTGATGCGCTCGAAAACCGGCTCGGTGCGCCGGGTGATCTATCGCAACCCGATCCGCTGAGGCAACGCT..ACGGGGGGGGGGCCGCcccccccc                                                                          >  1:3322962/1‑98 (MQ=255)
                      aaCGATCCTGATGCGCTCGAAAACCGGCTCGGTGCGCCGGGTGATCTATCGCAACCCGATCCGCTGAGGCAACGCT..AC.GGGGGACTGCCGTCCCCCGCAcc                                                                        <  1:1777865/101‑1 (MQ=255)
                        cGATCCTGATGCGCTCGAAAACCGGCTCGGTGCGCCGGGTGATCTATCGCAACCCGATCCGCTGAGGCAACGCT..ACCGGGGGGCCGCCGCCCCCCCCCccc                                                                       >  1:4195625/1‑101 (MQ=255)
                        cGATCCTGATGCGCGCGAAAACCGGCTCGGTGCGCCGGGTGATCTATCGCAACCCGATCCGCTGAGGCAACGCT..AC.GGGGGACTGCCGTCCCCCGCAcccc                                                                      <  2:3844312/101‑1 (MQ=255)
                         gATCCTGATGCGCTCGAAAACCGGCTCGGTGCGCCGGGTGATCTATCGCAACCCGATCCGCTGAGGCAACGCT..ACCGGGGGGGGGCCGCCCCCCCCCcccc                                                                      >  2:253551/1‑101 (MQ=255)
                          aTCCTGATGCGCTCGAAAACCGGCTCGGTGCGCCGGGTGATCTATCGCAACCCGATCCGCTGAGGCAACGCT..AC.GGGGGACTGCCGTCCCCCGCAccccca                                                                    <  2:1123913/101‑2 (MQ=255)
                            ccTGATGCGCTCGAAAACCGGCTCGGTGCGCCGGGTGATCTATCGCAACCCGATCCGCTGAGGCAACGCT..ACGGGGGG.CTGGCGGCCCCCGCCcccccccg                                                                  >  2:1662006/1‑98 (MQ=255)
                              tGATGCGCTCGAAAACCGGCTCGGTGCGCCGGGTGATCTATCGCAACCCGATCCGCTGAGGCAACGCT..ACGGGGGG.CCGCCGCCCCCCCCCCCCCCCgccg                                                                >  2:3321893/1‑99 (MQ=255)
                               gATGCGCTCGAAAACCGGCTCGGTGCGCCGGGTGATCTATCGCAACCCGATCCGCTGAGGCAACGCT..AC.GGGGGACTGCCGTCCCCCGCACCCCCTGcgcg                                                               <  1:2971603/101‑1 (MQ=255)
                                 tGCGCTCGAAAACCGGCTCGGTGCGCCGGGTGATCTATCGCAACCCGATCCGCTGAGGCAACGCT..AC.GGGGGACTGCCGTCCCCCGCACCCCCTGCGCGTa                                                             <  2:3363307/101‑1 (MQ=255)
                                  gcgcTCGAAAACCGGCTCGGTGCGCCGGGTGATCTATCGCAACCCGATCCGCTGAGGCAACGCT..ACGGGGGG.CTGCCGTCCCCCCCCCCCCCCCCGCGGAt                                                            >  2:1886523/1‑101 (MQ=255)
                                  gcgcTCGAAAACCGGCTCGGTGCGCCGGGTGATCTATCGCAACCCGATCCGCTGAGGCAACGCT..AC.GGGGGACTGCCGTCCCCCGCACCCCCTGCGCGTAt                                                            <  1:2740884/101‑1 (MQ=255)
                                  gcgcTCGAAAACCGGCTCGGTGCGCCGGGTGATCTATCGCAACCCGATCCGCTGAGGCAACGCT..AC.GGGGGACTGCCGTCCCCCGCACCCCCTGCGCGTAt                                                            <  1:3083259/101‑1 (MQ=255)
                                  gcgcTCGAAAACCGGCTCGGTGCGCCGGGTGATCTATCGCAACCCGATCCGCTGAGGCAACGCT..AC.GGGGGACTGCCGTCCCCCGCACCCCCTGCGCGTAt                                                            <  1:405319/101‑1 (MQ=255)
                                   cgcTCGAAAACCGGCTCGGTGCGCCGGGTGATCTATCGCAACCCGATCCGCTGAGGCAACGCT..AC.GGGGGACTGCCGTCCCCCGCACCCCCTGCGCGTAtt                                                           <  2:4458734/101‑1 (MQ=255)
                                    gcTCGAAAACCGGCTCGGTGCGCCGGGTGATCTATCGCAACCCGATCCGCTGAGGCAACGCT..ACGGGGGGGGGGCCGCCCCCCCCCcccccccccccgtct                                                           >  1:384967/1‑91 (MQ=255)
                                    gcTCGAAAACCGGCTCGGTGCGCCGGGTGATCTATCGCAACCCGATCCGCTGAGGCAACGCT..AC.GGGGGACTGCCGTCCCCCGCACCCCCTGCGCGTAttt                                                          <  2:5025911/101‑1 (MQ=255)
                                    gcTCGAAAACCGGCTCGGTGCGCCGGGCGATCTATCGCAACCCGATCCGCTGAGGCAACGCT..ACGGGGGG.CGGCCGTCCCCCCCCCCCCCCgccccgactt                                                          >  2:3652988/1‑93 (MQ=255)
                                     cTCGAAAACCGGCTCGGTGCGCCGGGTGATCTATCGCAACCCGATCCGCTGAGGCAACGCT..AC.GGGGGACTGCCGTCCCCCGCACCCCCTGCGCGTAtttt                                                         <  2:3420066/101‑1 (MQ=255)
                                     cTCGAAAACCGGCTCGGTGCGCCGGGTGATCTATCGCAACCCGATCCGCTGAGGCAACGCT..AC.GGGGGACTGCCGTCCCCCGCACCCCCTGCGCGTAtttt                                                         <  1:2370579/101‑1 (MQ=255)
                                     cTCGAAAACCGGCTCGGTGCGCCGGGTGATCTATCGCAACCCGATCCGCTGAGGCAACGCT..AC.GGGGGACTGCCGTCCCCCGCACCCCCTGCGCGTAtttt                                                         <  2:5220841/101‑1 (MQ=255)
                                     cTCGAAAACCGGCTCGGTGCGCCGGGTGATCTATCGCAACCCGATCCGCTGAGGCAACGCT..AC.GGGGGACTGCCGTCCCCCGCACCCCCTGCGCGTAtttt                                                         <  2:548766/101‑1 (MQ=255)
                                       cGAAAACCGGCTCGGTGCGCCGGGTGATCTATCGCAACCCGATCCGCTGAGGCAACGCT..AC.GGGGGACTGCCGTCCCCCGCACCCCCTGCGCGTATTTTGc                                                       <  2:4754089/101‑1 (MQ=255)
                                       cGAAAACCGGCTCGGTGCGCCGGGTGATCTATCGCAACCCGATCCGCTGAGGCAACGCT..AC.GGGGGACTGCCGTCCCCCGCACCCCCTGCGCGTATTTTGc                                                       <  2:5431382/101‑1 (MQ=255)
                                        gAAAACCGGCTCGGTGCGCCGGGTGATCTATCGCAACCCGATCCGCTGAGGCAACGCT..AC.GGGGGACTGCCGTCCCCCGCACCCCCTGCGCGTATTTTGcc                                                      >  1:2622808/1‑101 (MQ=255)
                                        gAAAACCGGCTCGGTGCGCCGGGTGATCTATCGCAACCCGATCCGCTGAGGCAACGCT..AC.GGGGGACTGCCGTCCCCCGCACCCCCTGCGCGTATTTTGcc                                                      <  1:2796618/101‑1 (MQ=255)
                                          aaaCCGGCTCGGTGCGCCGGGTGATCTATCGCAACCCGATCCGCTGAGGCAACGCT..AC.GGGGGACTGCCGTCCCCCGCACCCCCTGCGCGTATTTTGCCaa                                                    <  2:3484991/101‑1 (MQ=255)
                                            aCCGGCTCGGTGCGCCGGGTGATCTATCGCAACCCGATCCGCTGAGGCAACGCT..AC.GGGGGACTGCCGCCCCCCGCACCCCCTGCGCGTATTTTGCCAAgg                                                  <  2:1715738/101‑1 (MQ=38)
                                             ccGGCTCGGTGCGCCGGGTGATCTATCGCAACCCGATCCGCTGAGGCAACGCT..AC.GGGGGACTGCCGTCCCCCGCACCCCCTGCGCGTATTTTGCCAAGGa                                                 <  2:3006333/101‑1 (MQ=39)
                                              cGGCTCGGTGCGCCGGGTGATCTATCGCAACCCGATCCGCTGAGGCAACGCT..AC.GGGGGACTGCCGTCCCCCGCACCCCCTGCGCGTATTTTGCCAAGGaa                                                <  2:614838/101‑1 (MQ=38)
                                              cGGCTCGGTGCGCCGGGTGATCTATCGCAACCCGATCCGCTGAGGCAACGCT..AC.GGGGGACTGCCGTCCCCCGCACCCCCTGCGCGTATTTTGCCAAGGaa                                                <  2:5045547/101‑1 (MQ=38)
                                              cGGCTCGGTGCGCCGGGTGATCTATCGCAACCCGATCCGCTGAGGCAACGCT..AC.GGGGGACTGCCGTCCCCCGCACCCCCTGCGCGTATTTTGCCAAGGaa                                                <  1:337685/101‑1 (MQ=38)
                                               ggCTCGGTGCGCCGGGTGATCTATCGCAACCCGATCCGCTGAGGCAACGCT..AC.GGGGGACTGCCGTCCCCCGCACCCCCTGCGCGTATTTTGCCAAGGaaa                                               <  2:476193/101‑1 (MQ=38)
                                                gCTCGGTGCGCCGGGTGATCTATCGCAACCCGATCCGCTGAGGCAACGCT..AC.GGGGGACTGCCGTCCCCCGCACCCCCTGCGCGTATTTTGCCAAGGaaaa                                              <  2:232075/101‑1 (MQ=38)
                                                  gcGGTGCGCCGGTTGATCTACCGCAACCCGATCCGCTGAGGCAACGCT..AC.GGGGGACTGCCGTCCCCCGCACCCCCTGCGCGTATTTTGCCAAGGAAAACg                                            <  1:2153588/100‑1 (MQ=37)
                                                    ggTGCGCCGGGTGATCTATCGCAACCCGATCCGCTGAGGCAACGCT..AC.GGGGGACTGCCGTCCCCCGCACCCCCTGCGCGTATTTTGCCAAGGAAAACggg                                          <  2:2843371/101‑1 (MQ=38)
                                                       gcgcCGGGTGATCTATCGCAACCCGATCCGCTGAGGCAACGCT..AC.GGGGGACTGCCGTCCCCCGCACCCCCTGCGCGTATTTTGCCAAGGAAAACGGGGcc                                       <  1:4494400/101‑1 (MQ=38)
                                                       gcgcCGGGTGATCTATCGCAACCCGATCCGCTGAGGCAACGCT..AC.GGGGGACTGCCGTCCCCCGCACCCCCTGCGCGTATTTTGCCAAGGAAAACGGGGcc                                       <  2:3115378/101‑1 (MQ=38)
                                                       gcgcCGGGTGATCTATCGCAACCCGATCCGCTGAGGCAACGCT..AC.GGGGGACTGCCGTCCCCCGCACCCCCTGCGCGTATTTTGCCAAGGAAAACGGGGcc                                       <  2:733453/101‑1 (MQ=38)
                                                           cGGGTGATCTATCGCAACCCGATCCGCTGAGGCAACGCT..AC.GGGGGACTGCCGTCCCCCGCACCCCCTGCGCGTATTTTGCCAAGGAAAACGGGGCCGttt                                   <  1:4065254/101‑1 (MQ=38)
                                                            gggTGATCTATCGCAACCCGATCCGCTGAGGCAACGCT..AC.GGGGGACTGCCGTCCCCCGCACCCCCTGCGCGTATTTTGCCAAGGAAAACGGGGCCGTTTc                                  <  1:2470363/101‑1 (MQ=38)
                                                            gggTGATCTATCCCAACCCGATGCGCTGAGGCAACGCT..AC.GGGGGACTGCCGTCCCCCGCACCCCCTGCGCGTATTTTGCCAAGGAAAACGGGGCCGTTTc                                  <  2:4380048/101‑1 (MQ=37)
                                                             ggTGATCTATCGCAACCCGATCCGCTGAGGCAACGCT..AC.GGGGGACTGCCGTCCCCCGCACCCCCTGCGCGTATTTTGCCAAGGAAAACGGGGCCGTTTCa                                 <  1:141205/101‑1 (MQ=38)
                                                                 aTCTATCGCAACCCGATCCGCTGAGGCAACGCT..AC.GGGGGACTGCCGTCCCCCGCACCCCCTGCGCGTATTTTGCCAAGGAAAACGGGGCCGTTTCACCtt                             <  2:1448245/101‑1 (MQ=38)
                                                                 aTCTATCGCAACCCGATCCGCTGAGGCAACGCT..AC.GGGGGACTGCCGTCCCCCGCACCCCCTGCGCGTATTTTGCCAAGGAAAACGGGGCCGTTTCACCtt                             <  1:997170/101‑1 (MQ=38)
                                                                 aTCTATCGCAACCCGATCCGCTGAGGCAACGCT..AC.GGGGGACTGCCGTCCCCCGCACCCCCTGCGCGTATTTTGCCAAGGAAAACGGGGCCGTTTCACCtt                             <  1:880784/101‑1 (MQ=38)
                                                                 aTCTATCGCAACCCGATCCGCTGAGGCAACGCT..AC.GGGGGACTGCCGTCCCCCGCACCCCCTGCGCGTATTTTGCCAAGGAAAACGGGGCCGTTTCACCtt                             <  2:3732761/101‑1 (MQ=38)
                                                                  tCTATCGCAACCCGATCCGCTGAGGCAACGCT..AC.GGGGGACTGCCGTCCCCCGCACCCCCTGCGCGTATTTTGCCAAGGAAAACGGGGCCGTTTCACCTTg                            <  2:222343/101‑1 (MQ=38)
                                                                   cTATCGCAACCCGATCCGCTGAGGCAACGCT..AC.GGGGGACTGCCGTCCCCCGCACCCCCTGCGCGTATTTTGCCAAGGAAAACGGGGCCGTTTCACCTTgg                           <  1:930778/101‑1 (MQ=38)
                                                                     aTCGCAACCCGATCCGCTGAGGCAACGCT..AC.GGGGGACTGCCGTCCCCCGCACCCCCTGCGCGTATTTTGCCAAGGAAAACGGGGCCGTTTCACCTTGGcc                         <  2:5359993/101‑1 (MQ=38)
                                                                          aaCCCGATCCGCTGAGGCAACGCT..AC.GGGGGACTGCCGTCCCCCGCACCCCCTGCGCGTATTTTGCCAAGGAAAACGGGGCCGTTTCACCTTGGCCGAAat                    <  1:2410252/101‑1 (MQ=38)
                                                                           aCCCGATCCGCTGAGGCAACGCT..AC.GGGGGACTGCCGTCCCCCGCACCCCCTGCGCGTATTTTGCCAAGGAAAACGGGGCCGTTTCACCTTGGCCGAAata                   <  1:571765/101‑1 (MQ=38)
                                                                           aCCCGATCCGCTGAGGCAACGCT..AC.GGGGGACTGCCGTCCCCCGCACCCCCTGCGCGTATTTTGCCAAGGAAAACGGGGCCGTTTCACCTTGGCCGAAata                   <  2:828378/101‑1 (MQ=38)
                                                                           aCCCGATCCGCTGAGGCAACGCT..AC.GGGGGACTGCCGTCCCCCGCACCCCCTGCGCGTATTTTGCCAAGGAAAACGGGGCCGTTTCACCTTGGCCGAAata                   <  1:4789119/101‑1 (MQ=38)
                                                                           aCCCGATCCGCTGAGGCAACGCT..AC.GGGGGACTGCCGTCCCCCGCACCCCCTGCGCGTATTTTGCCAAGGAAAACGGGGCCGTTTCACCTTGGCCGAAata                   <  2:5106764/101‑1 (MQ=38)
                                                                            cccGATCCGCTGAGGCAACGCT..AC.GGGTGACTGCCGTCCCCCGCACCCCCTGCGCGTATTTTGCCAAGGAAAACGGGGCCGTTTCACCTTGGCCGAAatat                  <  1:1690472/101‑1 (MQ=38)
                                                                            cccGATCCGCTGAGGCAACGCT..AC.GGGGGACTGCCGTCCCCCGCACCCCCTGCGCGTATTTTGCCAAGGAAAACGGGGCCGTTTCACCTTGGCCGAAatat                  <  1:5204979/101‑1 (MQ=38)
                                                                              cGATCCGCTGAGGCAACGCT..AC.GGGGGACTGCCGTCCCCCGCACCCCCTGCGCGTATTTTGCCAAGGAAAACGGGGCCGTTTCACCTTGGCCGAAATATcc                <  1:4302186/101‑1 (MQ=38)
                                                                               gATCCGCTGAGGCAACGCT..ACGGGGGGGGGGCGGCCCCCCCCCCCCCCCCCCCCACTTTTTCCAAGAAAAAAAGGGCGGGTTCCCCTTGCCGGAAAAATcc                >  2:1195259/1‑101 (MQ=255)
                                                                                 tCCGCTGAGGCAACGCT..AC.GGGGGACTGCCGTCCCCCGCACCCCCTGCGCGTATTTTGCCAAGGAAAACGGGGCCGTTTCACCTTGGCCGAAATATCCCgg             <  1:849669/101‑1 (MQ=38)
                                                                                  ccGCTGAGGCAACGCT..ACGGGGGG.CTGCCGTCCCCCGCACCCCCCGCGCGTATTTTGCCAAGGAAAACGGGGCCGGTTCACCTTGGCCGAAATATTCCggg            >  2:4143249/1‑101 (MQ=25)
                                                                                  ccGCTGAGGCAACGCT..AC.GGGGGACTGCCGTCCCCCGCACCCCCTGCGCGTATTTTGCCAAGGAAAACGGGGCCGTTTCACCTTGGCCGAAATATCCCggg            <  1:3137156/101‑1 (MQ=38)
                                                                                     cTGAGGCAACGCT..AC.GGGGGACTGCCGTCCCCCGCACCCCCTGCGCGTATTTTGCCAAGGAAAACGGGGCCGTTTCACCTTGGCCGAAATATCCCgggggg         <  2:3803774/101‑1 (MQ=38)
                                                                                     cTGAGGCAACGCT..AC.GGGGGACTGCCGTCCCCCGCACCCCCTGCGCGTATTTTGCCAAGGAAAACGGGGCCGTTTCACCTTGGCCGAAATATCCCgggggg         <  1:1862143/101‑1 (MQ=38)
                                                                                       gAGGCAACGCT..AC.GGGGGACTGCCGTCCCCCGCACCCCCTGCGCGTATTTTGCCAAGGAAAACGGGGCCGTTTCACCTTGGCCGAAATATCCCGGGGGGGa       <  2:2660640/101‑1 (MQ=38)
                                                                                       gAGGCAACGCT..AC.GGGGGACTGCCGTCCCCCGCACCCCCTGCGCGTATTTTGCCAAGGAAAACGGGGCCGTTTCACCTTGGCCGAAATATCCCGGGGGGGa       <  1:4962442/101‑1 (MQ=38)
                                                                                         ggCAACGCT..AC.GGGGGACTGCCGTCCCCCGCACCCCCTGCGCGTATTTTGCCAAGGAAAACGGGGCCGTTTCACCTTGGCCGAAATATCCCGGGGGGGAAt     <  1:2066951/101‑1 (MQ=38)
                                                                                         ggCAACGCT..AC.GGGGGACTGCCGTCCCCCGCACCCCCTGCGCGTATTTTGCCAAGGAAAACGGGGCCGTTTCACCTTGGCCGAAATATCCCGGGGGGGAAt     <  2:1350410/101‑1 (MQ=38)
                                                                                          gCAACGCT..AC.GGGGGACTGCCGTCCCCCGCACCCCCTGCGCGTATTTTGCCAAGGAAAACGGGGCCGTTTCACCTTGGCCGAAATATCCCGGGGGGGAAtt    <  2:5632654/101‑1 (MQ=38)
                                                                                          gCAACGCT..AC.GGGGGACTGCCGTCCCCCGCACCCCCTGCGCGTATTTTGCCAAGGAAAACGGGGCCGTTTCACCTTGGCCGAAATATCCCGGGGGGGAAtt    <  1:1322052/101‑1 (MQ=38)
                                                                                           cAACGCT..AC.GGGGGACTGCCGTCCCCCGCACCCCCTGCGCGTATTTTGCCAAGGAAAACGGGGCCGTTTCACCTTGGCCAAAATATCCCGGGGGGGAAttt   <  2:2000562/101‑1 (MQ=37)
                                                                                                      |                                                                                             
ACCCGAAATTCATCGAGACCGAAACGATCCTGATGCGCTCGAAAACCGGCTCGGTGCGCCGGGTGATCTATCGCAACCCGATCCGCTGAGGCAACGCT..CCGGGGGGTATGCCGTCCCCCGCACCCCCTGCGCGTATTTTGCCAAGGAAAACGGGGCCGTTTCACCTTGGCCGAAATATCCCGGGGGGGAATTTG  >  NC_014034/3254587‑3254780

Base quality scores:  ATCG  < 3 ≤  ATCG  < 20 ≤  ATCG  < 27 ≤  ATCG  < 33 ≤  ATCG  < 40 ≤  ATCG