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% 49.9 118intergenic (‑211/+59)ispZ/ftsYintracellular septation protein A/cell division protein FtsY
Reads supporting (aligned to +/- strand):  new base (89/12):  ref base (1/16):  total (90/28)
Fisher's exact test strand distribution p-value = 1.97e-11
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 5.46e-06

Only 100 of 176 total aligned reads displayed.

GGTTTCCTGACAGGATTGTGTCCTTCGCGCCGCAGCAGGCCCGCTTTTCCTTGGGCGAAA.C.GCCGGGGGGTTGCGCGTCAGCGCAAGGGGGGGGCAGCGCCCCCCC.GAAGCGCCTGAGGCGCGGCGGGCTTGCCCGCCCAACGAAGGCGCAACCACCTTCAGTCGGCTTCGATCCCCACCAGCGCGGCGGCGAATTCCTC  >  NC_014034/2559836‑2560035
                                                                                                     |                                                                                                     
ggTTTCCTGACAGGATTGTGTCCTTCGCGCCGCAGCAGGCCCGCTTTTCCTTGGGCGAAA.C.GCCGGGGGGTTGCGCGTCAGCGCAAGCGGGGGGCAGCGcc                                                                                                      <  2:5027845/101‑1 (MQ=255)
ggTTTCCTGACAGGATTGTGTCCTTCGCGCCGCAGCAGGCCCGCTTTTCCTTGGGCGAAA.C.GCCGGGGGGGTGCGCGTCAGCGCAAGCGGGGGGCAGAGcc                                                                                                      >  2:2806276/1‑101 (MQ=255)
 gTTTCCTGACAGGATTGTGTCCTTCGCGCCGCAGCAGGCCCGCTTTTCCTTGGGCGAAA.C.GCCGGGGGGTTGCGCGTCAGCGCAAGCGGGGGGCAGCGccc                                                                                                     >  2:4477573/1‑101 (MQ=255)
   ttCCTGACAGGATTGTGTCCTTCGCGCCGCAGCAGGCCCGCTTTTCCTTGGGCGAAA.C.GCCGGGGGGTTGCGCGTCAGCGCAAGCGGGGGGCAGCGccccc                                                                                                   <  1:736028/101‑1 (MQ=255)
   ttCCTGACAGGATTGTGTCCTTCGCGCCGCAGCAGGCCCGCTTTTCCTTGGGCGAAA.C.GCCGGGGGGTTGCGCGTCAGCGCAAGCGGGGGGCAGCGccccc                                                                                                   >  1:977499/1‑101 (MQ=255)
   ttCCTGACAGGATTGTGTCCTTCGCGCCGCAGCAGGCCCGCTTTTCCTTGGGCGAAA.C.GCCGGGGGGTTGCGCGTCAGCGCAAGCGGGGGGCAGCGccccc                                                                                                   >  1:1518147/1‑101 (MQ=255)
   ttCCTGACAGGATTGTGTCCTTCGCGCCGCAGCAGGCCCGCTTTTCCTTGGGCGAAA.C.GCCGGGGGGTTGCGCGTCAGCGCAAGCGGGGGGCAGCGccccc                                                                                                   >  2:1208084/1‑101 (MQ=255)
   ttCCTGACAGGATTGTGTCCTTCGCGCCGCAGCAGGCCCGCTTTTCCTTGGGCGAAA.C.GCCGGGGGGTTGCGCGTCAGCGCAAGCGGGGGGCAGCGccccc                                                                                                   <  1:3398114/101‑1 (MQ=255)
   ttCCAGACAGGATTGTGTCCTTCGCGCCGCAGCAGGCCCGCTTTTCCTTGGGCGAAA.C.GCCGGGGGGTTGCGCGTCAGCGCAAGCGGGGGGCAGCGccccc                                                                                                   <  2:1872984/101‑1 (MQ=255)
      cTGACAGGATTGTGTCCTTCGCGCCGCAGCAGGCCCGCTTTTCCTTGGGCGAAA.C.GCCGGGGGGTTGCGCGTCAGCGCAAGCGGGGGGCAGCGccccccga                                                                                                >  1:1404127/1‑99 (MQ=255)
        gACAGGATTGTGTCCTTCGCGCCGCAGCAGGCCCGCTTTTCCTTGGGCGAAA.C.GCCGGGGGGTTGCGCGTCAGCGCAAGCGGGGGGCAGCGccccccgaag                                                                                              >  2:566707/1‑97 (MQ=255)
        gACAGGATTGTGTCCTTCGCGCCGCAGCAGGCCCGCTTTTCCTTGGGCGAAA.C.GCCGGGGGGTTGCGCGTCAGCGCAAGCGGGGGGCAGCGccccccgaag                                                                                              >  1:5004939/1‑97 (MQ=255)
         aCAGGATTGTGTCCTTCGCGCCGCAGCAGGCCCGCTTTTCCTTGGGCGAAA.C.GCCGGGGGGTTGCGCGTCAGCGCAAGCGGGGGGCAGCG.CCCCCC.GAAgc                                                                                           >  1:70303/1‑101 (MQ=255)
         aCAGGATTGTGTCCTTCGCGCCGCAGCAGGCCCGCTTTTCCTTGGGCGAAA.C.GCCGGGGGGTTGCGCGTCAGCGCAAGCGGGGGGCAGCG.CCCCCC.GAAgc                                                                                           >  2:3605131/1‑101 (MQ=255)
          cAGGATTGTGTCCTTCGCGCCGCAGCAGGCCCGCTTTTCCTTGGGCGAAA.C.GCCGGGGGGTTGCGCGTCAGCGCAAGCGGGGGGCAGCG.CCCCCC.GAAgcg                                                                                          >  1:4850989/1‑101 (MQ=255)
           aGGATTGTGTCCTTCGCGCCGCAGCAGGCCCGCTTTTCCTTGGGCGAAA.C.GCCGGGGGGTTGCGCGTCAGCGCAAGCGGGGGGCAGCG.CCCCCC.GAAgcgc                                                                                         >  1:3054635/1‑101 (MQ=255)
           aGGATTGTGTCCTTCGCGCCGCAGCAGGCCCGCTTTTCCTTGGGCGAAA.C.GCCGGGGGGTTGCGCGTCAGCGCAAGCGGGGGGCAGCG.CCCCCC.GAAgcgc                                                                                         >  2:3076008/1‑101 (MQ=255)
           aGGATTGTGTCCTTCGCGCCGCAGCAGGCCCGCTTTTCCTTGGGCGAAA.C.GCCGGGGGGTTGCGCGTCAGCGCAAGCGGGGGGCAGCG.CCCCCC.GAAgcgc                                                                                         >  1:100287/1‑101 (MQ=255)
             gATTGTGTCCTTCGCGCCGCAGCAGGCCCGCTTTTCCTTGGGCGAAA.C.GCCGGGGGGTTGCGCGTCAGCGCAAGCGGGGGGCAGCG.CCCCCC.GAAGCGCCt                                                                                       >  1:1560507/1‑101 (MQ=255)
                  tgtCCTTCGCTCCGCAGCAGGCCCGCTTTTCCTTGGGCGAAA.C.GCCGGGGGGTTGCGCGTCAGCGCAAGCGGGGGGCAGCG.CCCCCC.GAAGCGCCTGAGgc                                                                                  >  1:3123975/1‑101 (MQ=255)
                  tgtCCTTCGCGCCGCAGCAGGCCCGCTTTTCCTTGGGCGAAA.C.GCCGGGGGGTTGCGCGTCAGCGCAAGCGGGGGGCAGCG.CCCCCC.GAAGCGCCTGAGgc                                                                                  >  2:769773/1‑101 (MQ=255)
                    tCCTTCGCGCCGCAGCAGGCCCGCTTTTCCTTGGGCGAAA.C.GCCGGGGGGTTGCGCGTCAGCGCAAGCGGGGGGCAGCG.CCCCCC.GAAGCGCCTGAGgcgc                                                                                >  2:1999556/1‑101 (MQ=255)
                     ccTTCGCGCCGCAGCAGGCCCGCTTTTCCTTGGGCGAAA.C.GCCGGGGGGTTGCGCGTCAGCGCAAGCGGGGGGCAGCG.CCCCCC.GAAGCGCCTGAGGcgcg                                                                               >  1:1474774/1‑101 (MQ=255)
                     ccTTCGCGCCGCAGCAGGCCCGCTTTTCCTTGGGCGAAA.C.GCCGGGGGGTTGCGCGTCAGCGCAAGCGGGGGGCAACG.CCCCCC.GAAGCGCCGGAGGcgcg                                                                               >  2:2412604/1‑101 (MQ=255)
                      cTTCGCGCCGCAGCAGGCCCGCTTTTCCTTGGGCGAAA.C.GCCGGGGGGTTGCGCGTCAGCGCAAGCGGGGGGCAGCG.CCCCCC.GAAGCGCCTGAGGCgcgg                                                                              >  1:35279/1‑101 (MQ=255)
                         cgcgCCGCAGCAGGCCCGCTTTTCCTTGGGCGAAA.C.GCCGGGGGGTTGCGCGTCAGCGCAAGCGGGGGGCAGCGCCCCCCC..AAGCGCCTGCGGCGcggcgg                                                                           >  1:1794232/1‑101 (MQ=255)
                           cgccgcAGCGGGCCCGCTTTTCCTTGGGCGAAA.C.GCCGGGGGGTGGCGCGTCAGCGCAAGCGGGGGGCAGCG.CCCCCC.GAAGCGCCTGAGGCGCGGCGGGc                                                                         <  2:681952/101‑1 (MQ=39)
                            gccgcAGCAGGCCCGCTTTTCCTTGGGCGAAA.C.GCCGGGGGGTGGCGAGTCAGCGCAAGCGGGGCGCAGCG.CCCCCC.GAAGCGCCTGAGGCGCGGGGGGcg                                                                        >  2:3185781/1‑100 (MQ=255)
                             ccgcAGCGGGCCCCTTTTTCCGGGGGGGAAA.CGGCGGGGGGGTGGGGCGCCGCCGAAGGGGGGGGGCGGC.CCCCCCC.GAAGCGCCTGAGGCGCGGCGGGCt                                                                        <  2:3727280/101‑1 (MQ=255)
                                cagcagGCCCGCTTTTCCTTGGGCGAAA.C.GCCGGGGGGTTGCGCGTCAGCGCAAGCGGGGGGCAGCG.CCCCCC.GAAGCGCCTGAGGCGCGGCGGGCTTggc                                                                    >  1:2166461/1‑99 (MQ=255)
                                   cagGCCCGCTTTTCCTTGGGCGAAA.C.GCCGGGGGGTTGCGCGTCAGCGCAAGCGGGGGGCAGCG.CCCCCC.GAAGCGCCTGAGGCGCGGCGGGCTTgcccgc                                                                 >  1:4300450/1‑101 (MQ=255)
                                    agGCCCGCTTTTCCTTGGGCGAAA.C.GCCGGGGGGTTGCGCGTCAGCGCAAGCGGGGGGCAGCG.CCCCCC.GAAGCGCCTGAGGCGCGGCGGGCTTgcccgcc                                                                >  1:281988/1‑101 (MQ=255)
                                    agGCCCGCTTTTCCTTGGGCGAAA.C.GCCGGGGGGTTGCGCGTCAGCGCAAGCGGGGGGCAGCG.CCCCCC.GAAGCGCCTGAGGCGCGGCGGGCTTgcccgcc                                                                >  2:3364748/1‑101 (MQ=255)
                                    agGCCCGCTTTTCCTTGGGCGAAA.C.GCCGGGGGGTTGCGCGTCAGCGCAAGCGGGGGGCAGCG.CCCCCC.GAAGCGCCTGAGGCGCGGCGGGCTTgcccgcc                                                                >  1:5037168/1‑101 (MQ=255)
                                     gcccgttttttCTTTGGGCAAAC.C.GCGGGGGGGTGGCGGGCCGGCG.AAGGGGGGGGCCGCCCCCCCCG.GAAGCGCCTGAGGCGCGGCGGGCTTgcccgccc                                                               <  2:4601188/94‑1 (MQ=16)
                                      gCCCGTTTTTCCTTGGGCGAAA.C.GCGGGGGGGTTGCGCGTCGGCGCAAGGGGGGGGCAGC.CCCCCCC.GAAGCGCCTGAGGCGCGGCGGGCTTGCCCGCCCa                                                              <  2:1540604/101‑1 (MQ=37)
                                        ccGTTTTTCCTTGGGCGAAC.C.GCGGGGGGGTTGCGCGTCAGCGCAAGCGGGGGGCAGCG.CCCCCC.GAAGCGCCTGAGGCGCGGCGGGCTTGCCCGCCCAAc                                                            <  1:1475666/101‑1 (MQ=25)
                                        ccGTTTTTCCTGGGGCGAAG.C.GCGGGGGGGTTGCGCGGCAGCGCAGGCGGGGGCCAGC.CCCCCCC.GAAGCGCCTGAGGCGCGGCGGGCTTGCCCGCCCAAc                                                            <  1:1629381/101‑1 (MQ=25)
                                         cGCTTTTCCTTGGGCGAAA.C.GCCGGGGGGTTGCGCGTCAGCGCAAGCGGGGGGCAGCG.CCCCCC.GAAGCGCCTGAGGCGCGGCGGGCTTGCCCGCCCAACg                                                           >  2:45247/1‑101 (MQ=255)
                                          gtgtttCCTTGGGCGAAC.C.GCCGGGGGTTTGCGCGTCAGCGAAGGCGGGGGGCACCG.CCCCCC.GAAGCGCCTGAGGCGCGGCGGGCTTGCCCGCCCAACGa                                                          <  2:1463443/98‑1 (MQ=25)
                                          gCTTTTCCTTGGGCGAAA.C.GCCGGGGGGTTGCGCGTCAGCGCAAGCGGGGGGCAGCG.CCCCCC.GAAGCGCCTGAGGCGCGGCGGGCTTGCCCGCCCAACGa                                                          >  1:1289550/1‑101 (MQ=255)
                                           cTTTTCCTTGGGCGAAG.C.GCGGGGGGGTTGCGCTTCAGCGCAAGCGGGGGGCAGCG.CCCCCC.GAAGCGCCTGAGGCGCGGCGGGCTTGCCCGCCCAACGaa                                                         <  2:4211037/101‑1 (MQ=25)
                                           cTTTTCCTTGGGCGAAA.C.GCCGGGGGGTTGCGCGTCAGCGCAAGCGGGGGGCAGCG.CCCCCC.GAAGCGCCTGAGGCGCGGCGGGGCTGCCCGCCCCACGaa                                                         >  2:3515088/1‑101 (MQ=255)
                                             tttCCTTGGGCGAAA.C.GCCGGGGGGTTGCGCGTCAGCGCAAGCGGGGGGCAGCG.CCCCCC.GAAGCGCCTGAGGCGCGGCGGGCTTGCCCGCCCCACGCAgg                                                       >  1:3820394/1‑101 (MQ=37)
                                             tttCCTTGGGCGAAA.C.GCCGGGGGGTTGCGCGTCAGCGCAAGCGGGGGGCAGCG.CCCCCC.GAAGCGCCTGAGGCGCGGCGGGCTTGCCCGCCCAACCaaag                                                       >  1:2680281/1‑99 (MQ=37)
                                              ttCCTTGGGCGAAA.C.GCCGGGGGGTTGCGCGTCAGCGCAAGCGGGGGGCAGCG.CCCCCC.GAAGCGCCTGAGGCGCGGCGGGGTTGACCGCCCAACGAAggg                                                      >  1:228711/1‑100 (MQ=255)
                                              ttCCTTGGGCGAAA.C.GCCGGGGGGTTGCGCGTCAGCGCAAGCGGGGGGCAGCG.CCCCCC.GAAGCGCCTGAGGCGCGGCGGGCTTGCGCGCCCCACGCAggg                                                      >  1:1086905/1‑100 (MQ=255)
                                              ttCCTTGGGCGAAA.C.GCCGGGGGGTTGCGCGTCAGCGCAAGCGGGGGGCAGCG.CCCCCC.GAAGCGCCTGAGGCGCGGCGGGCTTGCCCGCCCAACGAAGgc                                                      >  1:3623211/1‑101 (MQ=255)
                                              ttCCTTGGGCGAAA.C.GCCGGGGGGTTGCGCGTCAGCGCAAGCGGGGGGCAGCG.CCCCCC.GAAGCGCCTGAGGCGCGGCGGGCTTGCCCGCCCAACGAAGgc                                                      >  1:2063718/1‑101 (MQ=255)
                                              tccggggggAGAAC.C.GCGGGGGGGTGGGGGGCCAGCGAAGGGGGGGGGCCGC.CCCCCCC.GAAGCGCCTGAGGCGCGGCGGGCTTGCCCGCCCAACGAAGgc                                                      <  2:4684183/95‑1 (MQ=16)
                                               tCCTGGGGCAAAACC.GCCGGGGGGTGGCGCGTGAGCGCAGGGGGGGGGCAGCG.CCCCCC.GAAGCGCCTGAGGCGCGGCGGGCTTGCCCGCCCAACGAAGgc                                                      <  2:562188/101‑1 (MQ=21)
                                                ccTGGGGGGAAA.C.GCGGGGGGGTTGCGCGTCAGCGCAGGGGGGGGGCAGCG.CCCCCC.GAAGCGCCTGAGGCGCGGCGGGCTTGCCCGCCCAACGAAGgcgc                                                    <  2:794336/101‑1 (MQ=25)
                                                 cTTGGGCGAAA.C.GCCGGGGGGTTGTGCGTCAGCGCAAGCGGGGGGCAGCG.CCCCCC.GAAGCGCCTGAGGCGCGGCGGGCCTGCCCGCCCAACCAAAGCGCa                                                   >  1:1748331/1‑101 (MQ=255)
                                                 cTTGGGCGAAA.C.GCCGGGGGGTTGCGCGTCAGCGCAAGCGGGGGGCAGCG.CCCCCC.GAAGCGCCTGAGGCGCGGCGGGCTTGCCCGCCCAACGAAGGCGCa                                                   >  2:4894539/1‑101 (MQ=255)
                                                 cTTGGGCGAAA.C.GCCGGGGGGTTGCGCGTCAGCGCAAGCGGGGGGCAGCG.CCCCCC.GAAGCGCCTGAGGCGCGGCGGGCTTGCCCGCCCAACGAAGGCGCa                                                   >  1:2853944/1‑101 (MQ=255)
                                                 cTTGGGCGAAA.C.GCCGGGGGGTTGCGCGTCAGCGCAAGCGGGGGGCAGCG.CCCCCC.GAAGCGCCTGAGGCGCGGCGGGCTTGCCCCCCCCACGAAAGCGCa                                                   >  2:1062746/1‑101 (MQ=37)
                                                  ttGGGCGAAA.C.GCCGGGGGGTTGCGCGTCAGCGCAAGCGGGGGGCAGCG.CCCCCC.GAAGCGCCTGAGGCGCGGCGGGCTTGCCCGCCCAACGAAGgcgcca                                                  >  1:4236887/1‑99 (MQ=255)
                                                  ttGGGCGAAA.C.GCCGGGGGGTTGCGCGTCAGCGCAAGCGGGGGGCAGCG.CCCCCC.GAAGCGCCTGAGGCGCGGCGGGCTTGCCCGCCCAACGAAGGCGCaa                                                  >  2:164674/1‑101 (MQ=255)
                                                  ttGGGCGAAA.C.GCCGGGGGGTTGCGCGTCAGCGCAAGCGGGGGGCAGCG.CCCCCC.GAAGCGCCTGAGGCGCGGCGGGCGTGCCCGCCCCACGAAGgcgccc                                                  >  2:151014/1‑99 (MQ=37)
                                                   tGGGCGAAA.C.GCCGGGGGGTTGCGCGTCAGCGCAAGCGGGGGGCAGCG.CCCCCC.GAAGCGCCTGAGGCGCGGCGGGCTTGCCCGCCCAACGAAGGCGCAac                                                 >  1:4013873/1‑101 (MQ=255)
                                                    gggCGAAA.C.GCCGGGGGGTTGCGCGTCAGCGCAAGCGGGGGGCAGCG.CCCCCC.GAAGCGCCTGAGGCGCGGCGGGCTTGCCCGCCCAACGAAGGCGCAacc                                                >  1:2776926/1‑101 (MQ=255)
                                                        gAAA.C.GCCGGGGGGTTGCGCGTCAGCGCAAGCGGGGGGCAGCG.CCCCCC.GAAGCGCCTGAGGCGCGGCGGGCCTTTCCGCCCAACGAAGGCGCAacccact                                            >  2:3078467/1‑97 (MQ=25)
                                                          gaagcgggggggggTGGCGCGTCAGCGCGAGGGGGGGGCAGC.CCCCCCC.GAAGCGCCTGAGGCGCGGCGGGCTTGCCCGCCCAACGAAGGCGCAACCACCt                                            <  1:4094638/93‑1 (MQ=25)
                                                          ga.C.GCCGGGGGGTTCCGCCTCAGCGCAAGCGGGCGCCAGCG.CGCCCC.GAAGCGCCTGAGGCGCGGCGGGCTTGCCCGCCCAACGAAGGCGCAACCACCTTc                                          <  1:5030259/100‑1 (MQ=21)
                                                          aa.C.GCCGGGGGGTTGCGGGCAAGCGCAAGGGGGGGGCAGC.CCCCCCC.GAAGCGCCTGAGGCGCGGCGGGCTTGCCCGCCCAACGAAGGCGCAACCACCTTc                                          <  2:1560377/101‑1 (MQ=25)
                                                          aa.C.GCCGGGGGGTTGCGCGTCAGCGCAAGCGGGGGGCAGCG.CCCCCC.GAAGCGCCTGAGGCGCGGCGGGCTTGCCCGCCCAACGAAGGCGCAACCACCTTc                                          >  2:2102239/1‑101 (MQ=255)
                                                           a.C.GCCGGGGGGTTGCGCGTCAGCGCAAGCGGGGGGCAGCG.CCCCCC.GAAGCGCCTGAGGCGCGGCGGGCTTGCCCGCCCAACGAAGGCGCAACCACCTTCa                                         >  1:365492/1‑101 (MQ=255)
                                                             c.GCGGGGGGGTTGCGCGCCGGCGCAGGGGGGGGGCAGC.CCCCCCC.GAAGCGCCTGAGGCGCGGCGGGCTTGCCCGCCCAACGAAGGCGCAACCACCTTCAg                                        <  2:281954/101‑1 (MQ=25)
                                                              ggggggggtgTGGCGCGCCGGCG.AAGGGGGGGGGCCCCCCCCCCC.GAAGCGCCTGAGGCGCGGCGGGCTTGCCCGCCCATCGAAGGCGCAACCACCTTCAg                                        <  1:2872338/97‑1 (MQ=18)
                                                                    ggggTTGCGCGTCAGCGCAAGCGGGGGGCAGCG.CCCCCC.GAAGCGCCTGAGGCGCGGCGGGCTTGCCCGCCCAACGAAGGCGCAACCCCCCTCAGGCGGCt                                  >  2:41754/1‑101 (MQ=25)
                                                                    ggggTTGCGCGTCAGCGCAAGCGGGGGGCAGCG.CCCCCC.GAAGCGCCTGAGGCGCGGCGGGCTTGCCCGCCCAACGAAGGCGCAACCACCTTCAGTCGGCt                                  >  1:323272/1‑101 (MQ=255)
                                                                     gggTTGCGCGTCGGCGCAGGGGGGGGGCGGCG.CCCCCC.GAAGCGCCTGAGGCGCGGCGGGCTTGCCCGCCCAACGAAGGCGCAACCACCTTCAGTCGGCtt                                 <  1:3674185/101‑1 (MQ=25)
                                                                     gggTTGCGCGTCAGCGCAAGCGGGGGGCAGCG.CCCCCC.GAAGCGCCTGAGGCGCGGCGGGGTTTCCCGCCCCACGAAGGGGCAACCCCCTTTCGTCGGCtt                                 >  1:2631561/1‑101 (MQ=255)
                                                                     gggTTGCGCGCCAGCGCAAGCGGGGGGCAGC.CCCCCCC.GAAGCGCCTGAGGCGCGGCGGGCTTGCCCGCCCAACGAAGGCGCAACCACCTTCAGTCGGCtt                                 <  2:4122084/101‑1 (MQ=255)
                                                                       gTTGCGCGTCAGCGCAAGCGGGGGGCAGCG.CCCCCC.GAAGCGCCTGAGGCGCGGCGGGCTTGCCCGCCCAACGAAGGCGCAACCACCTTCAGTCGGCTTCg                               >  1:4870274/1‑101 (MQ=255)
                                                                            cgggccgggcAAGGCGGGGGCCCCCCCCCCCC.GAAGCGCCTGAGGCGCGGCGGGCTTGCCCGCCCAACGAAGGCGCAACCACCTTCAGTCGGCTTCGATcc                           <  2:3995331/93‑1 (MQ=25)
                                                                             cgTCAGCGCAAGCGGGGGGCAGCG.CCCCCC.GAAGCGCCTGAGGCGCGGCGGGCTTGCCCGCCCAACGAAGGCGCAACCACCTTCAGTCGGCTTCGATcccc                         >  2:3952195/1‑101 (MQ=255)
                                                                              ggcgcgcgggggggggggCCCCCCCCCCCCGGAAGCGCCTGAGGCGCGGCGGGCTTGCCCGCCCAACGAAGGCGCAACCACCTTCAGTCGGCTTCGATccc                          <  2:4416104/91‑1 (MQ=21)
                                                                              gccagggaaaGGGGGGGGCGGC.CCCCCCC.GAAGCGCCTGAGGCGCGGCGGGCTTGCCCGCCCAACGAAGGCGCAACCACCTTCAGTCGGCTTCGATCCcca                        <  2:2346692/93‑1 (MQ=25)
                                                                              gTCAGCGCAAGCGGGGGGCAGCG.CCCCCC.GAAGCGCCTGAGGCGCGGCGGGCTTGCCCGCCCAACGAAGGCGCAACCACCTTCAGTCGGCTTCGATCCcca                        >  2:236061/1‑101 (MQ=255)
                                                                              gTCAGCGCAAGCGGGGGGCAGCG.CCCCCC.GAAGCGCCTGAGGCGCGGCGGGCTTGCCCGCCCAACGAAGGCGCAACCACCTTCAGTCCGCCTCGATCCcca                        >  2:426122/1‑101 (MQ=25)
                                                                              gTCAGCGCAAGCGGGGGGCAGCG.CCCCCC.GAAGCGCCTGAGGCGCGGCGGGCTTGCCCGCCCAACGAAAGCGCAACCACCTTCACCCGGCTTAGTTCCcca                        >  2:3652926/1‑101 (MQ=21)
                                                                                gcgcgaAAGGGGGGGCCA.CCCCCCCCC.GAAGCGCCTGAGGCGCGGCGGGCTTGCCCGCCCAACGAAGGCGCAACCACCTTCAGTCGGCTTCGATCCccacc                      <  2:4827617/99‑1 (MQ=25)
                                                                                cggggcAAGGGGGGGGCCCGCCCCCCCC.GAAGCGCCTGAGGCGCGGCGGGCTTGCCCGCCCAACGAAGGCGCAACCACCTTCAGTCGGCTTCGATCCccac                       <  2:191899/97‑1 (MQ=25)
                                                                                cAGCGCAAGCGGGGGGCAGCG.CCCCCC.GAAGCGCCTGAGGCGCGGCGGGCTTGCCCGCCCAACGAAGGCGCAACCACCTTCAGTCGGCTTCGATCCccacc                      >  1:3922865/1‑101 (MQ=255)
                                                                                acgggaaggcggggggCCGC.CCCCCCC.GAAGCGCCTGAGGCGCGGCGGGCTTGCCCGCCCAACGAAGGCGCAACCACCTTCAGTCGGCTTCGATCCccacc                      <  2:3231158/91‑1 (MQ=25)
                                                                                 cggcgggggggggggGCGCCCCCCCCC.GAAGCGCCTGAGGCGCGGCGGGCTTGCCCGCCCAACGAAGGCGCAACCACCTTCAGTCGGCTTCGATCCccacc                      <  1:774751/94‑1 (MQ=25)
                                                                                 cgcgtggggggggggCACCCCCCCCAG.GAAGCGCCTGAGGCGCGGCGGGCTTGCCCGCCCAACGAAGGCGCAACCACCTTCAGTCGGCTTCGATCCccacc                      <  2:228677/94‑1 (MQ=25)
                                                                                  agcccaGGGGGGGGCCCCCCCCCCCC.GAAGCGCCTGAGGCGCGGCGGGCTTGCCCGCCCAACGAAGGCGCAACCACCTTCAGTCGGCTTCGATCCccacca                     <  1:539393/96‑1 (MQ=37)
                                                                                   cgcAAGGGGGGGCCAGC.CCCCCCC.GAAGCGCCTGAGGCGCGGCGGGCTTGCCCGCCCAACGAAGGCGCAACCACCTTCAGTCGGCTTCGATCCCCACCAgc                   <  2:1480460/101‑1 (MQ=255)
                                                                                    gcAAGCGGGGGGCAGCG.CCCCCC.GAAGCGCCTGAGGCGCGGCGGGCTTGCCCGCCCAACGAAGGCGCAACCACCTTCAGTCGGCTTCGATCCCCACCAgcg                  >  2:3071031/1‑101 (MQ=255)
                                                                                    gcAAGCGGGGGGCACCG.CCCCCC.GAAGCGCCTGAGGCGCGGCGGGCTTGCCCGCCCAACGAAGGCGCAAACACCTTCAGTCCGNTTCGATTCCCACCAgcg                  >  2:2139689/1‑101 (MQ=25)
                                                                                     cAAGCGGGGGGCAGCG.CCCCCC.GAAGCGCCTGAGGCGCGGCGGGCTTGCCCGCCCAACGAAGGCGCAACCACCTTCAGTCGGCTTCGATCCCCACCAgcgc                 >  2:2861921/1‑101 (MQ=255)
                                                                                     cAAGCGGGGGGCAGCG.CCCCCC.GAAGCGCCTGAGGCGCGGCGGGCCTTACCGCCCAACGAAAGGGCCACCACCTTCAGACCGGTTCGTTCACCACCAacgc                 >  2:3411984/1‑101 (MQ=16)
                                                                                           gggggCAGCG.CCCCCC.GAAGCGCCTGAGGCGCGGCGGGCTTGCCCGCCCAACGAAGGCGCAACCACCTTCAGTCGGCTTCGATCCCCACCAGCgcggcggc           >  1:2416176/1‑101 (MQ=255)
                                                                                            ggggCAGCG.CCCCCC.GAAGCGCCTGAGGCGCGGCGGGCTTGCCCGCCCAACGAAGGCGCAACCACCTTCAGTCGGCTTCGATCCCCACCAGCgcggcggcg          >  2:1753660/1‑101 (MQ=255)
                                                                                            ggggCAGCG.CCCCCC.GAAGCGCCTGAGGCGCGGCGGGCTTGCCCGCCCAACGAAGGCGCAACCACCTTCAGTCGGCTTCGATCCCCACCAGCgcggcggcg          <  1:821013/101‑1 (MQ=255)
                                                                                             gggCAGCG.CCCCCC.GAAGCGCCTGAGGCGCGGCGGGCTTGCCCGCCCAACGAAGGCGCAACCACCTTCAGTCGGCTTCGATCCCCCCCAGCGCGGCGGCGa         >  2:3713135/1‑101 (MQ=255)
                                                                                                 ctgcccccccc.AAAGCGCCTGAGGCGCGGCGGGCTTGCCCGCCCAACGAAGGCGCAACCCCCTTCAGTCGGCTTCGATCCCCACCAGCGCGGCGGCGAAtt      <  1:3320491/97‑1 (MQ=255)
                                                                                                    cccccccc.GAAGCGCCTGAGGCGCGGCGGGCTTGCCCGCCCAACGAAGGCGCAACCACCTTCAGTCGGCTTCGATCCCCACCAGCGCGGCGGCGAATTCCt   <  1:1969926/100‑1 (MQ=255)
                                                                                                     |                                                                                                     
GGTTTCCTGACAGGATTGTGTCCTTCGCGCCGCAGCAGGCCCGCTTTTCCTTGGGCGAAA.C.GCCGGGGGGTTGCGCGTCAGCGCAAGGGGGGGGCAGCGCCCCCCC.GAAGCGCCTGAGGCGCGGCGGGCTTGCCCGCCCAACGAAGGCGCAACCACCTTCAGTCGGCTTCGATCCCCACCAGCGCGGCGGCGAATTCCTC  >  NC_014034/2559836‑2560035

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