BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_014034	2,494,014	A→G	6.2%	intergenic (+62/‑180)	RCAP_rcc02319 → / → RCAP_rcc02320	hypothetical protein/hypothetical protein

Read alignment evidence...
	seq id	position		change	freq	score	reads	annotation	genes	product
*	NC_014034	2,494,014	0	A→G	6.2%	15.3 ‑1.7 ‑1.1	145	intergenic (+62/‑180)	RCAP_rcc02319/RCAP_rcc02320	hypothetical protein/hypothetical protein
Reads supporting (aligned to +/- strand): new base (6/3): ref base (57/79): total (63/82)
Fisher's exact test strand distribution p-value = 1.77e-01
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 8.04e-02

Only 100 of 154 total aligned reads displayed.


AGCGCCGCGGCGGGCGCGCGGCACGACAGGCCGAGTAA.GGGGG.G.GCGCTGCCCCCGTTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGC..G.AG.CAGA.AGGGGGGCAAGCG.CCCCCTGCCCCGGCGGTCCGTGGCGGCGCCT  >  NC_014034/2493915‑2494111


aGCGCGGCGGCGGGCGCGCGGCACGACAGGCCGAGTAAGGGGGG.G.GCGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGaaa                                                                                                          <  2:1850840/101‑1 (MQ=18)


aGCGCCGCGGCGGGCGCGCGGCACGACAGGCCGAGTAAGGGGGG.G.GCGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGaaa                                                                                                          <  2:1071748/101‑1 (MQ=21)


aGCGCCGCGGCGGGCGCGCGGCACGACAGGCCGAGTAAGGGGGG.G.GCGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGaaa                                                                                                          <  1:3263804/101‑1 (MQ=21)


  ctccgcGGCGGGCGCGCGGCACGACAGGCCGAGTAAGGGGGG.G.GCGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAAcc                                                                                                        <  2:720224/99‑1 (MQ=18)


     cgcgGCGGGCGCGCGGCACGACAGGCCGAGTAAGGGGGG.G.GCGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGa                                                                                                     <  1:3827145/101‑1 (MQ=21)


     cgcgGCGGGCGCGCGGCACGACAGGCCGAGTAAGGGGGG.G.GCGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGa                                                                                                     <  1:2415967/101‑1 (MQ=21)


     cgcgCCGGGCGCGCGCCCCGACAGGCCGAGTAAGGGGGG.G.GCGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGa                                                                                                     <  2:4138218/101‑1 (MQ=17)


      gcggcgGGCGCGCGGCACGACAGGCCGAGTAAGGGGGG.G.GCGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGAt                                                                                                    <  1:3873785/101‑1 (MQ=21)


      gcggcgGGCGCGCGGCACGACAGGCCGAGTAAGGGGGG.G.GCGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGAt                                                                                                    <  2:3799038/101‑1 (MQ=21)


        ggcgggcgCGCGGCACGACAGGCCGAGTAAGGGGGG.G.GCGCTGCCCCCCGTCCGCAAGCGGACTCCCACCGGGAGATTTGCAGCAAGATGAAACAGGCTGc                                                                                                  >  1:4013220/1‑101 (MQ=18)


         gcgggcgCGCGGCACGACAGGCCGAGTAAGGGGGG.G.GCGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGcc                                                                                                 <  2:3405661/101‑1 (MQ=21)


          cgggcgCGCGGCACGACAGGCCGAGTAAGGGGGG.G.GCGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCt                                                                                                <  1:576105/101‑1 (MQ=21)


           gggcgCGCGGCACGACAGGCCGAGTAAGGGGGG.G.GCGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCtt                                                                                               <  1:4577252/101‑1 (MQ=21)


           gggcgCGCGGCACGACAGGCCGAGTAAGGGGGG.G.GCGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCtt                                                                                               <  2:3907789/101‑1 (MQ=21)


           gggcgCGCGGCACGACAGGCCGAGTAAGGGGGG.G.GCGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCtt                                                                                               <  2:3691524/101‑1 (MQ=21)


            ggcgCGCGGCACGACAGGCCGAGTAAGGGGGG.G.GCGCTGCCCCGCGTCCGCAAGCGGACTCCCCCCGGGGTATTTGCAGCAAGATGAAACCGGATGCCttt                                                                                              <  1:1624367/101‑1 (MQ=21)


            ggcgCGCGGCACGACAGGCCGAGTAAGGGGGG.G.GCGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCttt                                                                                              <  2:399985/101‑1 (MQ=21)


                cgcgGCACGACAGGCCGAGTAAGGGGGG.G.GCGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATc                                                                                          >  2:2974540/1‑101 (MQ=21)


                  cgGCACGACAGGCCGAGTAAGGGGGG.G.GCGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCtt                                                                                        <  1:1966740/101‑1 (MQ=21)


                  cgGCACGACAGGCCGAGTAAGGGGGG.G.GCGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCtt                                                                                        <  1:2097046/101‑1 (MQ=21)


                   ggCACGACAGGCCGAGTAAGGGGGG.G.GCGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTg                                                                                       >  2:2736730/1‑101 (MQ=21)


                   ggCACGACAGGCCGAGTAAGGGGGG.G.GCGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTg                                                                                       >  2:1364262/1‑101 (MQ=21)


                     cACGACAGGCCGAGTAAGGGGGG.G.GCGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGcc                                                                                     >  2:1311081/1‑101 (MQ=21)


                      aCGACAGGCCGAGTAAGGGGGG.G.GCGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATTTTGccc                                                                                    >  2:3578950/1‑101 (MQ=18)


                       cGACAGGCCGAGTAAGGGGGG.G.GCGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCg                                                                                   <  1:2581620/101‑1 (MQ=21)


                       cGACAGGCCGAGTAAGGGGGG.G.GCGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCg                                                                                   <  1:2796169/101‑1 (MQ=21)


                       cGACAGGCCGAGTAAGGGGGG.G.GCGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCg                                                                                   <  2:4362146/101‑1 (MQ=21)


                       cGACAGGCCGAGTAAGGGGGG.G.GCGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCg                                                                                   <  1:1254986/101‑1 (MQ=21)


                        gACAGGCCGAGTAAGGGGGG.G.GCGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGa                                                                                  <  2:110926/101‑1 (MQ=21)


                        gACAGGCCGAGTAAGGGGGG.G.GCGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGa                                                                                  <  1:3525447/101‑1 (MQ=21)


                        gACAGGCCGAGTAAGGGGGG.G.GCGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGACGCCTTTCATCTTTGCCGa                                                                                  >  2:1674435/1‑101 (MQ=17)


                         aCAGGCCGAGTAAGGGGGG.G.GCGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGaa                                                                                 <  2:4427230/101‑1 (MQ=21)


                          cAGGCCGAGTAAGGGGGG.G.GCGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGaaa                                                                                <  2:1981220/101‑1 (MQ=21)


                           aGGCCGAGTAAGGGGGG.G.GCGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAat                                                                               <  2:647773/101‑1 (MQ=21)


                            ggCCGAGTAAGGGGGG.G.GCGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAata                                                                              <  2:3559894/101‑1 (MQ=21)


                            ggCCGAGTAA.GGGGG.GCGGTCGGCCCCGCGCCCGCAAGGGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCCTCTTGCCCGAAata                                                                              <  2:2306986/101‑1 (MQ=17)


                             gCCGAGTAAGGGGGG.G.GCGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAatat                                                                             <  1:3778516/101‑1 (MQ=21)


                              tcGAGTAAGGGGGG.G.GCGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATc                                                                            <  1:4054669/100‑1 (MQ=18)


                               cGAGTAAGGGGGG.G.GCGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATcc                                                                           >  2:31263/1‑101 (MQ=18)


                               cGAGTAAGGGGGG.G.GCGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATcc                                                                           <  1:1686562/101‑1 (MQ=18)


                               cGAGTAAGGGGGG.G.GCCCTGCCCCCCGTCCGCAAGCCGACTCCCCCCGGCATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATcc                                                                           <  2:3151033/101‑1 (MQ=18)


                               agAGAAA.GGGGGAG.TCGCTGCCCCCCGTCCACAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATcc                                                                           <  1:530034/100‑1 (MQ=18)


                                gAGTAAGGGGGG.G.GCGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATccc                                                                          <  2:3429093/101‑1 (MQ=18)


                                gAGTAAGGGGGG.G.GCGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATccc                                                                          >  1:3418496/1‑101 (MQ=18)


                                 aGTAAGGGGGG.G.GCGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCg                                                                         >  1:1000152/1‑101 (MQ=18)


                                  gtaagggggg.g.gCGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCgg                                                                        >  1:1486152/6‑101 (MQ=255)


                                  gtaagggggg.g.gCGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCgg                                                                        <  1:1350477/96‑1 (MQ=255)


                                   taagggggg.g.gCGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCCGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCggg                                                                       >  1:3180596/5‑101 (MQ=18)


                                    aagggggg.g.gGCCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCgggg                                                                      <  2:1000069/98‑1 (MQ=18)


                                    aagggggg.g.gCGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCgggg                                                                      >  1:1628290/4‑101 (MQ=255)


                                      gggggg.g.gCGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCgggggg                                                                    >  1:2629620/2‑101 (MQ=255)


                                       ggggg.g.gCGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCggggggg                                                                   >  2:3028533/1‑100 (MQ=255)


                                       ggggg.g.gCGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGt                                                                   <  2:960816/101‑1 (MQ=255)


                                        gggg.g.gCGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCgggggggg                                                                  >  2:1621830/1‑99 (MQ=255)


                                        gggg.g.gCGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTg                                                                  <  1:1668650/101‑1 (MQ=255)


                                        gggg.g.gCGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTg                                                                  <  1:4511415/101‑1 (MQ=255)


                                         ggg.g.gCGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGa                                                                 >  1:3823557/1‑101 (MQ=255)


                                         ggg.g.gCGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGa                                                                 <  2:2208328/101‑1 (MQ=255)


                                                 gcTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGC..g                                                        >  2:2159991/1‑101 (MQ=255)


                                                  cTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGGGAA.TGGG..G.a                                                      >  1:4349045/1‑101 (MQ=18)


                                                  cTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGAAGCAAGATGAAACCGGATGCCTTTCATTTTGACCGAAATATCCCGGGGGGTGAA.TGGC..G.a                                                      >  2:2703254/1‑101 (MQ=18)


                                                   tGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGC..G.aa                                                     >  1:14024/1‑100 (MQ=255)


                                                    gCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGC..G.aag                                                    <  1:704842/101‑3 (MQ=255)


                                                     cccccG.TCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGC..G.aagcc                                                  <  2:1568902/101‑5 (MQ=255)


                                                     cccccG.TCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGC..G.aagcc                                                  >  1:422015/1‑97 (MQ=255)


                                                     cccccG.TCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGC..G.aagcc                                                  <  1:434851/101‑5 (MQ=255)


                                                     cccccCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGC..G.aagc                                                   <  1:198496/101‑4 (MQ=255)


                                                      ccccG.TCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGC..G.aagcca                                                 <  2:3670665/101‑6 (MQ=255)


                                                         ccgtCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGC..G.aagccaga                                               <  2:432097/98‑8 (MQ=255)


                                                           gtCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGC..G.aagccagagg                                             >  2:4129599/2‑92 (MQ=255)


                                                                 aaGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGC..GAAG.CCAG.AGGGGGGCa                                      <  1:814009/101‑1 (MQ=18)


                                                                  aGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGC..GAAG.CCAG.AGGGGGGCag                                     <  2:1264520/101‑2 (MQ=18)


                                                                    cGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGC..GAAG.CCAG.AGGGGGGCagcg                                   >  1:2943247/1‑98 (MQ=18)


                                                                       aCTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGG..GAAG.CCAG.AGGGGGGC.AGCG.ccc                              >  2:4255159/1‑101 (MQ=17)


                                                                        cTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGC..GAAG.CCAG.AGGGGGGC.AGCG.cccc                             >  2:486759/1‑101 (MQ=18)


                                                                         tCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGC..GAAGCCAG..AGGGGGGC.AGCG.ccccc                            <  2:1714876/101‑1 (MQ=17)


                                                                         tCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGC..GAAG.CCAG.AGGGGGGC.AGCG.ccccc                            >  1:2173295/1‑101 (MQ=18)


                                                                           cccccGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGC..GAAG.CCAG.AGGGGGGC.AGCG.cccccct                          <  1:2160141/101‑3 (MQ=17)


                                                                           cccccGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGC..GAAG.CCAG.AGGGGGGC.AGCG.cccccct                          <  2:2243612/101‑3 (MQ=17)


                                                                               cGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGC..G.AA.GACAGAGGGGGGC.AGCG.cccccctgacc                      >  1:1814513/1‑95 (MQ=17)


                                                                                 ggATATTTGCGGCAAGATGAAGGGGCATG.ATTTCATCTTGCCCGAAATATCCC.AGGGGTGAATTGGCCAA.AG.GCCA.AGAGGGGGCAGCG.CCCCCTGCCCcg                     <  1:2338829/101‑1 (MQ=0)


                                                                                 ggATATTTGCGGCAAGATGAAGGGGCATG.ATTTCATCTTGCCCGAAATATCCC.AGGGGTGAATTGGCCAA.AG.GCCA.AGAGGGGGCAGCG.CCCCCTGCCCcg                     >  2:1903843/1‑101 (MQ=0)


                                                                                 ggATATTTGCGGCAAGATGAAGGGGCATG.ATTTCATCTTGCCCGAAATATCCC.AGGGGTGAATTGGCCAA.AG.GCCA.AGAGGGGGCAGCG.CCCCCTGCCCcg                     >  1:1282178/1‑101 (MQ=2)


                                                                                 ggATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGC..GAAG.CCAG.AGGGGGGC.AGCGCCCCCCTGCCCcg                     >  1:991533/1‑101 (MQ=18)


                                                                                  gATATTTGCGGCAAGATGAAGGGGCATG.ATTTCATCTTGCCCGAAATATCCC.AGGGGTGAATTGGCCAA.AG.GCCA.AGAGGGGGCAGCG.CCCCCTGCCCcgc                    >  2:4290240/1‑100 (MQ=2)


                                                                                   atatTTGCGGCAAGATGAAGGGGCATG.ATTTCATCTTGCCCGAAATATCCC.AGGGGTGAATTGGCCAA.AG.GCCA.AGAGGGGGCAGCG.CCCCCTGCCCcgcc                   <  1:1127028/101‑3 (MQ=2)


                                                                                   atatTTGCGGCAAGATGAAGGGGCATG.ATTTCATCTTGCCCGAAATATCCC.AGGGGTGAATTGGCCAA.AG.GCCA.AGAGGGGGCAGCG.CCCCCTGCCCcgcc                   <  2:2444885/101‑3 (MQ=2)


                                                                                    tatTTGCGGCAAGATGAAGGGGCATG.ATTTCATCTTGCCCGAAATATCCC.AGGGGTGAATTGGCCAA.AG.GCCA.AGAGGGGGCAGCG.CCCCCTGCCCcgccg                  <  2:2091493/101‑1 (MQ=2)


                                                                                     atTTGCGGCAAGATGAAGGGGCATG.ATTTCATCTTGCCCGAAATATCCC.AGGGGTGAATTGGCCAA.AG.GCCA.AGAGGGGGCAGCG.CCCCCTGCCCcgccgg                 >  1:4479250/1‑101 (MQ=2)


                                                                                       ttGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGC..GAAG.CCAG.AGGGGGGCAGCGG.CCCCCGGCCCCGGCGGTc               >  2:714330/1‑101 (MQ=18)


                                                                                       ttGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGC..GAAG.CCAG.AGGGGGGC.AGCGCCCCCCTGCCCCGGCGGTc               <  1:1467624/101‑1 (MQ=17)


                                                                                        tGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGC..GAAG.CCAG.AGGGGGGC.AGCGCCCCCCTGCCCCGGCGGTcc              >  1:3926664/1‑101 (MQ=18)


                                                                                         gcagcaAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGC..GAAG.CCAC.AGGGGGGC.AGCGCCCCCCTGCCCCGGCGGTCCg             >  1:4037011/1‑101 (MQ=18)


                                                                                          cagcaAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGC..GAAG.CCAG.AGGGGGGC.AGCG.CCCCCCGCCCCCGCGGGcccgg           >  1:4092509/1‑98 (MQ=17)


                                                                                          cagcaAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGCGAA.TGGC..GAAG.CCAC.AGGGGGGC.AGCGCCCCCCTGCCCCGGCGGTCCGt            >  1:2256617/1‑101 (MQ=18)


                                                                                            gcaAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAG.TGGC..GAAG.CCAG.AGGGGGGC.AGCGCCCCCCTGCCCCGGCGGTCCGTgg          >  2:294835/1‑101 (MQ=18)


                                                                                             caAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGC..GAAG.CCAG.AGGGGGGC.AGCGCCCCCCTGCCCCGGCGGTCCGTggc         >  2:2049065/1‑101 (MQ=18)


                                                                                             caAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGAGAA.TGGC..GAAG.CCAG.AGGGGGGCCGACC.CCCCCTGCCCCGGCGGTCCGCggg         >  1:677081/1‑100 (MQ=17)


                                                                                               aGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGC..GAAG.CCAG.AGGGGGGC.AGCGCCCCCCTGCCCCGGCGGTCCGTggcgg       <  2:3823272/101‑1 (MQ=17)


                                                                                                  tGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGC..GAAG.CCAG.AGGGGGGC.AGCGCCCCCCTGCCCCGGCGGTCCGTGGCGgcgc    >  1:4604971/1‑101 (MQ=17)


AGCGCCGCGGCGGGCGCGCGGCACGACAGGCCGAGTAA.GGGGG.G.GCGCTGCCCCCGTTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGC..G.AG.CAGA.AGGGGGGCAAGCG.CCCCCTGCCCCGGCGGTCCGTGGCGGCGCCT  >  NC_014034/2493915‑2494111

Base quality scores:


ATCG

< 3 ≤


ATCG

< 18 ≤


ATCG

< 25 ≤


ATCG

< 33 ≤


ATCG

< 41 ≤


ATCG