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% 20.9 129intergenic (‑211/+59)ispZ/ftsYintracellular septation protein A/cell division protein FtsY
Reads supporting (aligned to +/- strand):  new base (111/9):  ref base (0/9):  total (111/18)
Fisher's exact test strand distribution p-value = 2.37e-09
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 2.10e-07

Only 100 of 176 total aligned reads displayed.

CGGTTTCCTGACAGGATTGTGTCCTTCGCGCCGCAGCAGGCCCGCTTTTCCTTGGGCGAAACGCCGGGGGGTTGCGCGTCAGCGCAAGG.GGGGGGCAGCGCCCCCCCGAAGCGCCTGAGGCGCGGCGGGCTTGCCCGCCCAACGAAGGCGCAACCACCTTCAGTCGGCTTCGATCCCCACCAGCGCGGCGGCGAATTCCTC  >  NC_014034/2559835‑2560035
                                                                                                     |                                                                                                    
cGGTTGCCTGACATGATTTTGTCCTTTGCGACGCAGCAGGCCCTCTTTTCCTAGGGCGAAACGCTAGGGGGTTGCGCGTCAGCGCAAGG.GGTGGGCAgcgc                                                                                                      <  2:1911132/101‑1 (MQ=255)
 ggTTTCCTGACAGGATTGTGTCCTTCGCGCCGCAGCAGGCCCGCTTTTCCTTGGGCGAAACGCCGGGGGGTTGCGCGTCAGCGCAAGC.GGGGGGCAGCGcc                                                                                                     >  1:3917589/1‑101 (MQ=255)
    ttCCTGACCGGATTGTGTCCTTCGCGCCGCAGCAGGCCCGCTTTTCCTTGGGCGAAACGCCGGGGGGTTGCGCGTCAGCGCAAGC.GGGGGGCAGCGccccc                                                                                                  <  2:1660893/101‑1 (MQ=255)
       cTGACAGGATTGTGTCCTCCGCGCCCCAGCAGCCCCGCTTTTCCTGGGGCGAAACGCGGGGGGGTTGCGCGTCAGCGCAAGG.GGGGGGCAGC.CCCCCCCGa                                                                                              <  2:577936/101‑1 (MQ=255)
         gACAGGATTGTGTCCTTCGCGCCGCAGCAGGCCCGCTTTTCCTTGGGCGAAACGCCGGGGGGTTGCGCGTCAGCGCAAGC.GGGGGGCAGCGccccccgaag                                                                                             >  2:3971604/1‑97 (MQ=255)
         gACAGGATTGTGTCCTTCGCGCCGCAGCAGGCCCGCTTTTCCTTGGGCGAAACGCCGGGGGGCTGCGCGTCAGCGCAAGC.GGGGGGCAGCGccccccgaag                                                                                             >  2:246481/1‑97 (MQ=255)
            aGGATTGTGTCCTTCGCGCCGCAGCAGGCCCGCTTTTCCTTGGGCGAAACGCCGGGGGGTTGCGCGTCAGCGCAAGC.GGGGGGCAGCG.CCCCCCGAAgcgc                                                                                         >  1:5113397/1‑101 (MQ=255)
            aGGATTGTGTCCTACGCGCCGCAGCAGGCCCGCTTTTCCTTGGGCGAAACGCCGGGGGGTTGCGCGTCAGCGCAAGC.GGGGGGCAGCG.CCCCCCGAAgcgc                                                                                         >  1:347752/1‑101 (MQ=255)
             ggATTGTGTCCTTCGCGCCGCAGCAGGCCCGCTTTTCCTTGGGCGAAACGCCGGGGGGTTGCGCGTCAGCGCAAGC.GGGGGGCAGCG.CCCCCCGAAGCGcc                                                                                        >  2:4232947/1‑101 (MQ=255)
             ggATTGTGTCCTTCGCGCCGCAGCAGGCCCGCTTTTCCTTGGGCGAAACGCCGGGGGGTTGCGCGTCAGCGCAAGC.GGGGGGCAGCG.CCCCCCGAAGCGcc                                                                                        >  1:5476700/1‑101 (MQ=255)
                ttGTGTCCTTCGCGCCGCAGCAGGCCCGCTTTTCCTTGGGCGAAACGCCGGGGGGTTGCGCGTCAGCGCAAGC.GGGGGGCAGCG.CCCCCCGAAGCGCCTGa                                                                                     >  2:4716938/1‑101 (MQ=255)
                 tgtgTCCTTCGCGCCGCAGCAGGCCCGCTTTTCCTTGGGCGAAACGCCGGGGGGTTGCGCGTCAGCGCAAGC.GGGGGGCAGCG.CCCCCCGAAGCGCCTGAg                                                                                    >  2:543043/1‑101 (MQ=255)
                  gtgtCCTTCGCGCCGCAGCAGGCCCGCTTTTCCTTGGGCGAAACGCCGGGGGGTTGCGCGTCAGCGCAAGC.GGGGGGCAGCG.CCCCCCGAAGCGCCTGAgg                                                                                   >  2:1178635/1‑101 (MQ=255)
                   tgtCCTTCGCGCCGCAGCAGGCCCGCTTTTCCTTGGGCGAAACGCCGGGGGGTTGCGCGTCAGCGCAAGC.GGGGGGCAGCG.CCCCCCGAAGCGCCTGAGgc                                                                                  >  2:3267478/1‑101 (MQ=255)
                       cTTCGCGCCGCAGCAGGCCCGCTTTTCCTTGGGCGAAACGCCGGGGGGTTGCGCGTCAGCGCAAGC.GGGGGGCAGCG.CCCCCCGAAGCGCCTGAGGCgcgg                                                                              >  2:4950281/1‑101 (MQ=255)
                       cTTCGCGCCGCAGCAGGCCCGCTTTTCCTTGGGCGAAACGCCGGGGGGTTGCGCGTCAGCGCAAGC.GGGGGGCAGCG.CCCCCCGAAGCGCCTGAGGCgcgg                                                                              >  1:710708/1‑101 (MQ=255)
                        ttCGCGCCGCAGCAGGCCCGCTTTTCCTTGGGCGAAACGCCGGGGGGTTGCGCGTCAGCGCAAGC.GGGGGGCAGCG.CCCCCCGAAGCGCCTGAGGCgcggc                                                                             >  1:4334207/1‑101 (MQ=255)
                        ctCGCGCCCGGGCCGGCCCGTTTTTCCTTGGGGCAAACGGCGGGGGGTTGGCCCGCCGCCGAAGG.GGGGGGCCCCCCCCCCCCGAAGCGCCTGAGGCgcgg                                                                              <  2:1524440/100‑1 (MQ=255)
                            cgccgcAGCAGGCCCGCTTTTCCTTGGGCGAAACGCCGGGGGGTTGCGCGTCAGCGCAAGC.GGGGGGCAGCG.CCCCCCGCAGCGCCTGAGGCGCGGCGGGc                                                                         >  1:5422269/1‑101 (MQ=255)
                             gccgcAGCAGGCCCGCTTTTCCTTGGGCGAAACGCCGGGGGGTTGCGCGTCAGCGCAAGC.GGGGGGCAGCG.CCCCCCGAAGCGCCTGAGGCGCGGCGGGCt                                                                        >  2:3714896/1‑101 (MQ=255)
                               cgcAGCAGGCCCGCTTTTCCTTGGGCGAAACGCCGGGGGGTTGCGCGTCAGCGCAAGC.GGGGGGCAGCG.CCCCCCGAAGCGCCTGAGGCgcggcgagatcg                                                                      >  1:4686479/1‑95 (MQ=255)
                               cgcAGCAGGCCCGCTTTTCCTTGGGCGAAACGCCGGGGGGTTGCGCGTCAGCGCAAGC.GGGGGGCAGCG.CCCCCCGAAGCGCCTGAGGCGCGGCGGGCTTg                                                                      >  2:641641/1‑101 (MQ=255)
                                gcagcaGGCCCGCTTTTCCTTGGGCGAAACGCCGGGGGGTTGCGCGTCAGCGCAAGC.GGGGGGCAGCG.CCCCCCGAAGCGCCTGAGGCGCGGCGGGCTTgc                                                                     >  1:2859284/1‑101 (MQ=255)
                                caggcagCCCCTCTTTCCCTTGGGCGAACCGCGGGGGGGTGCCGCGTCAGCCCACGC.GGGGGGCAGC.CCCCCCCGAAGCGCCTGAGGCGCGGCGGGCTTgc                                                                     <  1:2632711/98‑1 (MQ=19)
                                 cagcagGCCCGCTTTTCCTTGGGCGAAACGCCGGGGGGTTGCGCGTCAGCGCAAGC.GGGGGGCAGCG.CCCCCCGAAGCGCCTGAGGCGCGGCGGGCTTgcc                                                                    >  2:306016/1‑101 (MQ=255)
                                    cagGCCCGCTTTTCCTTGGGCGAAACGCCGGGGGGTTGCGCGTCAGCGCAAGC.GGGGGGCAGCG.CCCCCCGAAGCGCCTGAGGCGCGGCGGGCTTgccccc                                                                 >  1:4140274/1‑99 (MQ=255)
                                       gCCCGCTTTTCCTTGGGCGAAACGCCGGGGGGTTGCGCGTCAGCGCAAGC.GGGGGGCAGCG.CCCCCCGAAGCGCCTGAGGCGCGGCGGGCTTGCCCGCCCa                                                              >  1:4933404/1‑101 (MQ=255)
                                         ccGTTTTTCCTGGGGGGAAGCGCGGGGGGGTGGCGCGCCAGCGCAGGG.GGGGGGCGGC.CCCCCCCGAAGCGCCTGAGGCGCGGCGGGCTTGCCCGCCCAAc                                                            <  1:2692387/101‑1 (MQ=21)
                                          cGCTTTTCCTTGGGCGAAACGCCGGGGGGTTGCGCGTCAGCGCAAGC.GGGGGGCAGCG.CCCCCCGAAGCGCCTGAGGCGCGGCGGGCTTGGCCGCCCAACg                                                           >  2:473247/1‑101 (MQ=255)
                                          cGCTTTTCCTTGGGCGAAACGCCGGGGGGTTGCGCGTCAGCGCAAGC.GGGGGGCAGCG.CCCCCCGAAGCGCCTGAGGCGCGGCGGGCTTGCCCGCCCCACg                                                           >  1:2019608/1‑101 (MQ=255)
                                          cGCTTTTCCTTGGGCGAAACGCCGGGGGGTTGCGCGTCAGCGCAAGC.GGGGGGCAGCG.CCCCCCGAAGCGCCTGAGGCGCGGCGGGCTTGCCCGCCCAACg                                                           >  2:33096/1‑101 (MQ=255)
                                            cTTTTCCTTGGGCGAAACGCCGGGGGGTTGCGCGTCAGCGCAAGC.GGGGGGCAGCG.CCCCCCGAAGCGCCTGAGGCGCGGCGGGCTTGCCCGCCCAACGaa                                                         >  1:426479/1‑101 (MQ=255)
                                             ttttCCTTGGGCGAAACGCCGGGGGGTTGCGCGTCAGCGCAAGC.GGGGGGCAGCG.CCCCCCGAAGCGCCTGAGGCGCGGCGGGCTTGCCCGCCCCACGAAg                                                        >  2:3505754/1‑101 (MQ=255)
                                             ttttCCTTGGGCGAAACGCCGGGGGGTTGCGCGTCAGCGCAAGC.GGGGGGCAGCG.CCCCCCGAAGCGCCTGAGGCGCGGCGGGCTTGCCCGCCCAACGAAg                                                        >  1:3861993/1‑101 (MQ=255)
                                             ttttCCTTGGGCGAAACGCCGGGGGGTTGCGCGTCAGCGCAAGC.GGGGGGCAGCG.CCCCCCGAAGCGCCTGAGGCGCGGCGGGCTTGCCCGCCCAACGAAg                                                        >  1:273065/1‑101 (MQ=255)
                                             ttttCCTTGGGCGAAACGCCGGGGGGTTGCGCGTCAGCGCAAGC.GGGGGGCAGCG.CCCCCCGAAGCGCCTGAGGCGCGGCGGGCTTGCCCGCCCAACGAAg                                                        >  2:4592652/1‑101 (MQ=255)
                                              tttCCTTGGGCGAAACGCCGGGGGGTTGCGCGTCAGCGCAAGC.GGGGGGCGGCG.CCCCCCGAAGCGCCTGAGGCGCGGCGGGCCTTCCCCCCCAAcccagg                                                       >  2:1985646/1‑96 (MQ=255)
                                              tttCCTTGGGCGAAACGCCGGGGGGTTGCGCGTCAGCGCAAGC.GGGGGGCAGCG.CCCCCCGAAGCGCCTGAGGCGCGGCGGGCTTGCCCGCCCAACGAAgg                                                       >  1:4688703/1‑101 (MQ=255)
                                              ttcgcgggggGGAAACGGGGGGGGGTTGGGGCGCCGCCGAAGG.GGGGGGCCCCCCCCCCCGGAAGCGCCTGAGGCGCGGCGGGCTTGCCCGCCCAACGAAg                                                        <  2:1234484/94‑1 (MQ=12)
                                               ttCCTTGGGCGAAACGCCGGGGGGTTGCGCGTCAGCGCAAGC.GGGGTGCGGGGACCCCCC.AAGAGCCAGGGGAGCGGCGGCCTTACCCGCCCAGCCAAGgc                                                      >  1:1821433/1‑101 (MQ=255)
                                                tCCTTGGGCGAAACGCCGGGGGGTTGCGCGTCAGCGCAAGC.GGGGGGCAGCG.CCCCCCGAAGCGCCTGAGGCGCGGCGGGCTTGCCCGCCCAACGAAAgcg                                                     >  1:618540/1‑101 (MQ=255)
                                                tCCTGGGGCAAACCGCGGGGGGTTTGCGCGCCAGCGCAGGG.GGGGGGCAGC.CCCCCCCGAAGCGCCTGAGGCGCGGCGGGCTTGCCCGCCCAACGAAGgcg                                                     <  1:2809807/101‑1 (MQ=21)
                                                 ccTTGGGCGAAACGCCGGGGGGTTGCGCGTCAGCGCAAGC.GGGGGGCAGCG.CCCCCCGAAGCGCCTGAGGCGCGGCGGGCTTGCCCGCCCAACGAAGgcgc                                                    >  1:1676393/1‑101 (MQ=255)
                                                    ggggCGAAGCGCGGGGGGGCGCCGCG.GAGCCCAAGGCGGGGGGCCGC.CCCCCCCGAAGCGCCTGAGGCGCGGCGGGCTTGCCCGCCCAACGAAGGCGCAac                                                 <  1:3568796/100‑1 (MQ=12)
                                                     gggCGAAACGCCGGGGGGTTGCGCGTCAGCGCAAGC.GGGGGGCAGCG.CCCCCCGAAGCGCCTGAGGCGCGGCGGGCTTGCCCGCCCAACGAAGGCGCAacc                                                >  2:5783993/1‑101 (MQ=255)
                                                         gAAACGCCGGGGGGTTGCGCGTCAGCGCAAGC.GGGGGGCAGCG.CCCCCCGAAGCGCCTGAGGCGCGGGGGGGTTGCCCCCCCCACCAAAGGGCCacccccc                                            >  1:1136625/1‑100 (MQ=21)
                                                          aaaCGCCGGGGGGTTGCGCGTCAGCGCAAGC.GGGGGGCAGCG.CCCCCCGAAGCGCCTGAGGCGCGGCGGGCTTTCCCGCCCCACGAAGGCGCAACCACCtt                                           >  2:5061798/1‑101 (MQ=255)
                                                           gcggcgggggggTGGCGCGCCGGCGCAGGG.GGGGGCCAGC.CCCCCCCGAAGCGCCTGAGGCGCGGCGGGCTTGCCCGCCCAACGAAGGCGCAACCACCTTc                                          <  2:330437/95‑1 (MQ=21)
                                                           aaCGCCGGGGGGTTGCGCGTCAGCGCAAGC.GGGGGGCAGCG.CCCCCCGAAGCGCCTGAGGCGCGGCGGGCTTGCCCGCCCAACGAAGGCGCAACCACCTTc                                          >  2:2809701/1‑101 (MQ=255)
                                                           aaCGCCGGGGGGTGGCGCGTCAGCGCAAGC.GGGGGGCAGCG.CCCCCCGAAGCGCCTGAGGCGCGGCGGGCTTGCCCGCCCAACGAAGGCGCAACCACCTTc                                          <  2:3861751/101‑1 (MQ=255)
                                                             cGCCGGGGGGTTGCGCGTCAGCGCAAGC.GGGGGGCAGCG.CCCCCCGAAGCGCCTGAGGCGCGGCGGGCTTGCCCCCCCAACGAAGGCGCAACCACCTTCAg                                        >  2:876274/1‑101 (MQ=255)
                                                              gCCGGGGGGTTGCGCGTCAGCGCAAGC.GGGGGGCAGCG.CCCCCCGAAGCGCCTGAGGCGCGGCGGGCTTGCCCGCCCAACGAAGGCGCAACCACCTTCAGt                                       >  2:5532521/1‑101 (MQ=255)
                                                               ccGGGGGGTTGCGCGTCAGCGCAAGC.GGGGGGCAGCG.CCCCCCGAAGCGCCTGAGGCGCGGCGGGCTTTCCCGCCCAACGAAGGCGCAACCCCCCTCAGTc                                      >  1:4407981/1‑101 (MQ=25)
                                                               ccGGGGGGTTGCGCGTCAGCGCAAGC.GGGGGGCAGCG.CCCCCCGAAGCGCCTGAGGCGCGGCGGGCTTGCCCGCCCAACGAAGGGGCAACCAACTTCAGTc                                      >  2:1378391/1‑101 (MQ=25)
                                                               ccGGGGGGTTGCGCGTCAGCGCAAGC.GGGGGGCAGCG.CCCCCCGAAGCGCCTGAGGCGCGGCGGGCTTGCCCGCCCAACGAAGGCGCAACCCCCTTCAGTc                                      >  2:2198293/1‑101 (MQ=255)
                                                               ccGGGGGGTTGCGCGTCAGCGCAAGC.GGGGGGCAGCG.CCCCCCGAAGCGCCTGAGGCGCGGCGGGCTTGCCCGCCCAACGAAGGCGCAACCACCTTCAGTc                                      <  2:1676294/101‑1 (MQ=255)
                                                                gggggggTGGCGCGTCAGCGCAAGG.GGGGGCCAGCG.CCCCCCGAAGCGCCTGAGGCGCGGCGGGCTTGCCCGCCCAACGAAGGCGCAACCACCTTCAGTCg                                     <  2:3156330/100‑1 (MQ=37)
                                                                cGGGGGGTTGCGCGTCAGCGCAAGC.GGGGGGCAGCG.CCCCCCGAAGCGCCTGAGGCGCGGCGGGCTTGCCCGCCCAACGAAGGCGCAACCACCTTCAGTCg                                     >  1:383079/1‑101 (MQ=255)
                                                                cGGGGGGTTGCGCGTCAGCGCAAGC.GGGGGGCAGCG.CCCCCCGAAGCGCCTGAGGCGCGGCGGGCTTGCCCGCCCAACGAAGGCGCAACCACCTTCAGTCg                                     >  2:2695955/1‑101 (MQ=255)
                                                                  gggggTTGCGCGTCAGCGCAAGC.GGGGGGCAGCG.CCCCCCGAAGCGCCTGAGGCGCGGCGGGCTTGCCCGCCCAACGAAGGCGCAACCCCCTTCAGTCGGc                                   >  1:465442/1‑101 (MQ=25)
                                                                   ggggTGGCGCGTCAGCGCAAGC.GGGGGGCAGC.CCCCCCCGAAGCGCCTGAGGCGCGGCGGGCTTGCCCGCCCAACGAAGGCGCAACCACCTTCAGTCGGCt                                  <  2:347708/101‑1 (MQ=255)
                                                                    gggTTGCGCGTCAGCGCAAGC.GGGGGGCAGCG.CCCCCCGAAGCGCCTGAGGCGCGGCGGGCTTGCCCGCCCAACGAAGGCGCAACCCCCTTCAGTCGGCtt                                 >  2:4624191/1‑101 (MQ=255)
                                                                    gggTTGCGCGTCAGCGCAAGC.GGGGGGCAGCG.CCCCCCGAAGCGCCTGAGGCGCGGCGGGCTTGCCCGCCCAACGAAGGCGCAACCACCTTCAGTCGGCtt                                 >  1:4564117/1‑101 (MQ=255)
                                                                     gtttGCGCGTCAGCCCAGGC.GGGGGGC.GCCCCCCCCCGAAGCGCCTGAGGCGCGGCGGGCTTGCCCGCCCAACGAAGGCGCAACCACCTTCAGTCGGCTTc                                <  1:3811236/99‑1 (MQ=25)
                                                                     ggTTGCCCGCCAGCGCAGGC.GGGGGCCAGC.CCCCCCCGAAGCGCCTGAGGCGCGGCGGGCTTGCCCGCCCAACGAAGGCGCAACCACCTTCAGTCGGCTTc                                <  2:5282864/101‑1 (MQ=25)
                                                                     ggTGGCGCGTCAGCGCAGGC.GGGGGGCAGCG.CCCCCCGAAGCGCCTGAGGCGCGGCGGGCTTGCCCGCCCAACGAAGGCGCAACCACCTTCAGTCGGCTTc                                <  2:3765942/101‑1 (MQ=25)
                                                                      gTTGCGCGTCAGCGCAAGC.GGGGGGCAGCG.CCCCCCGAAGCGCCTGAGGCGCGGCGGGCTTGCCCGCCCAACGAAGGCGCAACCACCTTCAGTCGGCTTCg                               >  1:2882894/1‑101 (MQ=255)
                                                                      gTTGCGCGGCAGCGAAAGC.GGGGGGCAGCG.CCCCCGGAAGCGCCTGAGGCGCGGCGGGCTTGCCCGCCCAACGAAGGCGCAACCACCTTCAGTCGGCTTCg                               <  2:710651/101‑1 (MQ=25)
                                                                       ttGCGCGTCAGCGCAAGC.GGGGGGCAGCG.CCCCCCGAAGCGCCTGAGGCGCGGCGGGCTTGCCCGCCCAACGAAGGCGCAACCACCTTCAGTCGGCTTCGa                              >  1:5053226/1‑101 (MQ=255)
                                                                        tGCGCGTCAGCGCAAGC.GGGGGGCAGCG.CCCCCCGAAGCGCCTGAGGCGCGGCGGGCTTGCCCGCCCAACGAAGGCGCAACCACCTTCAGTCGGCTTCGAt                             >  2:2380413/1‑101 (MQ=255)
                                                                         gcgcGTGCGCGCAGGC.GGGGCGCAGCG.CCCCGCGAAGCGCCTGAGGCGCGGCGGGCTTGCCCGCCCAACGAAGGCGCAACCACCTTCAGTCGGCTTCGATc                            <  1:1475168/101‑1 (MQ=21)
                                                                           gcgTCAGCGCAAGC.GGGGGGCAGCG.CCCCCCGAAGCGCCTGAGGCGCGGCGGGCTTGCCCGCCCAACGAAGGCGCAACCACCTTCAGGCGGCTTCGATccc                          >  1:1069817/1‑101 (MQ=25)
                                                                            cgTCAGCGCAAGC.GGGGGGCAGCG.CCCCCCGAAGCGCCTGAGGCGCGGCGGGCTTGCCCGCCCAACGAAGGCGCAACCACCTTCAGTCGGCTTCGATcccc                         >  1:2340910/1‑101 (MQ=255)
                                                                             gtccgggaaggg.ggggggGGCCCCCCCCCCGAAGCGCCTGAGGCGCGGCGGGCTTGCCCGCCCAACGAAGGCGCAACCACCTTCAGTCGGCTTCGATcccc                         <  1:2575527/91‑1 (MQ=21)
                                                                             gTCAGCGCAAGC.GGGGGGCAGCG.CCCCCCGAAGCGCCTGAGGCGCGGCGGGCTTGCCCGCCCCACGAAGGCGCAACCACCTTCAGCCCGCTTCGATccccc                        >  1:4392300/1‑100 (MQ=25)
                                                                             gTCAGCGCAAGC.GGGGGGCAGCG.CCCCCCGAAGCGCCTGAGGCGCGGCGGGCTTGCCCGCCCAACGAAGGCGCAACCACCTTCAGTCGGCTTCGATccccc                        >  1:2687429/1‑100 (MQ=255)
                                                                             gTCAGCGCAAGC.GGGGGGCAGCG.CCCCCCGAAGCGCCTGAGGCGCGGCGGGCTTGCCCGCCCAACGAAGGCGCAACCACCTTCAGTCGGCTTCGATCCcca                        >  2:3249893/1‑101 (MQ=255)
                                                                              tCAGCGCAAGC.GGGGGGCAGCG.CCCCCCGAAGCGCCTGAGGCGCGGCGGGCTTGCCCGCCCAACGAAGGCGCAACCACCTTCAGTCGGCTTCGATcccccc                       >  2:1276644/1‑99 (MQ=255)
                                                                              tCAGCGCAAGC.GGGGGGCAGCG.CCCCCCGAAGCGCCTGAGGCGCGGCGGGCTTGCCCGCCCAACGAAGGCGCAACCACCTTCAGTCGGCTTCGATCCccac                       >  2:5581116/1‑101 (MQ=255)
                                                                                ggcgcACGG.GGGGGGCAGC.CCCCCCCGAAGCGCCTGAGGCGCGGCGGGCTTGCCCGCCCAACGAAGGCGCAACCACCTTCAGTCGGCTTCGATCCccacca                     <  2:5699412/100‑1 (MQ=255)
                                                                                gccgcAGGG.GGGGGGCCCCCCCCCCCCGAAGCGCCTGAGGCGCGGCGGGCTTGCCCGCCCAACGAAGGCGCAACCACCTTCAGTCGGCTTCGATCCccacc                      <  1:2259456/99‑1 (MQ=37)
                                                                                 gcgcAAGC.GGGGGGCAGCG.CCCCCCGAAGCGCCTGAGGCGCGGCGGGGTTGCCCTCCCAACGAAGGCGGCACAACCTTTAGTCGGCCTCGATCCCCACCAg                    >  1:5792060/1‑101 (MQ=21)
                                                                                 gcgcAAGC.GGGGGGCAGCG.CCCCCCGAAGCGCCTGAGGCGCGGCGGGCTTGCCCGCCCAACGAAGGCGCAACCACCTTCAGTCGGCTTCGATCCCCACCAg                    >  1:380261/1‑101 (MQ=255)
                                                                                 agccgaGG.GGGGGGGGGCCCCCCCCCGAAGCGCCTGAGGCGCGGCGGGCTTGCCCGCCCAACGAAGGCGCAACCACCTTCAGTCGGCTTCGATCCccacca                     <  1:867329/96‑1 (MQ=25)
                                                                                  cgcAAGC.GGGGGGCAGC.CCCCCCCGAAGCGCCTGAGGCGCGGCGGGCTTGCCCGCCCAACGAAGGCGCAACCACCTTCAGTCGGCTTCGATCCCCACCAgc                   <  2:5321262/101‑1 (MQ=255)
                                                                                    cAAGG.GGGGGGCCGCCCCCCCCGGAAGCGCCTGAGGCGCGGCGGGCTTGCCCGCCCAACGAAGGCGCAACCACCTTCAGTCGGCTTCGATCCCCACCAgcg                  <  1:5155266/101‑1 (MQ=255)
                                                                                    cAAGC.GGGGGGCAGCG.CCCCCGGAAGCGCCTGAGGCGCGGCGGGCTTGCCCGCCCAACGAAGGCGCAACCACCTTCAGTCGGCTTCGATCCCCACCAgcgc                 >  1:3285136/1‑101 (MQ=255)
                                                                                     aaGC.GGGGGGCAGCG.CCCCCCGAAGCGCCTGAGGCGCGGCGGGCTTGCCCGCCCAACCAAGGCGCAACCACCCTCAGTCGGCTTCGATCCCCCCCACcgcg                >  1:1957812/1‑101 (MQ=25)
                                                                                       gg.ggggggCCGC.CCCCCCGGAAGCGCCTGAGGCGCGGCGGGCTTGCCCGCCCAACGAAGGCGCAACCACCTTCAGTCGGCTTCGATCCCCACCAGCgcggc              <  2:4306562/101‑1 (MQ=255)
                                                                                       agcggggggCAGCG.CCCCCCGAAGCGCCTGAGGCGCGGCGGGCTTGCCCGCCCAACGAAGGCGCAACCACCTTCAGTCGGCTTCGATCCCCACCAGCgcgg               >  2:1802175/4‑101 (MQ=255)
                                                                                         cggggggCAGCG.CCCCCCGAAGCGCCTGAGGCGCGGCGGGCTTGCCCGCCCAACGAAGGCGCAACCACCTTCAGTCGGCTTCGATCCCCACCAGCgcggcg             >  2:3407760/2‑101 (MQ=255)
                                                                                          ggggggCAGCG.CCCCCCGAAGCGCCTGAGGCGCGGCGGGCTTGCCCGCCCAACGAAGGCGCAACCACCTTCAGTCGGCTTCGATCCCCACCAGCgcggcgg            >  2:1974169/1‑101 (MQ=255)
                                                                                          ggggggCAGCG.CCCCCCGAAGCGCCTGAGGCGCGGCGGGCTTGCCCGCCCAACGAAGCCGCAACCACCTTCAGTCGGCTTCGATCCCCACCAGCgcggcgg            >  2:4793415/1‑101 (MQ=255)
                                                                                           gggggCAGCG.CCCCCCGAAGCGCCTGAGGCGCGGCGGGCTTGCCCGCCCAACGAAGGCGCAACCACCTTCAGTCGGCTTCGATCCCCACCAGCgcggcggc           >  1:339763/1‑101 (MQ=255)
                                                                                            ggggCAGC.CCCCCCCGAAGCGCCTGAGGCGCGGCGGGCTTGCCCGCCCAACGAAGGCGCAACCACCTTCAGTCGGCTTCGATCCCCACCAGCgcggcggcg          <  1:2786101/101‑1 (MQ=255)
                                                                                                cAGCG.CCCCCCGAAGCGCCTGAGGCGCGGCGGGCTTGCCCGCCCAACGAAGGCGCAACCACCTTCAGTCGGCTTCGATCCCCACCAGCGCGGCGGCGGAtt      >  1:1305728/1‑101 (MQ=255)
                                                                                                cAGCG.CCCCCCGAAGCGCCTGAGGCGCGGCGGGCTTGCCCGCCCAACGAAGGCGCAACCACCTTCAGTCGGCTTCGATCCCCACCAGCGCGGCGGCGAAtt      >  1:4351288/1‑101 (MQ=255)
                                                                                                  ggcccccccgGAAGCGCCTGAGGCGCGGCGGGCTTGCCCGCCCAACGAAGGCGCAACCACCTTCAGTCGGCTTCGATCCCCACCAGCGCGGCGGCGAATTc     <  2:5113112/98‑1 (MQ=255)
                                                                                                     ccccccgGAAGCGCCTGAGGCGCGGCGGGCTTGCCCGCCCAACGAAGGCGCAACCACCTTCAGTCGGCTTCGATCCCCACCAGCGCGGCGGCGAATTCCTc  <  1:3505971/101‑1 (MQ=255)
                                                                                                     cccccccGAAGCGCCTGAGGCGCGGCGGGCTTGCCCGCCCAACGAAGGCGCAACCACCTTCAGTCGGCTTCGATCCCCACCAGCGCGGCGGCGAATTCCTc  <  2:1694159/101‑1 (MQ=255)
                                                                                                     |                                                                                                    
CGGTTTCCTGACAGGATTGTGTCCTTCGCGCCGCAGCAGGCCCGCTTTTCCTTGGGCGAAACGCCGGGGGGTTGCGCGTCAGCGCAAGG.GGGGGGCAGCGCCCCCCCGAAGCGCCTGAGGCGCGGCGGGCTTGCCCGCCCAACGAAGGCGCAACCACCTTCAGTCGGCTTCGATCCCCACCAGCGCGGCGGCGAATTCCTC  >  NC_014034/2559835‑2560035

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