BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,129,541	C→T	pseudogene (603/1101 nt)	yghO ←	pseudogene, IS‑interrupted

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,129,541	0	C	T	100.0%	135.7 / NA	43	pseudogene (603/1101 nt)	yghO	pseudogene, IS‑interrupted
Reads supporting (aligned to +/- strand): ref base C (0/0); new base T (22/21); total (22/21)

CGTCCGGCACCAGATATTTAAGTTGATCGCCCATGGTCGCGAATTCATGTTCGGTAAACGGCACAAATCCCCAGTTGTGCTGCCAGCCAGAGTTGAAAATCTCACGCAGGATCTGCATCTCTTCGGCAAACCGCTGACGATTGATGCAGCGAATGGTCACCTTTTTGCGCACCTGATCCATCAGTTTTTTTAGCGCCGGAGAGAAAGTGAGATCGGTTCG  >  NC_000913/3129426‑3129645

cGTCCGGCACCAGATATTTAAGTTGATCGCCCATGGTCGCGAATTCATGTTCGGTAAACGGCACAAATCCCCAGTTGTGCTGCCAGCCAGAGTTGAAAATCTCACGCAGGATCTGTATCTCTTCGGCAAACCGCTGACg                                                                                   >  2:893599/1‑139 (MQ=255)

  tCCGGCACCAGATATTTAAGTTGATCGCCCATGGTCGCGAATTCATGTTCGGTAAACGGCACAAATCCCCAGTTGTGCTGCCAGCCAGAGTTGAAAATCTCACGCAGGATCTGTATCTCTTCGGCAAACCGCTGACGAt                                                                                 >  1:1049159/1‑139 (MQ=255)

            gATATTTAAGTTGATCGCCCATGGTCGCGAATTCATGTTCGGTAAACGGCACAAATCCCCAGTTGTGCTGCCAGCCAGAGTTGAAAATCTCACGCAGGATCTGTATCTCTTCGGCAAACCGCTGACGATTGATGCAGCg                                                                       >  2:220792/1‑139 (MQ=255)

             cagtttAAGTTGAGCGCCCATGGTCGCGAATTCATGTTCGGTAAACGGCAGAAATCCCCAGTTGTTCTGCCAGCCAGAGTTGAAGATCTCACGCAGGATCTGTCTCTCATCGGCAAACCGCTGACGATTGATGCAGCGa                                                                      <  1:796895/136‑1 (MQ=255)

                 ttAAGTTGATCGCCCATGGTCGCGAATTCATGTTCGGTAAACGGCACAAATCCCCAGTTGTGCTGCCAGCCAGAGTTGAAAATCTCACGCAGGATCTGTATCTCTTCGGCAAACCGCTGACGATTGATGCAGCGAATgg                                                                  >  1:416026/1‑139 (MQ=255)

                       tGATCGCCCATGGTCGCGAATTCATGTTCGGTAAACGGCACAAATCCCCAGTTGTGCTGCCAGCCAGAGTTGAAAATCTCACGCAGGATCTGTATCTCTTCGGCAAACCGCTGACGAt                                                                                 >  1:821609/1‑118 (MQ=255)

                       tGATCGCCCATGGTCGCGAATTCATGTTCGGTAAACGGCACAAATCCCCAGTTGTGCTGCCAGCCAGAGTTGAAAATCTCACGCAGGATCTGTATCTCTTCGGCAAACCGCTGACGAt                                                                                 <  2:821609/118‑1 (MQ=255)

                       tGATCGCCCATGGTCGCGAATTCATGTTCGGTAAACGGCACAAATCCCCAGTTGTGCTGCCAGCCAGAGTTGAAAATCTCACGCAGGATCTGTATCTCTTCGGCAAACCGCTGACGATTGATGCAGCGAATGGTCACCt                                                            >  1:109469/1‑139 (MQ=255)

                       tGATCGCCCATGGTCGCGAATTCATGTTCGGTAAACGGCACAAATCCCCAGTTGTGCTGCCAGCCAGAGTTGAAAATCTCACGCAGGATCTGTATCTCTTCGGCAAACCGCTGACGATTGATGCAGCGAATGGTCACCt                                                            >  1:246493/1‑139 (MQ=255)

                        gATCGCCCATGGTCGCGAATTCATGTTCGGTAAACGGCACAAATCCCCAGTTGTGCTGCCAGCCAGAGTTGAAAATCTCACGCAGGATCTGTATCTCTTCGGCAAACCGCTg                                                                                      <  2:848569/112‑1 (MQ=255)

                        gATCGCCCATGGTCGCGAATTCATGTTCGGTAAACGGCACAAATCCCCAGTTGTGCTGCCAGCCAGAGTTGAAAATCTCACGCAGGATCTGTATCTCTTCGGCAAACCGCTg                                                                                      >  1:848569/1‑112 (MQ=255)

                           cGCCCCTGGTCGCGAATTCATGTTCGGTAAACGGCACAAATCCCCAGTTGTGCTCCCAGCCAGAGTTGAAAACCTCACGCAGGATCTGTATCTCTTCGGCAAACCGCTGACGATTCATGCAGCGAATTGTCACCttctg                                                        >  2:796895/1‑136 (MQ=255)

                           cGCCCATGGTCGCGAATTCATGTTCGGTAAACGGCACAAATCCCCAGTTGTGCTGCCAGCCAGAGTTGAAAATCTCACGCAGGATCTGTATCTCTTCGGCAAACCGCTGACGATTGATGCAGCGAATGGTCACCttttt                                                        >  2:403577/1‑139 (MQ=255)

                            gCCCATGGTCGCGAATTCATGTTCGGTAAACGGCACAAATCCCCAGTTGTGCTGCCAGCCAGAGTTGAAAATCTCACGCAGGATCTGTATCTCTTCGGCAAACCGCTGACGATTGATGCAGCGAATGGTCACCTTTTTg                                                       <  1:220792/139‑1 (MQ=255)

                            gCCCATGGTCGCGAATTCATGTTCGGTAAACGGCACAAATCCCCAGTTGTGCTGCCAGCCAGAGTTGAAAATCTCACGCAGGATCTGTATCTCTTCGGCAAACCGCTGACGATTGATGCAGCGAATGGTCACCTTTTTg                                                       <  1:893599/139‑1 (MQ=255)

                              ccATGGTCGCGAATTCATGTTCGGTAAACGGCACAAATCCCCAGTTGTGCTGCCAGCCAGAGTTGAAAATCTCACGCAGGATCTGTATCTCTTCGGCAAACCGCTGACGATTGATGCAGCGAATGGTCACCTTTTTgcg                                                     <  2:246493/139‑1 (MQ=255)

                                   gTCGCGAATTCATGTTCGGTAAACGGCACAAATCCCCAGTTGTGCTGCCAGCCAGAGTTGAAGATCTCACGCAGGATCTGTATCTCTTCGGCAAACCGCTGACGATTGATGCAGCGAATGGTCACCTTTTTGCGCACCt                                                <  2:439841/139‑1 (MQ=255)

                                    tCGCGAATTCATGTTCGGTAAACGGCACAAATCCCCAGTTGTGCTGCCAGCCAGAGTTGAAAATCTCACGCAGGATCTGTATCTCTTCGGCAAACCGCTGACGATTGATGCAGCGAATGGTCACCTTTTTGCGCACCTg                                               >  2:143960/1‑139 (MQ=255)

                                        gAATTCATGTTCGGTAAACGGCACAAATCCCCAGTTGTGCTGCCAGCCAGAGTTGAAAATCTCACGCAGGATCTGTATCTCTTCGGCAAACCGCTGACGATTGATGCAGCGAATGGTCACCTTTTTGCGCAcc                                                 <  2:129630/133‑1 (MQ=255)

                                        gAATTCATGTTCGGTAAACGGCACAAATCCCCAGTTGTGCTGCCAGCCAGAGTTGAAAATCTCACGCAGGATCTGTATCTCTTCGGCAAACCGCTGACGATTGATGCAGCGAATGGTCACCTTTTTGCGCAcc                                                 >  1:129630/1‑133 (MQ=255)

                                                    ggTAAGCGGCACAAATCCCCAGTTGTGCTGCCAGCCAGAGTTGAAAATCTCACGCAGGATCTGTATCTCTTCGGCAAACCGCTGACGATTGATGc                                                                           <  1:609483/95‑1 (MQ=255)

                                                    ggTAAACGGCACAAATCCCCAGTTGTGCTGCCAGCCAGAGTTGAAAATCTCACGCAGGATCTGTATCTCTTCGGCAAACCGCTGACGATTGATGc                                                                           >  2:609483/1‑95 (MQ=255)

                                                      tAAACGGCACAAATCCCCAGTTGTGCTGCCAGCCAGAGTTGAAAATCTCACGCAGGATCTGTATCTCTTCGGCAAACCGCTGACGATTGATGc                                                                           <  2:354848/93‑1 (MQ=255)

                                                      tAAACGGCACAAATCCCCAGTTGTGCTGCCAGCCAGAGTTGAAAATCTCACGCAGGATCTGTATCTCTTCGGCAAACCGCTGACGATTGATGc                                                                           >  1:354848/1‑93 (MQ=255)

                                                        aaCGGCACAAATCCCCAGTTGTGCTGCCAGCCAGAGTTGAAAATCTCACGCAGGATCTGTATCTCTTCGGCAAACCGCTGACGATTGATGCAGCGAATGGTCACCTTTTTgc                                                      >  2:916967/1‑112 (MQ=255)

                                                        aaCGGCACAAATCCCCAGTTGTGCTGCCAGCCAGAGTTGAAAATCTCACGCAGGATCTGTATCTCTTCGGCAAACCGCTGACGATTGATGCAGCGAATGGTCACCTTTTTgc                                                      <  1:916967/112‑1 (MQ=255)

                                                        aaCGGCACAAATCCCCAGTTGTGCTGCCAGCCAGAGTTGAAAATCTCACGCAGGATCTGTATCTCTTCGGCAAACCGCTGACGATTGATGCAGCGAATGGTCACCTTTTTGCGCACCTGATCCATCAGTTTTTTTAgcg                           <  2:1049159/139‑1 (MQ=255)

                                                           ggCACAAATCCCCAGTTGTGCTGCCAGCCAGAGTTGAAAATCTCACGCAGGATCTGTATCTCTTCGGCAAACCGCTGACGATTGATGCAGCGAATGGTCACCTTTTTGCGCACCTGATCCATCAGTTTTTTTAGCGCCg                        >  2:866723/1‑139 (MQ=255)

                                                             cacaAATCCCCAGTTGTGCTGCCAGCCAGAGTTGAAAATCTCACGCAGGATCTGTAt                                                                                                        <  1:703044/57‑1 (MQ=255)

                                                             cacaAATCCCCAGTTGTGCTGCCAGCCAGAGTTGAAAATCTCACGCAGGATCTGTATCTCTTCGGCAAACCGCTGACGATTGATGCAGCGAATGGTCACCTTTTTGCGCACCTGATCCATCAGTTTTTTTAGCGCCGga                      >  1:57882/1‑139 (MQ=255)

                                                               caAATCCCCAGTTGTGCTGCCAGCCAGAGTTGAAAATCTCACGCAGGATCTGTATCTCTTCGGCAAACCGCTGACGATTGATGCAGCGAATGGTCACCTTTTTGCGCACCTGATCCATCAGTTTTTTTAGCGCCGgaga                    <  1:403577/139‑1 (MQ=255)

                                                                aaaTCCCCAGTTGTGCTGCCAGCCAGAGTTGAAAATCTCACGCAGGATCTGTATCTCTTCGGCAAACCGCTg                                                                                      >  1:737194/1‑72 (MQ=255)

                                                                aaaTCCCCAGTTGTGCTGCCAGCCAGAGTTGAAAATCTCACGCAGGATCTGTATCTCTTCGGCAAACCGCTg                                                                                      <  2:737194/72‑1 (MQ=255)

                                                                aaaTCCCCAGTTGTGCTGCCAGCCAGAGTTGAAAATCTCACGCAGGATCTGTATCTCTTCGGCAAACCGCTGACGATTGATGCAGCGAATGGTCACCTTTTTGCGCACCTGATCCATCAGTTTTTTTAGCGCCGgagag                   <  2:416026/139‑1 (MQ=255)

                                                                aaaTCCCCAGTTGTGCTGCCAGCCAGAGTTGAAAATCTCACGCAGGATCTGTATCTCTTCGGCAAACCGCTGACGATTGATGCAGCGAATGGTCACCTTTTTGCGCACCTGATCCATCAGTTTTTTTAGCGCCGgagag                   >  2:178247/1‑139 (MQ=255)

                                                                     cccAGTTGTGCTGCCAGCCAGAGTTGAAAATCTCACGCAGGATCTGTATCTCTTCGGCAAACCGCTGACGATTGATGCAGCGAATGGTCACCTTTTTGCGCACCTGATCCATCAGTTTTTTTAGCGCCGGAGAGAAAGt              <  1:770502/139‑1 (MQ=255)

                                                                        aGTTGTGCTGCCAGCCAGAGTTGAAAATCTCACGCAGGATCTGTATCTCTTCGGCAAACCGCTGACGATTGATGCAGCGAATGGTCACCTTTTTGCGCACCTGATCCATCAGTTTTTTTAGCGCCGGAGAGAAAGTgag           <  1:178247/139‑1 (MQ=255)

                                                                             tgctgcCAGCCAGAGTTGAAAATCTCACGCAGGATCTGTATCTCTTCGGCAAACCGCTGACGATTGATGCAGCGAATGGTCACCTTTTTGCGCACCTGATCCATCAGTTTTTTTAGCGCCGGAGAGAAAGTGAGATCgg      >  1:895607/1‑139 (MQ=255)

                                                                             tgctgcCAGCCAGAGTTGAAAATCTCACGCAGGATCTGTATCTCTTCGGCAAACCGCTGACGATTGATGCAGCGAATGGTCACCTTTTTGCGCACCTGATCCATCAGTTTTTTTAGCGCCGGAGAGAAAGTGAGATCgg      >  2:672426/1‑139 (MQ=255)

                                                                              gctgcCAGCCAGAGTTGAAAATCTCACGCAGGATCTGTATCTCTTCGGCAAACCGCTGACGATTGATGCAGCGAATGGTCACCTTTTTGCGCACCTGATCCATCAGTTTTTTTAGCGCCGGAGAGAAAGTGAGATCGGt     <  1:131414/139‑1 (MQ=255)

                                                                               ctgcCAGCCAGAGTTGAAAATCTCACGCAGGATCTGTATCTCTTCGGCAAACCGCTGACGATTGATGCAGCGAATGGTCACCTTTTTGCGCACCTGATCCATCAGTTTTTTTAGCGCCGGAGAGAAAGTGAGATCGGtt    >  1:561708/1‑139 (MQ=255)

                                                                                 gccagccaGAGTTGAAAATCTCACGCAGGATCTGTATCTCTTCGGCAAACCGCTGACGATTGATGCAGCGAATGGTCACCTTTTTGCGCACCTGATCCATCAGTTTTTTTAGCGCCGGAGAGAAAGTGAGATCGGTTCg  <  1:866723/139‑1 (MQ=255)

                                                                                 gccagccaGAGTTGAAAATCTCACGCAGGATCTGTATCTCTTCGGCAAACCGCTGACGATTGATGCAGCGAATGGTCACCTTTTTGCGCACCTGATCCATCAGTTTTTTTAGCGCCGGAGAGAAAGTGAGATCGGTTCg  <  2:109469/139‑1 (MQ=255)

CGTCCGGCACCAGATATTTAAGTTGATCGCCCATGGTCGCGAATTCATGTTCGGTAAACGGCACAAATCCCCAGTTGTGCTGCCAGCCAGAGTTGAAAATCTCACGCAGGATCTGCATCTCTTCGGCAAACCGCTGACGATTGATGCAGCGAATGGTCACCTTTTTGCGCACCTGATCCATCAGTTTTTTTAGCGCCGGAGAGAAAGTGAGATCGGTTCG  >  NC_000913/3129426‑3129645

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