Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_014034 2,559,935 Δ1 bp 100% intergenic (‑211/+59) ispZ ← / ← ftsY intracellular septation protein A/cell division protein FtsY

Read alignment evidence...
  seq id position change freq score reads annotation genes product
*NC_0140342,559,9350C→.100.0% 31.7 112intergenic (‑211/+59)ispZ/ftsYintracellular septation protein A/cell division protein FtsY
Reads supporting (aligned to +/- strand):  new base (93/7):  ref base (0/12):  total (93/19)
Fisher's exact test strand distribution p-value = 1.14e-11
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 4.25e-08

Only 100 of 145 total aligned reads displayed.

CGGTTTCCTGACAGGATTGTGTCCTTCGCGCCGCAGCAGGCCCGCTTTTCCTTGGGCGAAACGCCGGGGGGT.T.G.C.GCGTCAGC.GCAAGGGGGGGGCAGCGCCCCC.CC.G.AAGCGCCTGAGGCGCGGCGGGCTTGCCCGCCCAACGAAGGCGCAACCACCTTCAGTCGGCTTCGATCCCCACCAGCGCGGCGGCGAATTCC  >  NC_014034/2559835‑2560033
                                                                                                         |                                                                                                     
cGGTTTCCTGACAGGATTGTGTCCTTCGCGCCGCAGCAGGCCCGCTTTTCCTTGGGCGAAACGCCGGGGGGT.T.G.C.GCGTCAGC.GCAAGCGGGGGGCAgcgc                                                                                                       >  2:2687287/1‑101 (MQ=255)
cGGTTTCCTGACAGGATTGTGTCCTTCGCGCCGCAGCAGGCCCGCTTTTCCTTGGGCGAAACGCCGGGGGGT.T.G.C.GCGTCAGC.GCAAGCGGGGGGCAgcgc                                                                                                       >  1:3906194/1‑101 (MQ=255)
  gTTTCCTGACAGGATTGTGTCCTTCGCGCCGCAGCAGGCCCGCTTTTCCTTGGGCGAAACGCCGGGGGGT.T.G.C.GCGTCAGC.GCAAGCGGGGGGCAGCGccc                                                                                                     >  1:286445/1‑101 (MQ=255)
   tttCCTGACAGGATTGTGTCCTTCGCGCCGCAGCAGGCCCGCTTTTCCTTGGGCGAAACGCCGGGGGGT.T.G.C.GCGTCAGC.GCAAGCGGGGGGCAGCGcccc                                                                                                    >  2:3148481/1‑101 (MQ=255)
    ttCCTGACAGGATTGTGTCCTTCGCGCCGCAGCAGGCCCGCTTTTCCTTGGGCGAAACGCCGGGGGGT.T.G.C.GCGTCAGC.GCAAGCGGGGGGCAGCGccccc                                                                                                   >  1:1705889/1‑101 (MQ=255)
    ttCCTGACAGGATTGTGTCCTTCGCGCCGCAGCAGGCCCGCTTTTCCTTGGGCGAAACGCCGGGGGGT.T.G.C.GCGTCAGC.GCAAGCGGGGGGCAGCGccccc                                                                                                   >  2:220582/1‑101 (MQ=255)
     tCCTGACAGGATTGTGTCCTTCGCGCCGCAGCAGGCCCGCTTTTCCTTGGGCGAAACGCCGGGGGGT.T.G.C.GCGTCAGC.GCAAGCGGGGGGCAGCGccccc.c                                                                                                 >  1:1333639/1‑101 (MQ=255)
     tCCTGACAGGATTGTGTCCTTCGCGCCGCAGCAGGCCCGCTTTTCCTTGGGCGAAACGCCGGGGGGT.T.G.C.GCGTCAGC.GCAAGCGGGGGGCAGCGccccc.c                                                                                                 <  1:2723177/101‑1 (MQ=255)
     tCCTGACAGGATTGTGTCCTTCGCGCCGCAGCAGGCCCGCTTTTCCTTGGGCGAAACGCCCGGGGGT.T.G.C.GCGTCAGC.GCAAGCGGGGGGCAGCGccccc.c                                                                                                 >  1:2005552/1‑101 (MQ=255)
     tCCTGACAGGATTGTGTCCTCCGCGCCGCAGCAGCCCCGCTTTTCCTTGGGCGAACCGCCGGGGGGT.G.G.C.GCGTCAGC.GCAAGCGGGGGGCAGCGccccc.c                                                                                                 <  2:436103/101‑1 (MQ=255)
       cTGACAGGATTGTGTCCTTCGCGCCGCAGCAGGCCCGCTTTTCCTTGGGCGAAACGCCGGGGGGT.T.G.C.GCGTCAGC.GCAAGCGGGGGGCAGCGccccc.cga                                                                                               >  1:3692509/1‑99 (MQ=255)
        ggACAGGGTTGTGTCCTTCGGGCCGCAGCAGGCCCGCTTTTCCTGGGGCGAAACGCGGGGGGGT.T.G.C.GCGTCAGC.GCAAGCGGGGGGCAGCGccccc.cgaa                                                                                              <  2:3735634/100‑4 (MQ=255)
            aGGATTGTGTCCTTCGCGCCGCAGCAGGCCCGCTTTTCCTTGGGCGAAACGCCGGGGGGT.T.T.C.GCGTCAGC.GCAAGCGGGGGGCAGCG.CCCC.CC.G.AAgcgc                                                                                       >  1:3744546/1‑101 (MQ=255)
            aGGATTGTGTCCTTCGCGCCGCAGCAGGCCCGCTTTTCCTTGGGCGAAACGCCGGGGGGT.T.G.C.GCGTCAGC.GCAAGCGGGGGGCAGCG.CCCC.CC.G.AAgcgc                                                                                       >  1:118300/1‑101 (MQ=255)
             ggATTGTGTCCTTCGCGCCGCAGCAGGCCCGCTTTTCCTTGGGCGAAACGCCGGGGGGT.T.G.C.GCGTCAGC.GCAAGCGGGGGGCAGCG.CCCC.CC.G.AAGCGcc                                                                                      >  1:1359000/1‑101 (MQ=255)
             ggATTGTGTCCTTCGCGCCGCAGCAGGCCCGCTTTTCCTTGGGCGAAACGCCGGGGGGT.T.G.C.GCGTCAGC.GCAAGCGGGGGGCAGCG.CCCC.CC.G.AAGCGcc                                                                                      >  1:938693/1‑101 (MQ=255)
                 tgtgTCCTTCGCGCCGCAGCGGCCCCGTTTTTCCTGGGGCGAACCGCGGGGGGGT.T.G.C.GCGTCAGC.GCAAGCGGGGGGCAGCG.CCCC.CC.G.AAGCGCCTGAg                                                                                  <  2:2688947/101‑1 (MQ=255)
                   tgtCCTTCGCGCCGCAGCAGGCCCGCTTTTCCTTGGGCGAAACGCCGGGGGGT.T.G.C.GCGTCAGC.GCAAGCGGGGGGCAGCG.CCCC.CC.G.AAGCGCCTGAGgc                                                                                >  2:1723621/1‑101 (MQ=255)
                       cTTCGCGCCGCAGCAGGCCCGCTTTTCCTTGGGCGAAACGCCGGGGGGT.T.G.C.GCGTCAGC.GCAAGCGGGGGGCAGCG.CCCC.CC.G.AAGCGCCTGAGGCgcgg                                                                            >  1:360826/1‑101 (MQ=255)
                       cTTCGCGCCGCAGCAGGCCCGCTTTTCCTTGGGCGAAACGCCGGGGGGT.T.G.C.GCGTCAGC.GCAAGCGGGGGGCAGCG.CCCC.CC.G.AAGCGCCTGAGGCgcgg                                                                            >  2:2752758/1‑101 (MQ=255)
                          cgcgCGGCAGCGGCCCCTTTTTTCCTTGGGCGAACCGCGGGGGGTT.G.G.C.GCGTCAGC..CAAGGGGGGGGGCGGCCCCCC.CC.G.AAGCGCCTGAGGCGcggcgg                                                                         <  1:611808/101‑1 (MQ=255)
                           cccccgcAGCGGGCCCGTTTTCCCTGGGGCGAACCGCGGGGGGGT.T.G.C.GCGTCAGC.GCAGGGGGGGGGCGGC.CCCCC.CC.G.AAGCGCCTGAGGCGCGGCggg                                                                        <  2:2527254/98‑1 (MQ=255)
                             gccgcAGCAGGCCCGCTTTTCCTTGGGCGAAACGCCGGGGGGT.T.G.C.GCGTCAGC.GCAAGCGGGGGGCAGCG.CCCC.CC.G.AAGCGCCTGAGGCGCGGCGGGCt                                                                      >  2:3399678/1‑101 (MQ=255)
                             gccgcAGCAGGCCCGCTTTTCCTTGGGCGAAACGCCGGGGGGT.T.G.C.GCGTCAGC.GCAAGCGGGGGGCAGCG.CCCC.CC.G.AAGCGCCTGAGGCGCGGCGGGCt                                                                      >  1:1961460/1‑101 (MQ=255)
                             gccccAGCAGGCCCGCTTTTCCTTGGGCAAAACGCCGGGGGGT.T.G.C.GCGTCAGC.GCAAGCGGGGGGCAGCG.CCCC.CC.G.AAGCGCCTGAGGCGCGGCGGGCt                                                                      <  2:2833210/101‑1 (MQ=38)
                              ccggcgccggCCCCGTTTTCCTTGGGGTCAACCCCGGGGGGG.T.GCC.GCGCCGCC.CCACGGGGGGGGC.CCCCCCCC.CC.G.AAGCGCCTGAGGCGCGGCGGGCt                                                                      <  2:3530955/93‑1 (MQ=255)
                              ccgcGGCGGCCCCGTTTTTCCTGGGGCGAACCGCGGGGGGGT.T.G.C.GCGTCAGC.GCAGGGGGGGGGCAGC.CCCCC.CC.G.AAGCGCCTGAGGCGCGGCGGGCtt                                                                     <  2:3342982/101‑1 (MQ=255)
                              ccgcAGCAGGCCCGCTTTTCCTTGGGCGAAACGCCGGGGGGT.T.G.C.GCGTCAGC.GCAAGCGGGGGGCAGCG.CCCC.CC.G.AAGCGCCTGAGGCGCGGCGGGCtt                                                                     >  1:209665/1‑101 (MQ=255)
                              ccgcAGCAGCCCCGCTTTCCCTGGGGCGAACCGCCGGGGGGT.T.G.C.GCGTCAGC.GCAAGGGGGGGGCAGC.CCCCC.CC.G.AAGCGCCTGAGGCGCGGCGGGCtt                                                                     <  1:2687493/101‑1 (MQ=255)
                               ccgcggcggCCCCGTTTTCCTTGGGGCAAACCGGGGGGGGT.TGG.C.GCGCCGGC.AAAGGGGGGGGGCA.CCCCCCC.CC.G.AAGCGCCTGAGGCGCGGCGGGCtt                                                                     <  2:2332146/94‑1 (MQ=255)
                                gcagcaGGCCCGCTTTTCCTTGGGCGAAACGCCGGGGGGT.T.G.C.GCGTCAGC.GCAAGCGGGGGGCAGCG.CCCC.CC.G.AAGCGCCTGAGGCGCGGCgggggttc                                                                   >  1:2935076/1‑96 (MQ=255)
                                gcagcaGGCCCGCTTTTCCTTGGGCGAAACGCCGGGGGGT.T.G.C.GCGTCAGC.GCAAGCGGGGGGCAGCG.CCCC.CC.G.AAGCGCCTGAGGCGCGGCGGGCTTgc                                                                   >  1:809162/1‑101 (MQ=255)
                                gcagcaGGCCCGCTTTTCCTTGGGCGAAACGCCGGGGGGT.T.G.C.GCGTCAGC.GCAAGCGGGGGGCAGCG.CCCC.CC.G.AAGCGCCTGAGGCGCGGCGGGCTTgc                                                                   >  2:2439647/1‑101 (MQ=255)
                                 cagcagGCCCGCTTTCCCTGGGGCGAACCGCCGGGGGTT.G.G.C.GCGTCAGC.GAAAGGGGGGGGCAGC.CCCCC.CC.G.AAGCGCCTGAGGCGCGGCGGGCTTgcc                                                                  <  1:2527534/101‑1 (MQ=37)
                                   gcggccccttttttCCCTGGGC.AAACGCGGGGGGGT.T.C.G.GCGGGCGCGGGGAGGGGGGGGGCGCCCCCCC.CG.G.AAGCGCCTGAGGCGCGGCGGGCTTgccc                                                                 <  2:4029393/91‑1 (MQ=255)
                                    cagGCCCGCTTTTCCTTGGGCGAAACGCCGGGGGGT.T.G.C.GCGTCAGC.GCAAGCGGGGGGCAGCG.CCCC.CC.G.AAGCGCCTGAGGCGCGGCGGGCTTgcccgc                                                               >  2:2703489/1‑101 (MQ=255)
                                     agGCCCGCTTTTCCTTGGGCGAAACGCCGGGGGGT.T.G.C.GCGTCAGC.GCAAGCGGGGGGCAGCG.CCCC.CC.G.AAGCGCCTGAGGCGCGGCGGGCTTgcccgcc                                                              >  1:2074361/1‑101 (MQ=255)
                                      ggCCCGCTTTTCCTTGGGCGAAACGCCGGGGGGT.T.G.C.GCGTCAGC.GCAAGCGGGGGGCAGCG.CCCC.CC.G.AAGCGCCTGAGGCGCGGCGGGCTTgcccgccc                                                             >  1:2482394/1‑101 (MQ=255)
                                      ggCCCGCTTTTCATTGGGCGAAACGCCGGGGGGT.T.G.C.GCGTCAGC.GCAAGCGGGGGGCAGCG.CCCC.CC.G.AACCGCCTGCGGCGCGGCGGGGTTcccccccc                                                             >  1:2296599/1‑101 (MQ=255)
                                       gCCCGCTTTTCCTTGGGCGAAACGCCGGGGGGT.T.G.C.GCGTCAGC.GCAAGCGGGGGGCAGCG.CCCC.CC.G.AAGCGCCTGAGGCGCGGCGGGCTTGCCCGCCCa                                                            >  1:1049247/1‑101 (MQ=255)
                                        cccttttttCCTGGGGC.AAACGGGGGGGGGG.T.G.CGGCGCCGGC.GCACGGGGGGGGC.GCCCCCCC.CG.G.AAGCGCCTGAGGCGCGGCGGGCTTGCCCGCCCaa                                                           <  2:3770593/96‑1 (MQ=12)
                                          cGCTTTTCCTTGGGCGAAACGCCGGGGGGT.T.G.C.GCGTCAGC.GCAAGCGGGGGGCAGCG.CCCC.CC.G.AAGCGCCTGAGGCGCGGCGGGCTTGCCCGGCCAACg                                                         >  2:1126843/1‑101 (MQ=255)
                                            cTTTTCCTTGGGCGAAACGCCGGGGGGT.T.G.C.GCGTCAGC.GCAAGCGGGGGGCAGCG.CCCC.CC.G.AAGCGCCTGAGGCGCGGCGGGCTTGCCCGCCCAACGaa                                                       >  2:1894995/1‑101 (MQ=255)
                                              tttCCTGGGGCGAAGCGCGGGGGGGT.T.G.C.GCGTCCGC.GCAAGGGGGGGGCGCCCCCCCC.CC.G.AAGCGCCTGAGGCGCGGCGGGCTTGCCCGCCCAACGAAg                                                      <  2:1195769/101‑1 (MQ=37)
                                               ttCCTTGGGCGAAACGCCGGGGGGT.T.G.C.GCGTCAGC.GCAAGCGGGGGGCAGCG.CGCC.CC.G.AAGCGCCTGAGGCGTGGCGGGCTTGCCCGCCCAAGGAAGgc                                                    >  2:3494779/1‑101 (MQ=25)
                                               ttCCTTGGGCGAAACGCCGGGGGGT.T.G.C.GCGTCAGC.GCAAGCGGGGGGCAGCG.CCCC.CC.G.AAGCGCCTGAGGCGCGGCGGGCTTGCCCGCCCAACGAAggg                                                    >  2:191669/1‑100 (MQ=255)
                                               ttCCTTGGGCGAAACGCCGGGGGGT.T.G.C.GCGTCAGC.GCAAGCGGGGGGCAGCG.CCCC.CC.G.AAGCGCCTGAGGCGCGGCGGGCCTGCCCGCCCAACGAAggg                                                    >  1:3163682/1‑100 (MQ=255)
                                                   ttGGGCGAAACGCCGGGGGGT.T.G.C.GCGTCAGC.GCAAGCGGGGGGCAGCG.CCCC.CC.G.AAGCGCCTGAGGCGCGGCGGGGTTGCCCGCCCAACGAAGGCGCaa                                                >  2:1489145/1‑101 (MQ=255)
                                                    tGGGCGAAACGCCGGGGGGT.T.G.C.GCGTCAGC.GCAAGCGGGGGGCAGCG.CCCC.CC.G.AAGCGCCTGAGGCGCGGCGGGCTTGCCCGCCCAACGAAGGCGCAac                                               >  2:2501449/1‑101 (MQ=255)
                                                     gggCGAAACGCCGGGGGGT.T.G.C.GCGTCAGC.GCAAGCGGGGGGCAGCG.CCCC.CC.G.AAGCGCCTGAGGCGCGGCGGGCTTGCCCGCCCAACGAAGGCGCAacc                                              >  2:2226295/1‑101 (MQ=255)
                                                     gggCGAAACGCCGGGGGGT.T.G.C.GCGTCAGC.GCAAGCGGGGGGCAGCG.CCCC.CC.G.AAGCGCCTGAGGCGCGGCGGGCCTGCCCGCCCAACGAAGGCGCAacc                                              >  1:3426237/1‑101 (MQ=255)
                                                      ggggAAACGGCGGGGGGG.T.G.C.GCCGCCCC.CAGAGGGGGGGGCCGCCGCCCC.CC.G.AAGCGCCTGAGGCGCGGCGGGCTTGCCCGCCCAACGAAGGCGCAacc                                              <  2:3295713/98‑1 (MQ=12)
                                                       gCGAAACGCCGGGGGGT.T.G.C.GCGTCAGC.GCAAGCGGGGGGCAGCG.CCCC.CC.G.AAGCGCCTGAGGCGCGGCGGGCTTGCCCGCCCCACGAAGGCGCAacccc                                            >  1:1542790/1‑99 (MQ=25)
                                                        cGAAACGCCGGGGGGT.T.G.C.GCGTCAGC.GCAAGCGGGGGGCAGCG.CCCC.CC.G.AAGCGCCTGAGGCGCGGCGGGCTTGCCCGCCCAACGAAGGCGCAaacacc                                           >  1:4070473/1‑101 (MQ=255)
                                                         gAAACGCGGGGGGGT.G.G.C.GCGTCAGC.GAAAGCGGGGGGCCGC.CCCCC.CC.G.AAGCGCCTGAGGCGCGGCGGGCTTGCCCGCCCAACGAAGGCGCAACCACCt                                          <  2:4066625/101‑1 (MQ=25)
                                                         gAAACGCCGGGGGGT.T.G.C.GCGTCAGC.GCAAGCGGGGGGCAGCG.CCCC.CC.G.AAGCGCCTGAGGCGCGGCGGGCTTGCCCGCCCAACGAAGGCGCAACCACCt                                          >  1:3757323/1‑101 (MQ=255)
                                                          aaaCGCCGGGGGGT.T.G.C.GCGTCAGC.GCAAGCGGGGGGCAGCGCCCCC.CC...AAGCGCCTGCGGCGCGGCGGGCTTTCCCGCCCCACGAACGCCCCACCACCtc                                         >  2:3990590/1‑100 (MQ=255)
                                                          aaaCGCCGGGGGGT.T.G.C.GCGTCAGC.GCAAGCGGGGGGCAGCG.CCCC.CC.G.CAGCGCCTGAGGCGCGGCGGGCTTGCCCCCCCAACGAAGGCGCAacccccct                                         >  1:3181607/1‑99 (MQ=25)
                                                          aaaCGCCGGGGGGT.T.G.C.GCGTCAGC.GCAAGCGGGGGGCAGCG.CCCC.CC.G.AAGCGCCTGAGGCGCGGCGGGCTTGCCCGCCCAACGAAGGCGCAACCACCtt                                         <  2:2005426/101‑1 (MQ=255)
                                                               gcGGGGGGTGG.G.C.GCGTGAGC..CGAGGGGGGGGC.CCGCCCCC.CCGG.AAGCGCCTGAGGCGCGGCGGGCTTGCCCGCCCAACGAAGGCGCAACCACCTTCAGt                                     <  1:2729149/100‑1 (MQ=18)
                                                               ccGGGGGGT.T.G.C.GCGTCAGC.GCAAGCGGGGGGCAGCG.CCCC.CC.G.AAGCGCCTGAGGCGCGGCGGGCTTGCCCGCCCAACGAAGGCGCAACCCCCTTCAGTc                                    >  2:1930683/1‑101 (MQ=255)
                                                                cGGGGGGT.T.G.C.GCGTCAGC.GCAAGCGGGGGGCAGCG.CCCC.CC.G.AAGCGCCTGAGGCGCGGCGGGCTTGCCCGCCCAACGAAGGCGCAACCACCTTCAGTCg                                   >  2:2767215/1‑101 (MQ=255)
                                                                 ggggggT.T.G.C.GCGTCAGC.GCAAGCGGGGGGCAGCG.CCCC.CC.G.AAGCGCCTGAGGCGCGGCGGGCTTGCCCGCCCAACGAAGGCGCAACCACCTTCAGTCgg                                  >  1:623155/1‑101 (MQ=255)
                                                                 ggggggT.T.G.C.GCGTCAGC.GCAAGCGGGGGGCAGCG.CCCC.CC.G.AAGCGCCTGAGGCGCGGCGGGCTTGCCCGCCCAACGAAGGCGCAACCACCTTCAGTCgg                                  >  2:2359289/1‑101 (MQ=255)
                                                                  gggggT.T.G.C.GCGTCAGC.GCAAGCGGGGGGCAGCG.CCCC.CC.G.AAGCGCCTGAGGCGCGGCGGGCTTGCCCGCCCCACGAAGGCGCAACCACCCTCCGTCGGc                                 >  1:1467093/1‑101 (MQ=25)
                                                                  gggggT.T.G.C.GCGTCAGC.GCAAGCGGGGGGCAGCG.CCCC.CC.G.AAGCGCCTGAGGCGCGGCGGGCTTGCCCGCCCAACGAAGGCGCAACCACCTTCAGTCGGc                                 >  2:2066494/1‑101 (MQ=255)
                                                                   ggggT.T.G.C.GCGTCAGC.GCAAGCGGGGGGCAGCG.CCCC.CC.G.AAGCGCCTGAGGCGCGGCGGGCTTGCCCGCCCCACGAAGGCGCAACCACCTTCAGTCGGCt                                >  1:954171/1‑101 (MQ=25)
                                                                   ggggT.T.G.C.GCGTCAGC.GCAAGCGGGGGGCAGCG.CCCC.CC.G.AAGCGCCTGAGGCGCGGCGGGCTTGCCCGCCCAACGAAGGCGCAACCACCTTCAGTCGGCt                                >  2:3932002/1‑101 (MQ=255)
                                                                   ggggT.T.G.C.GCCTCAGC.GCAGGGGGGGGGCAGC.CCCCC.CC.G.AAGCGCCTGAGGCGCGGCGGGCTTGCCCGCCCAACGAAGGCGCAACCACCTTCAGTCGGCt                                <  2:2562354/101‑1 (MQ=255)
                                                                    gggT.T.G.C.GCGTCAGC.GCAAGCGGGGGGCAGCG.CCCC.CC.G.AAGCGCCTGAGGCGCGGCGGGCTTGCCCGCCCAACGAAGGCGCAACCACCTTCAGTCGGCtt                               >  2:944910/1‑101 (MQ=255)
                                                                    gggT.T.G.C.GCGTCAGC.GAAAGCGGGGGGCAGCG.CCCC.CC.G.AAGCGCCTGAGGCGCGGCGGGCTTGCCCGCCCAACGAAGGCGCAACCACCTTCAGTCGGCtt                               <  2:2822106/101‑1 (MQ=255)
                                                                      gT.T.G.C.GCGCCGGC.GCAAGGGGGGGGCACCG.CCCC.CC.G.AAGCGCCTGAGGCGCGGCGGGCTTGCCCGCCCAACGAAGGCGCAACCACCTTCAGTCGGCTTCg                             <  2:546626/101‑1 (MQ=25)
                                                                               gcgTCAGC.GCAAGGGGGGGGCCGCG.CCCC.CC.G.AAGCGCCTGAGGCGCGGCGGGCTTGCCCGCCCAACGAAGGCGCAACCACCTTCAGTCGGCTTCGATccc                        <  1:3999324/101‑1 (MQ=255)
                                                                               gcgTCAGC.GCAAGCGGGGGGCAGCG.CCCC.CC.G.AAGCGCCTGAGGCGCGGCGGGCTTGCCCGCCCAACGAAGGCGCAACCACCTTCAGTCGGCTTCGATccc                        >  1:1865637/1‑101 (MQ=255)
                                                                               gcgTCAGC.GCAAGCGGGGGGCAGCG.CCCC.CC.G.AAGCGCCTGAGGCGCGGCGGGCTTGCCCGCCCAACGAAGGCGCAACCACCTTCAGTCGGCTTCGATccc                        >  2:2416981/1‑101 (MQ=255)
                                                                                cgTCAGC.GCAAGCGGGGGGCAGCG.CCCC.CC.G.AAGCGCCTGAGGCGCGGCGGGCTTGCCCGCCCAACGAAGGCGCAACCACCTTCAGTCGGCTTCGATcccc                       >  2:2303807/1‑101 (MQ=255)
                                                                                 gTCAGC.GCAAGCGGGGGGCAGCG.CCCC.CC.G.AAGCGCCTGAGGCGCGGCGGGCTTGCCCGCCCAACGAAGGCGCAACCACCTTCAGTCGGCTTCGATCCcca                      >  2:3864902/1‑101 (MQ=255)
                                                                                  tCAGC.GCAAGCGGGGGGCAGCG.CCCC.CC.G.AAGCGCCTGAGGCGCGGCGGGCTTGCCCGCCCAACGAAGGCGCAACCACCTTCAGTCGGCTTCGATCCccac                     >  2:161456/1‑101 (MQ=255)
                                                                                   gtacgcgaaGGGGGGGGCCGCCCCCCC.CC.G.AAGCGCCTGAGGCGCGGCGGGCTTGCCCGCCCAACGAAGGCGCAACCACCTTCAGTCGGCTTCGATCCcca                      <  2:1845326/94‑1 (MQ=37)
                                                                                   cAGC.GCAAGCGGGGGGCAGCG.CCCC.CC.G.AAGCGCCTGAGGCGCGGCGGGCTTGCCCGCGCAACGAAGGCGGCACCACCTTCAGNCGGGTTCGGTCCccacc                    >  2:289214/1‑101 (MQ=25)
                                                                                    aGC.GCAAGCGGGGGGCAGCG.CCCC.CC.G.AAGCGCCTGAGGCGCGGCGGGCTTGCCCGCCCAACGAAGGCGCAACCACCTTCAGTCGGCTTCGATCCccacca                   >  1:4155397/1‑101 (MQ=255)
                                                                                     gcgcgggggggggggCCCCCCCCCC.CC.GAAAGCGCCTGAGGCGCGGCGGGCTTGCCCGCCCAACGAAGGCGCAACCACCTTCAGTCGGCTTCGATCCccac                     <  1:274570/94‑1 (MQ=25)
                                                                                        ggagggggggggCGCCCCCCCCGCG.G.AAGCGCCTGAGGCGCGGCGGGCTTGCCCGCCCAACGAAGGCGCAACCACCTTCAGTCGGCTTCGATCCCCACCAg                  <  2:944842/97‑1 (MQ=25)
                                                                                         cAAGCGGGGGGCTGCG.CCCC.CC.G.AAGCGCCTGAGGCGCGGCGGGCTTGCCCGCCCAACGAAGGCGCAACCACCTTCAGTCGGCTTCGATCCCCCCCAgcgc               >  1:2297861/1‑101 (MQ=25)
                                                                                         cAAGCGGGGGGCAGCG.CCCC.CC.G.AAGCGCCTGAGGCGCGGCGGGCTTGCCCGCCCAACGAAGGCGCAACCACCTTCAGTCGGCTTCGATCCCCACCAgcgc               >  1:1157509/1‑101 (MQ=255)
                                                                                          aaGCGGGGGGCAGCG.CCCC.CC.G.AAGCGCCTGAGGCGCGGCGGGCTTGCCCGCCCAACGAAGGCGCAACCACCTTCAGTCGGCTTCGATCCCCCCCAGcgcg              >  2:3094699/1‑101 (MQ=255)
                                                                                             cggggggCAGCG.CCCC.CC.G.AAGCGCCTGAGGCGCGGGGGGCCTGCCCCCCCAACGCAGGCGCAACCCCCCTCCGTCGGCTTCGCTCCCCCCCAGCgcgggg           >  2:788678/2‑99 (MQ=255)
                                                                                             cggggggCAGCG.CCCC.CC.G.AAGCGCCTGAGGCGCGGCGGGCTTGCCCGCCCCACGAAAGCGCAACCACCTTCAGTCGGCTTCGATCCCCACCAGCgcggcg           >  1:3676403/2‑101 (MQ=37)
                                                                                             cggggggCAGCG.CCCC.CC.G.AAGCGCCTGAGGCGCGGCGGGCTTGCCCGCCCAACGAAGGCGCAACCACCTTCAGTCGGCTTCGATCCCCACCAGCgcggcg           >  1:1736666/2‑101 (MQ=255)
                                                                                              ggggggGAGGCCCCCC.CC.G.AAGCGCCTGAGGCGCGGCGGGCTTGCCCGCCCAACGAAGGCGCAACCACCTTCAGTCGGCTTCGATCCCCACCAGCgcggcg           <  2:3280170/101‑1 (MQ=255)
                                                                                               gggggCAGCG.CCCC.CC.G.AAGCGCCTGAGGCGCGGCGGGGTTTGCCCCCCCACGAAAGCGCAACCACCCTTCGTCGGGCTCGGTTCCCCCCACCgcgggggg         >  2:743329/1‑100 (MQ=12)
                                                                                               gggggCAGCG.CCCC.CC.G.AAGCGCCTGAGGCGCGGCGGGCTTGCCCGCCCAACGAAGGCGCAACCACCTTCAGTCGGCTTCGATCCCCACCAGCgcggcggc         >  1:3735573/1‑101 (MQ=255)
                                                                                               gggggCAGCG.CCCC.CC.G.AAGCGCCTGAGGCGCGGCGGGCTTGCCCGCCCAACGAAGGCGCAACCACCTTCAGTCGGCTTCGATCCCCACCAGCgcggcggc         >  2:2321797/1‑101 (MQ=255)
                                                                                               gggggCAGCG.CCCC.CC.G.AAGCGCCTGAGGCGCGGCGGGCTTGCCCGCCCAACGAAGGCGCAACCACCTTCAGTCGGCTTCGATCCCCACCAGCgcggcggc         >  2:2906201/1‑101 (MQ=255)
                                                                                                ggggCAGCG.CCCC.CC.G.AAGCGCCTGAGGCGCGGCGGGCTTGCCCGCCCAACGAAGGCGCAACCACCTTCAGTCGGCTTCGATCCCCACCAGCgcggcggcg        >  2:189406/1‑101 (MQ=255)
                                                                                                ggggCAGCG.CCCC.CC.G.AAGCGCCTGAGGCGCGGCGGGCTTGCCCGCCCAAAGGAGGCGCAAACACCTTCAGTCGGCTTCGATCCCCACCAGCgggggggcc        >  2:1552760/1‑100 (MQ=25)
                                                                                                   gCAGCG.CCCC.CC.G.AAGCGCCTGAGGCGCGGCGGGCTTGCCCGCCCAACGAAGGCGCAACCACCTTCAGTCGGCTTCGATCCCCACCAGCGCGGCGGCGAAt     >  2:3784052/1‑101 (MQ=255)
                                                                                                   gCAGCG.CCCC.CC.G.AAGCGCCTGAGGCGCGGCGGGCTTGCCCGCCCAACGAAGGCGCAACCACCTTCAGTCGGCTTCGATCCCCACCAGCGCGGCGGCGAAt     >  1:2682185/1‑101 (MQ=255)
                                                                                                    cccccccccc.gc.GAAAGCGCCTGAGGCGCGGCGGGCTTGCCCGCCCAACGAAGGCGCAACCACCTTCAGTCGGCTTCGATCCCCACCAGCGCGGCGGCGaa      <  2:1447878/96‑1 (MQ=37)
                                                                                                    cAGCG.CCCC.CC.G.AAGCGCCTGAGGCGCGGCGGGCTTGCCCGCCCAACGAAGGCGCAACCACCTTCAGTCGGCTTCGATCCCCACCAGCGCGGCGGCGAAtt    >  2:1829051/1‑101 (MQ=255)
                                                                                                         |                                                                                                     
CGGTTTCCTGACAGGATTGTGTCCTTCGCGCCGCAGCAGGCCCGCTTTTCCTTGGGCGAAACGCCGGGGGGT.T.G.C.GCGTCAGC.GCAAGGGGGGGGCAGCGCCCCC.CC.G.AAGCGCCTGAGGCGCGGCGGGCTTGCCCGCCCAACGAAGGCGCAACCACCTTCAGTCGGCTTCGATCCCCACCAGCGCGGCGGCGAATTCC  >  NC_014034/2559835‑2560033

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