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% 21.7 114intergenic (‑211/+59)ispZ/ftsYintracellular septation protein A/cell division protein FtsY
Reads supporting (aligned to +/- strand):  new base (101/4):  ref base (0/9):  total (101/13)
Fisher's exact test strand distribution p-value = 1.10e-10
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 3.57e-06

Only 100 of 158 total aligned reads displayed.

CGGTTTCCTGACAGGATTGTGTCCTTCGCGCCGCAGCAGGCCCGCTTTTCC.TTGGGCGAAACGCCGGGGGGTT.G.CGCG.TCAGCGCAAGGGGGGGGCAGCGCCCCCCCGAAGCGCCTGAGGCGCGGCGGGCTTGCCCGCCCAACGAAGGCGCAACCACCTTCAGTCGGCTTCGATCCCCACCAGCGCGGCGGCGAATTCCTC  >  NC_014034/2559835‑2560035
                                                                                                        |                                                                                                    
cGGTTTCCTGACAGGATTGTGTCCTTCGCGCCGCAGCAGGCCCGCTTTTCC.TTGGGCGAAACGCCGGGGGGTT.G.CGCG.TCAGCGCAAGCGGGGGGCAgcgc                                                                                                      >  1:1788032/1‑101 (MQ=255)
 ggTTTCCTGACAGGATTGTGTCCTTCGCGCCGCAGCAGGCCCGCTTTTCC.TTGGGCGAAACGCCGGGGGGTT.G.CGCG.TCAGCGCAAGCGGGGGGCAGCGcc                                                                                                     >  1:2702325/1‑101 (MQ=255)
 ggTTTCCTGACAGGATGGTGTCCTTCGCGCCGCAGCAGGCCCGCTTTTCC.TTGGGCGAAACGCCGGGGGGTT.G.CGCG.TCAGCGCAAGCGGGGGGCAGCGcc                                                                                                     <  2:1071124/101‑1 (MQ=255)
    ttCCTGACAGGATTGTGTCCTTCGCGCCGCAGCAGGCCCGCTTTTCC.TTGGGCGAAACGCCGGGGGGTT.G.CGCG.TCAGCGCAAGCGGGGGGCAGCGccccc                                                                                                  >  1:4395736/1‑101 (MQ=255)
     tCCTGACAGGATTGTGTCCTTCGCGCCGCAGCAGGCCCGCTTTTCC.TTGGGCGAAACGCCGGGGGGTT.G.CGCG.TCAGCGCAAGCGGGGGGCAGCGcccccc                                                                                                 <  2:2135245/101‑1 (MQ=255)
          aCAGGATTGTGTCCTTCGCGCCGCAGCAGGCCCGCTTTTCC.TTGGGCGAAACGCCGGGGGGTT.G.CGCG.TCAGCGCAAGCGGGGGGCAGCG.CCCCCCGAAgc                                                                                           >  2:1414504/1‑101 (MQ=255)
          aCAGGATTGTGTCCTTCGCGCCGCAGCAGGCCCGCTTTTCC.TTGGGCGAAACGCCGGGGGGTT.G.CGCG.TCAGCGCAAGCGGGGGGCAGCG.CCCCCCGAAgc                                                                                           >  2:2585130/1‑101 (MQ=255)
          aCAGGATTGTGTCCTTCGCGCCGCAGCAGGCCCGCTTTTCC.TTGGGCGAAACGCCGGGGGGTT.G.CGCG.TCAGCGCAAGCGGGGGGCAGCG.CCCCCCGAAgc                                                                                           >  2:2062046/1‑101 (MQ=255)
           cAGGATTGTGTCCTTCGCGCCGCAGCAGGCCCGCTTTTCC.TTGGGCGAAACGCCGGGGGGTT.G.CGCG.TCAGCGCAAGCGGGGGGCAGCG.CCCCCCGAAgcg                                                                                          >  1:2020621/1‑101 (MQ=255)
           cAGGATTGTGTCCTTCGCGCCGCAGCAGGCCCGCTTTTCC.TTGGGCGAAACGCCGGGGGGTT.G.CGCG.TCAGCGCAAGCGGGGGGCAGCG.CCCCCCGAAgcg                                                                                          >  2:3408296/1‑101 (MQ=255)
            aGGATTGTGTCCTTCGCGCCGCAGCAGGCCCGCTTTTCC.TTGGGCGAAACGCCGGGGGGTT.G.CGCG.TCAGCGCAAGCGGGGGGCAGCG.CCCCCCGAAgcgc                                                                                         >  2:1405844/1‑101 (MQ=255)
            aGGATTGTGTCCTTCGCGCCGCAGCAGGCCCGCTTTTCC.TTGGGCGAAACGCCGGGGGGTT.G.CGCG.TCAGCGCAAGCGGGGGGCAGCG.CCCCCCGAAgcgc                                                                                         >  2:2018869/1‑101 (MQ=255)
             ggATTGTGTCCTTCGCGCCGCAGCAGGCCCGCTTTTCC.TTGGGCGAAACGCCGGGGGGTT.G.CGCG.TCAGCGCAAGCGGGGGGCAGCG.CCCCCCGAAGCGcc                                                                                        >  2:3282853/1‑101 (MQ=255)
              gATTGTGTCCTTCGCGCCGCAGCAGGCCCGCTTTTCC.TTGGGCGAAACGCCGGGGGGTT.G.CGCG.TCAGCGCAAGCGGGGGGCAGCG.CCCCCCGAAGCGCCt                                                                                       >  1:2480041/1‑101 (MQ=255)
                ttGTGTCCTTCGCGCCGCAGCAGGCCCGCTTTTCC.TTGGGCGAAACGCCGGGGGGTT.G.CGCG.TCAGCGCAAGCGGGGGGCAGCG.CCCCCCGAAGCGCCTGa                                                                                     >  1:2919459/1‑101 (MQ=255)
                 tgtgTCCTTCGCGCCGCAGCAGGCCCGCTTTTCC.TTGGGCGAAACGCCGGGGGGTT.G.CGCG.TCAGCGCAAGCGGGGGGCGGCG.CCCCCCGAAGCGCCGGAg                                                                                    >  1:4254386/1‑101 (MQ=255)
                 tgtgTCCTTCGCGCCGCAGCAGGCCCGCTTTTCC.TTGGGCGAAACGCCGGGGGGTT.G.CGCG.TCAGCGCAAGCGGGGGGCAGCG.CCCCCCGAAGCGCCTGAg                                                                                    >  2:1457468/1‑101 (MQ=255)
                  gtgtCCTTCGCGCCGCAGCAGGCCCGCTTTTCC.TTGGGCGAAACGCCGGGGGGTT.G.CGCG.TCAGCGCAAGCGGGGGGCAGCG.CCCCCCGAAGCGCCTGAgg                                                                                   >  2:3408338/1‑101 (MQ=255)
                  gtgtCCTTCGCGCCGCAGCAGGCCCGCTTTTCC.TTGGGCGAAACGCCGGGGGGTT.G.CGCG.TCAGCGCAAGCGGGGGGCAGCG.CCCCCCGAAGCGCCTGAgg                                                                                   >  2:863048/1‑101 (MQ=255)
                  gtgtCCTTCGCGCCGCAGCAGGCCCGCTTTTCC.TTGGGCGAAACGCCGGGGGGGT.G.CGCG.TCAGCGCAAGCGGGGGGCAGCG.CCCCCCTAAGCGCCTGAgg                                                                                   >  2:376368/1‑101 (MQ=255)
                  gtgtCCTTCGCGCCGCAGCAGGCCCGCTTTTCC.TTGGGCGAAACGCCAGGGGGTT.G.CGCG.TCAGCGCAAGCGGGGGGCAGCG.CCCCCCGAAGCGCCTGAgg                                                                                   >  2:3857034/1‑101 (MQ=255)
                   tgtCCTTCGCGCCGCAGCAGGCCCGCTTTTCC.TTGGGCGAAACGCCGGGGGGTT.G.CGCG.TCAGCGCAAGCGGGGGGCAGCG.CCCCCCGAAGCGCCTGAGgc                                                                                  >  1:3607436/1‑101 (MQ=255)
                      ccTTCGCGCCGCCGCGGGCCCGCTTTTCC.TTGGGCGAAACGCGGGGGGGTT.G.CGCG.TCAGCGCAAGGGGGGGGCAGCG.CCCCCCGAAGCGCCTGAGGcgcg                                                                               <  2:2017553/101‑1 (MQ=255)
                      ccTTCGCGCCGCAGCAGGCCCGCTTTTCC.TTGGGCGAAACGCCGGGGGGTT.G.CGCG.TCAGCGCAAGCGGGGGGCAGCG.CCCCCCGAAGCGCCTGAGGcgcg                                                                               >  1:2283246/1‑101 (MQ=255)
                       ctgcgcgCCGCCGCCGGCCCGCTTTTCT.TTGGGAGAACCCGGGGGGGGTG.C.CGCG.GCACCACAAGGGGGGGGCA.CCCCCCCCCGAAGCGCCTGAGGCgcgg                                                                              <  2:2605222/98‑1 (MQ=255)
                       cTTCGCGCCGCAGCAGGCCCGCTTTTCC.TTGGGCGAAACGCCGGGGGGTT.G.CGCG.TCAGCGCAAGCGGGGGGCAGCG.CCCCCCGAAGCGCCTGAGGCgcgg                                                                              >  1:1368060/1‑101 (MQ=255)
                          cgcgCCGCAGCAGGCCCGCTTTTCC.TTGGGCGAAACGCCGGGGGGTT.G.CGCG.TCAGCGCAAGCGGGGGGCAGCG.CCCCCCGAAGCGCCTGAGGCGcggcgg                                                                           >  2:2039006/1‑101 (MQ=255)
                           gcgcCGCAGCAGGCCCGCTTTTCC.TTGGGCGAAACGCCGGGGGGTT.G.CGCG.TCAGCGCAAGCGGGGGGCAGCG.CCCCCCGAAGCGCCTGAGGCGCGGCggg                                                                          >  2:2281332/1‑101 (MQ=255)
                            cgccgcAGCAGGCCCGCTTTTCC.TTGGGCGAAACGCCGGGGGGTT.G.CGCG.TCAGCGCAAGCGGGGGGCAGCG.CCCCCCGAAGCGCCTGAGGCGCGGCGGGc                                                                         >  1:2309492/1‑101 (MQ=255)
                            cgccgcAGCAGGCCCGCTTTTCC.TTGGGCGAAACGCCGGGGGGTT.G.CGCG.TCAGCGCAAGCGGGGGGCAGCG.CCCCCCGAACCGCCTGAGGCGCGGCGGGc                                                                         >  2:1301411/1‑101 (MQ=255)
                             gccgcAGCAGGCCCGCTTTTCC.TTGGGCGAAACGCCGGGGGGTT.G.CGCG.TCAGCGCAAGCGGGGGGCAGCG.CCCCCCGAAGCGCCTGAGGCGCGGCGGGCt                                                                        >  1:3527176/1‑101 (MQ=255)
                              ccgcAGCGGGCCCGTTTTTCC.TTGGGCGAAACGCGGGGGGGTG.G.CGCG.TCAGCGCAGGCGGGGGGCAGC.CCCCCCCGAAGCGCCTGAGGCGCGGCGGGCtt                                                                       <  2:2728147/101‑1 (MQ=255)
                              ccgcAGCAGGCCCGCTTTTCC.TTGGGCGAAACGCCGGGGGGTT.G.CGCG.TCAGCGCAAGCGGGGGGCAGCG.CCCCCCGAAGCGCCTGAGGCGCGGCGGGCtt                                                                       >  2:1752995/1‑101 (MQ=255)
                                 cagcagGCCCGCTTTTCC.TTGGGCGAAACGCCGGGGGGTT.G.CGCG.TCAGCGCAAGCGGGGGGCAGCG.CCCCCCGAAGCGCCTGAGGCGCGGCGGGCTTgcc                                                                    >  2:3281730/1‑101 (MQ=255)
                                  agcagGCCCGCTTTTCC.TTGGGCGAAACGCCGGGGGGTT.G.CGCG.TCAGCGCAAGCGGGGGGCAGCG.CCCCCCGAAGCGGCTGAGGCGCGGCGGGCGTgccc                                                                   >  2:2875900/1‑101 (MQ=255)
                                  agcagGCCCGCTTTTCC.TTGGGCGAAACGCCGGGGGGTT.G.CGCG.TCAGCGCAAGCGGGGGGCAGCG.CCCCCCGAAGCGCCTGAGGCGCGGCGGGCTTtccc                                                                   >  1:2078987/1‑101 (MQ=255)
                                  agcagGCCCGCTTTTCC.TTGGGCGAAACGCCGGGGGGTT.G.CGCG.TCAGCGCAAGCGGGGGGCAGCG.CCCCCCGAAGCGCCTGAGGCGCGGCGGGCCTgccc                                                                   >  2:619532/1‑101 (MQ=255)
                                   gcagGCCCGCTTTTCC.TTGGGCGAAACGCCGGGGGGTT.G.CGCG.TCAGCGCAAGCGGGGGGCAGCG.CCCCCCGAAGCGCCTGAGGCGCGGCGGGCTTgcccc                                                                  >  1:2870217/1‑100 (MQ=255)
                                    cagGCCCGCTTTTCC.TTGGGCGAAACGCCGGGGGGTT.G.CGCG.TCAGCGCAAGCGGGGGGCAGCG.CCCCCCGAAGCGCCTGAGGCGCGGCGGGCTTgcccgc                                                                 >  1:277030/1‑101 (MQ=255)
                                     agGCCCGCTTTTCC.TTGGGCGAAACGCCGGGGGGTT.G.CGCG.TCAGCGCAAGCGGGGGGCAGCG.CCCCCCGAAGCGCCTGAGGCGCGGCGGGCTTgcccgcc                                                                >  1:3419082/1‑101 (MQ=255)
                                       gCCCGCTTTTCC.TTGGGCGAAACGCCGGGGGGTT.G.CGCG.TCAGCGCAAGCGGGGGGCAGCG.CCCCCCGAAGCGCCTGAGGCGCGGCGGGCTTGCCCGCCCa                                                              >  2:2187736/1‑101 (MQ=255)
                                       cccggtttttCC.TGGGGCGAACCG.CGGGGGGTT.GCCGCG.TCAGC.CAAGGGGGGGGCCCGCCCCCCCCGAAGCGCCTGAGGCGCGGCGGGCTTGCCCGCCCa                                                              <  1:2125215/95‑1 (MQ=16)
                                         ccGTTTTTCCGTTGGCCGAACGGCGGGGGGGTG.G.CG.G.TGAGGCGAAGCGGGGGGACACCCCCCCCGGAAGCGCCTGAGGCGCGGCGGGCTTGCCCGCCCaa                                                             <  2:4240316/101‑1 (MQ=9)
                                          cccTTTTCC.TTGGGGCAAACCGGGGGGGGTTGG.CGGG.GCGGCG.AAGGGGGGGGGGGCCCCCCCCCGAAGCGCCTGAGGCGCGGCGGGCTTGCCCGCCCAAc                                                            <  1:4207462/99‑1 (MQ=12)
                                           tttttcCC.TGGGGCGAACCGCGGGGGGGTG.G.CCCG.CCCGCGCAGGGGGGGGGCCGC.CCCCCCCGAAGCGCCTGAGGCGCGGCGGGCTTGCCCGCCCAACGa                                                          <  1:661131/99‑1 (MQ=18)
                                           ctttttCC.CTGGGGAAAACGGCGGGGGGGT.G..GCGCTCCGCCCCAGGGGGGGGGCGGCCCCCCCCGAAGCGCCTGAGGCGCGGCGGGCTTGCCCGCCCAACg                                                           <  2:3526928/99‑1 (MQ=12)
                                               ttCC.TGGGGCGAAACGCGGGGGGGTG.G.CGCG.CCAGCGCAGGGGGGGGGCAGC.CCCCCCCGAAGCGCCTGAGGCGCGGCGGGCTTGCCCGCCCAACGAAGgc                                                      <  2:497770/101‑1 (MQ=25)
                                                tCC.TTGGGCGAAACGCCGGGGGGTT.G.CGCG.TCAGCGCAAGCGGGGGGCAGCG.CCCCCCGAAGCGCCTGAGGCGCGGCGGGCTTGCCCGCCCAACGAAGgcg                                                     >  2:3720160/1‑101 (MQ=255)
                                                 cc.TTGGGCGAAACGCCGGGGGGTT.G.CGCG.TCAGCGCAAGCGGGGGGCAGCG.CCCCCCGAAGCGCCTGAGGCGCGGCGGGCTTGGCCGCCCAACGAAGgcgc                                                    >  1:3983496/1‑101 (MQ=255)
                                                 cc.TTGGGCGAAACGCCGGGGGGTT.G.CGCG.TCAGCGCAAGCGGGGGGCAGCG.CCCCCCGAAGCGCCTGAGGCGCGGCGGGCTTGCCCGCCCAACGAAGgcgc                                                    >  1:3647158/1‑101 (MQ=255)
                                                 cc.TTGGGCGAAACGCCGGGGGGTT.G.CGCG.TCAGCGCAAGCGGGGGGCAGCG.CCCCCCGAAGCGCCTGAGGCGCGGCGGGCTTGCCCGCCCAACGAAGgcgc                                                    >  1:295467/1‑101 (MQ=255)
                                                  c.TTGGGCAGACCGCGGGGGTGTG.G.CCCG.CCGGGGCAGGTCGGGTGCGGCG.CCCCTCGAAGCGCCTGAGGCGCGGCGGGCTTGCCCGCCCAACGAAGGCGCa                                                   <  2:3670982/101‑1 (MQ=9)
                                                   ctggggCGAAACGCGGGGGGGTT.G.CGCG.TCAGCGCAAGCGGGGGGCAGCG.CCCCCCGAAGCGCCTGAGGCGCGGCGGGCTTGCCCGCCCAACGAAGGCGCa                                                   <  2:2660677/98‑1 (MQ=25)
                                                    gggggcagaccGCCGGGGGGGT.G.CGCG.CCCGGGCAGGGGGGGGGCGGCG.CCCCCCAAAGCGCCTGAGGCGCGGCGGGCTTGCCCGCCCAACGCAGGCGCaa                                                  <  2:1472950/91‑1 (MQ=16)
                                                      gtgCGAAACGCCGGGGGGTT.G.CGCG.TCAGCGCAAGCGGGGGGCAGCG.CCCCCCGAAGCGCCTGAGGCGCGGCGGGCTTGCCCGCCCAACGAAGGCGCAacc                                                >  1:2079446/3‑101 (MQ=255)
                                                      gggCGACCCGCGGGGGGGTG.G.GGCG.CCAGCG.AAGGGGGGGGCCGGCCCCCCCCGAAGCGCCTGAGGCGCGGCGGGCTTGCCCGCCCAACGAAGGCGCAacc                                                <  2:3494889/101‑1 (MQ=18)
                                                      gggCGAAACGCCGGGGGGTT.G.CGCG.TCAGCGCAAGCGGGGGGCAGCG.CCCCCCGAAGCGCCTGAGGCGCGGCGGGCTTGCCCGCCCAACGAAGGCGCAacc                                                >  2:1095020/1‑101 (MQ=255)
                                                        gCGAAACGCCGGGGGGTT.G.CGCG.TCAGCGCAAGCGGGGGGCAGCG.CCCCCCGAAGCGCCTGAGGCGCGGCGGGCTTGCCCGCCCAACGAAGGCGCAaccac                                              >  1:1011932/1‑101 (MQ=255)
                                                         cGAAACGCCGGGGGGTT.G.CGCG.TCAGCGCAAGCGGGGGGCAGCG.CCCCCCGAAGCGCCTGAGGCGCGGCGGGTATACCCACCCAAAGAAgggggcacaacc                                             >  2:3870442/1‑91 (MQ=255)
                                                          gAAACGCCGGGGGGTT.G.CGCG.TCAGCGCAAGCGGGGGGCAGCG.CCCCCCGAAGCGCCTGAGGCGCGGCGGGCTTGCCCGCCCAACGAAGGCGCAACCACCt                                            >  1:2092625/1‑101 (MQ=255)
                                                          gAAACGCCGGGGGGTT.G.CGCG.TCAGCGCAAGCGGGGGGCAGCG.CCCCCCGAAGCGCCTGAGGCGCGGCGGGCTTGCCCGCCCAACCAAGGCGCAacccccc                                            >  1:3161884/1‑100 (MQ=25)
                                                           aaaCGCCGGGGGGTT.G.CGCG.TCAGCGCAAGCGGGGGGCAGCG.CCCCCCGAAGCGCCTGAGGCGCGGCGGGCTTGCCCGCCCAACGAAAGCGCaaaacacct                                           >  2:2326134/1‑94 (MQ=25)
                                                            aaCGCCGGGGGGTT.G.CGCG.TCAGCGCAAGCGGGGGGCAGCG.CCCCCCAAAGCGACTGAGGCGCGGCGGGCTTGCCCCCCCCACGAAGGCGCAACCACCCTc                                          >  2:3037360/1‑101 (MQ=25)
                                                             ccGCGGGGGGGTT.G.CGCG.CCCGCGAAAGGGGGGGGCAGC.CCCCCCCGAAGCGCCTGAGGCGCGGCGGGCTTGCCCGCCCAACGAAGGCGCAACCACCTTCa                                         <  2:2639171/100‑1 (MQ=25)
                                                             aCGCGGGGGGGTT.G.CGCG.CCAGCGAAGGGGGGGGGCAGCG.CCCCCCGAAGCGCCTGAGGCGCGGCGGGCTTGCCCGCCCAACGAAGGCGCAACCACCTTCa                                         <  2:2507758/101‑1 (MQ=21)
                                                               gCCGGGGGGTT.G.CGCG.TCAGCGCAAGCGGGGGGCAGCG.CCCCCCGAAGCGCCTGAGGCGCGGCGGGCTTGCCCGCCCAACGAAGGCGCAACCACCTTCAGt                                       >  2:1049258/1‑101 (MQ=255)
                                                                 gggggggGG.G.GGCG.TCGGGGCAGGGGGGGGGCAGC.CCCCCCCGAAGCGCCTGAGGCGCGGCGGGCTTGCCCGCCCAACGAAGGCGCAACCACCTTCAGTCg                                     <  1:1240456/100‑1 (MQ=21)
                                                                 cGGGGGGTT.G.CGCG.TCAGCGCAAGCGGGGGGCAGCG.CCCCCCGAAGCGCCTGAGGCGCGGCGGGCTTGCCCGCCCAACGAAGGCGCAACCACCTTCAGTCg                                     >  2:3201957/1‑101 (MQ=255)
                                                                 cGGGGGGTT.G.CGCG.TCAGCGCAAGCGGGGGGCAGCG.CCCCCCGAAGCGCCTGAGGCGCGGCGGGCTTGCCCGCCCAACGAAGGCGCAACCACCTTCAGTCg                                     >  2:2223762/1‑101 (MQ=255)
                                                                  ggggggTT.G.CGCG.TCAGCGCAAGCGGGGGGCAGCG.CCCCCCGAAGCGCCTGAGGCGCGGCGGGCTTGCCCGCCCAACGAAGGCGCAACCACCTTCCGTCgg                                    >  1:690263/1‑101 (MQ=255)
                                                                  ggggggTT.G.CGCG.TCAGCGCAAGCGGGGGGCAGCG.CCCCCCGAAGCGCCTGAGGCGCGGCGGGCTTGCCCGCCCAACGAAGGCGCAACCACCTTCAGTCgg                                    >  2:585119/1‑101 (MQ=255)
                                                                  ggggggTT.G.CGCG.TCAGCGCAAGCGGGGGGCAGCG.CCCCCCGAAGCGCCTGAGGCGCGGCGGGCTTGCCCGCCCAACGAAGGCGCAACCACCCTCAGTCgg                                    >  2:3427144/1‑101 (MQ=255)
                                                                   gggggTT.G.CGCG.TCAGCGCAAGCGGGGGGCAGCG.CCCCCCGAAGCGCCTGAGGCGCGGCTGGCTTGCCCGCCCAACGAAGGCGCAACCACCTTCAGTCGGc                                   >  1:3503245/1‑101 (MQ=25)
                                                                    gggtTT.G.CGCG.CGGGCGCGGGGGGGGGCCCGC.CCCCCCCGAAGCGCCTGAGGCGCGGCGGGCTTGCCCGCCCAACGAAGGCGCAACCACCTTCAGTCGGCt                                  <  1:3408533/101‑1 (MQ=21)
                                                                     gggTC.G.CGCG.CCAGCGAAGGGGGGGGGCAGC.CCCCCCCGAAGCGCCTGAGGCGCGGCGGGCTTGCCCGCCCAACGAAGGCGCAACCACCTTCAGTCGGCtt                                 <  2:2805839/101‑1 (MQ=25)
                                                                       gTT.G.CGCG.TCAGCGCAAGCGGGGGGCAGCG.CCCCCCGAAGCGCCTGAGGCGCGGCGGGCTTGCCCGCCCAACGAAGGCGCAACCACCTTCAGTCGGCTTCg                               >  2:3594057/1‑101 (MQ=255)
                                                                        tt.G.CGCG.TCAGCGCAAGCGGGGGGCAGCG.CCCCCCGAAGCGCCTGAGGCGCGGCGGGCTTGCCCGCCCAACGAAGGCGCAACCACCTTCAGTCGGCTTCGa                              >  1:1873376/1‑101 (MQ=255)
                                                                         t.G.CGCG.TCAGCGCAAGCGGGGGGCAGCG.CCCCCCGAAGCGCCTGAGGCGCGGCGGGCTTGCCCGCCCAACGAAGGCGCAACCACCTTCAGTCGGCTTCGAt                             >  2:2888625/1‑101 (MQ=255)
                                                                           ggggcccgcAGCGCGGGGGGGGCGCCCCCCCCCCCGGAAGCGCCTGAGGCGCGGCGGGCTTGCCCGCCCAACGAAGGCGCAACCACCTTCAGTCGGCTTCg                               <  1:2594336/93‑1 (MQ=21)
                                                                              gcg.TCAGCGCAAGCGGGGGGCAGCG.CCCCCCGAAGCGCCTGAGGCGCGGCGGGCCTGGCCGCCCAACGAACGGGCAACCACCTTCAGTCGGCTTCGATccc                          >  2:3990581/1‑101 (MQ=37)
                                                                               cgcgccggggaaGGGGGGGGGCGCCCCCCCCCGAAGCGCCTGAGGCGCGGCGGGCTTGCCCGCCCAACGAAGGCGCAACCACCTTCAGTCGGCTTCGATcc                           <  1:3252239/91‑1 (MQ=25)
                                                                               cg.TCACGGCAGGCGGGGGGCAGCG.CCCCCCGAAGCGCCTGAGGCGCGGCGGGCTTGCCCGCCCAACGAAGGCGCAACCACCTTCAGTCGGCTTCGATcccc                         <  2:4395464/101‑1 (MQ=21)
                                                                                  tCAGCGCAAGCGGGGGGCAGCG.CCCCCCGAAGCGCCTGAGGCGCGGCGGGCTTGCCCGCCCAACGAAGGCGCAACCACCTTCAGTCGGCTTCGATCCccac                       >  2:1232615/1‑101 (MQ=255)
                                                                                   cAGCGCAAGGGGGGGGCA.CCCCCCCCCGAAGCGCCTGAGGCGCGGCGGGCTTGCCCGCCCAACGAAGGCGCAACCACCTTCAGTCGGCTTCGATCCccacc                      <  2:2939551/101‑1 (MQ=255)
                                                                                   cAGCGCAAGCGGGGGGCAGCG.CCCCCCGAAGCGCCTGAGGCGCGGCGGGCTTGCCCGCCCAACGAAGGCGCAACCACCTTCAGTCGGCTTCGATCCccacc                      >  2:2701657/1‑101 (MQ=255)
                                                                                    aGCGCAAGCGGGGGGCAGCG.CCCCCCGAAGCGCCTGAGGCGCGGCGGGCTTGCCCGCCCAACGAAGGCGCAACCACCTTCAGTCGGCTTCGATCCccacca                     >  2:2509464/1‑101 (MQ=255)
                                                                                       gcAAGCGGGGGGCAGCG.CCCCCCGAAGCGCCTGAGGCGCGGCGGGCTTGCCCGCCCAACGAAGGCGCAACCACCTTCAGTCGGCTTCGATCCCCACCAgcg                  >  2:413654/1‑101 (MQ=255)
                                                                                        gaaGGGGGGGGCAGC.CCCCCCCGAAGCGCCTGAGGCGCGGCGGGCTTGCCCGCCCAACGAAGGCGCAACCACCTTCAGTCGGCTTCGATCCCCACCAgcgc                 <  1:2151426/100‑1 (MQ=255)
                                                                                        cAAGCGGGGGGCAGCG.CCCCCCGAAGCGCCTGAGGCGCGGCGGGCTTGCCCGCCCAACGAAGGCGCAACCACCTTCAGTCGGCTTCGATCCCCACCAgcgc                 >  1:274309/1‑101 (MQ=255)
                                                                                         aaGCGGGGGGCAGCG.CCCCCCGAAGCGCCTGAGGCGCGGCGGGCTTGCCCGCCCAACGAAGGCGCAACCACCTTCCGTCGGCTTCGATCCCCACCAGcgcg                >  1:1953908/1‑101 (MQ=255)
                                                                                             ggggggCCGCCGCCACCGGAAGCGCCTGAGGCGCGGCGGGCTTGCCCGCCCCACGAACGCGCAACCCCCTTCAGTCGGCTTCGATCCCCACCAGCgcggcg             <  2:4313737/101‑1 (MQ=37)
                                                                                             ggggggCAGCG.CCCCCCGAAGCGCCTGAGGCGCGGCGGGCTTGCCCGCCCAACGAAGGCGCAACCACCTTCAGTCGGCTTCGATCCCCACCCGCgcggcgg            >  1:1041647/1‑101 (MQ=255)
                                                                                             ggggggCAGCG.CCCCCCGAAGCGCCTGAGGCGCGGCGGGCTTGCCCGCCCAACGAAGGCGCAACCACCTTCAGTCGGCTTCGATCCCCACCAGCgcggcgg            >  2:1718034/1‑101 (MQ=255)
                                                                                              gggggCAGCG.CCCCCCGAAGCGCCTGAGGCGCGGCGGGCTTGCCCGCCCAACGAAGGCGCAACCACCTTCAGTCGGCTTCGATCCCCACCAGCgcggcggc           >  1:3931368/1‑101 (MQ=255)
                                                                                              gggggCAGCG.CCCCCCGAAGCGCCTGAGGCGCGGCGGGCTTGCCCGCCCAACGAAGGCGCAACCACCTTCAGTCGGCTTCGATCCCCACCAGCgcggcggc           >  2:2604259/1‑101 (MQ=255)
                                                                                              gggggCAGC.CCCCCCCGAAGCGCCTGAGGCGCGGCGGGCTTGCCCGCCCAACGAAGGCGCAACCACCTTCAGTCGGCTTCGATCCCCACCAGCgcggcggc           <  1:2018998/101‑1 (MQ=255)
                                                                                                gggCAGCG.CCCCCCGAAGCGCCTGAGGCGCGGCGGGCTTGCCCGCCCAACGAAGGCGCAACCACCTTCAGTCGGCTTCGATCCCCACCAGCGCGGCGGCGa         >  1:905990/1‑101 (MQ=255)
                                                                                                   cccccccccccgGAAGCGCCTGAGGCGCGGCGGGCTTGCCCGCCCAACGAAGGCGCAACCACCTTCAGTCGGCTTCGATCCCCACCAGCGCGGCGGCGAAt       <  2:65085/96‑1 (MQ=255)
                                                                                                        cccccccGAAGCGCCTGAGGCGCGGCGGGCTTGCCCGCCCAACGAAGGCGCAACCACCTTCAGTCGGCTTCGATCCCCACCAGCGCGGCGGCGAATTCCTc  <  1:3565492/101‑1 (MQ=255)
                                                                                                        cccccccGAAGCGCCTGAGGCGCGGCGGGCTTGCCCGCCCAACGAAGGCGCAACCACCTTCAGTCGGCTTCGATCCCCACCAGCGCGGCGGCGAATTCCTc  <  2:658665/101‑1 (MQ=255)
                                                                                                        |                                                                                                    
CGGTTTCCTGACAGGATTGTGTCCTTCGCGCCGCAGCAGGCCCGCTTTTCC.TTGGGCGAAACGCCGGGGGGTT.G.CGCG.TCAGCGCAAGGGGGGGGCAGCGCCCCCCCGAAGCGCCTGAGGCGCGGCGGGCTTGCCCGCCCAACGAAGGCGCAACCACCTTCAGTCGGCTTCGATCCCCACCAGCGCGGCGGCGAATTCCTC  >  NC_014034/2559835‑2560035

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