BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_014034	1,395,760	T→G	16.0%	intergenic (+572/‑359)	RCAP_rcc01306 → / → RCAP_rcc01307	Cas2 family CRISPR‑associated protein/hypothetical protein

Read alignment evidence...
	seq id	position		change	freq	score	reads	annotation	genes	product
*	NC_014034	1,395,760	0	T→G	16.0%	125.8 ‑0.7 ‑0.2	288	intergenic (+572/‑359)	RCAP_rcc01306/RCAP_rcc01307	Cas2 family CRISPR‑associated protein/hypothetical protein
Reads supporting (aligned to +/- strand): new base (20/26): ref base (120/122): total (140/148)
Fisher's exact test strand distribution p-value = 5.21e-01
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 6.61e-01

Only 100 of 302 total aligned reads displayed.


GAAACATCTGCATGTAGTCGGCGGGCAGTTCGAGCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGGATCT.G..GCGCAAACGG..AAAGTCT.G.TCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCTGCGGCGTGAGTTGCCAGTGATAGACGCCGTGGACGTCGCTCCCCCCGCGGGGGCGTGGATCGA  >  NC_014034/1395660‑1395860


gAAACATCTGCATGTAGTCGGCGGGCAGTTCGAGCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGGATCT.G..GCGCAAACGG..AAAGTCt                                                                                                        >  2:4286436/1‑101 (MQ=255)


 aaaCATCTGCATGTAGTCGGCGGGCAGTTCGAGCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGGATCT.G..GCGCAAACGG..AAAGTCT.g                                                                                                      <  2:4403571/101‑1 (MQ=255)


 aaaCATCTGCATGTAGTCGGCGGGCAGTTCGAGCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGGATCT.G..GCGCAAACGG..AAAGTCT.g                                                                                                      <  1:89249/101‑1 (MQ=255)


    cATCTGCATGTAGTCGGCGGGCAGTTCGAGCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGGATCT.G..GCGCAAACGG..AAAGTCT.g.tcg                                                                                                  <  1:4909450/101‑1 (MQ=255)


    cATCTGCATGTAGTCGGCGGGCAGTTCGAGCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGGATCT.G..GCGCAAACGG..AAAGTCT.g.tcg                                                                                                  >  2:4155464/1‑101 (MQ=255)


     ctCTGCATGTGGCGGGCGGGCAGTTCGAGCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGGATCT.G..GCGCAAACGG..AAAGTCT.g.tcgt                                                                                                 <  2:1548907/100‑1 (MQ=255)


      tCTGCATGTAGTCGGCGGGCAGTTCGAGCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGGATCT.G..GCGCAAACGG..AAAGTCT.g.tcgtc                                                                                                >  2:2678765/1‑101 (MQ=255)


        tGCATGTAGTCGGCGGGCAGTTCGAGCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGGATCT.G..GCGCAAACGG..AAAGTCT.G.TCGTCGc                                                                                              >  2:2857642/1‑101 (MQ=255)


           aTGTAGTCGGCGGGCAGTTCGAGCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGGATCT.G..GCGCAAACGG..AAAGTCT.G.TCGTCGCTcc                                                                                           >  1:3816673/1‑101 (MQ=255)


             gTAGTCGGCGGGCAGTTCGAGCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGGATCT.G..GCGCAAACGG..AAAGTCT.G.TCGTCGCTcccc                                                                                         <  2:4599327/101‑1 (MQ=255)


               cggcggcggGCAGTTCGACCCGCTGTCGCTGCCCCCGCGGGGGCGTGGATCTAAACT.GGTTTGGATCT.G..GCGCAAACGG..AAAGTCT.G.TCGTCGCTcccccc                                                                                       <  1:798152/98‑1 (MQ=38)


                  cggcggGCAGTTCGAGCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGGATCT.G..GCGCAAACGG..AAAGTCT.G.TCGTCGCTCCCCCcgcg                                                                                    >  1:2091024/1‑101 (MQ=255)


                   ggcggGCAGTTCGAGCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGGATCT.G..GCGCAAACGG..AAAGTCT.G.TCGTCGCTCCCCCCGCgg                                                                                   >  2:2121547/1‑101 (MQ=255)


                      gggCAGTTCGAGCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACA.GGTATGGATCT.G..GCGCAAACGG..AAAGTCT.G.TCGTCGCTCCCCCCGCggggg                                                                                >  1:4274480/1‑101 (MQ=255)


                         cAGTTCGAGCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGGATCT.G..GCGCAAACGG..AAAGTCT.G.TCGTCGCTCCCCCCGCGGGGGCGt                                                                             <  1:2690571/101‑1 (MQ=255)


                          aGTTCGAGCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGGATCT.G..GCGCAAACGG..AAAGTCT.G.TCGTCGCTCCCCCCGCGGGGGCGTg                                                                            >  2:3345035/1‑101 (MQ=255)


                           gtgcGAGCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGGATCT.G..GCGCAAACGG..AAAGTCT.G.TCGTCGCTCCCCCCGCGGGGGCGTgg                                                                           <  1:2186852/98‑1 (MQ=255)


                           gTTCGAGCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGGATCT.G..GCGCAAACGG..AAAGTCT.G.TCGTCGCTCCCCCCGCGGGGGCGTgg                                                                           >  2:3049997/1‑101 (MQ=255)


                             tcaagattcAGTCGCGCCCCCCGCGGGGGCGTGGATCGAAACTCGGCAAGGACAT.C..GCCGATCTGGGTAAGGCCG.G.TGGTCGCTCCCCCCGCGGGGGCGTg                                                                            <  2:3251915/93‑1 (MQ=2)


                              cGAGCCGCAGTCGCCCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGGATCT.G..GCGCAAACGG..AAAGTCT.G.TCGTCGCTCCCCCCGCGGGGGCGTGGATc                                                                        <  2:4260367/101‑1 (MQ=255)


                               gAGCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGGATCT.G..GCGCAAACGG..AAAGTCT.G.TCGTCGCTCCCCCCGCGGGGGCGTGGATCg                                                                       <  1:4511956/101‑1 (MQ=255)


                               cagattcAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGGCAAGGACAT.C..GCCGATCTGGGTAAGGCCG.G.TGGTCGCTCCCCCCGCGGGGGCGTGGa                                                                          <  2:815732/95‑1 (MQ=2)


                                cgcgctggTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.TG.AT.GACCA.GATCCGCAATCGG..CCGGGCTCGATCGTCGCTCCCCCCGCGGGGGCGTGGATc                                                                        <  2:5018250/94‑1 (MQ=0)


                                 gcgctggTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.TG.AT.GACCA.GATCCGCAATCGG..CCGGGCTCGATCGTCGCTCCCCCCGCGGGGGCGTGGATCg                                                                       <  1:3977772/95‑1 (MQ=0)


                                 gattcAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGGCAAGGACAT.C..GCCGATCTGGGTAAGGCCG.G.TGGTCGCTCCCCCCGCGGGGGCGTGGATc                                                                        >  2:1027640/5‑101 (MQ=2)


                                  attcAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGGCAAGGACATCG..GCGATCTGGGT.AAGGCCG.G.TGGTCGCTCCCCCCGCGGGGGCGTGGATCg                                                                       >  1:2002408/4‑101 (MQ=2)


                                  attcAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGGCAAGGACAT.C..GCCGATCTGGGTAAGGCCG.G.TGGTCGCTCCCCCCGCGGGGGCGTGGATCg                                                                       >  1:212518/4‑101 (MQ=2)


                                    gCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGGATCT.G..GCGCAAACGG..AAAGTCT.G.TCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAAcc                                                                  >  2:273323/1‑101 (MQ=255)


                                    ctggTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.TG.AT.GACCA.GATCCGCAATCGG..CCGGGCTCGATCGTCGCTCCCCCCGCGGGGGCGTGGATCGaaa                                                                    <  2:4923984/98‑1 (MQ=2)


                                     cAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGGATCT.G..GCGCAAACGG..AAAGTCT.G.TCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCt                                                                 >  1:3220649/1‑101 (MQ=255)


                                     cAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGGATCT.G..GCGCAAACGG..AAAGTCT.G.TCGTCGCTCCCCCCGCGGGGGCGGGGATCGAAACCt                                                                 >  1:4532405/1‑101 (MQ=255)


                                      ggTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.TG.AT.GACCA.GATCCGCAATCGG..CCGGGCTCGATCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAAcc                                                                  >  1:130646/2‑101 (MQ=2)


                                      ggTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.TG.AT.GACCA.GATCCGCAATCGG..CCGGGCTCGATCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAAcc                                                                  >  2:2289756/2‑101 (MQ=2)


                                      aGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGGCAAGGACAT.C..GCCGATCTGGGTAAGGCCG.G.TGGTCGCTCCCTCCGCGGGGGCGTGGATCGAAAc                                                                   <  2:515778/101‑1 (MQ=2)


                                      aGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGGCAAGGACAT.C..GCCGATCTGGGTAAGGCCG.G.TGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAAc                                                                   <  2:4756481/101‑1 (MQ=2)


                                      aGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGGATCT.G..GCGCAAACGG..AAAGTCT.G.TCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCTg                                                                <  2:1488763/101‑1 (MQ=255)


                                       gTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGGATCT.G..GCGCAAACGG..AAAGTCT.G.TCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCTgc                                                               <  2:3744672/101‑1 (MQ=255)


                                       gTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGGATCT.G..GCGCAAACGG..AAAGTCT.G.TCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCTgc                                                               >  2:4007226/1‑101 (MQ=255)


                                       gTCGCTCCCCCCGCGGGGGAGTGGATCGAAACT.GGTATGGATCT.G..GCGTAAACGG..AAAGTCT.G.TCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCTgc                                                               <  1:2830236/101‑1 (MQ=255)


                                       gTCGCCCCCCCCGCGGGGGCGTGGATCGAAACTCGGCAAGGACAT.C..GCCGATCTGGGTAAGGCCG.G.TGGTCGCTCCCCCCGCGGGGGCGTGGATCGAGAca                                                                  <  1:1896147/101‑2 (MQ=2)


                                        tCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGGATCT.G..GCGCAAACGG..AAAGTCT.G.TCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCTgcg                                                              >  2:123139/1‑101 (MQ=255)


                                         cGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGGATCT.G..GCGCAAACGG..AAAGTCT.G.TCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCCgcgg                                                             >  1:2222896/1‑101 (MQ=255)


                                          gCTCCCCCCGCGGGGGCGTGGATCGAAACT.TG.AT.GACCA.GATCCGCAATCGG..CCGGGCTCGATCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCtatg                                                              >  1:4871950/1‑98 (MQ=2)


                                          gCTCCCCCCGCGGGGGCGTGGATCGAAACT.TG.AT.GACCA.GATCCGCAATCGG..CCGGGCTCGATCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCtatg                                                              >  1:378423/1‑98 (MQ=2)


                                          gCTCCCCCCGCGGGGGCGTGGATCGAAACT.TG.AT.GACCA.GATCCGCAATCGG..CCGGGCTCGATCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCtatg                                                              <  1:155601/101‑4 (MQ=2)


                                          gCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGGATCT.G..GCGCAAACGG..AAAGTCT.G.TCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCTgcggt                                                            <  2:4565843/101‑2 (MQ=255)


                                          gCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGGATCT.G..GCGCAAACGG..AAAGTCT.G.TCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCTgcggc                                                            <  1:618059/101‑1 (MQ=255)


                                          gCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGGATCT.G..GCGCAAACGG..AAAGTCT.G.TCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCTgcggc                                                            >  2:4992675/1‑101 (MQ=255)


                                           ccccccccGGGGGGGCGTGGATCGAAACTCGGCAAGGACAT.C..GCCGATCTGGGTAAGGCCG.G.TGGTCGCCCCCCCCGCGGGGGCGTGGATCGAAAcacccg                                                              <  1:2651617/99‑6 (MQ=2)


                                           cTCCCCCCGCGGGGGCGTGGATCGAAACTCGGCAAGGACAT.C..GCCGATCTGGGTAAGGCCG.G.TGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAAcacccg                                                              <  2:4893434/101‑6 (MQ=2)


                                             gatctcgcgGGGGCGTGGATCGAAACT.GGTATGGATCT.G..GCGCAAACGG..AAAGTCT.G.TCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCTGCGGCGTg                                                         <  1:5017625/96‑1 (MQ=255)


                                             ccccccGCGGGGGCGTGGATCGAAACT.GGTATGGATCT.G..GCGCAAACGG..AAAGTCT.G.TCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCTGCGGCGTg                                                         >  1:3330049/1‑101 (MQ=255)


                                                 ccGCGGGGGCGTGGATCGAAACT.GGTATGGATCT.G..GCGCAAACGG..AAAGTCT.G.TCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCTGCGGCGTGAGtt                                                     <  2:3792730/101‑1 (MQ=255)


                                                    cgGGGGCGTGGATCGAAACT.GGTATGGATCT.G..GCGCAAACGG..AAAGTCT.G.TCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCTGCGGCGTGAGTTGcc                                                  <  1:2679592/101‑1 (MQ=255)


                                                     gggggCGTGGATCGAAACT.GGTATGGATCT.G..GCGCAAACGG..AAAGTCT.G.TCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCTGCGGCGTGAGTTGCCa                                                 >  1:2319675/1‑101 (MQ=255)


                                                        ggCGTGGATCGAAACT.GGTATGGATCT.G..GCGCAAACGG..AAAGTCT.G.TCGTCGCTCCCCCCGCGGGGGCGCGGATCGAAACCTGCGGCGTGAGTTGCCAGTg                                              >  1:3621578/1‑101 (MQ=255)


                                                          cGTGGATCGAAACT.GGTATGGATCT.G..GCGCAAACGG..AAAGTCT.G.TCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCTGCGGCGTGAGTTGCCAGTGAt                                            >  2:2051248/1‑101 (MQ=255)


                                                          cGTGGATCGAAACT.GGTATGGATCT.G..GCGCAAACGG..AAAGTCT.G.TCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCTGCGGCGTGAGTTGCCAGTGAt                                            <  1:2384191/101‑1 (MQ=255)


                                                           gTGGATCGAAACT.GGTATGGATCT.G..GCGCAAACGG..AAAGTCT.G.TCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCTGCGGCGTGAGTTGCCAGTGATa                                           >  2:1974341/1‑101 (MQ=255)


                                                            gggATCGAAACT.GGTATGGATCT.G..GCGCAAACGG..AAAGTCT.G.TCGTCGCT.CCCCCGCGGGGGCGTGGATCGAAACCTGCGGCGTGAGTTGCCAGTGATAGa                                         <  1:3491890/100‑1 (MQ=255)


                                                             ggATCGAAACT.GGTATGGATCT.G..GCGCAAACGG..AAAGTCT.G.TCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCTGCGGCGTGAGTTGCCAGTGATAGa                                         >  1:4601131/1‑101 (MQ=255)


                                                                tCGAAACT.GGTATGGATCT.G..GCGCAAACGG..AAAGTCT.G.TCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCTGCGGCGTGAGTTGCCAGTGATAGACGc                                      <  1:4627997/101‑1 (MQ=255)


                                                                   aaaCT.GGTATGGATCT.G..GCGCAAACGG..AAAGTCT.G.TCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCTGCGGCGTGAGTTGCCAGTGATAGACGCCgg                                   >  1:1093676/1‑100 (MQ=255)


                                                                   aaaCT.GGTATGGATCT.G..GCGCAAACGG..AAAGTCT.G.TCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCTGCGGCGTGAGTTGCCAGTGATAGACGCCGt                                   >  2:1840929/1‑101 (MQ=255)


                                                                     aCT.GGTATGGATCT.G..GCGCAAACGG..AAAGTCT.G.TCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCTGCGGCGTGAGTTGCCAGTGATAGACGCCGTgg                                 >  1:1677645/1‑101 (MQ=255)


                                                                     aCT.GGTATGGATCT.G..GCGCAAACGG..AAAGTCT.G.TCGTCGCTCCCCCCGCGGGGGCGCGGATCGAAACCTGCGGCGTGAGTTGCCAGTGATAGACGCCGTgg                                 >  1:2971400/1‑101 (MQ=255)


                                                                      cT.GGTATGGATCT.G..GCGCAAACGG..AAAGTCT.G.TCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCTGCGGCGTGAGTTGCCAGTGATAGACGCCGTGGa                                >  2:2219080/1‑101 (MQ=255)


                                                                      cT.GGTATGGATCT.G..GCGCAAACGG..AAAGTCT.G.TCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCTGCGGCGTGAGTTGCCAGTGATAGACGCCGTGGa                                >  2:4311946/1‑101 (MQ=255)


                                                                       t.GGTATGGATCT.G..GCGCAAACGG..AAAGTCT.G.TCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCTGCGGCGTGAGTTGCCAGTGATAGACGCCGTGGAc                               <  1:3282583/101‑1 (MQ=255)


                                                                       t.GGTATGGATCT.G..GCGCAAACGG..AAAGTCT.G.TCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCAGCGGCGTGAGTTGCCAGTGATAGACGCCGTGGAc                               >  1:4332930/1‑101 (MQ=255)


                                                                          gTATGGATCT.G..GCGCAAACGG..AAAGTCT.G.TCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCTGCGGCGTGAGTTGCCAGTGATAGACGCCGTTGACGt                             <  2:1928784/101‑1 (MQ=255)


                                                                          gTATGGATCT.G..GCGCAAACGG..AAAGTCT.G.TCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCTGCGGCGTGAGTTGCCAGTGATAGACGCCGTGGACGt                             <  2:3621307/101‑1 (MQ=255)


                                                                           tATGGATCT.G..GCGCAAACGG..AAAGTCT.G.TCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCTGCGGCGTGAGTTGCCAGTGATAGACGCCGTGGACGTc                            >  1:3443522/1‑101 (MQ=255)


                                                                           tATGGATCT.G..GCGCAAACGG..AAAGTCT.G.TCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCTGCGGCGTGAGTTGCCAGTGATAGACGCCGTGGACGTc                            >  1:2558902/1‑101 (MQ=255)


                                                                            aTGGATCT.G..GCGCAAACGG..AAAGTCT.G.TCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCTGCGGCGTGAGTTGCCAGTGATAGACGCCGTGGACgctg                           >  1:2485892/1‑98 (MQ=255)


                                                                              ggATCT.G..GCGCAAACGG..AAAGTCT.G.TCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCTGCGGCGTGAGTTGCCAGTGATAGACGCCGTGGTCGTCGCt                         >  2:4501529/1‑101 (MQ=255)


                                                                              ggATCT.G..GCGCAAACGG..AAAGTCT.G.TCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCTGCGGCGTGAGTTGCCAGTCATAGACGCCGTGGACGTCGCt                         >  2:3399631/1‑101 (MQ=255)


                                                                              ggATCT.G..GCGCAAACGG..AAAGTCT.G.TCGTCGCTCCCCCCGCGGGGGCGTGGATCCAAACCTGCGGCGTGAGTTGCCAGTGATAGACGCCGTGGACGTCGCt                         <  2:4362772/101‑1 (MQ=255)


                                                                                     g..gCGCAAACGG..AAAGTCT.G.TCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCTGCGGCGTGAGTTGCCAGTGATAGACGCCGTGGACGTCGCTcccccc                   >  2:386676/1‑101 (MQ=255)


                                                                                     g..gCGCAAACGG..AAAGTCT.G.TCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCTGCGGCGTGAGTTGCCAGTGATAGACGCCGTGGACGTCGCTcccccc                   <  1:3050235/101‑1 (MQ=255)


                                                                                     g..gCGCAAACGG..AAAGTCT.G.TCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCTGCGGCGTGAGTTGCCAGTGATAGACGCCGTGGACGTCGCTcccccc                   >  2:4673375/1‑101 (MQ=255)


                                                                                     g..gCGCAAACGG..AAAGTCT.G.TCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCTGCGGCGTGAGTTGCCAGTGATAGACGCCGTGGACGTCGCTcccccc                   <  2:2233526/101‑1 (MQ=255)


                                                                                     g..gCGCAAACGG..AAAGTCT.G.TCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCTGCGGCGTGAGTTGCCAGTGATAGACGCCGTGGACGTCGCTcccccc                   <  1:929284/101‑1 (MQ=255)


                                                                                     g..gCGCAAACGG..AAAGTCT.G.TCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCTGCGGCGTGAGTTGCCAGTGATAGACGCCGTGGACGTCGCTcccccc                   <  1:4949833/101‑1 (MQ=255)


                                                                                        gcgcAAACGG..AAAGTCT.G.TCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCTGCGGCGTGAGTTGCCAGTGATAGACGCCGTGGACGTCGCTccccccc                  >  2:352936/1‑100 (MQ=255)


                                                                                        gcgcAAACGG..AAAGTCT.G.TCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCTGCGGCGTGAGTTGCCAGTGATAGACGCCGTGGACGTCGCTCCCCCcg                  <  2:369860/101‑1 (MQ=255)


                                                                                             aaCGG..AAAGTCT.G.TCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCTGCGGCGTAAGTTGCCAGTGATAGACGCCGTGGACGTCGCTCCCCCCGCgggg             <  1:4029958/101‑1 (MQ=255)


                                                                                               cGG..AAAGTCT.G.TCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCTGCGGCGTGAGTTGCCAGTGATAGACGCCGTGGACGTCGCTCCCCCCGCggggcc           >  2:2355307/1‑99 (MQ=255)


                                                                                               cGG..AAAGTCT.G.TCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCTGCGGCGTGAGTTGCCAGTGATAGACGCCGTGGACGTCGCTCCCCCCGCGGGGGc           <  1:1008604/101‑1 (MQ=255)


                                                                                               cGG..AAAGTCT.G.TCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCTGCGGCGTGAGTTGCCAGTGATAGACGCCGTGGACGTCGCTCCCCCCGCGGGGGc           >  1:1230956/1‑101 (MQ=255)


                                                                                               cGG..AAAGTCT.G.TCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCTGCGGCGTGAGTTGCCAGTGATAGACGCCGTGGACGTCGCTCCCCCCGCGGGGGc           >  1:4783394/1‑101 (MQ=255)


                                                                                                gg..AAAGTCT.G.TCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCTGCGGCGTGAGTTGCCAGTGATAGACGCCGTGGACGTCGCTCCCCCCGCGGGGGCg          >  1:3101266/1‑101 (MQ=255)


                                                                                                 g..AAAGTCT.G.TCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCTGCGGCGTGAGTTGCCAGTGATAGACGCCGTGGACGTCGCTCCCCCCGCGGGGGCgg         >  1:1446326/1‑100 (MQ=255)


                                                                                                    aaaGTCT.G.TCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCTGCGGCGTGAGTTGCCAGTGATAGACGCCGTGGACGTCGCTCCCCCCGCGGGGGCGTg        <  1:2752851/101‑1 (MQ=255)


                                                                                                     aaGTCT.G.TCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCTGCGGCGTGAGTTGCCAGTGATAGACGCCGTGGACGTCGCTCCCCCCGCGGGGGCGTgg       <  1:3665525/101‑1 (MQ=255)


                                                                                                      aGTCT.G.TCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCTGCGGCGTGAGTTGCCAGTGATAGCCGCCGTGGACGCCGCTCCCCCCCGGGCGGCGTGGa      >  2:4203924/1‑101 (MQ=37)


                                                                                                        tCT.G.TCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCTGCGGCGTGAGTTGCCAGTGATAGACGCCGTGGACGTCGCTCCCCCCGCGGGGGCGTGGATc    >  2:647316/1‑101 (MQ=255)


                                                                                                        tCT.G.TCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCTGCGGCGTGAGTTGCCAGGGATAGACGCCGTGGACGTCGCCCCCCCCCCGGGGGCGTGGATc    >  1:1781594/1‑101 (MQ=255)


                                                                                                         cT.G.TCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCTGCGGCGTGAGTTGCCAGTGATAGACGCCGTGGACGTCGCTCCCCCCGCGGGGGCGTGGATCg   <  2:253847/101‑1 (MQ=255)


                                                                                                         cT.G.TCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCTGCGGCGTGAGTTGCCAGTGATAGACGCCGTGGACGTCGCTCCCCCCGCGGGGGCGTGGATCg   <  1:4212407/101‑1 (MQ=255)


GAAACATCTGCATGTAGTCGGCGGGCAGTTCGAGCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGGATCT.G..GCGCAAACGG..AAAGTCT.G.TCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCTGCGGCGTGAGTTGCCAGTGATAGACGCCGTGGACGTCGCTCCCCCCGCGGGGGCGTGGATCGA  >  NC_014034/1395660‑1395860

Base quality scores:


ATCG

< 3 ≤


ATCG

< 18 ≤


ATCG

< 27 ≤


ATCG

< 33 ≤


ATCG

< 40 ≤


ATCG