BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,765,234	G→A	G129G (GGG→GGA)	yibJ →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,765,234	0	G	A	100.0%	122.2 / NA	42	G129G (GGG→GGA)	yibJ	hypothetical protein
Reads supporting (aligned to +/- strand): ref base G (0/0); new base A (18/24); total (18/24)

GCAGGGGCGGTATATCACTCAGGATCCGATTGGGCTGAAGGGGGGATGGAATTTTTATCAGTATCCGTTGAATCCGGTCATAAATGTAGATCCGCAAGGTTTGGT  >  NC_000913/3765201‑3765305

gctggggCGGTATATCACTCAGGATCCGATTGGACTGAAGGGGGGATGGAATTTTTATCAGTATCCGTTGAATCCg                               <  2:2587249/73‑1 (MQ=255)

gctggggCGGTATATCACTCAGGATCCGATTGGACTGAAGGGGGGATGGAATTTTTATCAGTATCCGTTGAATCCg                               >  1:5537101/4‑76 (MQ=255)

 ctggggCGGTATATCACTCAGGATCCGATTGGACTGGAGGGGGGGTGGGATTTTTATCAGTATCCGTTGAATCCgg                              >  1:3365677/3‑76 (MQ=11)

 ctggggCGGTATATCACTCAGGATCCGATTGGACTGAAGGGGGGATGGAATTTTTATCAGTATCCGTTGAATCCgg                              <  2:2499961/74‑1 (MQ=255)

  tggggCGGTATATCACTCAGGATCCGATTGGACTGAAGGGGGGATGGAATTTTTATCAGTATCCGTTGAATCCGGt                             >  2:4188765/2‑76 (MQ=255)

  tggggCGGTATATCACTCAGGATCCGATTGGACTGAAGGGGGGATGGAATTTTTATCAGTATCCGTTGAATCCGGt                             >  1:867635/2‑76 (MQ=255)

  tggggCGGTATATCACTCAGGATCCGATTGGACTGAAGGGGGGATGGAATTTTTATCAGTATCCGTTGAATCCGGt                             >  1:5708488/2‑76 (MQ=255)

    gggCGGTATATCACTCAGGATCCGATTGGACTGAAGGGGGGATGGAATTTTTATCAGTATCCGTTGAATCCGGTCa                           >  2:3662146/1‑76 (MQ=255)

      gCGGTATATCACTCAGGATCCGATTGGACTGAAGGGGGGATGGAATTTTTATCAGTATCCGTTGAATCCGGTCATa                         >  2:1320019/1‑76 (MQ=255)

       cGGTATATCACTCAGGATCCGATTGGACTGAAGGGGGGATGGAATTTTTATCAGTCTCCGTTGAATCCGGTCATaa                        <  1:4735134/76‑1 (MQ=14)

       cGGTATATCACTCAGGATCCGATTGGACTGAAGGGGGGATGGAATTTTTATCAGTATCCGTTGAATCCGGTCATaa                        >  2:444744/1‑76 (MQ=255)

        ggTATATCACTCAGGATCCGATTGGACTGAAGGGGGGATGGAATTTTTATCAGTATCCGTTGAATCCGGTCATaaa                       >  1:3719085/1‑76 (MQ=255)

         gtttatCATGCAGGATCCGATTGGACTGACGGGGGGATGGAATTTTTATCAGTATCCGTTGAATCCGGTCATAAAt                      <  2:3863352/73‑1 (MQ=14)

         gTATATCACTCAGGATCCGATTGGACTGAAGGGGGGATGGAATTTTTATCAGTATCCGTTGAATCCGGTCATAAAt                      <  2:830861/76‑1 (MQ=255)

         gTATATCACTCAGGATCCGATTGGACTGAAGGGGGGATGGAATTTTTATCAGTATCCGTTGAATCCGGTCATAAAt                      >  1:476782/1‑76 (MQ=255)

         gTATATCACTCAGGATCCGATTGGACTGAAGGGGCGATGGAATTTTTATCAGTATCCGTTGAATCCGGTCATAAAt                      >  1:4589135/1‑76 (MQ=14)

         gTATAACACTCAGGATCCGATTGGACTGAAGGGGGGATGGAATTTTTATCAGTATCCGTTGAATCCGGTCATAAAt                      <  2:969126/76‑1 (MQ=255)

          tataTCACTCAGGATCCGATTGGACTGAAGGGGGGATGGAATTTTTATCAGTATCCGTTGAATCCGGTCATAAATg                     >  2:3126223/1‑76 (MQ=255)

          tataTCACTCAGGATCCGATTGGACTGAAGGGGGGATGGAATTTTTATCAGTATCCGTTGAATCCGGTCATAAATg                     <  1:4342179/76‑1 (MQ=255)

          tataTCACTCAGGATCCGATTGGACTGAAGGGGGGATGGAATTTTTATCAGTATCCGTTGAATCCGGTCATAAATg                     >  2:1873170/1‑76 (MQ=255)

           atatCACTCAGGATCCGATTGGACTGAAGGGGGGATGGAATTTTTATCAGTATCCGTTGAATCCGGTCATAAATGt                    <  2:2297284/76‑1 (MQ=255)

              tctcTCAGGATCCGATTGGACTGCAGGGGGGATGGAATTTTTATCAGTATCCGTTGAATCCGGTCATAAATGTAGa                 <  1:629672/73‑1 (MQ=17)

              tCACTCAGGATCCGATTGGACTGAAGGGGGGATGGAATTTTTATCAGTATCCGTTGAATCCGGTCATAAATGTAGa                 >  1:979030/1‑76 (MQ=255)

              tCACTCAGGATCCGATTGGACTGAAGGGGGGATGGAATTTTTATCAGTATCCGTTGAATCCGGTCATAAATGTAGa                 <  2:2305659/76‑1 (MQ=255)

               cACTCAGGATCCGATTGGACTGAAGGGGGGATGGAATTTTTATCAGTATCCGTTGAATCCGGTCATAAATGTAGAt                <  2:5412638/76‑1 (MQ=255)

                aCTCAGGATCCGATTGGACTGAAGGGGGGATGGAATTTTTATCAGTATCCGTTGAATCCGGTCATAAATGTAGATc               <  2:1385932/76‑1 (MQ=255)

                aCTCAGGATCCGATTGGACTGAAGGGGGGATGGAATTTTTATCAGTATCCGTTGAATCCGGTCATAAATGTAGATc               <  1:4065059/76‑1 (MQ=255)

                 cTCAGGATCCGATTGGACTGAAGGGGGGATGGAATTTTTATCAGTATCCGTTGAATCCGGTCATAAATGTAGATcc              <  2:2147451/76‑1 (MQ=255)

                   cAGGATCCGATTGGACTGAAGGGGGGATGGAATTTTTATCAGTATCCGTTGAATCCGGTCATAAATGTAGATCCGc            <  2:979030/76‑1 (MQ=255)

                    aGGATCCGATTGGACTGAAGGGGGGATGGAATTTTTATCAGTATCCGTTGAATCCGGTCATAAATGTAGATCCGCa           >  2:1651287/1‑76 (MQ=255)

                     ggATCCGATTGGACTGAAGGGGGGATGGAATTTTTATCAGTATACGTTGAATCCGGTCATAAATGTAGATCCGCaa          >  1:2614220/1‑76 (MQ=255)

                      gATCCGATTGGACTGAAGGGGGGATGGAATTTTTATCAGTATCCGTTGAATCCGGTCATAAATGTAGATCCGCAAg         <  2:3792459/76‑1 (MQ=255)

                      gATCCGATTGGACTGAAGGGGGGATGGAATTTTTATCAGTATCCGTTGAATCCGGTCATAAATGTAGATCCGCAAg         <  1:3819729/76‑1 (MQ=255)

                        tCCGATTGGACTGAAGGGGGGATGGAATTTTTATCAGTATCCGTTGAATCCGGTCATAAATGTAGATCCGCAAGGt       >  1:5712886/1‑76 (MQ=255)

                        tCCGATTGGACTGAAGGGGGGATGGAATTTTTATCAGTATCCGTTGAATCCGGTCATAAATGTAGATCCGCAAGGt       <  1:435598/76‑1 (MQ=255)

                          cGATTGGACTGAAGGGGGGATGGAATTTTTATCAGTATCCGTTGAATCCGGTCATAAATGTAGATCCGCAAGGttt     <  2:1249261/76‑1 (MQ=255)

                          cGATTGGACTGAAGGGGGGATGGAATTTTTATCAGTATCCGTTGAATCCGGTCATAAATGTAGATCCGCAAGGttt     <  2:3719085/76‑1 (MQ=255)

                           gATTGGACTGAAGGGGGGATGGAATTTTTATCAGTATCCGTTGAATCCGGTCATAAATGTAGATCCGCAAGGTTTg    >  2:1563639/1‑76 (MQ=255)

                           gATTGGACTGAAGGGGGGATGGAATTTTTATCAGTATCCGTTGAATCCGGTCATAAATGTAGATCCGCAAGGTTTg    <  1:3583561/76‑1 (MQ=255)

                            aTTGGACTGAAGGGGGGATGGAATTTTTATCAGTATCCGTTGAATCCGGTCATAAATGTAGATCCGCAAGGTTTgg   <  2:3067474/76‑1 (MQ=255)

                             ttGGACTGAAGGGGGGATGGAATTTTTATCAGTATCCGTTGAATCCGGTCATAAATGTAGATCCGCAAGGTTTGGt  <  1:4545500/76‑1 (MQ=255)

                             ttGGACTGAAGGGGGGATGGAATTTTTATCAGTATCCGTTGAATCCGGTCATAAATGTAGATCCGCAAGGTTTGGt  <  1:4314990/76‑1 (MQ=255)

GCAGGGGCGGTATATCACTCAGGATCCGATTGGGCTGAAGGGGGGATGGAATTTTTATCAGTATCCGTTGAATCCGGTCATAAATGTAGATCCGCAAGGTTTGGT  >  NC_000913/3765201‑3765305

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