Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_014034	1,301,094	Δ1 bp	92.7%	intergenic (+51/‑58)	urtA → / → urtC	urea ABC transporter permease UrtA/urea ABC transporter permease UrtC

Read alignment evidence...
	seq id	position		change	freq	score	reads	annotation	genes	product
*	NC_014034	1,301,094	0	C→.	92.7%	8.6 ‑0.5 ‑0.4	68	intergenic (+51/‑58)	urtA/urtC	urea ABC transporter permease UrtA/urea ABC transporter permease UrtC
Reads supporting (aligned to +/- strand):  new base (44/19):  ref base (2/2):  total (46/22)
Fisher's exact test strand distribution p-value = 5.84e-01
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 3.67e-01

CCGCCCGCGCGGCATCATCGCGCTGCGCGGCCGCGCCGCCGGAGATTGAGCCGATG.CCCCCCGGCG.CCTTTCATCTTGCCCAAAATATCCCGC.GGGGGTCCCGGGGG.CGCGAAG.CCCCCGGCGCCGCCAGATCCCGGAGACCCCGATGACACAGCCGATGACCCAAGGGTTCCTGTCCAAAAACCCCTCGGTGCTG  >  NC_014034/1300995‑1301190                                                                                                       |                                                                                                   ccgcccgcGCGGCATCATCGCGCTGCGCGGCCGCGCCGCCGGAGATTGAGCCGATG.CCCCCCGCCGCCCTTTCATCTTGCCCAAAATATCCCGC.GGGGGTc                                                                                                    <  2:2295561/101‑1 (MQ=255) ccgcccgcGCGGCATCATCGCGCTGCGCGGCCGCGCCGCCGGAGATTGAGCCGATG.CCCCCCGCCGCCCTTTCATCTTGCCCAAAATATCCCGC.GGGGGTc                                                                                                    >  2:5420/1‑101 (MQ=255) ccgcccgcGCGGCATCATCGCGCTGCGCGGCCGCGCCGCCGGAGATTGAGCCGATG.CCCCCCGCCGCCCTTTCATCTTGCCCAAAATATCCCGC.GGGGGTc                                                                                                    <  1:4139198/101‑1 (MQ=255) ccgcccgcGCGGCATCATCGCGCTGCGCGGCCGCGCCGCCGGAGATTGAGCCGATG.CCCCCCGCCGCCCTTTCATCTTGCCCAAAATATCCCGC.GGGGGTc                                                                                                    <  2:4085626/101‑1 (MQ=255)  cgcccgcGCGGCATCATCGCGCTGCGCGGCCGCGCCGCCGGAGATTGAGCCGATG.CCCCCCGCCGCCCTTTCATCTTGCCCAAAATATCCCGC.GGGGGTcc                                                                                                   <  2:11419/101‑1 (MQ=255)   gcccgcGCGGCATCATCGCGCTGCGCGGCCGCGCCGCCGGAGATTGAGCCGATG.CCCCCCGCCGCCCTTTCATCTTGCCCAAAATATCCCGC.GGGGGTccg                                                                                                  <  1:1558052/101‑2 (MQ=255)    cccgcGCGGCATCATCGCGCTGCGCGGCCGCGCCGCCGGAGATTGAGCCGATG.CCCCCCGCCGCCCTTTCATCTTGCCCAAAATATCCCGC.GGGGGTccgg                                                                                                 >  1:56914/1‑99 (MQ=255)    cccgcGCGGCATCATCGCGCTGCGCGGCCGCGCCGCCGGAGATTGAGCCGATG.CCCCCCGCCGCCCTTTCATCTTGCCCAAAATATCCCGC.GGGGGTccgg                                                                                                 >  2:3911699/1‑99 (MQ=255)     ccgcGCGGCATCATCGCGCTGCGCGGCCGCGCCGCCGGAGATTGATCCGATG.CCCCCCGCCGCCCTTTCATCTTGCCCAAAATATCCCGC.GGGGGTCCGgg                                                                                                <  1:3927758/101‑1 (MQ=37)      cgcgcgGCATCATCGCGCTGCGCGGCCGCGCCGCCGGAGATTGAGCCGATG.CCCCCCGCCGCCCTTTCATCTTGCCCAAAATATCCCGC.GGGGGTCCGggg                                                                                               >  1:2493041/1‑101 (MQ=37)        cgcgGCATCATCGCGCTGCGCGGCCGCGCCGCCGGAGATTGAGCCGATG.CCCCCCGCCGCCCTTTCATCTTGCCCAAAATATCCCGC.GGGGGTCCGggggc                                                                                             >  1:3311159/1‑100 (MQ=37)        cgcgGCATCATCGCGCTGCGCGGCCGCGCCGCCGGAGATTGAGCCGATG.CCCCCCGCCGCCCTTTCATCTTGCCCAAAATATCCCGC.GGGGGTCCGggggc                                                                                             >  1:2843440/1‑100 (MQ=37)         gcgGCATCATCGCGCTGCGCGGCCGCGCCGCCGGAGATTGAGCCGATG.CCCCCCGCCGCCCTTTCATCTTGCCCAAAATATCCCGC.GGGGGT.CCGGGGg.cg                                                                                          >  2:43392/1‑101 (MQ=25)          cgGCATCATCGCGCTGCGCGGCCGCGCCGCCGGAGATTGAGCCGATG.CCCCCCGCCGCCCTTTCATCTTGCCCAAAATATCCCGC.GGGGGT.CCGGGGg.cgc                                                                                         <  1:2810712/101‑1 (MQ=25)           ggCATCATCGCGCTGCGCGGCCGCGCCGCCGGAGATTGAGCCGATG.CCCCCCGCCGCCCTTTCATCTTGCCCAAAATATCCCGC.GGGGGT.CCGGGGG.cgcg                                                                                        <  2:3936243/101‑1 (MQ=25)            gCATCATCGCGCTGCGCGGCCGCGCCGCCGGAGATTGAGCCGATG.CCCCCCGCCGCCCTTTCATCTTGCCCAAAATATCCCGC.GGGGGT.CCGGGGG.CGCGa                                                                                       >  1:1042902/1‑101 (MQ=25)             catcatCGCGCTGCGCGGCCGCGCCGCCGGAGATTGAGCCGATG.CCCCCCGCCGCCCTTTCATCTTGCCCAAAATATCCCGC.GGGGGT.CCGGGGG.CGCGaa                                                                                      <  2:3521178/101‑1 (MQ=25)             catcatCGCGCTGCGCGGCCGCGCCGCCGGAGATTGAGCCGATG.CCCCCCGCCGCCCTTTCATCTTGCCCAAAATATCCCGC.GGGGGT.CCGGGGG.CGCGaa                                                                                      <  1:760017/101‑1 (MQ=25)              atcatcGCGCTGCGCGGCCGCGCCGCCGGAGATTGAGCCGATG.CCCCCCGCCGCCCTTTCATCTTGCCCAAAATATCCCGC.GGGGGT.CCGGGGG.CGCGAAg                                                                                     <  2:1830844/101‑1 (MQ=25)              atcatcGCGCTGCGCGGCCGCGCCGCCGGAGATTGAGCCGATG.CCCCCCGCCGCCCTTTCATCTTGCCCAAAATATCCCGC.GGGGGT.CCGGGGG.CGCGAAg                                                                                     >  2:3512955/1‑101 (MQ=25)                catcGCGCTGCGCGGCCGCGCCGCCGGAGATTGAGCCGATG.CCCCCCGCCGCCCTTTCATCTTGCCCAAAATATCCCGC.GGGGGT.CCGGGGG.CGCGAAG.cc                                                                                  >  2:3073784/1‑101 (MQ=25)                catcGCGCTGCGCGGCCGCGCCGCCGGAGATTGAGCCGATG.CCCCCCGCCGCCCTTTCATCTTGCCCAAAATATCCCGC.GGGGGT.CCGGGGG.CGCGAAG.cc                                                                                  <  2:881746/101‑1 (MQ=25)                 atcGCGCTGCGCGGCCGCGCCGCCGGAGATTGAGCCGATG.CCCCCCGCCGCCCTTTCATCTTGCCCAAAATATCCCGC.GGGGGT.CCGGGGG.CGCGAAG.ccc                                                                                 <  1:3928336/101‑1 (MQ=21)                 atcGCGCTGCGCGGCCGCGCCGCCGGAGATTGAGCCGATG.CCCCCCGCCGCCCTTTCATCTTGCCCAAAATATCCCGC.GGGGGT.CCGGGGG.CGCGAAG.ccc                                                                                 <  1:1643343/101‑1 (MQ=21)                  tcGCGCTGCGCGGCCGCGCCGCCGGAGATTGAGCCGATG.CCCCCCGCCGCCCTTTCATCTTGCCCAAAATATCCCGC.GGGGGT.CCGGGGG.CGCGAAG.cccc                                                                                <  1:3314130/101‑1 (MQ=21)                  tcGCGCTGCGCGGCCGCGCCGCCGGAGATTGAGCCGATG.CCCCCCGCCGCCCTTTCATCTTGCCCAAAATATCCCGC.GGGGGT.CCGGGGG.CGCGAAG.cccc                                                                                <  2:1477606/101‑1 (MQ=21)                  tcGCGCTGCGCGGCCGCGCCGCCGGAGATTGAGCCGATG.CCCCCCGCCGCCCTTTCATCTTGCCCAAAATATCCCGC.GGGGGT.CCGGGGG.CGCGAAG.cccc                                                                                <  1:2013903/101‑1 (MQ=21)                  tcGCGCTGCGCGGCCGCGCCGCCGGAGATTGAGCCGATG.CCCCCCGCCGCCCTTTCATCTTGCCCAAAATATCCCGC.GGGGGT.CCGGGGG.CGCGAAG.cccc                                                                                <  2:3189798/101‑1 (MQ=21)                  tcGCGCTGCGCGGCCGCGCCGCCGGAGATTGAGCCGATG.CCCCCCGCCGCCCTTTCATCTTGCCCAAAATATCCCGC.GGGGGT.CCGGGGG.CGCGAAG.cccc                                                                                <  1:1007470/101‑1 (MQ=21)                   cgcgCTGCGCGGCCGCGCCGCCGGAGATTGAGCCGATT.CCCCCCGCCGCCCTTTCATCTTGCCCAAAATATCCCGC.GGGGGT.CCGGGGG.CGCGAAG.ccccc                                                                               <  1:3105191/101‑1 (MQ=21)                   cgcgCTGCGCGGCCGCGCCGCCGGAGATTGAGCCGATG.CCCCCCGCCGCCCTTTCATCTTGCCCAAAATATCCCGC.GGGGGT.CCGGGGG.CGCGAAG.ccccc                                                                               >  2:2310576/1‑101 (MQ=21)                    gcgcTGCGCGGCCGCGCCGCCGGAGATTGAGCCGATG.CCCCCCGCCGCCCTTTCATCTTGCCCAAAATATCCCGC.GGGGGT.CCGGGGG.CGCGAAG.CCCCCg                                                                              >  2:1467941/1‑101 (MQ=21)                    gcgcTGCGCGGCCCCGCGGCCGGAGATGGAGCCGAGG.CCCCCCGCCGCCCTTTCATCTTGCCCAAAATCTCCCGC..GGGGTTCCGGGGG.CGCGAAG.CCCCCg                                                                              <  2:541088/101‑1 (MQ=17)                      gcTGCGCGGCCGCGCCGCGGGAGATTGAGCCGACG.CCCCCCGCCGCCCTTTCATCTTGCCCAAAATATCCCGG.GGGGGT.CCGGGGG.CGCGAAG.CCCCCGgc                                                                            <  1:1465913/101‑1 (MQ=18)                           gcgGCCGCGCCGCCGGAGATTGAGCCGATG.CCCCCCGCCGCCCTTTCATCTTGCCCAAAATATCCCGC.GGGGGT.CCGGGGG.CGCGAAG.CCCCCGGcgccgc                                                                       >  2:2659376/1‑101 (MQ=18)                           gcgGCCGCGCCGCCGGAGATTGAGCCGATG.CCCCCCGCCGCCCTTTCATCTTGCCCAAAATATCCCGC.GGGGGT.CCGGGGG.CGCGAAG.CCCCCGGcgccac                                                                       >  1:1676430/1‑99 (MQ=18)                            cgGCCGCGCCGCCGGAGATTGAGCCGATGCCCCCCCCCCC.CCTTTCATCTTGCCCAAAATATCCCCC.GGGGG.CCCGGGGG.CGCGAGG.CCCCGGGCgccccc                                                                      >  1:2290954/1‑101 (MQ=18)                             ggCCGCGCCGCCGGAGATTGAGCCGATG.CCCCCCGCCGCCCTTTCATCTTGCCCAAAATATCCCGC.GGGGGT.CCGGAGG.CGCGACG.CCGCCGGCGCCGCCa                                                                     >  2:971359/1‑101 (MQ=18)                                    ccgccgGAGATTGAGCCGATG.CCCCCCGCCGCCCTTTCATCTTGCCCAAAATATCCCGC.GGGGGT.CCGGGGG.CGCGAAG.CCCCCGGCGCCGCCAGATCCCg                                                              >  2:4371213/1‑101 (MQ=18)                                         ggAGATTGAGCCGATG.CCCCCCGCCGCCCTTTCATCTTGCCCAAAATATCCCGC.GGGGGT.CCGGGGG.CGCGAAG.CCCCCGGCGCCGCCAGATCCCGGAGAc                                                         >  1:51340/1‑101 (MQ=18)                                             atgaacccGATC.CCCCCCGCCGCCCTTTCATCTTGCCCAAAATACCCCGC.GGGGGT.CGGGGGG.CGCGAAG.CCCCCGGCGCCGCCAGATCCCGGAGACCCCg                                                     <  2:448595/95‑1 (MQ=17)                                               tGAGCCGTCC.CCCCCCCCCC.CCTTTTCTCCTGCCCAAAAGCCCCCGG.GGGGGGCCGGGGGG.CGCGAGC.CCCCCGGCGCCGCCAGATCCCGGAGACCCCGAt                                                   <  2:1059855/101‑1 (MQ=255)                                               tGAGCCGATG.CCCCCCGCCGCCCTTTCATCTTGCCCAAAATATCCCGC.GGGGGT.CCGGGGG.CGCGAAG.CCCCCGGCGCCGCCAGATCCCGGAGACCCCGAt                                                   >  1:3717390/1‑101 (MQ=18)                                                     gATG.CCCCCCGCCGCCCTTTCATCTTGCCCAAAATATCCCGC.GGGGGT.CCGGGGG.CGCGAAG.CCCCCGGCGCCGCCAGATCCCGGAGACCCCGATGAcaca                                             >  1:1016026/1‑101 (MQ=18)                                                     gATG.CCCCCCGCCGCCCTTTCATCTTGCCCAAAATATCCCGC.GGGGGT.CCGGGGG.CGCGAAG.CCCCCGGCGCCGCCAGATCCCGGAGACCCCGATGAcaca                                             >  1:1605708/1‑101 (MQ=18)                                                     gATA.CCCCCCGCCGCCCTTTCATCTTGCCCAAAATATCCCGC.GGGGGT.CCGGGGG.CGCGAAG.CCCCCGGCGCCGCCAGATCCCGGAGACCCCGATGAcaca                                             >  2:1917592/1‑101 (MQ=18)                                                      aTG.CCCCCCGCCGCCCTTTCATCTTGCCCAAAATATCCCGC.GGGGGT.CCGGGGG.CGCGAAG.CCCCCGGCGCCGCCAGATCCCGGAGACCCCGATGACACAg                                            >  1:796546/1‑101 (MQ=18)                                                       tG.CCCCCCGCCGCCCTTTCATCTTGCCCAAAATATCCCGC.GGGGGT.CCGGGGG.CGCGAAG.CCCCCGGCGCCGCCAGATCCCGGAGACCCCGCTGACACAGc                                           >  1:1661465/1‑101 (MQ=18)                                                        g.CCCCCCGCCGCCCTTTCATCTTGCCCAAAATACCCCGG.GGGGGT.CGGGGGG.CGCGAAC.CCCCCGGCGCCGCCAGATCCCGGAGACCCCGATGACACAGcc                                          <  1:2891683/101‑1 (MQ=18)                                                        ccccccccCCCC.CCTTTCATCGCG.CCAAAAAGCCCCGG.GGGGGGGGGGGGGG.CGCCCCC.CCCCCGGCGCCGCCAGATCCCGGAGACCCCGATGACACAGcc                                          <  2:2281223/99‑1 (MQ=255)                                                           cccccGCCGCCCTTTCATCTTGCCCAAAATATCCCGC.GGGGGT.CCGGGGG.CGCGAAG.CCCCCGGCGCCGCCAGATCCCGGAGACCCCGATGACACAGCCGa                                        >  1:1971852/1‑101 (MQ=18)                                                             ccccccgcccTTCCATCTTGCCCAAAATCCCCCGC.GGGTG.CCCGGGGG.CGCGAAC.CCCCCGGCGCCGCCAGATCCCGGAGACCCCGATGACACAGCCGATg                                      <  1:1468046/93‑1 (MQ=255)                                                             ccccccccccTTTCCTC.TGCCCAAAA.ATCCCCCGGGGGGGCCCGGGGG.CGCGACC.CCCCCGGCGCCGCCAGATCCCGGAGACCCCGATGACACAGCCGATg                                      <  2:3310927/93‑1 (MQ=255)                                                                 ccgcccTTTCATCTTGCCAAAAATATCCCGG.GGGGGT.CGGGGGG.CGCGAAG.CCCCCGGCGCCGCCAGATCCCGGAGACCCCGATGACACAGCCGATGAccc                                  <  2:431478/97‑1 (MQ=21)                                                                  ccccttttCTTCTTCCCAAAAATA.CCCGG.GGGGGGCCGGGGGC.CGCGACC.CCCCCGGCGCCGCCAGATCCCGGAGACCCCGATGACACAGCCGATGACCCa                                 <  1:1931271/96‑1 (MQ=255)                                                                   gcccTTTCATCTTGCCCAAAATATCCCGC.GGGGGT.CCGGGGG.CGCGAAG.CCCCCGGCGCCGCCAGATCCCGGAGACCCCGATGACACAGCCGATGACCCaa                                >  1:2629247/3‑101 (MQ=255)                                                                   gcccTTTCATCTTGCCCAAAATATCCCGC.GGGGGT.CCGGGGG.CGCGAAG.CCCCCGGCGCCGCCAGATCCCGGAGACCCCGATGACACAGCCGATGACCCaa                                >  2:3278811/3‑101 (MQ=255)                                                                    cccTTTCATCTCGCCCAAAATCCCCCGG.GGGGGT.CGGGGGG.CGCGAAG.CCCCCGGCGCCGCCAGATCCCGGAGACCCCGATGACACAGCCGATGACCCAAg                               <  1:2931797/100‑1 (MQ=25)                                                                     ccTTTCATCTTGCCCAAAATATCCCGC.GGGGGT.CCGGGGG.CGCGAAG.CCCCCGGCGCCGCCAGATCCCGGAGACCCCGATGACACAGCCGCTGACCCAAgg                              >  2:3843704/1‑101 (MQ=255)                                                                     ccTTTCATCTTGCCCAAAATATCCCGC.GGGGGT.CCGGGGG.CGCGAAG.CCCCCGGCGCCGCCAGATCCCGGAGACCCCGATGACACAGCCGATGACCCAAgg                              >  1:1798554/1‑101 (MQ=255)                                                                     ccTTTCATCTTGCCCAAAATATCCCGC.GGGGGT.CCGGGGG.CGCGAAG.CCCCCGGCGCCGCCAGATCCCGGAGACCCCGATGACACAGCCGATGACCCAAgg                              <  2:4434993/101‑1 (MQ=255)                                                                       tttCATCTTGCCCAAAATATCCCGC.GGGGGT.CCGGGGG.CGCGAAG.CCCCCGGCGCCGCCAGATCCCGGAGACCCCGATGACACAGCCGATGACCCAAGGGt                            >  1:3345353/1‑101 (MQ=255)                                                                        ttCATCTTGCCCAAAATATCCCGC.GGGGGT.CCGGGGG.CGCGAAG.CCCCCGGCGCCGCCAGATCCCGGAGACCCCGATGACACAGCCGATGACCCAAGGGtt                           >  1:1099406/1‑101 (MQ=255)                                                                        ttCATCTCACCAAAAATA.CCCGG.GGGGGGCCGGGGGG.CGGAAAC.CCCCGGGCGCCGCCAGATCCCGGAGACCCCGATGACACAGCCGATGACCCAAGGGtt                           <  1:2659548/101‑1 (MQ=255)                                                                         tCATCTTGCCCAAAATATCCCGC.GGGGGT.CCGGGGG.CGCGAAG.CCCCCGGCGCCGCCCGATCCCGGAGAGCCCGATGACACAGCCGATGACCCAAAGGTTc                          >  1:824755/1‑101 (MQ=21)                                                                          tcatctccccAAAAATCCCCGG.GGGGGGCCGGGGGG.CGCGAAC.CCCCCGGCGCCGCCAGATCCCGGAGACCCCGATGACACAGCCGATGACCCAAGGGTTc                          <  1:1150275/94‑1 (MQ=255)                                                                          taTCTTTCCCAAAATATCCCGG.GGGGGGTCGGGGGG.C.CGTAGCCCCCCGGCGCCGCCAGATCCCGGAGACCCCGATGACACAGCCGATGACCCAAGGGTTc                          <  2:2613662/100‑1 (MQ=255)                                                                           ctcctggccAAAATACCCCGG.GGGGGGCCGGGGGC.CCCGAAC.CCCCGGGCGCCGCCAGATCCCGGAGACCCCGATGACACAGCCGATGACCCAAGGGTTcc                         <  2:4424363/94‑1 (MQ=255)                                                                           aTCTTGCCCAAAATATCCCGC.GGGGGT.CCGGGGG.CGCGAAG.CCCCCGGCGCCGCCAGATCCCGGAGACCCCGATGACACAGCCGATGACCCAAGGGTTCCt                        >  1:1304787/1‑101 (MQ=255)                                                                           aTCTTGCCCAAAATATCCCGC.GGGGGT.CCGGGGG.CGCGAAG.CCCCCGGCGCCGCCAGATCCCGGAGACCCCGATGACACAGCCGATGACCCAAGGGTTCCt                        >  2:1139423/1‑101 (MQ=255)                                                                           aTCTTGCCAAAAATATCCCGC.GGGGGT.CCGGGGG.CGCGAAG.CCCCCGGCGCCGCCAGATCCCGGAGACCCCGATGACACAGCCGATGACCCAAGGGTTCCt                        <  2:3391165/101‑1 (MQ=255)                                                                            tcctccccAAAAACCCCCGG.GGGGGGCCGGGGGC.CGCCACC.CCCCGGGCGCCGCCAGATCCCGGAGACCCCGATGACACAGCCGATGACCCAAGGGTTCCt                        <  1:4092007/96‑1 (MQ=255)                                                                             cTTGCCCAAAATATCCCGC.GGGGGT.CCGGGGG.CGCGAAG.CCCCCGGCGCCGCCAGATCCCGGAGACCCCGATGACACAGCCGATGACCCAAGGGTTCCTGt                      >  2:817659/1‑101 (MQ=255)                                                                              ttGCCCAAAATATCCCGC.GGGGGT.CCGGGGG.CGCGAAG.CCCCCGGCGCCGCCAGATCCCGGAGACCCCGATGACACAGCCGATGACCCAAGGGTTCCTGTc                     >  1:441357/1‑101 (MQ=255)                                                                              cctcccCAAAAAACCCGC.GGGGGGGCCGGGGGCCGCGAAC.CCCCGGGCGCCGCCAGATCCCGGAGACCCCGATGACACAGCCGATGACCCAAGGGTTCCTg                       <  1:3977723/98‑1 (MQ=255)                                                                               tGCCCAAAATATCCCGC.GGGGGT.CCGGGGG.CGCGAAG.CCCCCGGCGCCGCCAGATCCCGGAGACCCCGATGACACAGCCGATGACCCAAGGGTTCCTGTcc                    >  1:1184866/1‑101 (MQ=255)                                                                               ctcccAAAACTCCCCGC.GGGGG.GCCGGGGC.CGCGAAG.CCCCCGGCGCCGCCAGATCCCGGAGACCCCGATGACACAGCCGATGACCCAAGGGTTCCTGTcc                    <  1:3580280/99‑1 (MQ=37)                                                                                gCCCAAAATATCCCGC.GGGGGT.CCGGGGG.CGCGAAG.CCCCCGGCGCCGCCAGATCCCGGAGACCCCGATGACACAGCCGATGACCCAAGGGTTCCTGTCCa                   >  2:1898299/1‑101 (MQ=255)                                                                                 cccAAAATTCCCCGG.GGGGGT.CCGGGGG.CGCGAAG.CCCCCGGCGCCGCCAGATCCCGGAGACCCCGATGACACAGCCGATGACCCAAGGGTTCCTGTCCaa                  <  2:4385362/101‑1 (MQ=38)                                                                                  ccAAAATATCCCGC.GGGGGT.CCGGGGG.CGCGAAG.CCCCCGGCGCCGCCAGATCCCGGAGACCCCGATGACACAGCCGATGACCCAAGGGTTCCTGTCCaaa                 >  1:1567740/1‑101 (MQ=255)                                                                                  ccAAAATATCCCGC.GGGGGT.CCGGGGG.CGCGAAG.CCCCCGGCGCCGCCAGATCCCGGAGACCCCGATGACACAGCCGATGACCCAAGGGTTCCTGTCCaaa                 >  1:1308960/1‑101 (MQ=255)                                                                                    aaaaTATCCCGC.GGGGGT.CCGGGGG.CGCGAAG.CCCCCGGCGCCGCCAGATCCCGGAGACCCCGATGACACAGCCGATGACCCAAGGGTTCCTGTCCaaaaa               >  2:3502987/1‑101 (MQ=255)                                                                                      aaTATCCCGC.GGGGGT.CCGGGGG.CGCGAAG.CCCCCGGCGCCGCCAGATCCCGGAGACCCCGATGACACAGCCGATGACCCAAGGGTTCCTGTCCAAAAAcc             >  2:1618880/1‑101 (MQ=255)                                                                                       atatCCCGC.GGGGGT.CCGGGGG.CGCGAAG.CCCCCGGCGCCGCCAGATCCCGGAGACCCCGATGACACAGCCGATGACCCAAGGGTTCCTGTCCAAAAAccc            >  1:328934/1‑101 (MQ=255)                                                                                       aaagcgcccggggggGCCGGGGGG.CGCAACC.CCCCGGGCGCCGCCAGATCCCGGAGACCCCGATGACACAGCCGATGACCCAAGGGTTCCTGTCCAAAAAc              <  2:1178596/91‑1 (MQ=255)                                                                                        tatCCCGC.GGGGGT.CCGGGGG.CGCGAAG.CCCCCGGCGCCGCCAGATCCCGGAGACCCCGATGACACAGCCGATGACCCAAGGGTTCCTGTCCAAAAAcccc           >  1:1191724/1‑101 (MQ=255)                                                                                         atcccccggggggGGCCGGGGG.CGCGACC.CCCCCGGCGCCGCCAGATCCCGGAGACCCCGATGACACAGCCGATGACCCAAGGGTTCCTGTCCAAAAAccc            <  1:3498580/93‑1 (MQ=255)                                                                                          tCCCGC.GGGGGT.CCGGGGG.CGCGAAG.CCCCCGGCGCCGCCAGATCCCGGAGACCCCGATGACACAGCCGATGACCCAAGGGTTCCTGTCCAAAAACCCCTc         >  1:3545993/1‑101 (MQ=255)                                                                                           cccGC.GGGGGT.CCGGGGG.CGCGAAG.CCCCCGGCGCCGCCAGATCCCGGAGACCCCGATGACACAGCCGATGACCCAAGGGTTCCTGTCCAAAAACCCCTCg        >  1:2591254/1‑101 (MQ=255)                                                                                           cccGC.GGGGGT.CCGGGGG.CGCGAAG.CCCCCGGCGCCGCCAGATCCCGGAGACCCCGATGACACAGCCGATGACCCAAGGGTTCCTGTCCAAAAACCCCTCg        >  1:1156345/1‑101 (MQ=255)                                                                                           cccGC.GGGGGT.CCGGGGG.CGCGAAG.CCCCCGGCGCCGCCAGATCCCGGAGACCCCGATGACACAGCCGATGACCCAAGGGTTCCTGTCCAAAAACCCCTCg        >  2:886349/1‑101 (MQ=255)                                                                                              gc.gGGGGT.CGGGGGC.CGCGAAG.CCCCCGGCGCCGCCAGATCCCGGAGACCCCGATGACACAGCCGATGACCCAAGGGTTCCTGTCCAAAAACCCCTCGGTg     <  2:2843233/101‑1 (MQ=255)                                                                                              gc.gGGGGT.CCGGGGG.CGCGAAG.CCCCCGGCGCCGCCAGATCCCGGAGACCCCGATGACACAGCCGATGACCCAAGGGTTCCTGTCCAAAAACCCCTCGGTg     >  1:3428951/1‑101 (MQ=255)                                                                                                 gggggT.CCGGGGG.CGCGAAG.CCCCCGGCGCCGCCAGATCCCGGAGACCCCGATGACACAGCCGATGACCCAAGGGTTCCTGTCCAAAAACCCCTCGGTGCt   >  2:2458813/1‑101 (MQ=255)                                                                                                 gggggT.CCGGGGG.CGCGAAG.CCCCCGGCGCCGCCAGATCCCGGAGACCCCGATGACACAGCCGATGACCCAAGGGTTCCTGTCCAAAAACCCCTCGGTGCt   >  1:3768930/1‑101 (MQ=255)                                                                                                  ggggT.CCGGGGG.CGCGAAG.CCCCCGGCGCCGCCAGATCCCGGAGACCCCGATGACACAGCCGATGACCCAAGGGTTCCTGTCCAAAAACCCCTCGGTGCTg  >  1:4211836/1‑101 (MQ=255)                                                                                                       |                                                                                                   CCGCCCGCGCGGCATCATCGCGCTGCGCGGCCGCGCCGCCGGAGATTGAGCCGATG.CCCCCCGGCG.CCTTTCATCTTGCCCAAAATATCCCGC.GGGGGTCCCGGGGG.CGCGAAG.CCCCCGGCGCCGCCAGATCCCGGAGACCCCGATGACACAGCCGATGACCCAAGGGTTCCTGTCCAAAAACCCCTCGGTGCTG  >  NC_014034/1300995‑1301190 Base quality scores:  ATCG  < 3 ≤  ATCG  < 17 ≤  ATCG  < 26 ≤  ATCG  < 33 ≤  ATCG  < 40 ≤  ATCG