Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_014034 3,258,963 +C 100% coding (778/1971 nt) efeU → ferrous iron permease EfeU

Read alignment evidence...
  seq id position change freq score reads annotation genes product
*NC_0140343,258,9631.→C100.0% 111.3 317coding (778/1971 nt)efeUferrous iron permease EfeU
Reads supporting (aligned to +/- strand):  new base (152/129):  ref base (0/36):  total (152/165)
Fisher's exact test strand distribution p-value = 8.16e-12
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 3.33e-03

Only 100 of 347 total aligned reads displayed.

CAATGCCTCGGTCGCGGGCACCGTGCTTTCGGCGCTCTTTGGCTATCACGACAGCCCGACCCAGGGCGAGGTCGCCCTCTGGGCCGCGTTCCTGGCCGT.G.A.CCTGCCGCTCTTCCTGCGCCCCGCGCCCGCCCCCCTTCCGCAAAAGGCCTGATCCGATGACCCCTTCCGCACAAGCGCCGTCCGGGCGCCTCCTT  >  NC_014034/3258863‑3259058
                                                                                                       |                                                                                               
 aaTGCCTCGGTCGCGGGCACCGTGCTTTCGGCGCTCTTTGGCTATCACGACAGCCCGACCCAGGGGGCGGTCGCCCTCTGGGCCCGCTTCCTGGCCTT.G.A.c                                                                                                >  1:2290831/1‑101 (MQ=255)
  aTGCCTCGGTCGCGGGCACCGTGCTTTCGGCGCTCTTTGGCTATCACGACAGCCCGACCCAGGGCGAGGTCGCCCTCTGGGCCGCGTTCCTGGCCGT.G.A.cc                                                                                               <  1:5914282/101‑1 (MQ=255)
    gCCTCGGTCGCGGGCACCGTGCTTTCGGCGCTCTTTGGCTATCACGACAGCCCGACCCAGGGCGAGGTCGCCCTCTGGGCCGCGTTCCTGGCCGT.G.A.ccct                                                                                             <  1:3809787/101‑3 (MQ=255)
       tCGGTCGCGGGCACCGTGCTTTCGGCGCTCTTTGGCTATCACGACAGCCCGACCCAGGGCGAGGTCGCCCTCTGGGCCGCGTTCCTGGCCGT.G.A.ccctgcc                                                                                          >  1:3591189/1‑96 (MQ=255)
       tCGGTCGCGGGCACCGTGCTTTCGGCGCTCTTTGGCTATCACGACAGCCCGACCCAGGGCGAGGTCGCCCTCTGGGCCGCGTTCCTGGCCGT.G.A.ccctgcc                                                                                          >  1:2227665/1‑96 (MQ=255)
        cGGTCGCGGGCACCGTGCTTTCGGCGCTCTTTGGCTATCACGACAGCCCGACCCAGGGCGAGGTCGCCCTCTGGGCCGCGTTCCTGGCCGT.G.ACCCTGCCg                                                                                          >  2:1966881/1‑101 (MQ=255)
         ggTCGCGGGCACCGTGCTTTCGGCGCTCTTTGGCTATCACGACAGCCCGACCCAGGGCGAGGTCGCCCTCTGGGCCGCGTTCCTGGCCGT.G.ACCCTGCCGc                                                                                         >  2:6504064/1‑101 (MQ=255)
         ggTCGCGGGCACCGTGCTTTCGGCGCTCTTTGGCTATCACGACAGCCCGACCCAGGGCGAGGGCGCACTCTGGGCCGCGTTCCTGGCCCT.C.A.ccctgcccc                                                                                        >  1:262500/1‑94 (MQ=255)
           tCGAGGGCACCTTGCCTCCGCCGCGGTTTGGCTAGAACGACAGCGCGTCCCAGGGGGAGGT..CCCTGTGGGCCGCGTTCCTGGCCGTGG.C.CCTGCCGCTCtt                                                                                     <  2:205493/101‑1 (MQ=255)
            cgcgGGCACCGTACTTTCGGCGCTCTTTGGCTATCACGACAGCCCGACCCAGGGCGAGGTCGCCCTCTGGGCCGCGTTCCTGGCCGT.G.ACCCTGCCGctct                                                                                      >  2:6366713/1‑101 (MQ=255)
             gcgGGCACCGTGCTTTCGGCGCTCTTTGGCTATCACGACAGCCCGACCCAGGGCGAGGTCGCCCTCTGGGCCGCGTTCCTGGCCGT.G.ACCCTGCCGctctc                                                                                     >  2:2270001/1‑100 (MQ=255)
             gcgGGCACCGTGCTTTCGGCGCTCTTTGGCTATCACGACAGCCCGACCCAGGGCGAGGTCGCCCTCTGGGCCGCGTTCCTGGCCGT.G.ACCCTGCCGCTCtt                                                                                     >  2:5913698/1‑101 (MQ=255)
              cgGGCACCGTGCTTTCGGCGCTCTTTGGCTATCACGACAGCCCGACCCAGGGCGAGGTCGCCCTCTGGGCCGCGTTCCTGGCCGT.G.ACCCTGCCGCTCTTc                                                                                    >  2:4410806/1‑101 (MQ=255)
              cgGGCACCGTGCTTTCGGCGCTCTTTGGCTATCACGACAGCCCGACCCAGGGCGAGGTCGCCCTCTGGGCCGCGTTCCTGGCCGT.G.ACCCTGCCGCTCTTc                                                                                    <  2:1702337/101‑1 (MQ=255)
              cgGGCACCGTGCTTTCGGCGCTCTTTGGCTATCACGACAGCCCGACCCAGGGCGAGGTCGCCCTCTGGGCCGCGTTCCTGGCCGT.G.ACCCTGCCGCTCTTc                                                                                    >  1:900502/1‑101 (MQ=255)
               gggCACCGTGCTTTCGGCGCTCTTTGGCTATCACGACAGCCCGACCCAGGGCGAGGTCGCCCTCTGGGCCGCGTTCCTGGCCGT.G.ACCCTGCCGCTCTTcc                                                                                   >  1:6611260/1‑101 (MQ=255)
               gggCACCGTGCTTTCGGCGCTCTTTGGCTATCACGACAGCCCGACCCAGGGCGAGGTCGCCCCCTGGGCCGCGTTCCTGGCCGT.G.ACCCTGCCGCTCTTcc                                                                                   >  2:3228989/1‑101 (MQ=255)
                ggCCCCGTGCTTTCGGCGCTCTTTGGCTATCCCGACAGCCCGACCCAGGGCGAGGTCGCCCTCTGGGCCGCGTTCCTGGCCGT.G.ACCCTGCCGCTCTTCCt                                                                                  <  2:1098014/101‑1 (MQ=255)
                ggCCCCGTGCTTTCGGCGCTCTTTGGCTATCACGACAGCCCGACCCAGGGCGAGGTCGCCCTCGGGGCCGCGTTCCTGGCCGT.G.ACCCTGCCGCTCTTCCt                                                                                  <  2:3818229/101‑1 (MQ=255)
                ggCACCGTGCTTTCGGCGCTCTTTGGCTATCACGACAGCCCGACCCAGGGCGAGGTCGCCCTCTGGGCCGCGTTCCTGGCCGT.G.ACCCTGCCGCTCTTCCt                                                                                  <  2:4723434/101‑1 (MQ=255)
                 gccccGTTCTTTCGGCTCTCTTTGGCTATCACGACAACCCGAGCCGGGGCGAGGTCGCCCTCTGGGCCGCGTTCCTGGCCGT.G.ACCCTGCCGCTCTTCCTg                                                                                 <  2:3902388/98‑1 (MQ=255)
                   aCCGTGCTTTCGGCGCTCTTTGGCTATCACGACAGCCCGACCCAGGGCGAGGTCGCCCTCTGGGCCGCGTTCCTGGCCGT.G.ACCCTGCCGCTCTTCCTgcg                                                                               >  1:2437719/1‑101 (MQ=255)
                    ccGTGCTTTCGGCGCTCTTTGGCTATCACGACAGCCCGACCCAGGGCGAGGTCGCCCTCTGGGCCGCGTTCCTGGCCGT.G.ACCCTGCCGCTCTTCCTgcgc                                                                              <  1:5686867/101‑1 (MQ=255)
                      gTGCTTTCGGCGCTCTTTGGCTATCACGACAGCCCGACCCAGGGCGAGGTCGCCCTCTGGGCCGCGTTCCTGGCCGT.G.ACCCTGCCGCTCTTCCTGCGccc                                                                            <  1:5718610/101‑1 (MQ=255)
                       tGCTTTCGGCGCTCTTTGGCTATCACGACAGCCCGACCCAGGGCGAGGTCGCCCTCTGGGCCGCGTTCCTGGCCGT.G.ACCCTGCCGCTCTTCCTGCGcccc                                                                           <  1:6593999/101‑1 (MQ=255)
                        gCTTTCGGCGCTCTTTGGCTATCACGACAGCCCGACCCAGGGCGAGGTCGCCCTCTGGGCCGCGTTCCTGGCCGT.G.ACCCTGCCGCTCTTCCTGCGCCCcg                                                                          <  2:928889/101‑1 (MQ=255)
                             cGGCGCTCTTTGGCTATCACGACAGCCCGACCCAGGGCGAGGTCGCCCTCTGGGCCGCGTTCCTGGCCGT.G.ACCCTGCCGCTCTTCCTGCGCCCCGcgccc                                                                     <  1:6816963/101‑1 (MQ=255)
                             cGGCGCTCTTTGGCTATCACGACAGCCCGACCCAGGGCGAGGTCGCCCTCTGGGCCGCGTTCCTGGCCGT.G.ACCCTGCCGCTCTTCCTGCGCCCCGcgccc                                                                     <  1:6254874/101‑1 (MQ=255)
                              ggCGCTCTTTGGCTATCACGACAGCCCGACCCAGGGCGAGGTCGCCCTCTGGGCCGCGTTCCTGGCCGT.G.ACCCTGCCGCTCTTCCTGCGCCCCGcgcccg                                                                    <  1:4453732/101‑1 (MQ=255)
                               gcgcTCTTTGGCTATCACGACAGCCCGACCCAGGGCGAGGTCGCCCTCTGGGCCGCGTTCCTGGCCGT.G.ACCCTGCCGCTCTTCCTGCGCCCCGcgcccgc                                                                   <  1:125461/101‑1 (MQ=255)
                                 gcTCTTTGGCTATCACGACAGCCCGACCCAGGGCGAGGTCGCCCTCTGGGCCGCGTTCCTGGCCGT.G.ACCCTGCCGCTCTTCCTGCGCCCCGCgcccgccc                                                                 >  2:6581308/1‑101 (MQ=255)
                                  ctctTTGGCTACCACGACAGCCCGACCCAGGGCGAGGTCGCCCTCTGGGCCGCGTTCCTGGCCGT.G.ACCCTGCCGCTCTTCCTGCGCCCCGCGCCCGcccc                                                                <  2:4991066/101‑1 (MQ=255)
                                    ctTTGGCTATCACGACAGCCCGACCCAGGGCGAGGTCGCCCTCTGGGCCGCGTTCCTGGCCGT.G.ACCCTGCCGCTCTTCCTGCGCCCCGCGCCCGcccccc                                                              >  1:1150571/1‑101 (MQ=255)
                                     tttGGCTATCACGACAGCCCGACCCAGGGCGAGGTCGCCCTCTGGGCCGCGTTCCTGGCCGT.G.ACCCTGCCGCTCTTCCTGCGCCCCGCGCCCGCCCCCCt                                                             <  1:3214752/101‑1 (MQ=255)
                                      ttGGCTATCACGACAGCCCGACCCAGGGCGAGGTCGCCCTCTGGGCCGCGTTCCTGGCCGT.G.ACCCTGCCGCTCTTCCTGCGCCCCGCGCCCGCCCCCCtt                                                            >  1:3734320/1‑101 (MQ=255)
                                      ttGGCTATCACGACAGCCCGACCCAGGGCGAGGTCGCCCTCTGGGCCGCGTTCCTGGCCGT.G.ACCCTGCCGCTCTTCCTGCGCCCCGCGCCCGCCCCCCtt                                                            <  1:4278303/101‑1 (MQ=255)
                                      ttGGCTACCACGACAGCCCGACCCAGGGCGAGGTCGCCCTCTGGGCCGCGTTCCTGGCCGT.G.ACCCTGCCGCTCTTCCTGCGCCCCGCGCCCGCCCCCCtt                                                            <  2:2614845/101‑1 (MQ=255)
                                       tGGCTATCACGACAGCCCGCCCCAGGGCGAGGTCGCCCTCTGGGCCGCGTTCCTGGCCGT.GCC.CCTGCCGCTCTTCCTGCGCCCCGCGCCCGCCCCCCTTc                                                           <  1:3330921/101‑1 (MQ=255)
                                       tGGCTATCACGACAGCCCGACCCAGGGCGAGGTCGCCCTCTGGGCCGCGTTCCTGGCCGT.G.ACCCTGCCGCTCTTCCTGCGCCCCGCGCCCGCCCCCCTTc                                                           <  2:231696/101‑1 (MQ=255)
                                       tGGCTATCACGACAGCCCGACCCAGGGCGAGGTCGCCCTCTGGGCCGCGTTCCTGGCCGT.G.ACCCTGCCGCTCTTCCTGCGCCCCGCGCCCGCCCCCCTTc                                                           <  1:7045364/101‑1 (MQ=255)
                                       tGGCTATCACGACAGCCCGACCCAGGGCGAGGTCGCCCTCTGGGCCGCGTTCCTGGCCGT.G.ACCCTGCCGCTCTTCCTGCGCCCCGCGCCCGCCCCCCTTc                                                           >  2:4097906/1‑101 (MQ=255)
                                       tGGCTATCACGACAGCCCGACCCAGGGCGAGGTCGCCCTCTGGGCCGCGTTCCTGGCCGT.G.ACCCTGCCGCTCTTCCTGCGCCCCGCGCCCGCCCCCCTTc                                                           <  2:4027010/101‑1 (MQ=255)
                                        ggCTATCACGACAGCCCGACCCAGGGCGAGGTCGCCCTCTGGGCCGCGTTCCTGGCCGT.G.ACCCTGCCGCTCTTCCTGCGCCCCGCGCCCGCCCCCCTTcc                                                          >  1:1581519/1‑101 (MQ=255)
                                         gCTATCACGACAGCCCGACCCAGGGCGAGGTCGCCCTCTGGGCCGCGTTCCTGGCCGT.G.ACCCTGCCGCTCTTCCTGCGCCCCGCGCCCGCCCCCCTTCCg                                                         <  2:680325/101‑1 (MQ=255)
                                          cTATCACGACAGCCCGACCCAGGGCGAGGTCGCCCTCTGGGCCGCGTTCCTGGCCGT.G.ACCCTGCCGCTCTTCCTGCGCCCCGCGCCCGCCCCCCTTCCGc                                                        >  2:5153507/1‑101 (MQ=255)
                                                cGACAGCCCGACCCAGGGCGAGGTCGCCCTCTGGGCCGCGTTCCTGGCCGT.G.ACCCTGCCGCACTTCCTGCGCCCCGCGCCCGCCCCCCTTCCGCAAAAgg                                                  >  2:1018268/1‑101 (MQ=255)
                                                  aCAGCCCGACCCAGGGCGAGGTCGCCCTCTGGGCCGCGTTCCTGGCCGT.G.ACCCTGCCGCTCTTCCTGCGCCCCGCGCCCGCCCCCCTTCCGCAAAAggac                                                >  1:4812896/1‑99 (MQ=255)
                                                  aCAGCCCGACCCAGGGCGAGGTCGCCCTCTGGGCCGCGTTCCTGGCCGT.G.ACCCTGCCGCTCTTCCTGCGCCCCGCGCCCGCCCCCCTTCCGCAAAAGGcc                                                <  1:1901490/101‑1 (MQ=255)
                                                  aCAGCCCGACCCAGGGCGAGGTCGCCCTCTGGGCCGCGTTCCTGGCCGT.G.ACCCTGCCGCTCTTCCTGCGCCCCGCGCCCGCCCCCCTTCCGCAAAAGGcc                                                <  2:6992443/101‑1 (MQ=255)
                                                   cAGCCCGACCCAGGGCGAGGTCGCCCTCTGGGCCGCGTTCCTGGCCGT.G.ACCCTGCCGCTCTTCCTGCGCCCCGCGCCCGCCCCCCTTCCGCAAAAGGCCt                                               <  2:616462/101‑1 (MQ=255)
                                                    aGCCCGACCCAGGGCGAGGTCGCCCTCTGGGCCGCGTTCCTGGCCGT.G.ACCCTGCCGCTCTTCCTGCGCCCCGCGCCCGCCCCCCTTCCGCAAAAGGCCTg                                              <  2:5824208/101‑1 (MQ=255)
                                                      cccGACCCAGGGCGAGGTCGCCCTCTGGGCCGCGTTCCTGGCCGT.G.ACCCTGCCGCTCTTCCTGCGCCCCGCGCCCGCCCCCCTTCCGCAAAAGGCCTGAt                                            <  2:3306452/101‑1 (MQ=255)
                                                      cccGACCCAGGGCGAGGTCGCCCTCTGGGCCGCGTTCCTGGCCGT.G.ACCCTGCCGCTCTTCCTGCGCCCCGCGCCCGCCCCCCTTCCGCAAAAGGCCTGAt                                            >  1:1009814/1‑101 (MQ=255)
                                                      cccGACCCAGGGCGAGGTCGCCCTCTGGGCCGCGTTCCTGGCCGT.G.ACCCTGCCGCTCTTCCTGCGCCCCGCGCCCGCCCCCCTTCCGCAAAAGGCCTGAt                                            <  1:4120643/101‑1 (MQ=255)
                                                       ccGACCCAGGGCGAGGTCGCCCTCTGGGCCGCGTTCCTGGCCGT.GCC.CCTGCCGCTCTTCCTGCGCCCCGCGCCCGCCCCCCTTCCGCAAAAGGCCTGATc                                           <  2:2945897/101‑1 (MQ=255)
                                                       ccGACCCAGGGCGAGGTCACCCTCTGGGCCGCGTTCCTGGCCGT.G.ACCCTGCCGCTCTTCCTGCGCCCCGCGCCCGCCCCCCTTCCGCAAAAGGCCTGATc                                           >  1:2897220/1‑101 (MQ=255)
                                                        cGACCCAGGGCGAGGTCGCCCTCTGGGCCGCGTTCCTGGCCGT.G.ACCCTGCCGCTCTTCCTGCGCCCCGCGCCCGCCCCCCTTCCGCAAAAGGCCTGATcc                                          <  1:5527355/101‑1 (MQ=255)
                                                         gACCCAGGGCGAGGTCGCCCTCTGGGCCGCGTTCCTGGCCGT.GCC.CCTGCCGCTCTTCCTGCGCCCCGCGCCCGCCCCCCTTCCGCAAAAGGCCTGATCCg                                         <  1:5315609/101‑1 (MQ=255)
                                                             cggggCGAGGTCGCCCTCTGGGCCGCGTTCCTGGCCGT.G.ACCCTGCCGCTCTTCCTGCGCCCCGCGCCCGCCCCCCTTCCGCAAAAGGCCTGATCCGATGa                                     <  2:1170415/99‑1 (MQ=255)
                                                              aGGGCGAGGTCGCCCTCTGGGCCGCGTTCCTGGCCGT.G.ACCCTGCCGCTCTTCCTGCGCCCCGCGCCCGCCCCCCTTCCGCAAAAGGCCTGATCCGATGAc                                    >  1:5246479/1‑101 (MQ=255)
                                                               gggCGAGGTCGCCCTCTGGGCCGCGTTCCTGGCCGT.G.ACCCTGCCGCTCTTCCTGCGCCCCGCGCCCGCCCCCCTTCCGCAAAAGGCCTGATCCGATGAcc                                   <  1:6104454/101‑1 (MQ=255)
                                                               gggCGAGGTCGCCCTCGGGGCCGCGTTCCTGGCCGT.G.ACCCTGCCGCTCTTCCTGCGCCCCGCGCCCGCCCCCCTTCCGCAAAAGGCCTGATCCGATGAcc                                   <  2:1563527/101‑1 (MQ=255)
                                                                ggCGAGCTCCCCCTCTGGGCCGCGTTCCTGGCCGT.GCC.CCTGCCGCTCTTCCTGCGCCCCGCGCCCGTCCCCCTTCCGCAAAAGGCCTGATCCGATGAccc                                  <  1:6375042/101‑1 (MQ=255)
                                                                 gCGAGGTCGCCCTCTGGGCCGCGTTCCTGGCCGT.G.ACCCTGCCGCTCTTCCTGCGCCCCGCGCCCGCCCCCCTTCCGCAAAAGGCCTGATCCGATGAcccc                                 >  1:2869856/1‑101 (MQ=255)
                                                                 gCGAGGTCGCCCTCTGGGCCGCGTGCCTGGCCGT.G.ACCCTGCCGCTCTTCCTGCGCCCCGCGCCCGCCCCCCTTCCGCAAAAGGCCTGATCCGATGAcccc                                 <  2:5645389/101‑1 (MQ=255)
                                                                  cGAGGTCGCCCTCTGGGCCGCGTTCCTGGCCGT.G.ACCCTGCCGCTCTTCCTGCGCCCCGCGCCCGCCCCCCTTCCGCAAAAGGCCTGATCCGATGACCCCt                                <  2:1362125/101‑1 (MQ=255)
                                                                  cGAGGTCGCCCTCTGGGCCGCGTTCCTGGCCGT.G.ACCCTGCCGCTCTTCCTGCGCCCCGCGCCCGCCCCCCTTCCGCAAAAGGCCTGATCCGATGACCCCt                                <  2:2890290/101‑1 (MQ=255)
                                                                   gAGGTCGCCCTCTGGGCCGCGTTCCTGGCCGT.G.ACCCTGCCGCTCTTCCTGCGCCCCGCGCCCGCCCCCCTTCCGCAAAAGGCCTGATCCGATGACCCCtt                               >  2:5034729/1‑101 (MQ=255)
                                                                    aGGTCGCCCTCTGGGCCGCGTTCCTGGCCGT.G.ACCCTGCCGCTCTTCCTGCGCCCCGCGCCCGCCCCCCTTCCGCAAAAGGCCTGATCCGATGACCCCTTc                              >  2:4259139/1‑101 (MQ=255)
                                                                    aGGTCGCCCTCTGGGCCGCGTTCCTGGCCGT.G.ACCCTGCCGCTCTTCCTGCGCCCCGCGCCCGCCCCCCTTCCGCAAAAGGCCTGATCCGATGACCCCTTc                              <  1:2669595/101‑1 (MQ=255)
                                                                    aGGTCGCCCTCTGGGCCGCGTCCCTGGCCGT.G.ACCCTGCCGCTCTTCCTGCGCCCCGCGCCCGCCCCCCTTCCGCAAAAGGCCTGATCCGATGACCCCTTc                              <  2:4425238/101‑1 (MQ=255)
                                                                    aGGTCGCCCTCGGGGCCGCGTTCCTGGCCGT.G.ACCCTGCCGCTCTTCCTGCGCCCCGCGCCCGCCCCCCTTCCGCAAAAGGCCTGATCCGATGACCCCTTc                              <  2:5100487/101‑1 (MQ=255)
                                                                     ggTCGCCCTCTGGGCCGCGTTCCTGGCCGT.G.ACCCTGCCGCTCTTCCTGCGCCCCGCGCCCGCCCCCCTTCCGCAAAAGGCCTGATCCGATGACCCCTTcc                             <  1:1064400/101‑1 (MQ=255)
                                                                      gTCGCCCTCTGGGCCGCGTTCCTGGCCGT.G.ACCCTGCCGCTCTTCCTGCGCCCCGCGCCCGCCCCCCTTCCGCAAAAGGCCTGATCCGATGACCCCTTCCg                            <  2:1826697/101‑1 (MQ=255)
                                                                       tCGCCCTCTGGGCCGCGTTCCTGGCCGT.GCC.CCTGCCGCTCTTCCTGCGCCCCGCGCCCGCCCCCCTTCCGCAAAAGGCCTGATCCGATGACCCCTTCCGc                           <  1:2966881/101‑1 (MQ=255)
                                                                       tCGCCCTCTGGGCCGCGTTCCTGGCCGT.G.ACCCTGCCGCTCTTCCTGCGCCCCGCGCCCGCCCCCCTTCCGCAAAAGGCCTGATCCGATGACCCCTTCCGc                           >  2:3689958/1‑101 (MQ=255)
                                                                        cGCCCTCTGGGCCGCGTTCCTGGCCGT.G.ACCCTGCCGCTCTTCCTGCGCCCCGCGCCCGCCCCCCTTCCGCAAAAGGCCTGATCCGATGACCCCTTCCGcc                          >  1:4380664/1‑100 (MQ=255)
                                                                         gCCCTCTGGGCCGCGTTCCTGGCCGT.G.ACCCTGCCGCTCTTCCTGCGCCCCGCGCCCGCCCCCCTTCCGCAAAAGGCCTGATCCGATGACCCCTTCCGcac                         >  2:4227541/1‑101 (MQ=255)
                                                                         gCCCTCTGGGCCGCGTTCCTGGCCGT.G.ACCCTGCCGCTCTTCCTGCGCCCCGCGCCCGCCCCCCTTCCGCAAAAGGCCTGATCCGATGACCCCTTCCGcac                         >  2:2178997/1‑101 (MQ=255)
                                                                           ccTCTGGGCCGCGTTCCTGGCCGT.GCC.CCTGCCGCTCTTCCTGCGCCCCGCGCCCGCCCCCCTTCCGCAAAAGGCCTGATCCGATGACCCCTTCCGCACaa                       <  1:793697/101‑1 (MQ=255)
                                                                           ccTCTGGGCCGCGTTCCTGGCCGT.G.ACCCTGCCGCTCTTCCTGCGCCCCGCGCCCGCCCCCCTTCCGCAAAAGGCCTGATCCGATGACCCCTTCCGCACaa                       >  2:545487/1‑101 (MQ=255)
                                                                            ctctGGGCCGCGTTCCTGGCCGT.G.ACCCTGCCGCTCTTCCTGCGCCCCGCGCCCGCCCCCCTTCCGCAAAAGGCCTGATCCGATGACCCCTTCCGCACAAg                      <  1:2359414/101‑1 (MQ=255)
                                                                             tctGGGCCGCGTTCCTGGCCGT.GCC.CCTGCCGCTGTTCCTGCGCCCCGCGCCCGCCCCCCTTCCGCAAAAGGCCTGATCCGATGACCCCTTCCGCACAAgc                     <  2:4583400/101‑1 (MQ=255)
                                                                             ccggggCCGCGTTCCTGGCCGT.G.ACCCTGCCGCTCTTCCTGCGCCCCGCGCCCGCCCCCCTTCCGCAAAAGGCCTGATCCGATGACCCCTTCCGCACAAgc                     <  2:3108864/98‑1 (MQ=255)
                                                                              ctGGGCCGCGTTCCTGCCCGT.GCC.CCTGCCGCTCTTCCTGCGCCCCGCGCCCGCCCCCCTTCCGCAAAAGGCCTGATCCGATGACCCCTTCCGCACAAgcg                    <  1:6900697/101‑1 (MQ=255)
                                                                               tGGGCCGCGTTCCTGGCCGT.G.ACCCTGCCGCTCTTCCTGCGCCCCGCGCCCGCCCCCCTTCCGCAAAAGGCCTGATCCGATGACCCCTTCCGCACAAgcgc                   >  1:1376359/1‑101 (MQ=255)
                                                                                gggCCGCGTTCCTGGCCGT.G.ACCCTGCCGCTCTTCCTGCGCCCCGCGCCCGCCCCCCTTCCGCAAAAGGCCTGATCCGATGACCCCTTCCGCACAAGCGcc                  <  2:6047956/101‑1 (MQ=255)
                                                                                 ggCCGCGTTCCTGGCCGT.G.ACCCTGCCGCTCTTCCTGCGCCCCGCGCCCGCCCCCCTTCCGCAAAAGGCCTGATCCGATGACCCCTTCCGCACAAGCGCCg                 >  1:4169435/1‑101 (MQ=255)
                                                                                   ccGCGTTCCTGGCCGT.G.ACCCTGCCGCTCTTCCTGCGCCCCGCGCCCGCCCCCCTTCCGCAAAAGGCCTGATCCGATGACCCCTTCCGCACAAGCGCCGTc               <  1:1966953/101‑1 (MQ=255)
                                                                                   ccGCGTTCCTGCCCGT.G.ACCCTGCCGCTCTTCCTGCGCCCCGCGCCCGCCCCCCTTCCGCAAAAGGCCTGATCCGATGACCCCTTCCGCACAAGCGCCGTc               <  2:2721950/101‑1 (MQ=255)
                                                                                    cgcgTTCCTGGCCGT.G.ACCCTGCCGCTCTTCCTGCGCCCCGCGCCCGCCCCCCTTCCGCAAAAGGCCTGATCCGATGACCCCTTCCGCACAAGCGCCGTcc              <  1:6206357/101‑1 (MQ=255)
                                                                                      cgatcttGGCCGT.GCC.CCTGCCGCTCTTCCTGCGCCCCTCGCCCGCCCCCCTTCCGCAAAAGGCCTGATCCGATGACCCCTTCCGCACAAGCGCCGTCCgg            <  1:5094743/95‑1 (MQ=255)
                                                                                       gTTCCTGGCCGT.G.ACCCTGCCGCTCTTCCTGCGCCCCGCGCCCGCCCCCCTTCCGCAAAAGGCCTGATCCGATGACCCCTTCCGCACAAGCGCCGTCCggg           >  1:4535451/1‑101 (MQ=255)
                                                                                         tCCTGGCCGT.G.ACCCTGCCGCTCTTCCTGCGCCCCGCGCCCGCCCCCCTTCCGCAAAAGGCCTGATCCGATGACCCCTTCCGCACAAGCGCCGTCCGGgcg         >  1:304709/1‑101 (MQ=255)
                                                                                           cTGGCCGT.G.ACCCTGCCGCTCTTCCTGCGCCCCGCGCCCGCCCCCGTTCCGCAAAAGGCCTGATCCGATGACCCCTTCCGCACAAGCGCCGTCCGGGCGcc       >  2:4864062/1‑101 (MQ=255)
                                                                                           cTGGCCGT.G.ACCCTGCCGCTCTTCCTGCGCCCCGCGCCCGCCCCCCTTCCGCAAAAGGCCTGATCCGATGACCCCTTCCGCACAAGCGCCGTCCGGGCGcc       >  1:4874016/1‑101 (MQ=255)
                                                                                              gCCGT.G.ACCCTGCCGCTCTTCCTGCGCCCCGCGCCCGCCCCCCTTCCGCAAAAGGCCTGATCCGATGACCCCTTCCGCACAAGCGCCGTCCGGGCGcctcc    >  1:6463040/1‑101 (MQ=255)
                                                                                              gCCGT.G.ACCCTGCCGCTCTTCCTGCGCCCCGCGCCCGCCCCCCTTCCGCAAAAGGCCTGATCCGATGACCCCTTCCGCACAAGCGCCGTCCGGGCGcctcc    >  2:3615533/1‑101 (MQ=255)
                                                                                              gCCGT.G.ACCCTGCCGCTCTTCCTGCGCCCCGCGCCCGCCCCCCTTCCGCAAAAGGCCTGATCCGATGACCCCTTCCGCACAAGCGCCGTCCGGGCGcctcc    >  1:1698189/1‑101 (MQ=255)
                                                                                                cGT.G.ACCCTGCCGCTCTTCCTGCGCCCCGCGCCCGCCCCCCTTCCGCAAAAGGCCTGATCCGATGACCCCTTCCGCACAAGCGCCGTCCGGGCGCCTCCtt  >  1:5046205/1‑101 (MQ=255)
                                                                                                       |                                                                                               
CAATGCCTCGGTCGCGGGCACCGTGCTTTCGGCGCTCTTTGGCTATCACGACAGCCCGACCCAGGGCGAGGTCGCCCTCTGGGCCGCGTTCCTGGCCGT.G.A.CCTGCCGCTCTTCCTGCGCCCCGCGCCCGCCCCCCTTCCGCAAAAGGCCTGATCCGATGACCCCTTCCGCACAAGCGCCGTCCGGGCGCCTCCTT  >  NC_014034/3258863‑3259058

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