BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,146,011	C→T	195 (TAG→TAA)	yceF ←	m(7)GTP pyrophosphatase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,146,011	0	C	T	100.0%	118.2 / NA	37	195 (TAG→TAA)	yceF	m(7)GTP pyrophosphatase
Reads supporting (aligned to +/- strand): ref base C (0/0); new base T (21/16); total (21/16)

CTGCGTAACGCGCGCTAATACAAGCATATAAACCTGATAATGGGCGGTTGGCCTGATAAGGCGTTTACGCCGCATCAGGCCGCCAGCACCGATTGCCGGATGCGACGTAACCACATCCGGCGCAAAAAACTATCCCATCAGCGGGTTTTTCCCTTCCCGACGTAACATCTGACACAGCGCGATTAGCGGTAAACCAACCAGCGTGTTAGGATCACGCCCCTCTAAACGCTCAAACAGCGTAATGCCAAATCCTTCACTC  >  NC_000913/1145882‑1146140

cTGCGTAACGCGCGCTAATACAAGCATATAAACCTGATAATGGGCGGTTGGCCTGATAAGGCGTTTACGCCGCATCAGGCCGCCAGCACCGATTGCCGGATGCGACGTAACCACATCCGGCGCAAAAAATTATCCCATc                                                                                                                          <  1:532564/139‑1 (MQ=255)

   cGTAACGCGCGCTAATACAAGCATATAAACCTGATAATGGGCGGTTGGCCTGATAAGGCGTTTACGCCGCATCAGGCCGCCAGCACCGATTGCCGGATGCGACGTAACCACATCCGGCGCAAAAAATTATCCCATCAGc                                                                                                                       >  1:344492/1‑139 (MQ=255)

         gcgcgcTAATACAAGCATATAAACCTGATAATGGGCGGTTGGCCTGATAAGGCGTTTACGCCGCATCAGGCCGCCAGCACCGATTGCCGGATGCGACGTAACCACATCCGGCGCAAAAAATTATCCCATCAGCGGGttt                                                                                                                 >  1:136557/1‑139 (MQ=255)

           gcgcTAATACAAGCATATAAACCTGATAATGGGCGGTTGGCCTGATAAGGCGTTTACGCCGCATCAGGCCGCCAGCACCGATTGCCGGATGCGACGTAACCACATCCGGCGCAAAAAATTATCCCATCAGCGGGttttt                                                                                                               >  1:273563/1‑139 (MQ=255)

                aaTACAAGCATATAAACCTGATAATGGGCGGTTGGCCTGATAAGGCGTTTACGCCGCATCAGGCCGCCAGCACCGATTGCCGGATGCGACGTAACCACATCCGGCGCAAAAAATTATCCCATCAGCGGGTTTttccctt                                                                                                          <  1:150573/139‑1 (MQ=255)

                  tACAAGCATATAAACCTGATAATGGGCGGTTGGCCTGATAAGGCGTTTACGCCGCATCAGGCCGCCAGCACCGATTGCCGGATGCGACGTAACCACATCCGGCGCAAAAAATTATCCCATCAGCGGGTTTttcccttcc                                                                                                        <  1:683017/139‑1 (MQ=255)

                         atatAAACCTGATAATGGGCGGTTGGCCTGATAAGGCGTTTACGCCGCATCAGGCCGCCAGCACCGATTGCCGGATGCGACGTAACCACATCCGGCGCAAAAAATTATCCCATCAGCGGGTTTTTCCCTTCCCGACGTa                                                                                                 >  2:460106/1‑139 (MQ=255)

                         atatAAACCTGATAATGGGCGGTTGGCCTGATAAGGCGTTTACGCCGCATCAGGCCGCCAGCACCGATTGCCGGATGCGACGTAACCACATCCGGCGCAAAAAATTATCCCATCAGCGGGTTTTTCCCTTCCCGACGTa                                                                                                 <  2:456466/139‑1 (MQ=255)

                           ataAACCTGATAATGGGCGGTTGGCCTGATAAGGCGTTTACGCCGCATCAGGCCGCCAGCACCGATTGCCGGATGCGACGTAACCACATCCGGCGCAAAAAATTATCCCATCAGCGGGTTTTTCCCTTCCCGACGTa                                                                                                 >  2:756264/1‑137 (MQ=255)

                           ataAACCTGATAATGGGCGGTTGGCCTGATAAGGCGTTTACGCCGCATCAGGCCGCCAGCACCGATTGCCGGATGCGACGTAACCACATCCGGCGCAAAAAATTATCCCATCAGCGGGTTTTTCCCTTCCCGACGTa                                                                                                 <  1:756264/137‑1 (MQ=255)

                           ataAACCTGATAATGGGCGGTTGGCCTGATAAGGCGTTTACGCCGCATCAGGCCGCCAGCACCGATTGCCGGATGCGACGTAACCACATCCGGCGCAAAAAATTATCCCATCAGCGGGTTTTTCCCTTCCCGACGTAAc                                                                                               >  2:577515/1‑139 (MQ=255)

                           ataAACCTGATAATGGGCGGTTGGCCTGATAAGGCGTTTACGCCGCATCAGGCCGCCAGCACCGATTGCCGGATGCGACGTAACCACATCCGGCGCAAAAAATTATCCCATCAGCGGGTTTTTCCCTTCCCGACGTAAc                                                                                               >  1:14642/1‑139 (MQ=255)

                                 cTGATAATGGGCGGTTGGCCTGATAAGGCGTTTACGCCGCATCAGGCCGCCAGCACCGATTGCCGGATGCGACGTAACCACATCCGGCGCAAAAAATTATCCCATCAGCGGGTTTTTCCCTTCCCGACGTAACATCTGa                                                                                         >  1:269818/1‑139 (MQ=255)

                                    aTAATGGGCGGTTGGCCTGATAAGGCGTTTACGCCGCATCAGGCCGCCAGCACCGATTGCCGGATGCGACGTAACCACATCCGGCGCAAAAAATTATCCCATCAGCGGGTTTTTCCCTTCCCGACGTa                                                                                                 >  1:492305/1‑128 (MQ=255)

                                    aTAATGGGCGGTTGGCCTGATAAGGCGTTTACGCCGCATCAGGCCGCCAGCACCGATTGCCGGATGCGACGTAACCACATCCGGCGCAAAAAATTATCCCATCAGCGGGTTTTTCCCTTCCCGACGTa                                                                                                 <  2:492305/128‑1 (MQ=255)

                                          ggCGGTTGGCCTGATAAGGCGTTTACGCCGCATCAGGCCGCCAGCACCGATTGCCGGATGCGACGTAACCACATCCGGCGCAAAAAATTATCCCATCAGCGGGTTTTTCCCTTCCCGACGTAACATCTGACACAGcgcg                                                                                >  1:560504/1‑139 (MQ=255)

                                                         aaGGCGTTTACGCCGCATCAGGCCGCCAGCACCGATTGCCGGATGCGACGTAACCACATCCGGCGCAAAAAATTATCCCATCAGCGGGTTTTTCCCTTCCCGACGTAACATCTGACACAGCGCGATTAGCGGTAaacca                                                                 <  1:577515/139‑1 (MQ=255)

                                                            gCGTTTACGCCGCATCAGGCCGCCAGCACCGATTGCCGGATGCGACGTAACCACATCCGGCGCAAAAAATTATCCCATCAGCGGGTTTTTCCCTTCCCGACGTAACATCTGACACAGCGCGATTAGCGGTAaaccaacc                                                              >  2:580046/1‑139 (MQ=255)

                                                              gTTTACGCCGCATCAGGCCGCCAGCACCGATTGCCGGATGCGACGTAACCACATCCGGCGCAAAAAATTATCCCATCAGCGGGTTTTTCCCTTCCCGACGTAACATCTGACACAGCGCGATTAGCGGTAAACCAACCAg                                                            <  2:636628/139‑1 (MQ=255)

                                                              gTTTACGCCGCATCAGGCCGCCAGCACCGATTGCCGGATGCGACGTAACCACATCCGGCGCAAAAAATTATCCCATCAGCGGGTTTTTCCCTTCCCGACGTAACATCTGACACAGCGCGATTAGCGGTAAACCAACCAg                                                            <  1:580046/139‑1 (MQ=255)

                                                              gTTTACGCCGCATCAGGCCGCCAGCACCGATTGCCGGATGCGACGTAACCACATCCGGCGCAAAAAATTATCCCATCAGCGGGTTTTTCCCTTCCCGACGTAACATCTGACACAGCGCGATTAGCGGTAAACCAACCAg                                                            <  2:14642/139‑1 (MQ=255)

                                                                 tACGCCGCATCAGGCCGCCAGCACCGATTGCCGGATGCGACGTAACCACATCCGGCGCAAAAAATTATCCCATCAGCGGGTTTTTCCCTTCCCGACGTAACATCTGACACAGCGCGATTAGCGGTAAACCAACCAGc                                                           <  1:13795/137‑1 (MQ=255)

                                                                 tACGCCGCATCAGGCCGCCAGCACCGATTGCCGGATGCGACGTAACCACATCCGGCGCAAAAAATTATCCCATCAGCGGGTTTTTCCCTTCCCGACGTAACATCTGACACAGCGCGATTAGCGGTAAACCAACCAGc                                                           >  2:13795/1‑137 (MQ=255)

                                                                 tACGCCGCATCAGGCCGCCAGCACCGATTGCCGGATGCGACGTAACCACATCCGGCGCAAAAAATTATCCCATCAGCGGGTTTTTCCCTTCCCGACGTAACATCTGACACAGCGCGATTAGCGGTAAACCAACCAGCgt                                                         >  1:183742/1‑139 (MQ=255)

                                                                         aTCAGGCCGCCAGCACCGATTGCCGGATGCGACGTAACCACATCCGGCGCATAAAATTATCCCATCAGCGGGTTTttcccttcc                                                                                                        >  2:475958/1‑84 (MQ=255)

                                                                         aTCAGGCCGCCAGCACCGATTGCCGGATGCGACGTAACCACATCCGGCGCAAAAAATTATCCCATCAGCGGGTTTttcccttcc                                                                                                        <  1:475958/84‑1 (MQ=255)

                                                                              gccgccAGCACCGATTGCCGGATGCGACGTAACCACATCCGGCGCAAAAAATTATCCCATCAGCGGGTTTTTCCCTTCCCGACGTAACATCTGACACAGCGCGATTAGCGGTAAACCAACCAGCGTGTTAGGATTACGc                                            >  1:664831/1‑139 (MQ=255)

                                                                                          gATTGCCGGGTGCGACGTAACCACATCCGGCGCAAAAAATTTTCCCATCAGCGGGTTTTTCCCTTCCCGACGTAACATCTGACACAGCGCGATTAGCGGTAAACCAACCAGCGTGTTAGGATCACGCCCCTCTAAACGc                                >  1:354711/1‑139 (MQ=255)

                                                                                          gATTGCCGGATGCGACGTAACCACATCCGGCGCAAAAAATTATCCCATCAGCGGGTTTTTCCCTTCCCGACGTAACATCTGACACAGCGCGATTAGCGGTAAACCAACCAGCGTGTTAGGATCACGCCCCTCTAAACGc                                <  2:136557/139‑1 (MQ=255)

                                                                                                                  aTCCGGCGCAAAAAATTATCCCATCAGCGGGTTTTTCCCTTCCCGACGTAACATCTGACACAGCGCGATTAGCGGTAAACCAACCAGCGTGTTAGGATCACGCCCCTCTAAACGCTCAAACAGCGTAATGCCAAATCCt        >  1:709407/1‑139 (MQ=255)

                                                                                                                    ccGGCGCAAAAAATTATCCCATCAGCGGGTTTTTCCCTTCCCGACGTAACATCTGACACAGCGCGATTAGCGGTAAACCAACCAGCGTGTTAGGATCACGCCCCTCTAAACGCTCAAACAGCGTAATGCCAAATCCTTc      <  2:344492/139‑1 (MQ=255)

                                                                                                                    ccGGCGCAAAAAATTATCCCATCAGCGGGTTTTTCCCTTCCCGACGTAACATCTGACACAGCGCGATTAGCGGTAAACCAACCAGCGTGTTAGGATCACGCCCCTCTAAACGCTCAAACAGCGTAATGCCAAATCCTTc      <  2:431765/139‑1 (MQ=255)

                                                                                                                    ccGGCGCAAAAAATTATCCCATCAGCGGGTTTTTCCCTTCCCGACGTAACATCTGACACAGCGCGATTAGCGGTAAACCAACCAGCGTGTTAGGATCACGCCCCTCTAAACGCTCAAACAGCGTAATGCCAAATCCTTc      <  1:460106/139‑1 (MQ=255)

                                                                                                                       gcgcAAAAAATTATCCCATCAGCGGGTTTTTCCCTTCCCGACGTAACATCTGACACAGCGCGATTAGCGGTAAACCAACCAGCGTGTTAGGATCACGCCCCTCTAAACGCTCAAACAGCGTAATGCCAAATCCTTCAct   >  2:108117/1‑139 (MQ=255)

                                                                                                                       gcgcAAAAAATTATCCCATCAGCGGGTTTTTCCCTTCCCGACGTAACATCTGACACAGCGCGATTAGCGGTAAACCAACCAGCGTGTTAGGATCACGCCCCTCTAAACGCTCAAACAGCGTAATGCCAAATCCTTCAct   >  1:746131/1‑139 (MQ=255)

                                                                                                                        cgcAAAAAATTATCCCATCAGCGGGTTTTTCCCTTCCCGACGTAACATCTGACACAGCGCGATTAGCGGTAAACCAACCAGCGTGTTAGGATCACGCCCCTCTAAACGCTCAAACAGCGTAATGCCAAATCCTTCActc  >  1:713321/1‑139 (MQ=255)

                                                                                                                        cgcAAAAAATTATCCCATCAGCGGGTTTTTCCCTTCCCGACGTAACATCTGACACAGCGCGATTAGCGGTAAACCAACCAGCGTGTTAGGATCACGCCCCTCTAAACGCTCAAACAGCGTAATGCCAAATCCTTCActc  >  1:416268/1‑139 (MQ=255)

CTGCGTAACGCGCGCTAATACAAGCATATAAACCTGATAATGGGCGGTTGGCCTGATAAGGCGTTTACGCCGCATCAGGCCGCCAGCACCGATTGCCGGATGCGACGTAACCACATCCGGCGCAAAAAACTATCCCATCAGCGGGTTTTTCCCTTCCCGACGTAACATCTGACACAGCGCGATTAGCGGTAAACCAACCAGCGTGTTAGGATCACGCCCCTCTAAACGCTCAAACAGCGTAATGCCAAATCCTTCACTC  >  NC_000913/1145882‑1146140

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