Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	384,059	C→T	*186* (TAG→TAA)	yaiS ←	putative PIG‑L family deacetylase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	384,059	0	C	T	100.0%	150.7 / NA	47	*186* (TAG→TAA)	yaiS	putative PIG‑L family deacetylase
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (24/23);  total (24/23)

CCAGGTTTTCATGCAATTTGCTCAAAGAATCATTTTATGAATTACAAAGCCTTCACCCAGATCGCTATTGACCTGTTGTCCGCGAAATTGTGCAACTGCACGCAGGCGATCATGACGCATATAATCGCGTCGTTCCTAGCTTTTATGTTTTTTTAATGCAGCAAGTTTGACTGTGAAATATTCTTCTTTAACGGATTCAAAAACCTGAGGCATAAATGAAAGCCAGGTACTTGGGGTTTCGTAGCCGAGAATTTGCGGAATAGTTCGG  >  NC_000913/383924‑384191                                                                                                                                        |                                                                                                                                     ccAGGTTTTCATGCAATTTGCTCAAAGAATCATTTTATGAATTACAAAGCCTTCACCCAGATCGCTATTGACCTGTTGTCCGCGAAATTGTGCAACTGCACGCAGGCGATCATGACGCATATAATCGCGTCGTTCTTAg                                                                                                                                   >  1:709853/1‑139 (MQ=255)                        aaaGAATCATTTTATGAATTACAAAGCCTTCACCCAGATCGCTATTGACCTGTTGTCCGCGAAATTGTGCAACTGCACGCAGGCGATCATGACGCATATAATCGCGTCGTTCTTAGCTTTTATGTTTTTTTAATgcagc                                                                                                            <  2:388969/139‑1 (MQ=255)                        aaaGAATCATTTTATGAATTACAAAGCCTTCACCCAGATCGCTATTGACCTGTTGTCCGCGAAATTGTGCAACTGCACGCAGGCGATCATGACGCATATAATCGCGTCGTTCTTAGCTTTTATGTTTTTTTAATgcagc                                                                                                            <  2:824470/139‑1 (MQ=255)                                     aTGAATTACAAAGCCTTCACCCAGATCGCTATTGACCTGTTGTCCGCGAAATTGTGCAACTGCACGCAGGCGATCATGACGCATATAATCGCGTCGTTCTTAGCTTTTATGTTTTTTTAATGCAGCAAGTTTGACtgtg                                                                                               <  2:709853/139‑1 (MQ=255)                                                                   aTTGACCTGTTGTCCGCGAAATTGTGCAACTGCACGCAGGCGATCATGACGCATATAATCGCGTCGTTCTTAGCTTTTATGTTTTTTTAATg                                                                                                                <  1:34900/92‑1 (MQ=255)                                                                   aTTGACCTGTTGTCCGCGAAATTGTGCAACTGCACGCAGGCGATCATGACGCATATAATCGCGTCGTTCTTAGCTTTTATGTTTTTTTAATg                                                                                                                >  2:34900/1‑92 (MQ=255)                                                                       aCCTGTTGTCCGCGAAATTGTGCAACTGCACGCAGGCGATCATGACGCATATAATCGCGTCGTTCTTAGCTTTTATGTTTTTTTAATGCAGCAAGttt                                                                                                      <  2:528158/98‑1 (MQ=255)                                                                       aCCTGTTGTCCGCGAAATTGTGCAACTGCACGCAGGCGATCATGACGCATATAATCGCGTCGTTCTTAGCTTTTATGTTTTTTTAATGCAGCAAGttt                                                                                                      >  1:528158/1‑98 (MQ=255)                                                                            ttgtCCGCGAAATTGTGCAACTGCACGCAGGCGATCATGACGCATATAATCGCGTCGTTCTTAGCTTTTATGt                                                                                                                          >  1:477198/1‑73 (MQ=255)                                                                            ttgtCCGCGAAATTGTGCAACTGCACGCAGGCGATCATGACGCATATAATCGCGTCGTTCTTAGCTTTTATGt                                                                                                                          <  2:477198/73‑1 (MQ=255)                                                                            ttgtCCGCGAAATTGTGCAACTGCACGCAGGCGATCATGACGCATATAATCGCGTCGTTCTTAGCTTTTATGTTTTTTTAATGCAGCAAGTTTGACtgag                                                                                               >  1:254152/1‑98 (MQ=255)                                                                            ttgtCCGCGAAATTGTGCAACTGCACGCAGGCGATCATGACGCATATAATCGCGTCGTTCTTAGCTTTTATGTTTTTTTAATGCAGCAAGTTTGACtgag                                                                                               <  2:254152/100‑3 (MQ=255)                                                                            ttgtCCGCGAAATTGTGCAACTGCACGCAGGCGATCATGACGCATATAATCGCGTCGTTCTTAGCTTTTATGTTTTTTTAATGCAGCAAGTTTGACTGTGAAATAttc                                                                                       <  2:183364/108‑1 (MQ=255)                                                                            ttgtCCGCGAAATTGTGCAACTGCACGCAGGCGATCATGACGCATATAATCGCGTCGTTCTTAGCTTTTATGTTTTTTTAATGCAGCAAGTTTGACTGTGAAATAttc                                                                                       >  1:183364/1‑108 (MQ=255)                                                                            ttgtCCGCGAAATTGTGCAACTGCACGCAGGCGATCATGACGCATATAATCGCGTCGTTCTTAGCTTTTATGTTTTTTTAATGCAGCAAGTTTGACTGTGAAATATTCTTCTTTAACGGATTCAAAAACCTGAGGCATa                                                        >  1:935129/1‑139 (MQ=255)                                                                                 cgcgAAATTGTGCAACTGCACGCAGGCGATCATGACGCATATAATCGCGTCGTTCTTAGCTTTTATGTTTTTTTAATGCAGCAAGTTTGACTGTGAAATATTCTTCTTTAACGGATTCAAAAACCTGAGGCATAAATGa                                                   <  2:935129/139‑1 (MQ=255)                                                                                      aaTTGTGCAACTGCACGCAGGCGATCATGACGCATATAATCGCGTCGTTCTTAGCTTTTATGTTTTTTTAATGCAGCAAGTTTGACTGTGAAATATTCTTCTTTAACGGATTCAAAAACCTGAGGCATAAATGAAAGcc                                              <  1:1315735/139‑1 (MQ=255)                                                                                           tgCAACTGCACGCAGGCGATCATGACGCATATAATCGCGTCGTTCTTAGCTTTTATGTTTTTTTaa                                                                                                                  >  2:332918/1‑66 (MQ=255)                                                                                           tgCAACTGCACGCAGGCGATCATGACGCATATAATCGCGTCGTTCTTAGCTTTTATGTTTTTTTaa                                                                                                                  <  1:332918/66‑1 (MQ=255)                                                                                             cAACTGCACGCAGGCGATCATGACGCATATAATCGCGTCGTTCTTAGCTTTTATGttttttt                                                                                                                    <  1:753686/62‑1 (MQ=255)                                                                                             cAACTGCACGCAGGCGATCATGACGCATATAATCGCGTCGTTCTTAGCTTTTATGttttttt                                                                                                                    >  2:753686/1‑62 (MQ=255)                                                                                             cAACTGCACGCAGGCGATCATGACGCATATAATCGCGTCGTTCTTAGCTTTTATGTTTTTTTAATGCAGCAAGTTTGACTGTGAAatat                                                                                         <  1:492236/89‑1 (MQ=255)                                                                                             cAACTGCACGCAGGCGATCATGACGCATATAATCGCGTCGTTCTTAGCTTTTATGTTTTTTTAATGCAGCAAGTTTGACTGTGAAatat                                                                                         >  2:492236/1‑89 (MQ=255)                                                                                                     cGCAGGCGATCATGACGCATATAATCGCGTCGTTCTTAGCTTTTATGTTTTTTTAATGCAGCAAGTTTGACTGTGAAATATTCTTCTTTAACGGATTCAAAAACCTGAGGCATAAATGAAAGCCAGGTACTTGGGGttt                               >  1:246460/1‑139 (MQ=255)                                                                                                          gCGATCATGACGCATATAATCGCGTCGTTCTTAGCTTTTATGTTTTTTTAATGCAGCAAGTTTGACTGTGAAATATTCTTCTTTAACGGATTCAAAAACCTGAGGCATAAATGAAAGCCAGGTACTTGGGGTTTCGTAg                          >  2:549140/1‑139 (MQ=255)                                                                                                                    cGCATATAATCGCGTCGTTCTTAGCTTTTATGTTTTTTTAATGCAGCAAGTTTGACTGTGAAATATTCTTCTTTAACGGATTCAAAAACCTGAGGCATAAATGAAAGCCAGGTACt                                       >  2:578526/1‑116 (MQ=255)                                                                                                                    cGCATATAATCGCGTCGTTCTTAGCTTTTATGTTTTTTTAATGCAGCAAGTTTGACTGTGAAATATTCTTCTTTAACGGATTCAAAAACCTGAGGCATAAATGAAAGCCAGGTACt                                       <  1:578526/116‑1 (MQ=255)                                                                                                                       atatAATCGCGTCGTTCTTAGCTTTTATGTTTTTTTAATGCAGCAAGTTTGACTGTGAAATATTCTTCTTTAACGGATTCa                                                                       >  2:596613/1‑81 (MQ=255)                                                                                                                       atatAATCGCGTCGTTCTTAGCTTTTATGTTTTTTTAATGCAGCAAGTTTGACTGTGAAATATTCTTCTTTAACGGATTCa                                                                       <  1:596613/81‑1 (MQ=255)                                                                                                                        tataATCGCGTCGTTCTTAGCTTTTATGTTTTTTTAATGCAGCAAGTTTGACtgtg                                                                                               <  2:246553/56‑1 (MQ=255)                                                                                                                        tataATCGCGTCGTTCTTAGCTTTTATGTTTTTTTAATGCAGCAAGTTTGACtgtg                                                                                               >  1:246553/1‑56 (MQ=255)                                                                                                                        tataATCGCGTCGTTCTTAGCTTTTATGTTTTTTTAATGCAGCAAGTTTGACt                                                                                                  <  2:1058695/53‑1 (MQ=255)                                                                                                                        tataATCGCGTCGTTCTTAGCTTTTATGTTTTTTTAATGCAGCAAGTTTGACt                                                                                                  >  1:1058695/1‑53 (MQ=255)                                                                                                                        tataATCGCGTCGTTCTTAGCTTTTATGTTTTTTTAATGCAGCAAGTTTGACTGTGaa                                                                                             <  1:651854/58‑1 (MQ=255)                                                                                                                        tataATCGCGTCGTTCTTAGCTTTTATGTTTTTTTAATGCAGCAAGTTTGACTGTGaa                                                                                             >  2:651854/1‑58 (MQ=255)                                                                                                                             tCGCGTCGTTCTTAGCTTTTATGTTTTTTTAATGCAGCAAGTTTGACTGTGAAATATTCTTCTTTAACGGATTCAAAAACCt                                                                >  1:675696/1‑82 (MQ=255)                                                                                                                             tCGCGTCGTTCTTAGCTTTTATGTTTTTTTAATGCAGCAAGTTTGACTGTGAAATATTCTTCTTTAACGGATTCAAAAACCt                                                                <  2:675696/82‑1 (MQ=255)                                                                                                                             tCGCGTCGTTCTTAGCTTTTATGTTTTTTTAATGCAGCAAGTTTGACTGTGAAATATTCTTCTTTAACGGATTCAAAAACCTGAGGCATAAATGAAAGCCAGGTACTTGGGGTTTCGTAGCCGAGAATTTGCGGAATAg       >  2:705860/1‑139 (MQ=255)                                                                                                                                cgtcgtTCTTAGCTTTTATGTTTTTTTAATgcagc                                                                                                            >  1:77230/1‑35 (MQ=255)                                                                                                                                cgtcgtTCTTAGCTTTTATGTTTTTTTAATgcagc                                                                                                            <  2:77230/35‑1 (MQ=255)                                                                                                                                cgtcgtTCTTAGCTTTTATGTTTTTTTAATGCAGCAAGTTTGACTGTGAAatat                                                                                         >  2:618358/1‑54 (MQ=255)                                                                                                                                cgtcgtTCTTAGCTTTTATGTTTTTTTAATGCAGCAAGTTTGACTGTGAAatat                                                                                         <  1:618358/54‑1 (MQ=255)                                                                                                                                cgtcgtTCTTAGCTTTTATGTTTTTTTAATGCAGCAAGTTTGACTGTGAAATATTCTTCTTTAACGGATTCAAAAACCTGAGGCATAAATGAAAGCCAg                                            <  2:1245161/99‑1 (MQ=255)                                                                                                                                cgtcgtTCTTAGCTTTTATGTTTTTTTAATGCAGCAAGTTTGACTGTGAAATATTCTTCTTTAACGGATTCAAAAACCTGAGGCATAAATGAAAGCCAg                                            >  1:1245161/1‑99 (MQ=255)                                                                                                                                 gtcgtTCTTAGCTTTTATGTTTTTTTAATGCAGCAAGTTTGACTGTGAAATATTCTTCTTTAACGGATTCAAAAACCTGAGGCATAAATGAAAg                                                <  2:281678/94‑1 (MQ=255)                                                                                                                                 gtcgtTCTTAGCTTTTATGTTTTTTTAATGCAGCAAGTTTGACTGTGAAATATTCTTCTTTAACGGATTCAAAAACCTGAGGCATAAATGAAAg                                                >  1:281678/1‑94 (MQ=255)                                                                                                                                  tcgtTCTTAGCTTTTATGTTTTTTTAATGCAGCAAGTTTGACTGTGAAATATTCTTCTTTAACGGATTCAAAAACCTGAGGCATAAATGAAAGCCAGGTACTTGGGGTTTCGTAGCCGAGAATTTGCGGAATAGTTCgg  >  1:312396/1‑139 (MQ=255)                                                                                                                                        |                                                                                                                                     CCAGGTTTTCATGCAATTTGCTCAAAGAATCATTTTATGAATTACAAAGCCTTCACCCAGATCGCTATTGACCTGTTGTCCGCGAAATTGTGCAACTGCACGCAGGCGATCATGACGCATATAATCGCGTCGTTCCTAGCTTTTATGTTTTTTTAATGCAGCAAGTTTGACTGTGAAATATTCTTCTTTAACGGATTCAAAAACCTGAGGCATAAATGAAAGCCAGGTACTTGGGGTTTCGTAGCCGAGAATTTGCGGAATAGTTCGG  >  NC_000913/383924‑384191

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