BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_014034	1,395,732	G→T	8.4%	intergenic (+544/‑387)	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,732	0	G→T	8.4%	50.9 ‑0.4 ‑0.1	261	intergenic (+544/‑387)	RCAP_rcc01306/RCAP_rcc01307	Cas2 family CRISPR‑associated protein/hypothetical protein
Reads supporting (aligned to +/- strand): new base (10/12): ref base (123/116): total (133/128)
Fisher's exact test strand distribution p-value = 6.59e-01
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 8.85e-01

Only 100 of 274 total aligned reads displayed.


CGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCGGGCAGTTCGAGCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGGATCT.G..GCGCAAACGG..AAAGTCT.G.TCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCTGCGGCGTGAGTTGCCAGTGATAGACGCCGTGGACG  >  NC_014034/1395632‑1395832


ggccGCGCCCCCGGGGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCGGGCAGTTCGAGCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.g                                                                                                             <  1:598615/98‑1 (MQ=37)


cggcGCCCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCGGGCAGTTCGAGCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.g                                                                                                             <  1:3325631/98‑1 (MQ=255)


cGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCGGGCAGTTCGAGCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.g                                                                                                             >  1:14741/1‑101 (MQ=255)


 gTCGCCCCCCCCGCGGGGGCGTGGATCAAAACATCTGCATGTAGTGGGCGGGCAGTCCGAGCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.gg                                                                                                            <  1:848499/101‑1 (MQ=37)


  ccGCCCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCGGGCAGTTCGAGCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGt                                                                                                           <  1:1816451/100‑1 (MQ=255)


   cGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCGGGCAGTTCGAGCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.ggga                                                                                                          >  2:3414626/1‑99 (MQ=255)


          cccGCGGGGGGGTGGATCGAAACATCTGCATGTAGTCGGCGGGCAGTTCGAGCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGGATCt                                                                                                   <  1:2698946/101‑1 (MQ=255)


          cccGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCGGGCAGTTCGAGCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGGATCt                                                                                                   <  1:3439648/101‑1 (MQ=255)


            cgcgGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCGGGCAGTTCGAGCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGGATCT.g..g                                                                                              >  1:2612570/1‑101 (MQ=255)


            cgcgGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCGGGCAGTTCGAGCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGGATCT.g..g                                                                                              <  2:2499199/101‑1 (MQ=255)


               gggggCGTGGATCGAAACATCTGCATGTAGTCGGCGGGCAGTTCGAGCCTCAGTCGATTCCCCCGCGGGGGCGTGGATCGAAACT.GTTATGGATGT.gactc                                                                                             >  2:2072258/1‑97 (MQ=37)


               gggggCGTGGATCGAAACATCTGCATGTAGTCGGCGGGCAGTTCGAGCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGGATCT.G..gcgc                                                                                           >  2:3200016/1‑101 (MQ=255)


                    cGTGGATCGAAACATCTGCATGTAGTCGGCGGGCAGTTCGAGCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGGATCT.G..GCGCAAACg                                                                                      >  2:3742661/1‑101 (MQ=255)


                    cGTGGATCGAAACATCTGCATGTAGTCGGCGGGCAGTTCGAGCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGGATCT.G..GCGCAAACg                                                                                      >  1:1783396/1‑101 (MQ=255)


                      tGGATCGAAACATCTGCATGTAGTCGGCGGGCAGTTCGAGCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGGATCT.G..GCGCAAACGG..a                                                                                  >  2:1629877/1‑101 (MQ=255)


                      gggATCGAACCATCTGCACGTAGTCGGCGGGCAGTTCGAGCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGGATCT.G..GCGCAAACGG..a                                                                                  <  2:1913152/100‑1 (MQ=255)


                       ggATCGAAACATCTGCATGTAGTCGGCGGGCAGTTCGAGCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGGATCT.G..GCGCAAACGG..aa                                                                                 <  1:3194024/101‑1 (MQ=255)


                        gATCGAAACATCTGCATGTAGTCGGCGGGCAGTTCGAGCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGGATCT.G..GCGCAAACGG..aaa                                                                                >  2:1187649/1‑101 (MQ=255)


                        gATCGAAACATCTGCATGTAGTCGGCGGGCAGTTCGAGCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGGATCT.G..GCGCAAACGG..aaa                                                                                <  2:3930694/101‑1 (MQ=255)


                        gATCGAAACATCTGCATGTAGTCGGCGGGCAGTTCGAGCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAAAT.GGTATCGGTCT.G..GCGCAAACGG..aag                                                                                >  2:3092418/1‑100 (MQ=255)


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


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


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


                              aaCATCTGCATGTAGTCGGCGGGCAGTTCGAGCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGGATCT.G..GCGCAAACGG..AAAGTCT.g.t                                                                        <  1:2043358/101‑1 (MQ=255)


                                 aTCTGCATGTAGTCGGCGGGCAGTTCGAGCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGGATCT.G..GCGCAAACGG..AAAGTCT.g.tcgt                                                                     >  1:3550556/1‑101 (MQ=255)


                                 aTCTGCATGTAGTCGGCGGGCAGTTCGAGCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGGATCT.G..GCGCAAACGG..AAAGTCT.g.tcgt                                                                     <  2:591493/101‑1 (MQ=255)


                                  tCTGCATGTAGTCGGCGGGCAGTTCGAGCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGGATCT.G..GCGCAAACGG..AAAGTCT.g.gcgtc                                                                    >  1:919872/1‑101 (MQ=255)


                                     gCATGTAGTCGGCGGGCAGTTCGAGCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGGATCT.G..GCGCAAACGG..AAAGTCT.G.TCGTCGCt                                                                 >  2:1465308/1‑101 (MQ=255)


                                     gCATGTAGCCGGCGGGCAGTTCGAGCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGGATCT.G..GCGCAAACGG..AAAGTCT.G.TCGTCGCt                                                                 <  2:1876006/101‑1 (MQ=255)


                                      cATGTAGTCGGCGGGCAGTTCGAGCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGGATCT.G..GCGCAAACGG..AAAGTCT.G.TCGTCGCTc                                                                <  1:443112/101‑1 (MQ=255)


                                       aggTAGTCGGGGGGCAGTTCGAGCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGGATCT.G..GCGCAAACGG..AAAGTCT.G.TCGTCGCTcc                                                               <  1:979659/99‑1 (MQ=255)


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


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


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


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


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


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


                                                gcggGCAGTTCGAGCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGGATCT.G..GCGCAAACGG..AAAGTCT.G.TCGTCGCTCCCCCCGCggg                                                      >  1:3469124/1‑101 (MQ=255)


                                                  gggCAGTTCGAGCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGGATCT.G..GCGCAAACGG..AAAGTCT.G.TCGTCGCTCCCCCCGCggggg                                                    >  2:1402118/1‑101 (MQ=255)


                                                  gggCAGTTCGAGCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGGATCT.G..GCGCAAACGG..AAAGTCT.G.TCGTCGCTCCCCCCGCggggg                                                    <  1:4122783/101‑1 (MQ=255)


                                                      aGTTCGAGCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGGATCT.G..GCGCAAACGG..AAAGTCT.G.TCGTCGCTCCCCCCGCGGGGGCGTg                                                <  1:1015688/101‑1 (MQ=255)


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


                                                       ggtCGAGCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGGATCT.G..GCGCAAACGG..AAAGTCT.G.TCGTCGCTCCCCCCGCGGGGGCGTgg                                               <  2:3451297/99‑1 (MQ=255)


                                                        ttCGAGCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGGATCT.G..GCGCAAACGG..AAAGTCT.G.TCGTCGCTCCCCCCGCGGGGGCGTGGa                                              >  2:2350300/1‑101 (MQ=255)


                                                        ttCGAGCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGGATCT.G..GCGCAAACGG..AAAGTCT.G.TCGTCGCTCCCCCCGCGGGGGCGTGGa                                              >  1:1496272/1‑101 (MQ=255)


                                                         tcaagattcAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGGCAAGGACAT.C..GCCGATCTGGGTAAGGCCG.G.TGGTCGCTCCCCCCGCGGGGGCGTg                                                >  1:3764702/9‑101 (MQ=2)


                                                          tgcgcgctggTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.TG.AT.GACCA.GATCCGCAATCGG..CCGGGCTCGATCGTCGCTCCCCCCGCGGGGGCGTGGa                                              >  2:1027987/10‑101 (MQ=0)


                                                          cGAGCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGGATCT.G..GCGCAAACGG..AAAGTCT.G.TCGTCGCTCCCCCCGCGGGGGCGTGGATc                                            >  2:896153/1‑101 (MQ=255)


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


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


                                                            agattcAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGGCAAGGACAT.C..GCCGATCGGGGTAAGGCCG.G.TGGTCGCTCCCCCCGCGGGGGCGTGGAt                                             <  2:3673078/96‑1 (MQ=2)


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


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


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


                                                               ttcAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGGCAAGGACAT.C..GCCGATCTGGGTAAGGCCG.G.TGGTCGCTCCCCCCGCGGGGGCGTGGATCGa                                          <  1:2451083/99‑1 (MQ=2)


                                                               gctggTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.TG.AT.GACCA.GATCCGCAATCGG..CCGGGCTCGATCGTCGCTCCCCCCGCGGGGGCGTGGATCGaa                                         <  1:4030793/97‑1 (MQ=2)


                                                                tcAGTCCCGCCCCCCGCGGGGGCGTGGATCGAAACTCGGCAAGGACAT.C..GCCGATCGGGGTAAGGCCG.G.TGGTCGCTCCCCCCGCGGGGGCGTGGATCGaa                                         <  2:1266001/100‑1 (MQ=2)


                                                                gCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGGATCT.G..GCGCAAACGG..AAAGTCT.G.TCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAAcc                                      <  1:3056906/101‑1 (MQ=255)


                                                                 cAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGGCAAGGACAT.C..GCCGATCTGGGTAAGGCCG.G.TGGTCGCTCCCCCCGCGGGGGCGTGGATCGaaa                                        <  1:1148354/101‑1 (MQ=2)


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


                                                                  aGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGGATCT.G..GCGCAAACGG..AAAGTCT.G.TCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCta                                    >  2:2195220/1‑100 (MQ=255)


                                                                  aGTCGCTCCCACCGCGGGGGCGTGGATCGAAACT.GGTATGGATCT.G..GCGCAAACGG..AAAGTCT.G.TCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCTg                                    >  1:3659981/1‑101 (MQ=255)


                                                                   gTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGGCAAGGACAT.C..GCCGATCTGGGTAAGGCCG.G.TGGGCGCTCCCCCCGCGGGGGCGTGGATCGAAAca                                      >  1:1970189/1‑100 (MQ=2)


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


                                                                     cGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGGCAAGGACAT.C..GCCGATCTGGGTAAGGCCG.G.TGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAAcacc                                    >  2:1769359/1‑98 (MQ=2)


                                                                     cGCTCCCCCCGCGGGGGCGTGGATCGAAACT.TG.AT.GACCA.GATCCGCAATCGG..CCGGGCTCGATCGTCGCTCCCCCCGCGGGGGGGTGGATCGAAACCtat                                   >  1:3529207/1‑99 (MQ=2)


                                                                     cGCTCCCCCCGCGGGGGCGTGGATCGAAACT.TG.AT.GACCA.GATCCGCAATCGG..CCGGGCTCGATCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCtat                                   <  1:1116532/101‑3 (MQ=2)


                                                                      gCTCCCCCCGCGGGGGCGTGGATCGAAACTCGGCAAGGACATCG..CCGCTCTGGGT.AAGGCCG.G.TGGTCGCTCCCCCCGCGGGGGCGGGGATCGAGAcaccc                                   >  1:693912/1‑97 (MQ=2)


                                                                      gCTCCCCCCGCGGGGGCGTGGATCGAAACTCGGCAAGGACAT.C..GCCGATCTGGGTAAGGCCG.G.TGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAAcaccc                                   <  1:2735186/101‑5 (MQ=2)


                                                                      gCTCCCCCCGCGGGGGCGTGGATCGAAACTCGGCAAGGACAT.C..GCCGATCTGGGTAAGGCCG.G.TGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAAcaccc                                   <  1:566754/101‑5 (MQ=2)


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


                                                                      gCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGGATCT.G..GCGCAAACGG..AAAGTCT.G.TCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCCgcggg                                >  2:2319886/1‑100 (MQ=255)


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


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


                                                                       cTCCCCCCGCGGGGGCGTGGATCGAAACTCGGCAAGGACAT.C..GCCGATCTGGGTAAGGCCG.G.TGGTCGCTCCCCCCGCGGGGGCCTGGATCGAAAcacccg                                  <  1:3695981/101‑6 (MQ=2)


                                                                       cTCCCCCCGCGGGGGCGTGGATCGAAACT.TG.AT.GACCA.GATCCGCAATCGG..CCGGGCTCGATCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCtatga                                 <  2:21273/101‑5 (MQ=2)


                                                                        tCCCCCCGCGGGGGCGTGGATCGAAACTCGGCAAGGACAT.C..GCCGATCTGGGTAAGGCCG.G.TGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAAcacccgg                                 >  1:3393258/1‑95 (MQ=2)


                                                                        tCCCCCCGCGGGGGCGTGGATCGAAACT.TG.AT.GACCA.GATCCGCAATCGG..CCGGGCTCGATCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCtatgag                                >  1:1389987/1‑96 (MQ=2)


                                                                         ccccccGCGGGGGCGTGGATCGAAACTCGGCAAGGACAT.C..GCCGATCTGGGTAAGGCCG.G.TGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAAcacccggc                                <  1:3149155/101‑8 (MQ=0)


                                                                          cccccGCGGGGGCGTGGATCGAAACTCGGCAAGGACAT.C..GCCGATCTGGGTAAGGCCG.G.TGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAAcacccggcg                               <  1:4177721/101‑9 (MQ=2)


                                                                           ccccGCGGGGGCGTGGATCGAAACTCGGCAAGGACAT.C..GCCGATCTGGGTAAGGCCG.G.TGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAAcacccggcgc                              <  1:1621800/101‑10 (MQ=2)


                                                                            cccGCGGGGGCGTGGATCGAAACTCGGCAAGGACAT.C..GCCGATCTGGGTAAGGCCG.G.TGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAAcacccggcgcc                             <  1:4109695/101‑11 (MQ=2)


                                                                               gcgGGGGCGTGGATCGAAACT.GGTATGGATCT.G..GCGCAAACGG..AAAGTCT.G.TCGTCGCTCCCACCGCGGGGGCGTGGATCGAAACCTGCGGCGTGAGTTGc                       <  1:1909009/101‑1 (MQ=255)


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


                                                                                  ggggCGTGGATCGAAACT.GGTATGGATCT.G..GCGCAAACGG..AAAGTCT.G.TCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCTGCGGCGTGAGTTGCCAg                    <  1:1430911/101‑1 (MQ=255)


                                                                                    ggCGTGGATCGAAACT.GGTATGGATCT.G..GCGCAAACGG..AAAGTCT.G.TCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCTACGGCGTGAGTTGCCAGTg                  <  1:2195379/101‑1 (MQ=255)


                                                                                     gCGTGGATCGAAACT.GGTATGGATCT.G..GCGCAAACGG..AAAGTCT.G.TCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCTGCGGCGTGAGTTGCCAGTGa                 <  2:982545/101‑1 (MQ=255)


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


                                                                                       gTGGATCGAAACT.GGTATGGATCT.G..GCGCAAACGG..AAAGTCT.G.TCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCTGCGGCGTGAGTTGCCAGTGATa               <  1:1885997/101‑1 (MQ=255)


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


                                                                                          gATCGAAACT.GGTATGGATCT.G..GCGCAAACGG..AAAGTCT.G.TCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCTGCGGCGTGAGTTGCCAGTGATAGAc            >  2:806233/1‑101 (MQ=255)


                                                                                          gATCGAAACT.GGTATGGATCT.G..GCGCAAACGG..AAAGTCT.G.TCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCTGCGGCGTGAGTTGCCAGTGATAGAc            >  1:2993949/1‑101 (MQ=255)


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


                                                                                             cGAAACT.GGTATGGATCT.G..GCGCAAACGG..AAAGTCT.G.TCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCTGCGGCGTGAGTTGCCAGTGATAGACGcc         >  1:1645150/1‑101 (MQ=255)


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


                                                                                               aaaCT.GGTATGGATCT.G..GCGCAAACGG..AAAGTCT.G.TCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCTGCGGCGTGAGTTGCCAGTGATAGACGCCGt       <  1:1737144/101‑1 (MQ=255)


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


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


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


                                                                                                     ggTATGGATCT.G..GCGCAAACGG..AAAGTCT.G.TCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCTGCGGCGTGAGTTGCCAGTGATAGACGCCGTGGACg  <  1:1561851/101‑1 (MQ=255)


CGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCGGGCAGTTCGAGCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGGATCT.G..GCGCAAACGG..AAAGTCT.G.TCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCTGCGGCGTGAGTTGCCAGTGATAGACGCCGTGGACG  >  NC_014034/1395632‑1395832

Base quality scores:


ATCG

< 3 ≤


ATCG

< 18 ≤


ATCG

< 26 ≤


ATCG

< 33 ≤


ATCG

< 40 ≤


ATCG