BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,761,850	C→T	A36T (GCG→ACG)	yibF ←	glutathione S‑transferase homolog

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,761,850	0	C	T	100.0%	107.3 / NA	34	A36T (GCG→ACG)	yibF	glutathione S‑transferase homolog
Reads supporting (aligned to +/- strand): ref base C (0/0); new base T (20/14); total (20/14)

ATTCATTAATTCAATATATTCAGCGATGATCGGCGAATCAAACCAGCATTCGCCCTCTTCGGTCACCAGCACCGGCACTTTTCCTAACGGGTTAAATTGCGCCACGCCGTTGTCCGCGTTATAGGGCAGTTCATTAATAAATTCGAAAGTTATGCCCTTTTCTAACAACAGAATAGAAAGTTTGCGTACAAACGGGCTGGTGTAGCTACCGACGAGTTTCATGCCGAGT  >  NC_000913/3761734‑3761962

attcattAATTCAATATATTCAGCGATGATCGGCGAATCAAACCAGCATTCGCCCTCTTCGGTCACCAGCACCGGCACTTTTCCTAACGGGTTAAATTGCGCCACGCCGTTGTCCGTGTTATAGGGCAGTTCATtaata                                                                                            <  1:119964/139‑1 (MQ=255)

       aaTTCAATATATTCAGCGATGATCGGCGAATCAAACCAGCATTCGCCCTCTTCGGTCACCAGCACCGGCACTTTTCCTATCGGGTTAAATTGCGCCACGCCGTTGTCCGTGTTATAGGGCAGTTCATTAATAAATTCGa                                                                                     <  2:118113/139‑1 (MQ=255)

                      gCGATGATCGGCGAATCAAACCAGCATTCGCCCTCTTCGGTCACCAGCACCGGCACTTTTCCTAACGGGTTAAATTGCGCCACGCCGTTGTCCGTGTTATAGGGCAGTTCATTAATAAATTCGa                                                                                     <  2:200460/124‑1 (MQ=255)

                      gCGATGATCGGCGAATCAAACCAGCATTCGCCCTCTTCGGTCACCAGCACCGGCACTTTTCCTAACGGGTTAAATTGCGCCACGCCGTTGTCCGTGTTATAGGGCAGTTCATTAATAAATTCGa                                                                                     >  1:200460/1‑124 (MQ=255)

                       cGATGATCGGCGAATCAAACCAGCATTCGCCCTCTTCGGTCACCAGCACCGGCACTTTTCCTAACGGGTTAAATTGCGCCACGCCGTTGTCCGTGTTATAGGGCAGTTCATTAATAAATTCGAAAGTTATGCCCTTTTc                                                                     >  2:458268/1‑139 (MQ=255)

                        gatgatCGGCGAATCAAACCAGCATTCGCCCTCTTCGGTCACCAGCACCGGCACTTTTCCTAACGGGTTAAATTGCGCCACGCCGTTGTCCGTGTTATAGGGCAGTTCATTAATAAATTCGAAAGTTATGCCCTTTTCt                                                                    >  1:793288/1‑139 (MQ=255)

                            atCGGCGAATCAAACCAGCATTCGCCCTCTTCGGTCACCAGCACCGGCACTTTTCCTAACGGGTTAAATTGCGCCACGCCGTTGTCCGTGTTATAGGGCAGTTCATTAATAAATTCGAAAGTTATGCCCTTTTCTaaca                                                                >  2:262118/1‑139 (MQ=255)

                            atCGGCGAATCAAACCAGCATTCGCCCTCTTCGGTCACCAGCACCGGCACTTTTCCTAACGGGTTAAATTGCGCCACGCCGTTGTCCGTGTTATAGGGCAGTTCATTAATAAATTCGAAAGTTATGCCCTTTTCTaaca                                                                >  1:830293/1‑139 (MQ=255)

                            atCGGCGAATCAAACCAGCATTCGCCCTCTTCGGTCACCAGCACCGGCACTTTTCCTAACGGGTTAAATTGCGCCACGCCGTTGTCCGTGTTATAGGGCAGTTCATTAATAAATTCGAAAGTTATGCCCTTTTCTaaca                                                                >  2:139026/1‑139 (MQ=255)

                               ggCGAATCAAACCCGCATTCGCCCTCTTCGGTCACCAGCACCGGCACTTTGCCTAACCGGTTGACTTGCGCCACGCCGTTGTCCGTGTTATAGGGCAGTTCATTAACAAATTCGAAAGTTATGCCCTTTTATAccaaca                                                             >  2:829734/1‑139 (MQ=255)

                                     tCAAACCAGCATTCGCCCTCTTCGGTCACCAGCACCGGCACTTTTCCTAACGGGTTAAATTGCGCCACGCCGTTGTCCGTGTTATAGGGCAGTTCATtaat                                                                                             <  1:868206/101‑1 (MQ=255)

                                     tCAAACCAGCATTCGCCCTCTTCGGTCACCAGCACCGGCACTTTTCCTAACGGGTTAAATTGCGCCACGCCGTTGTCCGTGTTATAGGGCAGTTCATtaat                                                                                             >  2:868206/1‑101 (MQ=255)

                                           cAGCATTCGCCCTCTTCGGTCACCAGCACCGGCACTTTTCCTAACGGGTTAAATTGCGCCACGCCGTTGTCCGTGTTATAGGGCAGTTCATTAATAAATTCGAAAGTTATGCCCt                                                                         <  1:451878/115‑1 (MQ=255)

                                           cAGCATTCGCCCTCTTCGGTCACCAGCACCGGCACTTTTCCTAACGGGTTAAATTGCGCCACGCCGTTGTCCGTGTTATAGGGCAGTTCATTAATAAATTCGAAAGTTATGCCCt                                                                         >  2:451878/1‑115 (MQ=255)

                                                  cGCCCTCTTCGGTCACCAGCACCGGCACTTTTCCTAACGGGTTAAATTGCGCCACGCCGTTGTCCGTGTTATATGGCAGTTCATTAATAAATTCGAAAGTTATGCCCTTTTCTAACAACAGAATAGAAAGTTTGCGTAc                                          >  1:206179/1‑139 (MQ=255)

                                                     ccTCTTCGGTCACCAGCACCGGCACTTTTCCTAACGGGTTAAATTGCGCCACGCCGTTGTCCGTGTTATAggg                                                                                                         >  2:479246/1‑73 (MQ=255)

                                                     ccTCTTCGGTCACCAGCACCGGCACTTTTCCTAACGGGTTAAATTGCGCCACGCCGTTGTCCGTGTTATAggg                                                                                                         <  1:479246/73‑1 (MQ=255)

                                                        ctTCGGTCACCAGCACCGGCACTTTTCCTAACGGGTTAAATTGCGCCACGCCGTTGTCCGTGTTATAGGGCAGTTCATTAATAAATTCGAAAGTTATGCCCTTTTCTAACAACAGAATAGAAAGTTTGCGTACAAACgg                                    >  2:305772/1‑139 (MQ=255)

                                                          tCGGTCACCAGCACCGGCACTTTTCCTAACGGGTTAAATTGCGCCACGCCGTTGTCCGTGTTATAGGGCAGTTCATTAATAAATTCGAAAGTTATGCCCTTTTCTAACAACAGAATAGAAAGTTTGCGTACAAACGGGc                                  <  1:262118/139‑1 (MQ=255)

                                                           cGGTCACCAGCACCGGCACTTTTCCTAACGGGTTAAATTGCGCCACGCCGTTGTCCGTGTTATAGGGCAGTTCATTAATAAATTCGAAAGTTATGCCCTTTTCTAACAACAGAATAGAAAGTTTGCGTACAAACGGGCt                                 <  1:139026/139‑1 (MQ=255)

                                                                 ccAGCACCGGCACTTTTCCTAACGGGTTAAATTGCGCCACGCCGTTGTCCGTGTTATAGGGc                                                                                                        >  2:632164/1‑62 (MQ=255)

                                                                 ccAGCACCGGCACTTTTCCTAACGGGTTAAATTGCGCCACGACGTTGTCCGTGTTATAGGGc                                                                                                        <  1:632164/62‑1 (MQ=255)

                                                                  cAGCACCGGCACTTTTCCTAACGGGTTAAATTGCGCCACGCCGTTGTCCGTGTTATAGGGCAGTTCATTAATAAATTCGAAAGTTATGCCCTTTTCTAACAACAGAATAGAAAGTTTGCGTACAAACGGGCTGGTGTAg                          <  2:892118/139‑1 (MQ=255)

                                                                          gCACTTTTCCTAACGGGTTAAATTGCGCCACGCCGTTGTCCGTGTTATAGGGCAGTTCATTAATAAATTCGAAAGTTATGCCCTTTTCTAACAACAGAATAGAAAGTTTGCGTACAAACGGGCTGGTGTAGCTACcgac                  <  1:305772/139‑1 (MQ=255)

                                                                          gCACTTTTCCTAACGGGTTAAATTGCGCCACGCCGTTGTCCGTGTTATAGGGCAGTTCATTAATAAATTCGAAAGTTATGCCCTTTTCTAACAACAGAATAGAAAGTTTGCGTACAAACGGGCTGGTGTAGCTACcgac                  >  1:293297/1‑139 (MQ=255)

                                                                          gCACTTTTCCTAACGGGTTAAATTGCGCCACGCCGTTGTCCGTGTTATAGGGCAGTTCATTAATAAATTCGAAAGTTATGCCCTTTTCTAACAACAGAATAGAAAGTTTGCGTACAAACGGGCTGGTGTAGCTACcgac                  >  1:161422/1‑139 (MQ=255)

                                                                          gCACTTTTCCTAACGGGTTAAATTGCGCCACGCCGTTGTCCGTGTTATAGGGCAGTTCATTAATAAATTCGAAAGTTATGCCCTTTTCTAACAACAGAATAGAAAGTTTGCGAACAAACGGGCTGGTGTAGTTACggac                  >  2:500386/1‑139 (MQ=255)

                                                                                          gTTAAATTGCGCCACGCCGTTGTCCGTGTTATAGGGCAGTTCATtaata                                                                                            <  2:166994/49‑1 (MQ=255)

                                                                                          gTTAAATTGCGCCACGCCGTTGTCCGTGTTATAGGGCAGTTCATtaata                                                                                            >  1:166994/1‑49 (MQ=255)

                                                                                          gTTAAATTGCGCCACGCCGTTGTCCGTGTTATAGGGCAGTTCATTAATAAATTCGAAAGTTATGCCCTTTTCTaacaac                                                              <  1:489892/79‑1 (MQ=255)

                                                                                          gTTAAATTGCGCCACGCCGTTGTCCGTGTTATAGGGCAGTTCATTAATAAATTCGAAAGTTATGCCCTTTTCTaacaac                                                              >  2:489892/1‑79 (MQ=255)

                                                                                          gTTAAATTGCGCCACGCCGTTGTCCGTGTTATAGGGCAGTTCATTAATAAATTCGAAAGTTATGCCCTTTTCTaa                                                                  >  1:797224/1‑75 (MQ=255)

                                                                                          gTTAAATTGCGCCACGCCGTTGTCCGTGTTATAGGGCAGTTCATTAATAAATTCGAAAGTTATGCCCTTTTCTaa                                                                  <  2:797224/75‑1 (MQ=255)

                                                                                          gTTAAATTGCGCCACGCCGTTGTCCGTGTTATAGGGCAGTTCATTAATAAATTCGAAAGTTATGCCCTTTTCTAACAACAGAATAGAAAGTTTGCGTACAAACGGGCTGGTGTAGCTACCGACGAGTTTCATGCCGAGt  >  2:662601/1‑139 (MQ=255)

ATTCATTAATTCAATATATTCAGCGATGATCGGCGAATCAAACCAGCATTCGCCCTCTTCGGTCACCAGCACCGGCACTTTTCCTAACGGGTTAAATTGCGCCACGCCGTTGTCCGCGTTATAGGGCAGTTCATTAATAAATTCGAAAGTTATGCCCTTTTCTAACAACAGAATAGAAAGTTTGCGTACAAACGGGCTGGTGTAGCTACCGACGAGTTTCATGCCGAGT  >  NC_000913/3761734‑3761962

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: `‑`