Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,605,159	C→T	Q61* (CAG→TAG)	yhhL →	DUF1145 domain‑containing protein YhhL

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,605,159	0	C	T	100.0%	114.7 / NA	32	Q61* (CAG→TAG)	yhhL	DUF1145 domain‑containing protein YhhL
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (15/17);  total (15/17)

GTGGCGCTGATTTTTACCGTGCTGATGCATGGTATGCAGCTGGCGCTATTGAAATCCACTTTACCGAAAGATGGCCCGCAGATGACCACCGCCGAAAAGGTACGGATTTTCCTTTTCGGCGTGTTTGAACTGCTGGCCTGGCAGAAGAAATTTAA  >  NC_000913/3605081‑3605235                                                                               |                                                                             gTGGCGCTGATTTTTACCGTGCTGATGCATGGTATGCAGCTGGCGCTATTGAAATCCACTTTACCGAAAGATGGCCCGTAg                                                                            <  1:630526/81‑1 (MQ=255)       cTGATTTTTACCGTGCTGATGCATGGTATGCAGCTGGCGCTATTGAAATCCACTTTACCGAAAGATGGCCCGTAGATGac                                                                       >  1:630577/1‑80 (MQ=255)          aTTTTTACCGTGCTGATGCATGGTATGCAGCTGGCGCTATTGAAATCCACTTTACCGAAAGATGGCCCGTAGATGaccacc                                                                   <  1:630527/81‑1 (MQ=255)             tttACCGTGCTGATGCATGGTATGCAGCTGGCGCTATTGAAATCCACTTTACCGAAAGATGGCCCGTAGATGACCAccgcc                                                                >  1:630584/1‑81 (MQ=255)              ttACCGTGCTGATGCATGGTATGCAGCTGGCGCTATTGAAATCCACTTTACCGAAAGATGGCCCGTAGATGACCAccgcc                                                                >  1:630565/1‑80 (MQ=255)                aCCGTGCTGATGCATGGTATGCAGCTGGCGCTATTGAAATCCACTTTACCGAAAGATGGCCCGTAGATGACCACCGCCGa                                                              >  2:630579/1‑80 (MQ=255)                   gTGCTGATGCATGGTATGCAGCTGGCGCTATTGAAATCCACTTTACCGAAAGATGGCCCGTAGATGACCACCGCCGAAAAg                                                          >  2:630582/1‑81 (MQ=255)                         aTGCATGGTATGCAGCTGGCGCTATTGAAATCCACTTTACCGAAAGATGGCCCGTAGATGACCACCGCCGAAAAGGTACg                                                     >  2:630575/1‑80 (MQ=255)                          tGCATGGTATGCAGCTGGCGCTATTGAAATCCACTTTACCGAAAGATGGCCCGTAGATGACCACCGCCGAAAAGGTACGGa                                                   <  1:630528/81‑1 (MQ=255)                            cATGGTATGCAGCTGGCGCTATTGAAATCCACTTTACCGAAAGATGGCCCGTAGATGACCACCGCCGAAAAGGTACGGAtt                                                 <  1:630529/81‑1 (MQ=255)                             aTGGTATGCAGCTGGCGCTATTGAAATCCACTTTACCGAAAGATGGCCCGTAGATGACCACCGCCGAAAAGGTACGGAttt                                                >  1:630568/1‑81 (MQ=255)                              tGGTATGCAGCTGGCGCTATTGAAATCCACTTTACCGAAAGATGGCCCGTAGATGACCACCGCCGAAAAGGTACGGAtttt                                               <  2:630529/81‑1 (MQ=255)                              tGGTATGCAGCTGGCGCTATTGAAATCCACTTTACCGAAAGATGGCCCGTAGATGACCACCGCCGAAAAGGTACGGAtttt                                               <  2:630530/81‑1 (MQ=255)                                gTATGCAGCTGGCGCTATTGAAATCCACTTTACCGAAAGATGGCCCGTAGATGACCACCGCCGAAAAGGTACGGATTTTcc                                             >  2:630593/1‑81 (MQ=255)                                 tATGCAGCTGGCGCTATTGAAATCCACTTTACCGAAAGATGGCCCGTAGATGACCACCGCCGAAAAGGTACGGATTTTCCt                                            <  2:630531/81‑1 (MQ=255)                                             gcTATTGAAATCCACTTTACCGAAAGATGGCCCGTAGATGACCACCGCCGAAAAGGTACGGATTTTCCTTTTCGGCGTGtt                                >  1:630578/1‑81 (MQ=255)                                                   gAAATCCACTTTACCGAAAGATGGCCCGTAGATGACCACCGCCGAAAAGGTACGGATTTTCCTTTTCGGCGTGTTTGAAc                           <  1:630533/80‑1 (MQ=255)                                                    aaaTCCACTTTACCGAAAGATGGCCCGTAGATGACCACCGCCGAAAAGGTACGGATTTTCCTTTTCGGCGTGTTTGAActg                         >  1:630588/1‑81 (MQ=255)                                                    aaaTCCACTTTACCGAAAGATGGCCCGTAGATGACCACCGCCGAAAAGGTACGGATTTTCCTTTTCGGCGTGTTTGAAct                          <  1:630534/80‑1 (MQ=255)                                                            tttACCGAAAGATGGCCCGTAGATGACCACCGCCGAAAAGGTACGGATTTTCCTTTTCGGCGTGTTTGAACTGctggcctg                 >  1:630586/1‑81 (MQ=255)                                                            tttACCGAAAGATGGCCCGTAGATGACCACCGCCGAAAAGGTACGGATTTTCCTTTTCGGCGTGTTTGAACTGctggcctg                 >  1:630579/1‑81 (MQ=255)                                                             ttACCGAAAGATGGCCCGTAGATGACCACCGCCGAAAAGGTACGGATTTTCCTTTTCGGCGTGTTTGAACTGctggcctgg                >  2:630590/1‑81 (MQ=255)                                                              tACCGAAAGATGGCCCGTAGATGACCACCGCCGAAAAGGTACGGATTTTCCTTTTCGGCGTGTTTGAACTGctggcctgg                <  2:630535/80‑1 (MQ=255)                                                              tACCGAAAGATGGCCCGTAGATGACCACCGCCGAAAAGGTACGGATTTTCCTTTTCGGCGTGTTTGAACTGctggcctg                 <  1:630535/79‑1 (MQ=255)                                                                ccGAAAGATGGCCCGTAGATGACCACCGCCGAAAAGGTACGGATTTTCCTTTTCGGCGTGTTTGAACTGCTGGCCTGGCa              <  1:630537/80‑1 (MQ=255)                                                                   aaaGATGGCCCGTAGATGACCACCGCCGAAAAGGTACGGATTTTCCTTTTCGGCGTGTTTGAACTGCTGGCCTGGCagaag          <  2:630537/81‑1 (MQ=255)                                                                     aGATGGCCCGTAGATGACCACCGCCGAAAAGGTACGGATTTTCCTTTTCGGCGTGTTTGAACTGCTGGCCTGGCAgaagaa        <  2:630538/81‑1 (MQ=255)                                                                      gATGGCCCGTAGATGACCACCGCCGAAAAGGTACGGATTTTCCTTTTCGGCGTGTTTGAACTGCTGGCCTGGCAGAAGaaa       >  1:630575/1‑81 (MQ=255)                                                                        tGGCCCGTAGATGACCACCGCCGAAAAGGTACGGATTTTCCTTTTCGGCGTGTTTGAACTGCTGGCCTGGCAGAAGAAAtt     <  2:630539/81‑1 (MQ=255)                                                                         ggCCCGTAGATGACCACCGCCGAAAAGGTACGGATTTTCCTTTTCGGCGTGTTTGAACTGCTGGCCTGGCAGAAGAAAtt     >  2:630600/1‑80 (MQ=255)                                                                          gCCCGTAGATGACCACCGCCGAAAAGGTACGGATTTTCCTTTTCGGCGTGTTTGAACTGCTGGCCTGGCAGAAGAAAttt    <  1:630541/80‑1 (MQ=255)                                                                           cccGTAGATGACCACCGCCGAAAAGGTACGGGTTTTCCTTTTCGGCGTGTTTGAACTGCTGGCCTGGCAGAAGAAATTTaa  <  2:630541/81‑1 (MQ=255)                                                                               |                                                                             GTGGCGCTGATTTTTACCGTGCTGATGCATGGTATGCAGCTGGCGCTATTGAAATCCACTTTACCGAAAGATGGCCCGCAGATGACCACCGCCGAAAAGGTACGGATTTTCCTTTTCGGCGTGTTTGAACTGCTGGCCTGGCAGAAGAAATTTAA  >  NC_000913/3605081‑3605235

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