Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,301,599	G→A	G39S (GGC→AGC)	yhbU →	U32 peptidase family protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,301,599	0	G	A	100.0%	96.0 / NA	31	G39S (GGC→AGC)	yhbU	U32 peptidase family protein
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (14/17);  total (14/17)

CAACCTTCAGGAAAAATTATGGAGCTGCTCTGCCCTGCCGGAAATCTCCCGGCGCTTAAGGCGGCCATCGAAAACGGCGCAGATGCTGTTTATATCGGGCTAAAAGATGATACCAATGCCCGTCACTTCGCCGGCCTTAACTTTACCGAGAAAAAATTGCAGGAAGCGGTGAGTTTTGTCCATCAACATCGCCGCAAACTTCACATCGCGATTAACACTTTTGCGCATCCGGACGGTTACGCCCGTTGGCAGCGCGCCGT  >  NC_000913/3301467‑3301726                                                                                                                                     |                                                                                                                                cAACCTTCAGGAAAAATTATGGAGCTGCTCTGCCCTGCCGGAAATCTCCCGGCGCTTAAGGCGGCCATCGAAAACGGCGCAGATGCTGTTTATATCGGGCTAAAAGATGATACCAATGCCCGTCACTTCGCCAGCCTTa                                                                                                                           <  2:1047654/139‑1 (MQ=255)   aCCTTCAGGAAAAATTATGGAGCTGCTCTGCCCTGCCGGAAATCTCCCGGCGCTTAAGGCGGCCATCGAAAACGGCGCAGATGCTGTTTATATCGGGCTAAAAGATGATACCAATGCCCGTCACTTCGCCAGCCTTAAc                                                                                                                         >  1:109140/1‑139 (MQ=255)              aaaTTATGGAGCTGCTCTGCCCTGCCGGAAATCTCCCGGCGCTTAAGGCGGCCATCGAAAACGGCGCAGATGCTGTTTATATCGGGCTAAAAGATGATACCAATGCCCGTCACTTCGCCAGCCTTAACTTTACCGAGaa                                                                                                              >  1:111073/1‑139 (MQ=255)                         ctgctCTGCCCTGCCGGAAATCTCCCGGCGCTTAAGGCGGCCATCGAAAACGGCGCAGATGCTGTTTATATCGGGCTAAAAGATGATACCAATGCCCGTCACTTCGCCAGCCTTAACTTTACCGAGAAAAAATTGCAgg                                                                                                   >  1:661594/1‑139 (MQ=255)                              ctgccctgccGGAAATCTCCCGGCGCTTAAGGCGGCCATCGAAAACGGCGCAGATGCTGTTTATATCGGGCTAAAAGATGATACCAATGCCCGTCACTTCGCCAGCCTTAACTTTACCGAGAAAAAATTGCAGGAAGCg                                                                                              <  2:111073/139‑1 (MQ=255)                                      ccGGAAATCTCCCGGCGCTTAAGGCGGCCATCGAAAACGGCGCAGATGCTGTTTATATCGGGCTAAAAGATGATACCAATGCCCGTCACTTCGCCAGCCTTAACTTTACCGAGAAAAAATTGCAGGAAGCGGTGAGttt                                                                                      >  1:241914/1‑139 (MQ=255)                                                       cTTAAGGCGGCCATCGAAAACGGCGCAGATGCTGTTTATATCGGGCTAAAAGATGATACCAATGCCCGTCACTTCGCCAGCCTTAACTTTACCGAGAAAAAATTGCAGGAAGCGGTGAGTTTTGTCCATCAACa                                                                          >  1:1115907/1‑134 (MQ=255)                                                       cTTAAGGCGGCCATCGAAAACGGCGCAGATGCTGTTTATATCGGGCTAAAAGATGATACCAATGCCCGTCACTTCGCCAGCCTTAACTTTACCGAGAAAAAATTGCAGGAAGCGGTGAGTTTTGTCCATCAACa                                                                          <  2:1115907/134‑1 (MQ=255)                                                         tAAGGCGGCCATCGAAAACGGCGCAGATGCTGTTTATATCGGGCTAAAAGATGATACCAATGCCCGTCACTTCGCCAGCCTTAACTTTACCGAGAAAAAATTGCAGGAAGCGGTGAGTTTTGTCCATCAACATcgccgc                                                                   <  1:882251/139‑1 (MQ=255)                                                                     cGAAAACGGCGCAGATGCTGTTTATATCGGGCTAAACGATGATCCCAATGCCCGTCACTTCGCCAGCCTTAACTTTACCGAGAAAAAATTGCAGGAAGCGGTGAGTTTTGTCCATCAACATCGCCGCAAACTTCACATc                                                       <  2:241914/139‑1 (MQ=255)                                                                             gcgcAGATGCTGTTTATATCGGGCTAAAAGATGATACCAATGCCCGTCACTTCGCCAGCCTTAACTTTACCGAGAAAAAATTGCAGGAAGCGGTGAGTTTTGTCCATCAACATCGCCGCAAACTTCACATCGCGATTaa                                               <  1:501373/139‑1 (MQ=255)                                                                              cgcAGATGCTGTTTATATCGGGCTAAAAGATGATACCAATGCCCGTCACTTCGCCAGCCTTAACTTTACCGAGAAAAAATTGCAGGAAGCGGTGAGTTTTGTCCATCAACATCGCCGCAAACTTCACATCGCGATTAac                                              >  1:451791/1‑139 (MQ=255)                                                                                        gTTTATATCGGGCTAAAAGATGATACCAATGCCCGTCACTTCGCCAGCCTTAACTTTACCGAGAAAAAATTGCAGGAAGCGGTGAGTTTTGTCCATCAACATCGCCGCAAACTTCACATCGCGATTAACACTTTTgcgc                                    >  2:39649/1‑139 (MQ=255)                                                                                              atCGGGCTAAAAGATGATACCAATGCCCGTCACTTCGCCAGCCTTAACTTTACCGAGAAAAAATTGCAGGAAGCGGTGAGTTTTGTCCATCAACATCGCCGCAAACTTCACATCGCGATTAACACTTTTGCGCATCCgg                              <  2:661594/139‑1 (MQ=255)                                                                                               tCGGGCTAAAAGATGATACCAATGCCCGTCACTTCGCCAGCCTTAACTTTACCGAGAAAAAATTGCAGGAAGCGGTGAGTTTTGTCCATCAACATCGCCGCAAACTTCACATCGCGATTAACACTTTTGCGCATCCGGa                             <  2:645435/139‑1 (MQ=255)                                                                                               tCGGGCTAAAAGATGATACCAATGCCCGTCACTTCGCCAGCCTTAACTTTACCGAGAAAAAATTGCAGGAAGCGGTGAGTTTTGTCCATCAACATCGCCGCAAACTTCACATCGCGATTAACACTTTTGCGCATCCGGa                             <  1:39649/139‑1 (MQ=255)                                                                                                   gCTAAAAGATGATACCAATGCCCGTCACTTCGCCAGCCTTAACTTTACCGAGAAAAAATTGCAGGAAGCGGTGAGTTTTg                                                                                    <  1:801826/80‑1 (MQ=255)                                                                                                   gCTAAAAGATGATACCAATGCCCGTCACTTCGCCAGCCTTAACTGTACCGAGAAAAAATTGCAGGAAGCGGTGAGTTTTg                                                                                    >  2:801826/1‑80 (MQ=255)                                                                                                       aaaGATGATACCAATGCCCGTCACTTCGCCAGCCTTAACTTTACCGAGAAAAAATTGCAGGAAGCGGTGAGTTTTGTCCATCAACATCGCCGCAAACTTCACATCGCGATTAACACTTTTGCGCATCCGGACGGTTACg                     <  2:109140/139‑1 (MQ=255)                                                                                                             gatACCAATGCCCGTCACTTCGCCAGCCTTAACTTTACCGAGAAAAAATTGCAGGAAGCGGTGAGTTTTGTCCATCAACATcgccg                                                                    >  1:304935/1‑86 (MQ=255)                                                                                                             gatACCAATGCCCGTCACTTCGCCAGCCTTAACTTTACCGAGAAAAAATTGCAGGAAGCGGTGAGTTTTGTCCATCAACATcgccg                                                                    <  2:304935/86‑1 (MQ=255)                                                                                                             gatACCAATGCCCGTCACTTCGCCAGCCTTAACTTTACCGAGAAAAAATTGCAGGAAGCGGTGAGTTTTGTCCATCAACATCGCCGCAAACTTCACATCGCGATTAACACTTTTGCGCATCCGGACgg                          <  2:793310/128‑1 (MQ=255)                                                                                                             gatACCAATGCCCGTCACTTCGCCAGCCTTAACTTTACCGAGAAAAAATTGCAGGAAGCGGTGAGTTTTGTCCATCAACATCGCCGCAAACTTCACATCGCGATTAACACTTTTGCGCATCCGGACgg                          >  1:793310/1‑128 (MQ=255)                                                                                                                aCCAATGCCCGTCACTTCGCCAGCCTTAACTTTACCGAGAAAAAATTGCAGGAAGCGGTGAGTTTTGTCCATCAACATCGCCGCAAACTTCACATCGCGATTAACACTTTTGCGCATCCgg                              >  1:260074/1‑121 (MQ=255)                                                                                                                aCCAATGCCCGTCACTTCGCCAGCCTTAACTTTACCGAGAAAAAATTGCAGGAAGCGGTGAGTTTTGTCCATCAACATCGCCGCAAACTTCACATCGCGATTAACACTTTTGCGCATCCgg                              <  2:260074/121‑1 (MQ=255)                                                                                                                    aTGCCCGTCACTTCGCCAGCCTTAACTTTACCGAGAAAAAATTGCAGGAAGCGGTGAGTTTTGTCCATCAACATCGCCGCAAACTTCACATCGCGATTAACACTTTTGCGCATCCGGACGGTTACGCCCGTTGGCAgcg        >  1:122430/1‑139 (MQ=255)                                                                                                                         cGTCACTTCGCCAGCCTTAACTTTACCGAGAAAAAATTGCAGGAAGCGGTGAGTTTTGTCCATCAACATCGCCGCAAACTTCACATc                                                       <  2:807857/87‑1 (MQ=255)                                                                                                                         cGTCACTTCGCCAGCCTTAACTTTACCGAGAAAAAATTGCAGGAAGCGGTGAGTTTTGTCCATCAACATCGCCGCAAACTTCACATc                                                       >  1:807857/1‑87 (MQ=255)                                                                                                                          gTCACTTCGCCAGCCTTAACTTTACCGAGAAAAAATTGCAGGAAGCGGTGAGTTTTGTCCATCAACATCGCCGCAAACTTCACATCGCGATTAACACTTTTGCGCATCCGGACGGTTACGCCCGTTGGCAGCGCGCCGt  <  2:122430/139‑1 (MQ=255)                                                                                                                               ttCGCCAGCCTTAACTTTACCGAGAAAAAATTGCAGGAAGCGGTGAGTTTTGTCCATCAACATCGCCGCAAACTTCACATCGCGATTAACACTTTTGCGCATCCGGACGGTTACg                     <  2:890550/115‑1 (MQ=255)                                                                                                                               ttCGCCAGCCTTAACTTTACCGAGAAAAAATTGCAGGAAGCGGTGAGTTTTGTCCATCAACATCGCCGCAAACTTCACATCGCGATTAACACTTTTGCGCATCCGGACGGTTACg                     >  1:890550/1‑115 (MQ=255)                                                                                                                                     |                                                                                                                                CAACCTTCAGGAAAAATTATGGAGCTGCTCTGCCCTGCCGGAAATCTCCCGGCGCTTAAGGCGGCCATCGAAAACGGCGCAGATGCTGTTTATATCGGGCTAAAAGATGATACCAATGCCCGTCACTTCGCCGGCCTTAACTTTACCGAGAAAAAATTGCAGGAAGCGGTGAGTTTTGTCCATCAACATCGCCGCAAACTTCACATCGCGATTAACACTTTTGCGCATCCGGACGGTTACGCCCGTTGGCAGCGCGCCGT  >  NC_000913/3301467‑3301726

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