BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,442,007	T→C	E57G (GAG→GGG)	msrA ←	methionine sulfoxide reductase A

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,442,007	0	T	C	100.0%	131.7 / NA	37	E57G (GAG→GGG)	msrA	methionine sulfoxide reductase A
Reads supporting (aligned to +/- strand): ref base T (0/0); new base C (18/19); total (18/19)

TCGTAAACAATGCGTACCGCTTCGGCATGACCCGTATCACCGGAGCACACTTCCCGATAAGTCGGATTTGGCGTATAGCCGCCGGTATAGCCTGCGGCGGTGCTGTAAACGCCGGGTAACTGCCAGAACAGACGCTCCACACCCCAGAAACAACCCATCGCAAAAATGGCAATCTCCATTCCGTCAGGTACATTGGTCATTGAGTGACCGTTGACCGCATGCAGCGTGGCTACGGGCATCGGGGTGTTACGTCCAGG  >  NC_000913/4441872‑4442128

tCGTAAACAATGCGTACCGCTTCGGCATGACCCGTATCACCGGAGCACACTTCCCGATAAGTCGGATTTGGCGTATAGCCGCCGGTATAGCCTGCGGCGGTGCTGTAAACGCCGGGTAACTGCCAGAACAGACGCCCca                                                                                                                        <  2:71309/139‑1 (MQ=255)

 cGTAAACAATGCGTACCGCTTCGGCATGACCCGTATCACCGGAGCACACTTCCCGATAAGTCGGATTTGGCGTATAGCTGCCGGTATAGCCTGCGGCGGTGCTGTAAACGCCGGGTAACTGCCAGAACAGACGCCCcac                                                                                                                       <  2:258485/139‑1 (MQ=255)

      aCAATGCGTACCGCTTCGGCATGACCCGTATCACCGGAGCACACTTCCCGATAAGTCGGATTTGGCGTATAGCCGCCGGTATAGCCTGCGGCGGTGCTGTAAACGCCGGGTAACTGCCAGAACAGACGCCCCACAcccc                                                                                                                  <  1:736885/139‑1 (MQ=255)

          tGCGTACGGCTTCGGCATGACCCGTATCACCGGAGCACACATCCCGATAAGTCGGATTTGGCGTATAGCCGCCGGTATAGCCTGCGGCGGTGCTGTAAACGCCGGGTAACTGCCAGAACAGACGCCCCACACCCCAGaa                                                                                                              >  2:488650/1‑139 (MQ=255)

                  gCTTCGGCATGACCCGTATCACCGGAGCACACTTCCCGATAAGTCGGATTTGGCGTATAGCCGCCGGTATAGCCTGCGGCGGTGCTGTAAACGCCGGGTAACTGCCAGAACAGACGCCCCACACCCCAGAAACAACCCa                                                                                                      >  2:301327/1‑139 (MQ=255)

                    ttCGGCATGACCCGTATCACCGGAGCACACTTCCCGATAAGTCGGATTTGGCGTATAGCCGCCGGTATAGCCTGCGGCGGTGCTGTAAACGCCGGGTAACTGCCAGAACAGACGCCCCACACCCCAGAAACAACCCATc                                                                                                    >  1:679178/1‑139 (MQ=255)

                        gCATGACCCGTATCACCGGAGCACACTTCCCGATAAGTCGGATTTGGCGTATAGCCGCCGGTATAGCGTGCGGCGGTGCTGTAAACGCCGGGTAACTGCCAGAACAGACGCCCCACACCCCAGAAACAACCCATCGCaa                                                                                                >  2:499420/1‑139 (MQ=255)

                              cccGTATCACCGGAGCACACTTCCCGATAAGTCGGATTTGGCGTATAGCCGCCGGTATAGCCTGCGGCGGTGCTGTAAACGCCGGGTAACTGCCAGAACAGACGCCCCACACCCCAGAAACAACCCATCGCAAAAATgg                                                                                          >  1:100119/1‑139 (MQ=255)

                                  tATCACCGGAGCACACTTCCCGATAAGTCGGATTTGGCGTATAGCCGCCGGTATAGCCTGCGGCGGTGCTGTAAACGCCGGGTAACTGCCAGAACAGACGCCCCACACCCCAGAAACAACCCATCGCAAAAATGGCAAt                                                                                      <  1:301327/139‑1 (MQ=255)

                                       ccGGAGCACACTTCCCGATAAGTCGGATTTGGCGTATAGCCGCCGGTATAGCCTGCGGCGGTGCTGTAAACGCCGGGTAACTGCCAGAACAGACGCCCCACACCCCAGAAACAACCCATCGCAAAAATGGCAATCTCCa                                                                                 <  2:201464/139‑1 (MQ=255)

                                        cGGAGCACACTTCCCGATAAGTCGGATTTGGCGTATAGCCGCCGGTATAGCCTGCGGCGGTGCTGTAAACGCCGGGTAACTGCCAGAACAGACGCCCCACACCCCAGAAACAACCCATCGCAAAAATGGCAATCTCCAt                                                                                >  2:566908/1‑139 (MQ=255)

                                        cGGAGCACACTTCCCGATAAGTCGGATTTGGCGTATAGCCGCCGGTATAGCCTGCGGCGGTGCTGTAAACGCCGGGTAACTGCCAGAACAGACGCCCCACACCCCAGAAACAACCCATCGCAAAAATGGCAATCTCCAt                                                                                <  1:230785/139‑1 (MQ=255)

                                           aGCACACTTCCCGATAAGTCGGATTTGGCGTATAGCCGCCGGTATAGCCTGCGGCGGTGCTGTAAACGCCGGGTAACTGCCAGAACAGACGCCCCACACCCCAGAAACAACCCATCGCAAAAATGGCAATCTCCATTcc                                                                             >  1:403296/1‑139 (MQ=255)

                                            gCACACTTCCCGATAAGTCGGATTTGGCGTATAGCCGCCGGTATAGCCTGCGGCGGTGCTGTAAACGCCGGGTAACTGCCAGAACAGACGCCCCACACCCCAGAAACAACCCATCGCAAAAATGGCAATCTCCATTCCg                                                                            >  2:587931/1‑139 (MQ=255)

                                             cacaCTTCCCGATAAGTCGGATTTGGCGTATAGCCGCCGGTATAGCCTGCGGCGGTGCTGTAAACGCCGGGTAACTGCCAGAACAGACGCCCCACACCCCAGAAACAACCCATCGCAAAAATGGCAATCTCCATTc                                                                              >  1:637331/1‑136 (MQ=255)

                                             cacaCTTCCCGATAAGTCGGATTTGGCGTATAGCCGCCGGTATAGCCTGCGGCGGTGCTGTAAACGCCGGGTAACTGCCAGAACAGACGCCCCACACCCCAGAAACAACCCATCGCAAAAATGGCAATCTCCATTc                                                                              <  2:637331/136‑1 (MQ=255)

                                              acacTTCCCGATAAGTCGGATTTGGCGTATAGCCGCCGGTATAGCCTGCGGCGGTGCTGTAAACGCCGGGTAACTGCCAGAACAGACGCCCCACACCCCAGAAACAACCCATCGCAAAAATGGCAATCTCCATTCCGTc                                                                          <  1:496878/139‑1 (MQ=255)

                                                    cccGATAAGTCGGATTTGGCGTATAGCCGCCGGTATAGCCTGCGGCGGTGCTGTAAACGCCGGGTAACTGCCAGAACAGACGCCCCACACCCCAGAAACAACCCATCGCAAAAATGGCAATCTCCATTCCGTCAGGTAc                                                                    <  1:792255/139‑1 (MQ=255)

                                                      cGATAAGTCGGATTTGGCGTATAGCCGCCGGTATAGCCTGCGGCGGTGCTGTAAACGCCGGGTAACTGCCAGAACAGACGCCCCACACCCCAGAAACAACCCATCGCAAAAATGGCAATCTCCATTCCGTCAGGTACAt                                                                  <  2:100119/139‑1 (MQ=255)

                                                        aTAAGTCGGATTTGGCGTATAGCCGCCGGTATAGCCTGCGGCGGTGCTGTAAACGCCGGGTAACTGCCAGAACAGACGCCCCACACCCCAGAAACAACCCATCGCAAAAATGGCAATCTCCATTCCGTCAGGTACATTg                                                                <  2:679178/139‑1 (MQ=255)

                                                          aaGTCGGATTTGGCGTATAGCCGCCGGTATAGCCTGCGGCGGTGCTGTAAACGCCGGGTAACTGCCAGAACAGACGCCCCACACCCCAGAAACAACCCATCGCAAAAATGGc                                                                                         >  2:325615/1‑112 (MQ=255)

                                                          aaGTCGGATTTGGCGTATAGCCGCCGGTATAGCCTGCGGCGGTGCTGTAAACGCCGGGTAACTGCCAGAACAGACGCCCCACACCCCAGAAACAACCCATCGCAAAAATGGc                                                                                         <  1:325615/112‑1 (MQ=255)

                                                            gTCGGATTTGGCGTATAGCCGCCGGTATAGCCTGCGGCGGTGCTGTAAACGCCGGGTAACTGCCAGAACAGACGCCCCACACCCCAGAAACAACCCATCGCAAAAATGGCAATCTCCATTCCGTCAGGTACATTGGTCa                                                            >  1:394914/1‑139 (MQ=255)

                                                            gTCGGATTTGGCGTATAGCCGCCGGTATAGCCTGCGGCGGTGCTGTAAACGCCGGGTAACTGCCAGAACAGACGCCCCACACCCCAGAAACAACCCATCGCAAAAATGGCAATCTCCATTCCGTCAGGTACATTGGTCa                                                            >  2:21816/1‑139 (MQ=255)

                                                                            aGCCGCCGGTATAGCCTGCGGCGGTGCTGTAAACGCCGGGTAACTGCCAGAACAGACGCCCCACACCCCAGAAACAACCCATCGCAAAAATGGCAATCTCCATTCCGTCAGGTACATTGGTCATTGAGTGACCGTTGAc                                            <  2:231217/139‑1 (MQ=255)

                                                                                gccgGTATAGCCTGCGGCGGTGCTGTAAACGCCGGGTAACTGCCAGAACAGACGCCCCACACCCCAGAAACAACCCATCGCAAAAATGGCAATCTCCATTCCGTCAGGTACATTGGTCATTGAGTGACCGTTGACCGCa                                        >  2:192248/1‑139 (MQ=255)

                                                                                gccgGTATAGCCTGCGGCGGTGCTGTAAACGCCGGGTAACTGCCAGAACAGACGCCCCACACCCCAGAAACAACCCATCGCAAAAATGGCAATCTCCATTCCGTCAGGTACATTGGTCATTGAGTGACCGTTGACCGCa                                        >  2:801496/1‑139 (MQ=255)

                                                                                         gCCTGCGGCGGTGCTGTAAACGCCGGGTAACTGCCAGAACAGACGCCCCACACCCCAGAAACAACCCATCGCAAAAATGGCAATCTCCATTCCGTCAGGTACATTGGTCATTGAGTGACCGTTGACCGCATGCAGCGTg                               <  1:587931/139‑1 (MQ=255)

                                                                                          ccTGCGGCGGTGCTGTAAACGCCGGGTAACTGCCAGAACAGACGCCCCACACCCCAGAAACAACCCATCGCAAAAATGGCAATCTCCATTCCGTCAGGTACATTGGTCATTGAGTGACCGTTGACCGCATGCAGCGTgg                              >  2:769854/1‑139 (MQ=255)

                                                                                                   gTGCTGTAAACGCCGGGTAACTGCCAGAACAGACGCCCCACACCCCAGAAACAACCCATCGCAAAAATGGCAATCTCCATTCCGTCAGGTACATTGGTCATTGAGTGACCGTTGACCGCATGCAGCGTGGCTACGGGCa                     >  2:167677/1‑139 (MQ=255)

                                                                                                      cTGTAAACGCCGGGTAACTGCCAGAACAGACGCCCCACACCCCAGAAACAACCCATCGCAAAAATGGCAATCTCCATTCCGTCAGGTACATTGGTCATTGAGTGACCGTTGACCGCATGCAGCGTGGCTACGGGCATCg                  <  1:21816/139‑1 (MQ=255)

                                                                                                       tGTAAACGCCGGGTAACTGCCAGAACAGACGCCCCACACCCCAGAAACAACCCATCGCAAAAATGGCAATCTCCATTCCGTCAGGTACATTGGTCATTGAGTGACCGTTGACCGCATGCAGCGTGGCTACGGGCATCgg                 <  1:577489/139‑1 (MQ=255)

                                                                                                            aCGCCGGGTAACTGCCAGAACAGACGCCCCACACCCCAGAAACAACCCATCGCAAAAATGGCAATCTCCATTCCGTCAGGTACATTGGTCATTGAGTGACCGTTGACCGCATGCAGCGTGGCTACGGGCATCGGGgtgt            <  2:39087/139‑1 (MQ=255)

                                                                                                                cGGGTAACTGCCAGAACAGACGCCCCACACCCCAGAAACAACCCATCGCAAAAATGGCAATCTCCATTCCGTCAGGTACATTGGTCATTGAGTGACCGTTGACCGCATGCAGCGTGGCTACGGGCATCGGGGTGTTACg        >  1:562309/1‑139 (MQ=255)

                                                                                                                      aCTGCCAGAACAGACGCCCCACACCCCAGAAACAACCCATCGCAAAAATGGCAATCTCCATTCCGTCAGGTACATTGGTCATTGAGTGACCGTTGACCGCATGCAGCGTGGCTACGGGCATCGGGGTGTTACGTCcagg  <  1:192248/139‑1 (MQ=255)

                                                                                                                                  gACGCCCCACACCCCAGAAACAACCCATCGCAAAAATGGCAATCTCCATTCCGTCAGGTACATTGGTCATTGAGTGACCGTTGACCGCATg                                      <  1:693924/91‑1 (MQ=255)

                                                                                                                                  gACGCCCCACACCCCAGAAACAACCCATCGCAAAAATGGCAATCTCCATTCCGTCAGGTACATTGGTCATTGAGTGACCGTTGACCGCATg                                      >  2:693924/1‑91 (MQ=255)

TCGTAAACAATGCGTACCGCTTCGGCATGACCCGTATCACCGGAGCACACTTCCCGATAAGTCGGATTTGGCGTATAGCCGCCGGTATAGCCTGCGGCGGTGCTGTAAACGCCGGGTAACTGCCAGAACAGACGCTCCACACCCCAGAAACAACCCATCGCAAAAATGGCAATCTCCATTCCGTCAGGTACATTGGTCATTGAGTGACCGTTGACCGCATGCAGCGTGGCTACGGGCATCGGGGTGTTACGTCCAGG  >  NC_000913/4441872‑4442128

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