Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_014034	900,793	+G	4.0%	intergenic (‑122/‑165)	ldh ← / → ftsW	L‑lactate dehydrogenase/cell division protein FtsW

Read alignment evidence...
	seq id	position		change	freq	score	reads	annotation	genes	product
*	NC_014034	900,793	3	.→G	4.0%	2.9 ‑3.9 ‑0.9	75	intergenic (‑122/‑165)	ldh/ftsW	L‑lactate dehydrogenase/cell division protein FtsW
Reads supporting (aligned to +/- strand):  new base (0/3):  ref base (54/18):  total (54/21)
Fisher's exact test strand distribution p-value = 1.97e-02
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 1.16e-01

GCTACCATGGCACGCTCACGCGTCCGCGCAACCCCGTTTCCG.C.GCGG.CGCCGCTTTGCCGCAGCCCCGGTTTTCCTTGGTGCAAATAC.GCCC...GGGGGT...GC...CGCGGA..G.....CCCGCGGGTCGCCCCGGCCCTGCCGGGGCGAAAGCCCCCGCGCGCCGCGGCCCGGCAAAAGCTCTGGTCAGAG  >  NC_014034/900696‑900875                                                                                                          |                                                                                               gcTACCATGGCACGCTCACGCGTCCGCGCAACCCCGTTTCCGCC.GCGGCCGCCGCTTTGCCGCAGCCCCGGTTTTCCTTGGTGCAAATAC.GCCC...GGGGGT...g                                                                                             <  2:3494840/101‑1 (MQ=255) gcTACCATGGCACGCTCACGCGTCCGCGCAACCCCGTTTCCGCC.GCGGCCGCCGCTTTGCCGCAGCCCCGGTTTTCCTTGGTGCAAATAC.GCCC...GGGGGT...g                                                                                             <  2:3904566/101‑1 (MQ=255) gcTACCATGGCACGCTCACGCGTCCGCGCAACCCCGTTTCCGCC.GCGGCCGCCGCTTTGCCGCAGCCCCGGTTTTCCTTGGTGCAAATAC.GCCC...GGGGGT...g                                                                                             <  1:3967322/101‑1 (MQ=255) gcTACCATGGCACGCTCACGCGTCCGCGCAACCCCGTTTCCGCC.GCGGCCGCCGCTTTGCCGCAGCCCCGGTTTTCCTTGGTGCAAATAC.GCCC...GGGGGT...g                                                                                             >  1:3597102/1‑101 (MQ=255)  cTACCATGGCACGCTCACGCGTCCGCGCAACCCCGTTTCCGCC.GCGGCCGCCGCTTTGCCGCAGCCCCGGTTTTCCTTGGTGCAAATAC.GCCC...GGGGGT...Gc                                                                                            <  2:615114/101‑1 (MQ=255)   tACCATGGCGCGCGCACGCGTCCGCGCAACCCCGTTTCCGCC.GCGGCCGCCGCTTTGCCGCAGCCCCGGTTTTCCTTGGTGCAAATAC.GCCC...GGGGGT...Gcg                                                                                           <  1:1112891/101‑2 (MQ=255)   tACCATGGCACGCTCACGCGTCCGCGCAACCCCGTTTCCGCC.GCGGCCGCCGCTTTGCCGCAGCCCCGGTTTTCCTTGGTGCAAATAC.GCCC...GGGGGG...Gcg                                                                                           >  1:2509012/1‑100 (MQ=38)      cATGGCACGCTCACGCGTCCGCGCAACCCCGTTTCCGCC.GCGGCCGCCGCTTTGCCGCAGCCCCGGTTTTCCTTGGTGCAAATAC.GCCC...GGGGGT...Gcgggg                                                                                        >  1:785274/1‑97 (MQ=38)      cATGGCACGCTCACGCGTCCGCGCAACCCCGTTTCCGCC.GCGGCCGCCGCTTTGCCGCAGCCCCGGTTTTCCTTGGTGCAAATAC.GCCC...GGGGGT...Gcgggg                                                                                        <  1:4393492/101‑5 (MQ=38)       aTGGCCCGCTCACGCGTCCGCGCAGCCCCGTTTCCGCC.GCGGCCGCCGCTTTGCCGCAGCCCCGGTTTTCCTTGGTGCAAATAC.GCCC...GGGGGT...Gcggggg                                                                                       <  2:1950040/101‑6 (MQ=37)       aTGGCACGCTCACGCGTCCGCGCAACCCCGTTTCCGCC.GCGGCCGCCGCTTTGCCGCAGCCCCGGTTTTCCTTGGTGCAAATAC.GCCC...GGGGGT...Gcggggg                                                                                       <  1:5582310/101‑6 (MQ=38)          gCACGCTCACGCGCCCGCGCAAACCCGTTCCCGAC.GCGGCCGCCGCTTTGCCGCAGCCCCGGTTTTCCTTGGTGCAAATAC.GCCC...GGGGGT...gagggggcgg                                                                                    <  2:416132/101‑10 (MQ=255)           cACGCTCCCCCGACCCCGCAGCACAGATTCCA.CAGCGGCCGCCGCTTTGCCCCAGCCCCGGTTTTCCTTCGTGCAAATAC.GCCC...GGGGGT...Gcgggggcgga                                                                                   <  2:2516957/101‑10 (MQ=255)           cACGCTCACGCGTCCGCGCAACCCCGTTTCCGCC.GCGGCCGCCGCTTTGCCGCAGCCCCGGTTTTCCTTGGTGCAAATAC.GCCC...GGGGGT...Gcgggggcgga                                                                                   <  1:2581360/101‑10 (MQ=38)           cACGCTCACGCGTCCGCGCAACCCCGTTTCCGCC.GCGGCCGCCGCTTTGCCGCAGCCCCGGTTTTCCTTGGTGCAAATAC.GCCC...GGGGGT...Gcgggggcgga                                                                                   >  1:1080507/1‑92 (MQ=38)            aCGCTCACGCGTCCGCGCAACCCCGTTTCCGCC.GCGGCCGCCGCTTTGCCGCAGCCCCGGTTTTCCTTGGTGCAAATAC.GCCC...GGGGGT...Gcgggggcggag                                                                                  >  1:4201983/1‑91 (MQ=37)                             cAACCCCGTTTCCGCC.GCGGCCGCCGCTTTGCCGCAGCCCCGGTTTTCCTTGGTGCAAATAC.GCCC...GGGGGT...GCGGGGGCGGA..GCC...CCCGCGGGTCGcccc                                                            >  1:94253/1‑101 (MQ=16)                              aaCCCCGTTTCCGCC.GCGGCCGCCGCTTTGCCGCAGCCCCGGTTTTCCTTGGTGCAAATAC.GCCC...GGGGGT...GCGGGGGCGGA..GCC...CCCGCGGGTCGCCCCg                                                           >  2:3229056/1‑101 (MQ=16)                               aCCCCGTTTCCGCC.GCGGCCGCCGCTTTGCCGCAGCCCCGGTTTTCCTTGGTGCAAATAC.GCCC...GGGGGT...GCGGGGGCGGA..GCC...CCCGCGGGTCGCCCCgg                                                          >  1:5214695/1‑101 (MQ=16)                               aCCCCGTTTCCGCC.GCGGCCGCCGCTTTGCCGCAGCCCCGGTTTTCCTTGGTGCAAATAC.GCCC...GGGGGT...GCGGGGGCGGA..GCC...CCCGCGGGTCGCCCCgg                                                          >  2:55171/1‑101 (MQ=16)                                ccccGTTTCCGCC.GCGGCCGCCGCTTTGCCGCAGCCCCGGTTTTCCTTGGTGCAAATAC.GCCC...GGGGGT...GCGGGGGCGGA..GCC...CCCGCGGGTCGCCCCGGc                                                         >  1:4023335/1‑101 (MQ=16)                                ccccGTTCCCGCC.GCGGCCGCCGCTTTGCCGCAGCCCCGGTTTTCCTTGGTGCAAATAC.GCCC...GGGGGT...GCGGGGGCGGA..GCC...CCCGCGGGTCGCCCCGGc                                                         <  2:5390899/101‑1 (MQ=12)                                  cctttcccgcc.gcgGCCGCCTCTTGGCCGCAGCCCCGGTTTTCCTTGGTGCAAATCC.GCCG...GGGGGTGGGGG...GGCGGA..GCC...CCCGCGGGTCGCCCCGGccc                                                       <  2:3286248/91‑1 (MQ=2)                                  ccGTTTCCGCC.GCGGCCGCCGCTTTGCCGCAGCCCCGGTTTTCCTTGGTGCAAATAC.GCCC...GGGGGT...GCGGGGGCGGA..GCC...CCCGCGGGTCGCCCCGGccc                                                       >  1:388452/1‑101 (MQ=16)                                   cgttccggcc.gcgGCCGCCGCTTTGCCGCAGCCCCGGTTTTCCTTGGTGCAAATAC.GCCC...GGGGGT...GCGGGGGCGGA..GCC...CCCGCGGGTCGCCCCGGCCCt                                                      <  2:5792561/92‑1 (MQ=12)                                    gTTTCCGCC.GCGGCCGCCGCTTTGCCGCAGCCCCGGTTTTCCTTGGTGCAAATAC.GCCC...GGGGGT...GCGGGGGCGGA..GCC...CCCGCGGGTCGCCCCGGCCCTg                                                     >  1:2986400/1‑101 (MQ=16)                                    gTTTCCGCC.GCGGCCGCCGCTTTGCCGCAGCCCCGGTTTTCCTTGGTGCAAATAC.GCCC...GGGGGT...GCGGGGGCGGA..GCC...CCCGCGGGTCGCCCCGGCCCTg                                                     >  2:1182790/1‑101 (MQ=16)                                      ttCCGCC.GCGGCCGCCGCTTTGCCGCAGCCCCGGTTTTCCTTGGTGCAAATAC.GCCC...GGGGGT...GCGGGGGCGGA..GCC...CCCGCGGGTCGCCCCGGCCCTGcc                                                   >  2:1641585/1‑101 (MQ=16)                                        ccgcc.gcgGCCGCCGCTTTGCCGCAGCCCCGGTTTTCCTTGGTGCAAATAC.GCCC...GGGGGT...GCGGGGGCGGA..GCC...CCCGCGGGTCGCCCCGGCCCTGCCgg                                                 >  2:2911644/5‑101 (MQ=18)                                         cgcc.gcgGCCGCCGCTTTGCCGCAGCCCCGGTTTTCCTTGGTGCAAATAC.GCCC...GGGGGT...GCGGGGGCGGA..GCC...CCCGCGGGTCGCCCCGGCCCTGCCggg                                                >  2:1494438/4‑101 (MQ=18)                                         cgcc.gcgGCCGCCGCTTTGCCGCAGCCCCGGTTTTCCTTGGTGCAAATAC.GCCC...GGGGGT...GCGGGGGCGGA..GCC...CCCGCGGGTCGCCCCGGCCCTGCCggg                                                >  2:495586/4‑101 (MQ=18)                                          gcc.gcgGCCGCCGCTTTGCCGCAGCCCCGGTTTTCCTTGGTGCAAATCC.GCCC...GGGGGT...GCGGGGGCGGA..GCC...CCCGCGGGTCGCCCCGGCCCTGCCgggg                                               <  1:2629547/99‑1 (MQ=16)                                           cc.gcgGCCGCCGCTTTGCCGCAGCCCCGGTTTTCCTTGGTGCAAATAC.GCCC...GGGGGT...GCGGGGGCGGA..GCC...CCCGCGGGTCGCCCCGGCCCTGCCGGGGc                                              >  2:4168776/2‑101 (MQ=18)                                           cc.gcgGCCGCCGCTTTGCCGCAGCCCCGGTTTTCCTTGGTGCAAATAC.GCCC...GGGGGT...GCGGGGGCGGA..GCC...CCCGCGGGTCGCCCCGGCCCTGCCGGGGc                                              >  1:3826685/2‑101 (MQ=18)                                                ggccgcgccTTTCCCGCAGCCCCGGTTTTCCTTGGTGCAAATCC.GCCG...GGGGGTGGGGG...GGCGGA..GCC...CCCGCGGGTCGCCCCGGCCCTGCCGGGGCGaaa                                          <  1:2713222/93‑1 (MQ=9)                                                ggccgcctcTTCGCCGCAGCCCCGGTTTCCCTTGGTG.AAATGCGGCCG...GGGGGT...GCGGGGGCGGA..GCC...CCCGCGGGTCGCCCCGGCCCTGCCGGGGCGaaa                                          <  1:3212459/98‑1 (MQ=2)                                                ggccgccgcTTTGCCGCAGCCCCGGTTTTCCTTGGTGCAAATCC.GCCC...GGGGGT...GCGGGGGCGGA..GCC...CCCGCGGGTCGCCCCGGCCCTGCCGGGGCGaaa                                          <  2:3596885/98‑1 (MQ=18)                                                ggccgccgcTTTGCCGCAGCCCCGGTTTTCCTTGGTGCAAATAC.GCCC...GGGGGT...GCGGGGGCGGA..GCC...CCCGCGGGTCGCCCCGGCCCTGCCGGGGCGaaa                                          <  1:803891/98‑1 (MQ=18)                                                 gccgccgcTTGGCCGCAGCCCCNGTTNNNCNTGGTGCAAATAC.GCCCGGGGGGGGG...GG...GGCGGA..GCC...CCCGCGGGTCGCCCCGGCCCTGCCGGGGCGAAAg                                         <  1:3229256/99‑1 (MQ=12)                                                    gcctcTTGGCCGCAGCCCGGGTTTTCCTTGGGCCAAATCC.GCCG...GGGGGT...GC...GGGGGC..GGCCCCCCCGCGGGTCGCCCCGGCCCTGCCGGGGCGAAAGccc                                      <  2:4380333/101‑1 (MQ=2)                                                    gccgcTTTGCCGCAGCCCCGGTTTTCCTTGGTGCAAATAC.GCCC...GGGGGT...GCGGGGGCGGA..GCC...CCCGCGGGTCGCCCCGGCCCTGCCGGGGCGAAAGgcc                                      >  2:3836517/1‑101 (MQ=18)                                                     ccgcTTTGCCGCAGCCCCGGTTTTCCTTGGTGCAAATAC.GCCC...GGGGGT...GCGGGGGCGGA..GCC...CCCGCGGGTCGCCCCGGCCCTGCCGGGGCGAAAGcccc                                     >  2:3505171/1‑101 (MQ=21)                                                      cgcTGTGGCGCAGCCCCGGTTTTCCTGGGTGCAAATGC.GCCC...GGGGGTGGGGG...GGCGGG..GCC...CCCGCGGGTCGCCCCGGCCCTGCCGGGGCGAAAGccccc                                    <  2:3718672/101‑1 (MQ=12)                                                       gcTTTGCCGCAGCCCCGGTTTTCCTTGGTGCAAATAC.GCCC...GGGGGT...GCGGGGGCGGA..GCC...CCCGCGGGTCGCCCCGGCCCTGCCGGGGCGAAAGCCCCcg                                   >  1:2636607/1‑101 (MQ=18)                                                       gcTTTGCCGCAGCCCCGGTTTTCCTTGGTGCAAATAC.GCCC...GGGGGT...GCGGGGGCGGA..GCC...CCCGCGGGTCGCCCCCGCCCTGCCGGGGCGAAAGCCCCcg                                   >  2:3433540/1‑101 (MQ=18)                                                        cTTTGCCGCAGCCCCGGTTTTCCTTGGTGCAAATAC.GCCC...GGGGGT...GCGGGGGCGGA..GCC...CCCGCGGGTCGCCCCGGCCCTGCCGGGGCGAAAGCCCCcgc                                  >  1:1838899/1‑101 (MQ=18)                                                        cTTTGCCGCAGCCCCGGTTTTCCTTGGTGCAAATAC.GCCC...GGGGGT...GCGGGGGCGGA..GCC...CCCGCGGGTCGCCCCGGCCCTGCCGGGGCGAAAGCCCCcgc                                  >  1:739470/1‑101 (MQ=18)                                                        cTTTGCCGCAGCCCCGGTTTTCCTTGGTGCAAATAC.GCCC...GGGGGT...GCGGGGGCGGA..GCC...CCCGCGGGTCGCCCCGGCCCTGCCGGGGCGAAAGCCCCcgc                                  >  2:2923777/1‑101 (MQ=18)                                                         tttGCCGCAGCCCCGGTTTTCCTTGGTGCAAATAC.GCCC...GGGGGT...GCGGGGGCGGA..GCC...CCCGCGGGTCGCCCCGGCCCTGCCGGGGCGAAAGCCCCcgcg                                 >  2:3445841/1‑101 (MQ=18)                                                            gCCGCAGCCCCGGTTTTCCTTGGTGCAAATAC.GCCC...GGGGGT...GCGGGGGCGGA..GCC...CCCGCGGGTCGCCCCGGCCCTGCCGGGGCGAAAGCCCCCgcgcgc                              >  1:2393452/1‑101 (MQ=18)                                                            gCCGCAGCCCCGGTTTTCCTTGGTGCAAATAC.GCCC...GGGGGT...GCGGGGGCGGA..GCC...CCCGCGGGTCGCCCCGGCCCTGCCGGGGCGAAAGCCCCCgcgcgc                              >  1:2901180/1‑101 (MQ=18)                                                             ccGCAGCCCCGGTTTTCCTTGGTGCAAATAC.GCCC...GGGGGT...GCGGGGGCGGA..GCC...CCCGCGGGTCGCCCCGGCCCTGCCGGGGCGAAAGCCCCCGCGcgcc                             >  1:3498262/1‑101 (MQ=18)                                                              cGCAGCCCCGGTTTTCCTTGGTGCAAATAC.GCCC...GGGGGT...GCGGGGGCGGA..GCC...CCCGCGGGTCGCCCCGGCCCTGCCGGGGCGAAAGCCCCCGCGcgccg                            <  2:5553473/101‑1 (MQ=18)                                                              cGCAGCCCCGGTTTTCCTTGGTGCAAATAC.GCCC...GGGGGT...GCGGGGGCGGA..GCC...CCCGCGGGGCGGCCCGGCCCTGCCGGGGTGAAAGCCCCCGCGcgccg                            >  1:1150346/1‑101 (MQ=18)                                                               gCAGCCCCGGTTTTCCTTGGTGCAAATAC.GCCC...GGGGGT...GCGGGGGCGGA..GCC...CCCGCGGGTCGCCCCGGCCCTGCCGGGGCGAAAGCCCCCGCGcgccgc                           >  1:2321305/1‑101 (MQ=18)                                                               gCAGCCCCGGTTTTCCTTGGTGCAAATAC.GCCC...GGGGGC...GCGGGGGCGGA..GCC...CCCGCGGGTCGCCCCGGCCCTGCCGGGGCGAAAGCCCCCGCGcgccgc                           >  1:5641172/1‑101 (MQ=18)                                                                cAGCCCCGGTTTTCCTTGGTGCAAATAC.GCCC...GGGGGT...GCGGGGGCGGA..GCC...CCCGCGGGTCGCCCCGGCCCTGCCGGGGCGAAAGCCCCCGCGCGCcgcg                          >  2:748226/1‑101 (MQ=18)                                                                cAGCCCCGGTTTTCCTTGGTGCAAATAC.GCCC...GGGGGT...GCGGGGGCGGA..GCC...CCCGCGGGTCGCCCCGGCCCTGCCGGGGCGAAAGCCCCCGCGCGCcgcg                          >  2:3726079/1‑101 (MQ=18)                                                                cAGCCCCGGTTTTCCTTGGTGCAAATAC.GCCC...GGGGGT...GCGGGGGCGGA..GCC...CCCGCGGGTCGCCCCGGCCCTGCCGGGGCGAAAGCCCCCGCGCGCcgcg                          >  2:4114933/1‑101 (MQ=18)                                                                cAGCCCCGGTTTTCCTTGGTGCAAATAC.GCCC...GGGGGT...GCGGGGGCGGA..GCC...CCCGCGGGTCGCCCCGGCCCTGCCGGGGCGAAAGCCCCCGCGCGCcgcg                          <  2:2303378/101‑1 (MQ=18)                                                                cAGCCCCGGTTTTCCTTGGTGCAAATAC.GCCC...GGGGGT...GCGGGGGCGGA..GCC...CCCGCGGGTCGCCCCGGCCCTGCCGGGGCGAAAGCCCCCGCGCGCcgcg                          >  1:2347478/1‑101 (MQ=18)                                                                    cccGGTTTTCCTTGGTGCAAATAC.GCCC...GGGGGT...GCGGGGGCGGA..GCC...CCCGCGGGTCGCCCCGGCCCTGCCGGGGCGAAAGCCCCCGCGCGCCGCGGccc                      <  2:2523380/101‑1 (MQ=18)                                                                    cccGGTTTTCCTTGGTGCAAATAC.GCCC...GGGGGT...GCGGGGGCGGA..GCC...CCCGCGGGCCGCCCCGGCCCTGCCGGGCCGAAAGCCCCCGCGCGCCGCGGccc                      >  1:707089/1‑101 (MQ=18)                                                                     ccGGTTTTCCTGGGTGCAAATCC.GCCC...GGGGGTGGGGG...GGCGGA..GCC...CCCGCGGGTCGCCCCGGCCCTGCCGGGGCGAAAGCCCCCGCGCGCCGCGGCCCg                     <  1:4779651/101‑1 (MQ=17)                                                                        gTTTTCCTTGGTGCAAATAC.GCCC...GGGGGT...GCGGGGGCGGA..GCC...CCCGCGGGTCGCCCCGGCCCTGCCGGGGCGAAAGCCCCCGCGCGCCGCGGCCCGGCa                  >  1:880238/1‑101 (MQ=18)                                                                        gTTTTCCTTGGTGCAAATAC.GCCC...GGGGGT...GCGGGGGCGGA..GCC...CCCGCGGGTCGCCCCGGCCCTGCCGGGGCGAAAGCCCCCGCGCGCCGCGGCCCGGCa                  >  1:4665722/1‑101 (MQ=18)                                                                         ttttCCTTGGTGCAAATAC.GCCC...GGGGGT...GCGGGGGCGGA..GCC...CCCGCGGGTCGCCCCGGCCCTGCCGGGGCGAAAGCCCCCGCGCGCCGCGGCCCGGCaa                 >  1:1434348/1‑101 (MQ=18)                                                                           ttCCTTGGTGCAAATAC.GCCC...GGGGGT...GCGGGGGCGGA..GCC...CCCGCGGGTCGCCCCGGCCCTGCCGGGGCGAAAGCCCCCGCGCGCCGCGGCCCGGCaaaa               >  2:5112403/1‑101 (MQ=18)                                                                            tCCTTGGTGCAAATAC.GCCC...GGGGGT...GCGGGGGCGGA..GCC...CCCGCGGGTCGCCCCGGCCCTGCCGGGGCGAAAGCCCCCGCGCGCCGCGGCCCGGCAAAAg              >  1:2320623/1‑101 (MQ=18)                                                                             ccTTGGTGCAAATCC.GCCC...GGGGGT...GCGGGGGCGGACCC.....CCCGCGGGTCGCCCCGGCCCTGCCGGGGCGAAAGCCCCCGCGCGCCGCGGCCCGGCAAAAGc             <  2:739402/101‑1 (MQ=18)                                                                             ccTTGGTGCAAATAC.GCCC...GGGGGT...GCGGGGGCGGA..GCC...CCCGCGGGTCGCCCCGGCCCTGCCGGGGCGAAAGCCCCCGCGCGCCGCGGCCCGGCAAAAGc             >  1:2422126/1‑101 (MQ=18)                                                                              cTTGGTGCAAATAC.GCCC...GGGGGT...GCGGGGGCGGA..GCC...CCCGCGGGTCGCCCCGGCCCTGCCGGGGCGAAAGCCCCCGCGCGCCGCGGCCCGGCAAAAGct            >  1:2739236/1‑101 (MQ=18)                                                                                tGGTGCAAATAC.GCCC...GGGGGT...GCGGGGGCGGA..GCC...CCCGCGGGCCGCCCCGGCCCTGCCGGGGCCAAAGCCCCCGCGCGCCGCGGCCCGGCAAAAGctct          >  2:2679156/1‑101 (MQ=18)                                                                                 ggTGCAAATAC.GCCC...GGGGGT...GCGGGGGCGGA..GCC...CCCGCGGGTCGCCCCGGCCCTGCCGGGGCGAAAGCCCCCGCGCGCCGCGGCCCGGCAAAAGCTCTg         >  1:3872166/1‑101 (MQ=18)                                                                                  gTGCAAATAC.GCCC...GGGGGT...GCGGGGGCGGA..GCC...CCCGCGGGTCGCCCCGGCCCTGCCGGGGCGAAAGCCCCCGCGCGCCGCGGCCCGGCAAAAGCTCTgg        >  1:2173718/1‑101 (MQ=18)                                                                                     cAAATAC.GCCC...GGGGGT...GCGGGGGCGGA..GCC...CCCGCGGGTCGCCCCGGCCCTGCCGGGGCGAAAGTCCCCGCGCGCCGCGGCCCGGCAAAAGCTCTGGTCa     >  2:3086343/1‑101 (MQ=18)                                                                                     cAAATAC.GCCC...GGGGGT...GCGGGGGCGGA..GCC...CCCGCGGGTCGCCCCGGCCCTGCCGGGGCGAAAGCCCCCGCGCGCCGCGGCCCGGCAAAAGCTCTGGTCa     >  1:4120059/1‑101 (MQ=18)                                                                                     cAAATAC.GCCC...GGGGGT...GCGGGGGCGGA..GCC...CCCGCGGGTCGCCCCGGCCCTGCCGGGGCGAAAGCCCCCGCGCGCCGCGGCCCGGCAAAAGCTCTGGTCa     >  1:226993/1‑101 (MQ=18)                                                                                      aaaTAC.GCCC...GGGGGT...GCGGGGGCGGA..GCC...CCCGCGGGTCGCCCCGGCCCTGCCGGGGCGAAAGCCCCCGCGCGCCGCGGCCCGGCAAAAGCTCTGGTCag    >  2:1528748/1‑101 (MQ=18)                                                                                      aaaTAC.GCCC...GGGGGT...GCGGGGGCGGA..GCC...CCCGCGGGACGCCCCGGCCCTGCCGGGGCGAAAGCCCCCGCGCGCCGCGGCCCGGCAAAAGCTCTGGTCag    >  1:770264/1‑101 (MQ=18)                                                                                        aTAC.GCCC...GGGGGT...GCGGGGGCGGA..GCC...CCCGCGGGTCGCCCCGGCCCTGCCGGGGCGAAAGCCCCCGCGCGCCGCGGCCCGGCAAAAGCTCTGGTCagag  >  1:2423900/1‑101 (MQ=18)                                                                                        aTAC.GCCC...GGGGGT...GCGGGGGCGGA..GCC...CCCGCGGGTCGCCCCGGCCCTGCCGGGGCGAAAGCCCCCGCGCGCCGCGGCCCGGCAAAAGCTCTGGTCagag  >  2:1333455/1‑101 (MQ=18)                                                                                        aTAC.GCCC...GGGGGT...GCGGGGGCGGA..GCC...CCCGCGGGTCGCCCCGGCCCTGCCGGGGCGAAAGCCCCCGCGCGCCGCGGCCCGGCAAAAGCTCTGGTCagag  >  1:1314039/1‑101 (MQ=18)                                                                                                          |                                                                                               GCTACCATGGCACGCTCACGCGTCCGCGCAACCCCGTTTCCG.C.GCGG.CGCCGCTTTGCCGCAGCCCCGGTTTTCCTTGGTGCAAATAC.GCCC...GGGGGT...GC...CGCGGA..G.....CCCGCGGGTCGCCCCGGCCCTGCCGGGGCGAAAGCCCCCGCGCGCCGCGGCCCGGCAAAAGCTCTGGTCAGAG  >  NC_014034/900696‑900875 Base quality scores:  ATCG  < 3 ≤  ATCG  < 15 ≤  ATCG  < 22 ≤  ATCG  < 32 ≤  ATCG  < 40 ≤  ATCG