BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	812,947	G→A	226 (TAG→TAA)	bioD →	dethiobiotin synthetase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	812,947	0	G	A	100.0%	123.8 / NA	39	226 (TAG→TAA)	bioD	dethiobiotin synthetase
Reads supporting (aligned to +/- strand): ref base G (0/0); new base A (19/20); total (19/20)

GAAAACGTCACGCTGAATATATGACCACGCTCACCCGCATGATTCCCGCGCCGCTGCTGGGAGAGATCCCCTGGCTTGCAGAAAATCCAGAAAATGCGGCAACCGGAAAGTACATAAACCTTGCCTTGTTGTAGCCATTCTGTATTTGGTTAAATTGCGAGCGAGATCGCGTCTTCGATTGACTGCAATTTAACCAATTAAATTCTAAAATAATCACGAAAAAAATTTTACTTCCGCCTCATGCGGCG  >  NC_000913/812814‑813061

gAAAACGTCACGCTGAATATATGACCACGCTCACCCGCATGATTCCCGCGCCGCTGCTGGGAGAGATCCCCTGGCTTGCAGAAAATCCAGAAAATGCGGCAACCGGAAAGTACATAAACCTTGCCTTGTTGTAACCAtt                                                                                                               <  2:157386/139‑1 (MQ=255)

gAAAACGTCACGCTGAATATATGACCACGCTCACCCGCATGATTCCCGCGCCGCTGCTGGGAGAGATCCCCTGGCTTGCAGAAAATCCAGAAAATGCGGCAACCGGAAAGTACATAAACCTTGCCTTGTTGTAACCAtt                                                                                                               <  2:1027122/139‑1 (MQ=255)

      gTCACGCTGAATATATGACCACGCTCACCCGCATGATTCCCGCGCCGCTGCTGGGAGAGATCCCCTGGCTTACAGAAAATCCAGAAAATGCGGCAACCGGAAAGTACATAAACCTTGCCTTGTTGTAACCATTCTGTAt                                                                                                         >  1:354802/1‑139 (MQ=255)

                         cACGCTCACCCGCATGATTCCCGCGCCGCTGCTGGGAGAGATCCCCTGGCTTGCAGAAAATCCAGAAAATGCGGCAACCGGAAAGTACATAAACCTTGCCTTGTTGTAACCATTCTGTATTTGGTTAAATTgcgagcga                                                                                      >  1:575696/1‑139 (MQ=255)

                           cGCTCACCCGCATGATTCCCGCGCCGCTGCTGGGAGAGATCCCCTGGCTTGCAGAAAATCCAGAAAATGCGGCAACCGGAAAGTACATAAACCTTGCCTTGTTGTAACCATTCTGTATTTGGTTAAATTGCGAGCgaga                                                                                    >  2:7038/1‑139 (MQ=255)

                           cGCTCACCCGCATGATTCCCGCGCCGCTGCTGGGAGAGATCCCCTGGCTTGCAGAAAATCCAGAAAATGCGGCAACCGGAAAGTACATAAACCTTGCCTTGTTGTAACCATTCTGTATTTGGTTAAATTGCGAGCgaga                                                                                    >  2:641091/1‑139 (MQ=255)

                           cGCTCACCCGCATGATTCCCGCGCCGCTGCTGGGAGAGATCCCCTGGCTTGCAGAAAATCCAGAAAATGCGGCAACCGGAAAGTACATAAACCTTGCCTTGTTGTAACCATTCTGTATTTGGTTAAATTGCGAGCgaga                                                                                    >  2:1134497/1‑139 (MQ=255)

                            gCTCACCCGCATGATTCCCGCGCCGCTGCTGGGAGAGATCCCCTGGCTTGCAGAAAATCCAGAAAATGCGGCAACCGGAAAGTACATAAACCTTGCCTTGTTGTAACCATTCTGTATTTGGTTAAATTGCGAGCGAGAt                                                                                   >  1:410754/1‑139 (MQ=255)

                            gCTCACCCGCATGATTCCCGCGCCGCTGCTGGGAGAGATCCCCTGGCTTGCAGAAAATCCAGAAAATGCGGCAACCGGAAAGTACATAAACCTTGCCTTGTTGTAACCATTCTGTATTTGGTTAAATTGCGAGCGAGAt                                                                                   >  1:758678/1‑139 (MQ=255)

                               cACCCGCATGATTCCCGCGCCGCTGCTGGGAGAGATCCCCTGGCTTGCAGAAAATCCAGAAAATGCGGCAACCGGAAAGTACATAAACCTTGCCTTGTTGTAACCATTCTGTATTTGGTTAAATTGCGAGCGAGATcgc                                                                                <  1:72180/139‑1 (MQ=255)

                                     cATGATTCCCGCGCCGCTGCTGGGAGAGATCCCCTGGCTTGCAGAAAATCCAGAAAATGCGGCAACCGGAAAGTACATAAACCTTGTCTTGTTGTAACCATTCTGTATTTGGTTAAATTGCGAGCGAGATCGCGTCTTc                                                                          >  2:945972/1‑139 (MQ=255)

                                        gATTCCCGCGCCGCTGCTGGGAGAGATCCCCTGGCTTGCAGAAAATCCAGAAAATGCGGCAACCGGAAAGTACATAAACCTTGCCTTGTTGTAACCATTCTGTATTTGGTTAAATTGCGAGCGAGATCGCGTCTTCGAt                                                                       <  2:758678/139‑1 (MQ=255)

                                         aTTCCCGCGCCGCTGCTGGGAGAGATCCCCTGGCTTGCAGAAAATCCAGAAAATGCGGCAACCGGAAAGTACATAAACCTTGTCTTGTTGTAACCATTCTGTATTTGGTTAAATTGCGAGCGAGATCGCGTCTTCGAtt                                                                      <  1:945972/139‑1 (MQ=255)

                                         aTTCCCGCGCCGCTGCTGGGAGAGATCCCCTGGCTTGCAGAAAATCCAGAAAATGCGGCAACCGGAAAGTACATAAACCTTGCCTTGTTGTAACCATTCTGTATTTGGTTAAATTGCGAGCGAGATCGCGTCTTCGAtt                                                                      <  1:1112168/139‑1 (MQ=255)

                                            cccGCGCCGCTGCTGGGAGAGATCCCCTGGCTTGCAGAAAATCCAGAAAATGCGGCAACCGGAAAGTACATAAACCTTGCCTTGTTGTAACCATTCTGTATTTGGTTAAATTGCGAGCGAGATCGCGTCTTCGATTGAc                                                                   <  2:457345/139‑1 (MQ=255)

                                                 gccgcTGCTGGGAGAGATCCCCTGGCTTGCAGAAAATCCAGAAAATGCGGCAACCGGAAAGTACATAAACCTTGCCTTGTTGTAACCATTCTGTATTTGGTTAAATTGCGAGCGAGATCGCGTCTTCGATTGACTGCaa                                                              <  1:1251394/139‑1 (MQ=255)

                                                  ccgcTGCTGGGAGAGATCCCCTGGCTTGCAGAAAATCCAGAAAATGCGGCAACCGGAAAGTACATAAACCTTGCCTTGTTGTAACCTTTCTGTATTTGGTTAAATTGCGAGCGAGATCGCGTCTTCGATTGACTGCAAt                                                             >  1:868542/1‑139 (MQ=255)

                                                             agagaTCCCCTGGCTTGCAGAAAATCCAGAAAATGCGGCAACCGGAAAGTACATAAACCTTGCCTTGTTGTAACCATTCTGTATTTGGTTAAATTGCGAGCGAGATCGCGTCTTCGATTGACTGCAATTTAACCAATTa                                                  <  2:410754/139‑1 (MQ=255)

                                                                gaTCCCCTGGCTTGCAGAAAATCCAGAAAATGCGGCAACCGGAAAGTACATAAACCTTGCCTTGTTGTAACCATTCTGTCTTTGGTTAAATTgc                                                                                            >  2:797355/1‑94 (MQ=255)

                                                                gaTCCCCTGGCTTGCAGAAAATCCAGAAAATGCGGCAACCGGAAAGTACATAAACCTTGCCTTGTTGTAACCATTCTGTATTTGGTTAAATTgc                                                                                            <  1:797355/94‑1 (MQ=255)

                                                                  tCCCCTGGCTTGCAGAAAATCCAGAAAATGCGGCAACCGGAAAGTACATAAACCTTGCCTTGTTGTAACCATTCTGTATTTGGTTAAATTGCGAGCGAGATCGCGTCTTCGATTGACTGCAATTTAACCAATTAAATTc                                             >  2:709456/1‑139 (MQ=255)

                                                                     ccTGGCTTGCAGAAAATCCAGAAAATGCGGCAACCGGAAAGTACATAAACCTTGCCTTGTTGTAACCATTCTGTATTTGGTTAAATTgcg                                                                                           >  2:390102/1‑90 (MQ=255)

                                                                     ccTGGCTTGCAGAAAATCCAGAAAATGCGGCAACCGGAAAGTACATAAACCTTGCCTTGTTGTAACCATTCTGTATTTGGTTAAATTgcg                                                                                           <  1:390102/90‑1 (MQ=255)

                                                                           ttGCAGAAAATCCAGAAAATGCGGCAACCGGAAAGTACATAAACCTTGCCTTGTTGTAACCATTCTGTATTTGGTTAAATTGCGAGCGAGATCGCGTCTTCGATTGACTGCAATTTAACCAATTAAATTCTAAaataat                                    >  1:924133/1‑139 (MQ=255)

                                                                             gCAGAAAATCCAGAAAATGCGGCAACCGGAAAGTACATAAACCTTGCCTTGTTGTAACCATTCTGTATTTGGTTAAATTGCGAGCGAGATCGCGTCTTCGATTGACTGCAATTTCACCAATTAAATTCTAAAATAATCa                                  >  1:732693/1‑139 (MQ=255)

                                                                              cAGAAAATCCAGAAAATGCGGCAACCGGAAAGTACATAAACCTTGCCTTGTTGTAACCATTCTGTa                                                                                                          <  1:1291967/66‑1 (MQ=255)

                                                                              cAGAAAATCCAGAAAATGCGGCAACCGGAAAGTACATAAACCTTGCCTTGTTGTAACCATTCTGTa                                                                                                          >  2:1291967/1‑66 (MQ=255)

                                                                              cAGAAAATCCAGAAAATGCGGCAACCGGAAAGTACATAAACCTTGCCTTGTTGTAACCATTCTGTATTTGGTTAAATTGCGAGCGAGATCGCGTCTTCGATTGACTGCAATTTAACCAATTAAATTCTAAAATAATCAc                                 >  1:818655/1‑139 (MQ=255)

                                                                                  aaaTCCAGAAAATGCGGCAACCGGAAAGTACATAAACCTTGCCTTGTTGTAACCATTCTGTATTTGGTTAAATTGCGAGCGAGATCGCGTCTTCGATTGACTGCAATTTAACCAATTAAATTCTAAAATAATCACGaaa                             <  2:818655/139‑1 (MQ=255)

                                                                                   aagccAGAAAATGCGGCAACCGGAAAGTACATAAACCTTGCCTTGTTGTAACCATTCTGTATTTGGTTAAATTGCGAGCGAGATCGCGTCTTCGATTGACTGCAATTTAACCAATTAAATTCTAAAATAATCACGaaaa                            <  2:104434/136‑1 (MQ=255)

                                                                                   aaTCCAGAAAATGCGGCAACCGGAAAGTACATAAACCTTGCCTTGTTGTAACCATTCTGTATTTGGTTAAATTGCGAGCGAGATCGCGTCTTCGATTGACTGCAATTTAACCAATTAAATTCTAAAATAATCACGaaaa                            <  2:924133/139‑1 (MQ=255)

                                                                                   aaTCCAGAAAATGCGGCAACCGGAAAGTACATAAACCTTGCCTTGTTGTAACCATTCTGTATTTGGTTAAATTGCGAGCGAGATCGCGTCTTCGATTGACTGCAATTTAACCAATTAAATTCTAAAATAATCACGaaaa                            <  1:1134497/139‑1 (MQ=255)

                                                                                               gCGGCAACCGGAAAGTACATAAACCTTGCCTTGTTGTAACCATTCTGTa                                                                                                          >  2:654120/1‑49 (MQ=255)

                                                                                               gCGGCAACCGGAAAGTACATAAACCTTGCCTTGTTGTAACCATTCTGTa                                                                                                          <  1:654120/49‑1 (MQ=255)

                                                                                                    aaCCGGAAAGTACATAAACCTTGCCTTGTTGTAACCATTCTGTATTTGGTTAAATTGCGAGCGAGATCGCGTCTTCGATTGACTGCAATTTAACCAATTAAATTCTAAAATAATCACGAAAAAAATTTTACTTCCGCCt           >  2:696237/1‑139 (MQ=255)

                                                                                                    aaCCGGAAAGTACATAAACCTTGCCTTGTTGTAACCATTCTGTATTTGGTTAAATTGCGAGCGAGATCGCGTCTTCGATTGACTGCAATTTAACCAATTAAATTCTAAAATAATCACGAAAAAAATTTTACTTCCGCCt           <  1:641091/139‑1 (MQ=255)

                                                                                                             gTACATAAACCTTGCCTTGTTGTAACCATTCTGTATTTGGTTAAATTGCGAGCGAGATCGCGTCTTCGATTGACTGCAATTTAACCAATTAAATTCTAAAATAATCACGAAAAAAATTTTACTTCCGCCTCATgcggcg  <  1:696237/139‑1 (MQ=255)

                                                                                                                          gccttgTTGTAACCATTCTGTATTTGGTTAAATTGCGAGCGAGATCGCGTCTTCGATTGACTGCAATTTAACCAATTAAATTCTAAAATAATCACGAAAAAAATTTTACTTCCGCCTCATgcggcg  >  1:216273/1‑126 (MQ=255)

                                                                                                                          gccttgTTGTAACCATTCTGTATTTGGTTAAATTGCGAGCGAGATCGCGTCTTCGATTGACTGCAATTTAACCAATTAAATTCTAAAATAATCACGAAAAAAATTTTACTTCCGCCTCATgcggcg  <  2:216273/126‑1 (MQ=255)

GAAAACGTCACGCTGAATATATGACCACGCTCACCCGCATGATTCCCGCGCCGCTGCTGGGAGAGATCCCCTGGCTTGCAGAAAATCCAGAAAATGCGGCAACCGGAAAGTACATAAACCTTGCCTTGTTGTAGCCATTCTGTATTTGGTTAAATTGCGAGCGAGATCGCGTCTTCGATTGACTGCAATTTAACCAATTAAATTCTAAAATAATCACGAAAAAAATTTTACTTCCGCCTCATGCGGCG  >  NC_000913/812814‑813061

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