BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,441,591	A→G	Y196H (TAT→CAT)	msrA ←	methionine sulfoxide reductase A

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,441,591	0	A	G	100.0%	125.9 / NA	36	Y196H (TAT→CAT)	msrA	methionine sulfoxide reductase A
Reads supporting (aligned to +/- strand): ref base A (0/0); new base G (18/18); total (18/18)

CCTGGTAAGTCGTTTACGCCGCATCAGGCATAGACAAAGAGCACTTTATCAACAAACTACATTTGTACCCGTAACGCTATGCTTCCGGCGGCAGACAGACGCCAATTCCGCCAATTCCACAGTAACCATACGGGTTTTTATGCAGATATTGCTGGTGGTCATCTTCGGCATAATAAAACGGTGTGGCGTTAGCGATTTCCGTGGTGATGTGACGATCGTCATCGGCGGCAAGCATCGCCGCCTGAAAACGTTCCAGACTGGCGC  >  NC_000913/4441462‑4441725

ccTGGTAAGTCGTTTACGCCGCATCAGGCACAGACAAAGAGCACTTTATCAACAAACTACATTTGTACCCGTAACGCTATGCTTCCGGCGGCAGACAGACGCCAATTCCGCCAATTCCACAGTAACCATGCGGGTTTTTa                                                                                                                              >  2:268358/1‑140 (MQ=255)

 cTGGTAAGTCGTTTACGCCGCATCAGGCATAGACAAAGAGCACTTTATCAACAAACTACATTTGTACCCGTAACGCTATGCTTCCGGCGGCAGACAGACGCCAATTCCGCCAATTCCACAGTAACCATGCGGGTTTTTAt                                                                                                                             <  2:558616/140‑1 (MQ=255)

 cTGGTAAGTCGTTTACGCCGCATCAGGCATAGACAAAGAGCACTTTATCAACAAACTACATTTGTACCCGTAACGCTATGCTTCCGGCGGCAGACAGACGCCAATTCCGCCAATTCCACAGTAACCATGCGGGTTTTTAt                                                                                                                             <  2:444235/140‑1 (MQ=255)

   ggTAAGTCGTTTACGCCGCATCAGGCATAGACAAAGAGCACTTTATCAACAAACTACATTTGTACCCGTAACGCTATGCTTCCGGCGGCAGACAGACGCCAATTCCGCCAATTCCACAGTAACCATGCGGGTTTTTATGc                                                                                                                           >  1:347019/1‑140 (MQ=255)

   ggTAAGTCGTTTACGCCGCATCAGGCATAGACAAAGAGCACTTTATCAACAAACTACATTTGTACCCGTAACGCTATGCTTCCGGCGGCAGACAGACGCCAATTCCGCCAATTCCACAGTAACCATGCGGGTTTTTATGc                                                                                                                           <  2:354964/140‑1 (MQ=255)

    gTAAGTCGTTTACGCCGCATCAGGCATAGACAAAGAGCACTTTATCAACAAACTACATTTGTACCCGTAACGCTATGCTTCCGGCGGCAGACAGACGCCAATTCCGCCAATTCCACAGTAACCATGCGGGTTTTTATGCa                                                                                                                          <  2:521825/140‑1 (MQ=255)

              tACGCCGCATCAGGCATAGACAAAGAGCACTTTATCAACAAACTACATTTGTACCCGTAACTCTATGCTTCCGGCGGCAGACAGACGCCAATTCCGCCAATTCCACAGTAACCATGCGGGTTTTTATGCAGATATTGCtg                                                                                                                >  2:169087/1‑140 (MQ=255)

              tACGCCGCATCAGGCATAGACAAAGAGCACTTTATCAACAAACTACATTTGTACCCGTAACGCTATGCTTCCGGCGGCAGACAGACGCCAATTCCGCCAATTCCACAGTAACCATGCGGGTTTTTATGCAGATATTGc                                                                                                                  <  1:344211/138‑1 (MQ=255)

              tACGCCGCATCAGGCATAGACAAAGAGCACTTTATCAACAAACTACATTTGTACCCGTAACGCTATGCTTCCGGCGGCAGACAGACGCCAATTCCGCCAATTCCACAGTAACCATGCGGGTTTTTATGCAGATATTGc                                                                                                                  >  2:344211/1‑138 (MQ=255)

              tACGCCGCATCAGGCATAGACAAAGAGCACTTTATCAACAAACTACATTTGTACCCGTAACGCTATGCTTCCGGCGGCAGACAGACGCCAATTCCGCCAATTCCACAGTAACCATGCGGGTTTTTATGCAGATATTGCtg                                                                                                                >  2:572496/1‑140 (MQ=255)

                   cgcATCAGGCATAGACAAAGAGCACTTTATCAACAAACTACATTTGTACCCGTAACGCTATGCTTCCGGCGGCAGACAGACGCCAATTCCGCCAATTCCACAGTAACCATGCGGGTTTTTATGCAGATATTGCTggtggt                                                                                                           <  1:502092/140‑1 (MQ=255)

                     cATCAGGCATAGACAAAGAGCACTTTATCAACAAACTACATTTGTACCCGTAACGCTATGCTTCCGGCGGCAGACAGACGCCAATTCCGCCAATTCCACAGTAACCATGCGGGTTTTTATGCAGATATTGCTGGTGGTCa                                                                                                         >  1:174617/1‑140 (MQ=255)

                            cATAGACAAAGAGCACTTTATCAACAAACTACATTTGTACCCGTAACGCTATGCTTCCGGCGGCAGACAGACGCCAATTCCGCCAATTCCACAGTAACCATGCGGGTTTTTATGCAGATATTGCTGGTGGTCATCTTCgg                                                                                                  <  1:8852/140‑1 (MQ=255)

                              tAGACAAAGAGCACTTTATCAACAAACTACATTTGTACCCGTAACGCTATGCTTCCGGCGGCGGACAGACGCCAATTCCGCCAATTCCACAGTAACCATGCGGGTTTTTATGCAGATATTGCTGGTGGTCATCTTCGGCa                                                                                                <  1:388999/140‑1 (MQ=255)

                                gACAAAGAGCACTTTATCAACAAACTACATTTGTACCCGTAACGCTATGCTTCCGGCGGCAGACAGACGCCAATTCCGCCAATTCCACAGTAACCATGCGGGTTTTTATGCAGATATTGCTGGTGGTCATCTTCGGCata                                                                                              >  2:496243/1‑140 (MQ=255)

                                     agagCACTTTATCAACAAACTACATTTGTACCCGTAACGCTATGCTTCCGGCGGCAGACAGACGCCAATTCCGCCAATTCCACAGTAACCATGCGGGTTTTTATGCAGATATTGCTGGTGGTCATCTTCGGCa                                                                                                >  2:89063/1‑133 (MQ=255)

                                     agagCACTTTATCAACAAACTACATTTGTACCCGTAACGCTATGCTTCCGGCGGCAGACAGACGCCAATTCCGCCAATTCCACAGTAACCATGCGGGTTTTTATGCAGATATTGCTGGTGGTCATCTTCGGCa                                                                                                <  1:89063/133‑1 (MQ=255)

                                     agagCACTTTATCAACAAACTACATTTGTACCCGTAACGCTATGCTTCCGGCGGCAGACAGACGCCAATTCCGCCAATTCCACAGTAACCATGCGGGTTTTTATGCAGATATTGCTGGTGGTCATCTTCGGCATAATaaa                                                                                         <  2:179321/140‑1 (MQ=255)

                                      gagCACTTTATCAACAAACTACATTTGTACCCGTAACGCTATGCTTCCGGCGGCAGACAGACGCCAATTCCGCCAATTCCACAGTAACCATGCGGGTTTTTATGCAGATATTGCTGGTGGTCATCTTCGGCATAATaaaa                                                                                        <  1:173468/140‑1 (MQ=255)

                                         cACTTTATCAACAAACTACATTTGTACCCGTAACGCTATGCTTCCGGCGGCAGACAGACGCCAATTCCGCCAATTCCACAGTAACCATGCGGGTTTTTATGCAGATATTGCTGGTGGTCATCTTCGGCATAATAAAACgg                                                                                     >  2:457567/1‑140 (MQ=255)

                                                tCAACAAACTACATTTGTACCCGTAACGCTATGCTTCCGGCGGCAGACAGACGCCAATTCCGCCAATTCCACAGTAACCATGCGGGTTTTTATGCAGATATTGCTGGTGGTCATCTTCGGCATAATAAAACGGTGTGGCg                                                                              <  1:286245/140‑1 (MQ=255)

                                                 caacaaACTACATTTGTACCCGTAACGCTATGCTTCCGGCGGCAGACAGACGCCAATTCCGCCAATTCCACAGTAACCATGCGGGTTTTTATGCAGATATTGCTGGTGGTCATCTTCGGCATAATAAAACGGTGTGGCGt                                                                             <  2:163308/140‑1 (MQ=255)

                                                             tttGTACCCGTAACGCTATGCTTCCGGCGGCAGACAGACGCCAATTCCGCCAATTCCACAGTAACCATGCGGGTTTTTATGCAGATATTGCTGGTGGTCATCTTCGGCATAATAAAACGGTGTGGCGTTAGCGATTTCCg                                                                 <  1:568116/140‑1 (MQ=255)

                                                             tttGTACCCGTAACGCTATGCTTCCGGCGGCAGACAGACGCCAATTCCGCCAATTCCACAGTAACCATGCGGGTTTTTATGCAGATATTGCTGGTGGTCATCTTCGGCATAATAAAACGGTGTGGCGTTAGCGATTTCCg                                                                 <  1:135397/140‑1 (MQ=255)

                                                               tGTACCCGTAACGCTATGCTTCCGGCGGCAGACAGACGCCAATTCCGCCAATTCCACAGTAACCATGCGGGTTTTTATGCAGATATTGCTGGTGGTCATCTTCGGCATAATAAAACGGTGTGGCGTTAGCGATTTCCgtg                                                               >  1:50017/1‑140 (MQ=255)

                                                               tGTACCCGTAACGCTATGCTTCCGGCGGCAGACAGACGCCAATTCCGCCAATTCCACAGTAACCATGCGGGTTTTTATGCAGATATTGCTGGTGGTCATCTTCGGCATAATAAAACGGTGTGGCGTTAGCGATTTCCgtg                                                               <  2:55749/140‑1 (MQ=255)

                                                                gTACCCGTAACGCTATGCTTCCGGCGGCAGACAGACGCCAATTCCGCCAATTCCACAGTAACCATGCGGGTTTTTATGCAGATATTGCTGGTGGTCATCTTCGGCATAATAAAACGGTGTGGCGTTAGCGATTTCCgtgg                                                              >  2:14881/1‑140 (MQ=255)

                                                                 tACCCGTAACGCTATGCTTCCGGCGGCAGACAGACGCCAATTCCGCCAATTCCACAGTAACCATGCGGGTTTTTATGCAGATATTGCTGGTGGTCATCTTCGGCATAATAAAACGGTGTGGCGTTAGCGATTTCCgtggt                                                             >  1:490915/1‑140 (MQ=255)

                                                                   cccGTAACGCTATGCTTCCGGCGGCAGACAGACGCCAATTCCGCCAATTCCACAGTAACCATGCGGGTTTTTATGCAGATATTGCTGGTGGTCATCTTCGGCATAATAAAACGGTGTGGCGTTAGCGATTTCCGTGGTGa                                                           >  2:539469/1‑140 (MQ=255)

                                                                        aaCGCTATGCTTCCGGCGGCAGACAGACGCCAATTCCGCCAATTCCACAGTAACCATGCGGGTTTTTATGCAGATATTGCTGGTGGTCATCTTCGGCATAATAAAACGGTGTGGCGTTAGCGATTTCCGTGGTGATGTGa                                                      <  1:524595/140‑1 (MQ=255)

                                                                        aaCGCTATGCTTCCGGCGGCAGACAGACGCCAATTCCGCCAATTCCACAGTAACCATGCGGGTTTTTATGCAGATATTGCTGGTGGTCATCTTCGGCATAATAAAACGGTGTGGCGTTAGCGATTTCCGTGGTGATGTGa                                                      <  2:347019/140‑1 (MQ=255)

                                                                                 ctcccGGCGGCAGACAGACGCCAATTCCGCCAATTCCACAGTAACCATGCGGGTTTTTATGCAGATATTGCTGGTGGTCATCTTCGGCATAATAAAACGGTGTGGCGTTAGCGATTTCCGTGGTGATGTGACGATCGTCa                                             >  2:555150/4‑140 (MQ=255)

                                                                                    ccGGCGGCAGACAGACGCCAATTCCGCCAATTCCACAGTAACCATGCGGGTTTTTATGCAGATATTGCTGGTGGTCATCTTCGGCATAATAAAACGGTGTGGCGTTAGCGATTTCCGTGGTGATGTGACGATCGTCATcg                                          >  1:370105/1‑140 (MQ=255)

                                                                                     cggcggCAGACAGACGCCAATTCCGCCAATTCCACAGTAACCATGCGGGTTTTTATGCAGATATTGCTGGTGGTCATCTTCGGCATAATAAAACGGTGTGGCGTTAGCGATTTCCGTGGTGATGTGACGATCGTCATcgg                                         >  2:576883/1‑140 (MQ=255)

                                                                                                                            aCCATGCGGGTTTTTATGCAGATATTGCTGGTGGTCATCTTCGGCATAATAAAACGGTGTGGCGTTAGCGATTTCCGTGGTGATGTGACGATCGTCATCGGCGGCAAGCATCGCCGCCTGAAAACGTTCCAGACTGgcgc  >  2:303738/1‑140 (MQ=255)

                                                                                                                            aCCATGCGGGTTTTTATGCAGATATTGCTGGTGGTCATCTTCGGCATAATAAAACGGTGTGGCGTTAGCGATTTCCGTGGTGATGTGACGATCGTCATCGGCGGCAAGCATCGCCGCCTGAAAACGTTCCAGACTGgcgc  >  1:113665/1‑140 (MQ=255)

CCTGGTAAGTCGTTTACGCCGCATCAGGCATAGACAAAGAGCACTTTATCAACAAACTACATTTGTACCCGTAACGCTATGCTTCCGGCGGCAGACAGACGCCAATTCCGCCAATTCCACAGTAACCATACGGGTTTTTATGCAGATATTGCTGGTGGTCATCTTCGGCATAATAAAACGGTGTGGCGTTAGCGATTTCCGTGGTGATGTGACGATCGTCATCGGCGGCAAGCATCGCCGCCTGAAAACGTTCCAGACTGGCGC  >  NC_000913/4441462‑4441725

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