Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,302,238	C→A	L252I (CTC→ATC)	yhbU →	putative peptidase YhbU

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,302,238	0	C	A	100.0%	163.7 / NA	45	L252I (CTC→ATC)	yhbU	putative peptidase YhbU
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base A (15/30);  total (15/30)

GGTTATCCGACGCTATGTAAAGGGCGTTATCTGGTGGACGGCGAGCGCTATCACGCGCTGGAAGAACCAACCAGTCTCAATACCCTGGAACTGCTGCCGGAGTTAATGGCGGCGAATATTGCTTCGGTGAAAATTGAAGGCCGCCAGCGTAG  >  NC_000913/3302163‑3302314                                                                            |                                                                             ggTTATCCGACGCTATGTAAAGGGCGTTATCTGGTGGACGGCGAGCGCTATCACGCGCTGGAAGAACCAACCAGTATCaa                                                                          >  1:556986/1‑80 (MQ=255)  gTTATCCGACGCTATGTAAAGGGCGTTATCTGGTGGACGGCGAGCGCTATCACGCGCTGGAAGAACCAACCAGTATCAATa                                                                        >  1:556988/1‑81 (MQ=255)     aTCCGACGCTATGTAAAGGGCGTTATCTGGTGGACGGCGAGCGCTATCACGCGCTGGAAGAACCAACCAGTATCAATAccc                                                                     >  2:556981/1‑81 (MQ=255)     aTCCGACGCTATGTAAAGGGCGTTATCTGGTGGACGGCGAGCGCTATCACGCGCTGGAAGAACCAACCAGTATCAATAccc                                                                     >  2:556978/1‑81 (MQ=255)      tCCGACGCTATGTAAAGGGCGTTATCTGGTGAACGGCGAGCGCTATCACGCGCTGGAAGAACCAACCAGTATCAATACCCt                                                                    <  2:556920/81‑1 (MQ=255)        cGACGCTATGTAAAGGGCGTTATCTGGTGGACGGCGAGCGCTATCACGCGCTGGAAGAACCAACCAGTATCAATACCCTgg                                                                  <  1:556923/81‑1 (MQ=255)        cGACGCTATGTAAAGGGCGTTATCTGGTGGACGGCGAGCGCTATCACGCGCTGGAAGAACCAACCAGTATCAATACCCTgg                                                                  <  1:556922/81‑1 (MQ=255)         gACGCTATGTAAAGGGCGTTATCTGGTGGACGGCGAGCGCTATCACGCGCTGGAAGAACCAACCAGTATCAATACCCTGGa                                                                 <  2:556925/81‑1 (MQ=255)         gACGCTATGTAAAGGGCGTTATCTGGTGGACGGCGAGCGCTATCACGCGCTGGAAGAACCAACCAGTATCAATACCCTGGa                                                                 <  1:556925/81‑1 (MQ=255)         gACGCTATGTAAAGGGCGTTATCTGGTGGACAGCGAGCGCTATCACGCGCTGGAAGAACCAACCAGTATCAATACCCTGGa                                                                 <  2:556923/81‑1 (MQ=255)                tGTAAAGGGCGTTATCTGGTGGACGGCGAGCGCTATCACGCGCTGGAAGAACCAACCAGTATCAATACCCTGGAActgct                                                           <  2:556927/80‑1 (MQ=255)                tGTAAAGGGCGTTATCTGGTGGACGGCGAGCGCTATCACGCGCTGGAAGAACCAACCAGTATCAATACCCTGGAActgc                                                            <  2:556926/79‑1 (MQ=255)                 gTAAAGGGCGTTATCTGGTGGACGGCGAGCGCTATCACGCGCTGGAAGAACCAACCAGTATCAATACCCTGGAACtgctgc                                                         >  2:556988/1‑81 (MQ=255)                  tAAAGGGCGTTATCTGGTGGACGGCGAGCGCTATCACGCGCTGGAAGAACCGACCAGTATCAATACCCTGGAACTGCTGcc                                                        <  1:556929/81‑1 (MQ=255)                  tAAAGGGCGTTATCTGGTGGACGGCGAGCGCTATCACGCGCTGGAAGAACCAACCAGTATCAATACCCTGGAACTGCTGcc                                                        <  1:556930/81‑1 (MQ=255)                    aaGGGCGTTATCTGGTGGACGGCGAGCGCTATCACGCGCTGGAAGAACCAACCAGTATCAATACCCTGGAACTGCTGCCgg                                                      <  2:556930/81‑1 (MQ=255)                    aaGGGCGTTATCTGGTGGACGGCGAGCGCTATCACGCGCTGGAAGAACCAACCAGTATCAATACCCTGGAACTGCTGCCgg                                                      <  1:556932/81‑1 (MQ=255)                          gTTATCTGGTGGACGGCGAGCGCTATCACGCGCTGGAAGAACCAACCAGTATCAATACCCTGGAACTGCTGCCGGAGTTaa                                                >  2:556972/1‑81 (MQ=255)                                      aCGGCGAGCGCTATCACGCGCTGGAAGAACCAACCAGTATCAATACCCTGGAACTGCTGCCGGAGTTAATGGCGGCGAata                                    >  2:556984/1‑81 (MQ=255)                                      aCGGCGAGCGCTATCACGCGCTGGAAGAACCAACCAGTATCAATACCCTGGAACTGCTGCCGGAGTTAATGGCGGCGAata                                    <  1:556933/81‑1 (MQ=255)                                       cGGCGAGCGCTATCACGCGCTGGAAGAACCAACCAGTATCAATACCCTGGAACTGCTGCCGGAGTTAATGGCGGCGAatat                                   >  2:556977/1‑81 (MQ=255)                                       cGGCGAGCGCTATCACGCGCTGGAAGAACCAACCAGTATCAATACCCTGGAACTGCTGCCGGAGTTAATGGCGGCGAatat                                   <  2:556933/81‑1 (MQ=255)                                         gCGAGCGCTATCACGCGCTGGAAGAACCAACCAGTATCAATACCCTGGAACTGCTGCCGGAGTTAATGGCGGCGAATATTg                                 >  2:557002/1‑81 (MQ=255)                                          cGAGCGCTATCACGCGCTGGAAGAACCAACCAGTATCAATACCCTGGAACTGCTGCCGGAGTTAATGGCGGCGAATATTGc                                <  1:556935/81‑1 (MQ=255)                                          cGAGCGCTATCACGCGCTGGAAGAACCAACCAGTATCAATACCCTGGAACTGCTGCCGGAGTTAATGGCGGCGAATATTGc                                <  2:556935/81‑1 (MQ=255)                                            aGCGCTATCACGCGCTGGAAGAACCAACCAGTATCAATACCCTGGAACTGCTGCCGGAGTTAATGGCGGCGAATATTGCtt                              <  1:556937/81‑1 (MQ=255)                                             gcgcTATCACGCGCTGGAAGAACCAACCAGTATCAATACCCTGGAACTGCTGCCGGAGTTAATGGCGGCGAATATTGCTTc                             <  1:556938/81‑1 (MQ=255)                                              cgcTATCACGCGCTGGAAGAACCAACCAGTATCAATACCCTGGAGCTGCTGCCGGAGTTAATAGCGGCGAATATTGCTTCg                            <  2:556938/81‑1 (MQ=255)                                              cgcTATCACGCGCTGGAAGAACCAACCAGTATCAATACCCTGGAACTGCTGCCGGAGTTAATGGCGGCGAATATTGCTTCg                            >  1:556991/1‑81 (MQ=255)                                                cTATCACGCGCTGGAAGAACCAACCAGTATCAATACCCTGGAACTGCTGCCGGAGTTAATGGCGGCGAATATTGCTTCGGt                          <  1:556940/81‑1 (MQ=255)                                                 tATCACGCGCTGGAAGAACCAACCAGTATCAATACTCTGGAACTGCTGCCGGAGTTAATGGCGGCGAATATTGTTTCGGt                          >  2:556994/1‑80 (MQ=255)                                                      cgcgCTGGAAGAACCAACCAGTATCAATACCCTGGAACTGCTGCCGGAGTTAATGGCGGCGAATATTGCTTCGGTGAAAAt                    <  1:556941/81‑1 (MQ=255)                                                       gcgcTGGAAGAACCAACCAGTATCAATACCCTGGAACTGCTGCCGGAGTTAATGGCGGCGAATATTGCTTCGGTGAAAAtt                   >  1:556987/1‑81 (MQ=255)                                                         gcTGGAAGAACCAACCAGTATCAATACCCTGGAACTGCTGCCGGAGTTAATGGCGGCGAATATTGCTTCGGTGAAAATTg                  >  2:556989/1‑80 (MQ=255)                                                           tGGAAGAACCAACCAGTATCAATACCCTGGAACTGCTGCCGGAGTTAATGGCGGCGAATATTGCTTCGGTGAAAATTGAAg               >  1:557002/1‑81 (MQ=255)                                                             gaagaaCCAACCAGTATCAATACCCTGGAACTGCTGCCGGAGTTAATGGCGGCGAATATTGCTTCGGTGAAAATTGAAGgc             <  2:556941/81‑1 (MQ=255)                                                              aagaaCCAACCAGTATCAATACCCTGGAACTGCTGCCGGAGTTAATGGCGGCGAATATTGCTTCGGTGAAAATTGAAGgcc            <  2:556942/81‑1 (MQ=255)                                                               agaaCCAACCAGTATCAATACCCTGGAACTGCTGCCGGAGTTAATGGCGGCGAATATTGCTTCGGTGAAAATTGAAGgccg           <  1:556944/81‑1 (MQ=255)                                                                gaaCCAACCAGTATCAATACCCTGGAACTGCTGCCGGAGTTAATGGCGGCGAATATTGCTTCGGTGAAAATTGAAGgccgc          >  2:556991/1‑81 (MQ=255)                                                                 aaccaaccAGTATCAATACCCTGGAACTGCTGCCGGAGTTAATGGCGGCGAATATTGCTTCGGTGAAAATTGAAGgccgcc         <  2:556946/81‑1 (MQ=255)                                                                 aaccaaccAGTATCAATACCCTGGAACTGCTGCCGGAGTTAATGGCGGCGAATATTGCTTCGGTGAAAATTGAAGgccgcc         <  1:556946/81‑1 (MQ=255)                                                                 aaccaaccAGTATCAATACCCTGGAACTGCTGCCGGAGTTAATGGCGGCGAATATTGCTTCGGTGAAAATTGAAGgccgcc         <  1:556945/81‑1 (MQ=255)                                                                     aaccaGTATCAATACCCTGGAACTGCTGCCGGAGTTAATGGCGGCGAATATTGCTTCGGTGAAAATTGAAGGCCGCCAGCg     <  2:556947/81‑1 (MQ=255)                                                                      accaGTATCAATACCCTGGAACTGCTGCCGGAGTTAATGGCGGCGAATATTGCTTCGGTGAAAATTGAAGGCCGCCAGCGt    <  2:556948/81‑1 (MQ=255)                                                                        caGTATCAATACCCTGGAACTGCTGCCGGAGTTAATGGCGGCGAATATTGCTTCGGTGAAAATTGAAGGCCGCCAGCGTAg  <  2:556949/81‑1 (MQ=255)                                                                            |                                                                             GGTTATCCGACGCTATGTAAAGGGCGTTATCTGGTGGACGGCGAGCGCTATCACGCGCTGGAAGAACCAACCAGTCTCAATACCCTGGAACTGCTGCCGGAGTTAATGGCGGCGAATATTGCTTCGGTGAAAATTGAAGGCCGCCAGCGTAG  >  NC_000913/3302163‑3302314

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