BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	W3110S.gb	1,282,409	A→G	P323P (CCA→CCG)	narG →	nitrate reductase 1, alpha subunit

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	W3110S.gb	1,282,409	0	A	G	100.0%	221.5 / NA	64	P323P (CCA→CCG)	narG	nitrate reductase 1, alpha subunit
Reads supporting (aligned to +/- strand): ref base A (0/0); new base G (29/35); total (29/35)

GGCACCGATGCGGCAATGGCGCTGGCGATGGGCCACGTAATGCTGCGTGAATTCCACCTCGACAACCCAAGCCAGTATTTCACCGACTATGTGCGTCGCTACACCGACATGCCGATGCTGGT  >  W3110S.gb/1282341‑1282462

ggCACCGATGCGGCAATGGCGCTGGCGATGGGCCACGTAATGCTGCGTGAATTCCACCTCGACAACCCGAg                                                     >  1:3856177/1‑71 (MQ=255)

 gCACCGATGCGGCAATGGCGCTGGCGATGGGCCACGTAATGCTGCGTGAATTCCACCTCGACAACCCGAGc                                                    >  1:69257/1‑71 (MQ=255)

    ccGATGCGGCAATGGCGCTGGCGATGGGCCACGTAATGCTGCGTGAATTCCACCTCGACAACCCGAGCCAg                                                 >  1:691968/1‑71 (MQ=255)

      gATGCGGCAATGGCGCTGGCGATGGGCCACGTAATGCTGCGTGAATTCCACCTCGACAACCCGAGCCAGTa                                               >  1:1287575/1‑71 (MQ=255)

      gATGCGGCAATGGCGCTGGCGATGGGCCACGTAATGCTGCGTGAATTCCACCTCGACAACCCGAGCCAGTa                                               >  1:1373675/1‑71 (MQ=255)

       aTGCGGCAATGGCGCTGGCGATGGGCCACGTAATGCTGCGTGAATTCCACCTCGACAACCCGAGCCAGTAt                                              <  1:4458182/71‑1 (MQ=255)

        tGCGGCAATGGCGCTGGCGATGGGCCACGTAATGCTGCGTGAATTCCACCTCGACAACCCGAGCCAGTAtt                                             >  1:4006132/1‑71 (MQ=255)

        tGCGGCAATGGCGCTGGCGATGGGCCACGTAATGCTGCGTGAATTCCACCTCGACAACCCGAGCCAGTAtt                                             >  1:156013/1‑71 (MQ=255)

        tGCGGCAATGGCGCTGGCGATGGGCCACGTAATGCTGCGTGAATTCCACCTCGACAACCCGAGCCAGTAtt                                             >  1:2025657/1‑71 (MQ=255)

         gCGGCAATGGCGCTGGCGATGGGCCAGGTAATGCTGCGTGAATTCCACCTCGACAACCCGAGCCAGTAttt                                            >  1:3989143/1‑71 (MQ=255)

         gCGGCAATGGCGCTGGCGATGGGCCACGTAATGCTGCGTGAATTCCACCTCGACAACCCGAGCCAGTa                                               >  1:4406519/1‑68 (MQ=255)

         gCGGCAATGGCGCTGGCGATGGGCCACGTAATGCTGCGTGAATTCCACCTCGACAACCCGAGCCAGTAttt                                            >  1:4381139/1‑71 (MQ=255)

               aTGGCGCTGGCGATGGGCCACGTAATGCTGCGTGAATTCCACCTCGACAACCCGAGCCAGTa                                               <  1:3834750/62‑1 (MQ=255)

               aTGGCGCTGGCGATGGGCCACGTAATGCTGCGTGAATTCCACCTCGACAACCCGAGCCAGTATTTCAcc                                        <  1:2182750/69‑1 (MQ=255)

                tGGCGCTGGCGATGGGCCACGTAATGCTGCGTGAATTCCACCTCGACAACCCGAGCCAGTATTTCACCGAc                                     <  1:4582094/71‑1 (MQ=255)

                 ggCGCTGGCGATGGGCCACGTAATGCTGCGTGAATTCCACCTCGACAACCCGAGCCAGTATTTCACCGAc                                     >  1:3833289/1‑70 (MQ=255)

                 ggCGCTGGCGATGGGCCACGTAATGCTGCGTGAATTCCACCTCGACAACCCGAGCCAGTATTTCACCGAc                                     >  1:870248/1‑70 (MQ=255)

                 ggCGCTGGCGATGGGCCACGTAATGCTGCGTGAATTCCACCTCGACAACCCGAGCCAGTATTTCACCGAc                                     >  1:348486/1‑70 (MQ=255)

                 ggCGCTGGCGATGGGCCACGTAATGCTGCGTGAATTCCACCTCGACAACCCGAGCCAGTATTTCACCGAc                                     >  1:3385118/1‑70 (MQ=255)

                           aTGGGCCACGTAATGCTGCGTGAATTCCACCTCGACAACCCGAGCCAg                                                 <  1:3332106/48‑1 (MQ=255)

                                   cGTAATGCTGCGTGAATTCCACCTCGACAACCCGAGCCAGTa                                               <  1:1860224/42‑1 (MQ=255)

                                   cGTAATGCTGCGTGAATTCCACCTCGACAACCCGAGCCAGTATTTCACCGACTATGTGCGTCGCTACAcc                   >  1:903622/1‑70 (MQ=255)

                                   cGTAATGCTGCGTGAATTCCACCTCGACAACCCGAGCCAGTATTTCACCGACTATGTGCGTCGCTACAcc                   >  1:990681/1‑70 (MQ=255)

                                    gTAATGCTGCGTGAATTCCACCTCGACAACCCGAGCCAGTATTTCACCGACTATGTGCGTCGCTACACCGa                 <  1:4153990/71‑1 (MQ=255)

                                    gTAATGCTGCGTGAATTCCACCTCGACAACCCGAGCCAGTATTTCACCGACTATGTGCGTCGCTACACCGa                 <  1:428211/71‑1 (MQ=255)

                                    gTAATGCTGCGTGAATTCCACCTCGACAACCCGAGCCAGTATTTCACCGACTATGTGCGTCGCTACACCGa                 <  1:489215/71‑1 (MQ=255)

                                    gTAATGCTGCGTGAATTCCACCTCGACAACCCGAGCCAGTATTTCACCGACTATGTGCGTCGCTACACCGa                 <  1:1280048/71‑1 (MQ=255)

                                    gTAATGCTGCGTGAATTCCACCTCGACAACCCGAGCCAGTATTTCACCGACTATGTGCGTCGCTACACCGa                 <  1:1064035/71‑1 (MQ=255)

                                    gTAATGCTGCGTGAATTCCACCTCGACAACCCGAGCCAGTATTTCACCGACTATGTGCGTCGCTACACCGa                 <  1:1479906/71‑1 (MQ=255)

                                    gTAATGCTGCGTGAATTCCACCTCGACAACCCGAGCCAGTATTTCACCGACTATGTGCGTCGCTACACCGa                 <  1:2264217/71‑1 (MQ=255)

                                    gTAATGCTGCGTGAATTCCACCTCGACAACCCGAGCCAGTATTTCACCGACTATGTGCGTCGCTACACCGa                 <  1:1677429/71‑1 (MQ=255)

                                    gTAATGCTGCGTGAATTCCACCTCGACAACCCGAGCCAGTATTTCACCGACTATGTGCGTCGCTACACCGa                 <  1:2600263/71‑1 (MQ=255)

                                     tAATGCTGCGTGAATTCCACCTCGACAACCCGAGCCAGTa                                               <  1:1897490/40‑1 (MQ=255)

                                     tAATGCTGCGTGAATTCCACCTCGACAACCCGAGCCAGTa                                               <  1:4086286/40‑1 (MQ=255)

                                     tAATGCTGCGTGAATTCCACCTCGACAACCCGAGCCAGTATTTCACCGACTAtgtg                               >  1:4028090/1‑56 (MQ=255)

                                     tAATGCTGCGTGAATTCCACCTCGACAACCCGAGCCAGTATTTCACCGACTATGTGCGTCGCTACACCGa                 >  1:135928/1‑70 (MQ=255)

                                      aaTGCTGCGTGAATTCCACCTCGACAACCCGAGCCAGTATTTCACCGACTATGTGCGTCGCTACACCGACa               >  1:4586089/1‑71 (MQ=255)

                                      aaTGCTGCGTGAATTCCACCTCGACAACCCGAGACAGTAtt                                             >  1:461481/1‑41 (MQ=38)

                                       aTGCTGCGTGAATTCCACCTCGACAACCCGAGCCAGTATTTCACCg                                       <  1:517640/46‑1 (MQ=255)

                                       aTGCTGCGTGAATTCCACCTCGACAACCCGAGCCAGTATTTCACCGACTa                                   <  1:504011/50‑1 (MQ=255)

                                       aTGCTGCGTGAATTCCACCTCGACAACCCGAGCCAGTATTTCACCGACTAtg                                 <  1:3966003/52‑1 (MQ=255)

                                       aTGCTGCGTGAATTCCACCTCGACAACCCGAGCCAGTATTTCACCGACTAtg                                 >  1:2399166/1‑52 (MQ=255)

                                       aTGCTGCGTGAATTCCACCTCGACAACCCGAGCCAGTATTTCACCGACTATGTGCGTCGCTACACCGACAt              >  1:3485599/1‑71 (MQ=255)

                                        tgctgcGTGAATTCCACCTCGACAACCCGAGCCAGTATTTCACCGACTATGTGCGTCGCTACACCGAc                <  1:1713004/68‑1 (MQ=255)

                                        tgctgcGTGAATTCCACCTCGACAACCCGAGCCAGTATTTCACCGACTATGTGCGTCGCTACACCGAc                <  1:1813936/68‑1 (MQ=255)

                                         gctgcGTGAATTCCACCTCGACAACCCGAGCCAGTATTTCACCGAc                                     >  1:2401528/1‑46 (MQ=255)

                                         gctgcGTGAATTCCACCTCGACAACCCGAGCCAGTATTTCACCGAc                                     >  1:2744423/1‑46 (MQ=255)

                                         gctgcGTGAATTCCACCTCGACAACCCGAGCCAGTATTTCACCGAc                                     >  1:1833812/1‑46 (MQ=255)

                                          ctgcGTGAATTCCACCTCGACAACCCGAGCCAGTATTTCACCGACTATGTGCGTCGCTa                       >  1:3914132/1‑59 (MQ=255)

                                          ctgcGTGAATTCCACCTCGACAACCCGAGCCAGTATTTCACCGACTATGTGCGTCGCTa                       >  1:2669144/1‑59 (MQ=255)

                                                 aaTTCCACCTCGACAACCCGAGCCAGTATTTCACCGACTATGTGCGTCGCTACACCGACa               >  1:3532425/1‑60 (MQ=255)

                                                 aaTTCCACCTCGACAACCCGAGCCAGTATTTCACCGACTATGTGCGTCGCTACACCGACATGCCGATGCTg    <  1:3304739/71‑1 (MQ=255)

                                                   ttCCACCTCGACAACCCGAGCCAGTATTTCACCGACTATGTGCGTc                           <  1:4676478/46‑1 (MQ=255)

                                                   ttCCACCTCGACAACCCGAGCCAGTATTTCACCGACTATGTGCGTc                           <  1:4472188/46‑1 (MQ=255)

                                                   ttCCACCTCGACAACCCGAGCCAGTATTTCACCGACTATGTGCGTc                           <  1:208695/46‑1 (MQ=255)

                                                   ttCCACCTCGACAACCCGAGCCAGTATTTCACCGACTATGTGCGTc                           <  1:170412/46‑1 (MQ=255)

                                                   ttCCACCTCGACAACCCGAGCCAGTATTTCACCGACTATGTGCGTc                           <  1:3331356/46‑1 (MQ=255)

                                                   ttCCACCTCGACAACCCGAGCCAGTATTTCACCGACTATGTGCGTCGCTACACCGa                 <  1:901447/56‑1 (MQ=255)

                                                   ttCCACCTCGACAACCCGAGCCAGTATTTCACCGACTATGTGCGTCGCTACACCGACATGCCGATGCTGGt  <  1:437874/71‑1 (MQ=255)

                                                   ttCCACCTCGACAACCCGAGCCAGTATTTCACCGACTATGTGCGTCGCTACACCGACATGCCGATGCTGGt  <  1:2739120/71‑1 (MQ=255)

                                                   ttCCACCTCGACAACCCGAGCCAGTATTTCACCGACTATGTGCGTCGCTACACCGACATGCCGATGCTGGt  <  1:2846643/71‑1 (MQ=255)

                                                   ttCCACCTCGACAACCCGAGCCAGTATTTCACCGACTATGTGCGTCGCTACACCGACATGCCGATGCTGGt  <  1:3074308/71‑1 (MQ=255)

                                                    tCCACCTCGACAACCCGAGCCAGTATTTCACCGACTATGTGCGTCGCTACACCGACATGCCGAt        <  1:676999/64‑1 (MQ=255)

                                                            gACAACCCGAGCCAGTATTTCACCGACTATGTGCGTCGCTACACCGACATGcc           <  1:496545/53‑1 (MQ=255)

GGCACCGATGCGGCAATGGCGCTGGCGATGGGCCACGTAATGCTGCGTGAATTCCACCTCGACAACCCAAGCCAGTATTTCACCGACTATGTGCGTCGCTACACCGACATGCCGATGCTGGT  >  W3110S.gb/1282341‑1282462

Alignment Legend
Aligned base mismatch/match (shaded by quality score): `ATCG/ATCG` `< 3 ≤` `ATCG/ATCG` `< 14 ≤` `ATCG/ATCG` `< 32 ≤` `ATCG/ATCG` `< 36 ≤` `ATCG/ATCG`
Unaligned base: `atcg` Masked matching base: `atcg` Alignment gap: `‑` Deleted base: `‑`