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%	96.6 / NA	31	*186* (TAG→TAA)	yaiS	putative PIG‑L family deacetylase
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (20/11);  total (20/11)

CCAGGTTTTCATGCAATTTGCTCAAAGAATCATTTTATGAATTACAAAGCCTTCACCCAGATCGCTATTGACCTGTTGTCCGCGAAATTGTGCAACTGCACGCAGGCGATCATGACGCATATAATCGCGTCGTTCCTAGCTTTTATGTTTTTTTAATGCAGCAAGTTTGACTGTGAAATATTCTTCTTTAACGGATTCAAAAACCTGAGGCATAAATGAAAGCCAGGTACTTGGGGTTTCGTAGCCGAGAATTTGCGGAATAGTTCGGCAG  >  NC_000913/383924‑384194                                                                                                                                        |                                                                                                                                        ccAGGTTTTCATGCAATTTGCTCAAAGAATCATTTTATGAATTACAAAGCCTTCACCCAGATCGCTATTGACCTGTTGTCCGCGAAATTGTGCAACTGCACGCAGGCGATCATGACGCATATAATCGCGTCGTTCTTAg                                                                                                                                      >  2:641064/1‑139 (MQ=255)   aGGTTTTCATGCAATTTGCTCAAAGAATCATTTTATGAATTACAAAGCCTTCACCCAGATCGCTATTGACCTGTTGTCCGCGAAATTGTGCAACTGCACGCAGGCGATCATGACGCATATAATCGCGTCGTTCTTAGCt                                                                                                                                    >  2:192430/1‑139 (MQ=255)          cATGCAATTTGCTCAAAGAATCATTTTATGAATTACAAAGCCTTCACCCAGATCGCTATTGACCTGTTGTCCGCGAAATTGTGCAACTGCACGCAGGCGATCATGACGCATATAATCGCGTCGTTCTTAGCTTTTATGt                                                                                                                             >  2:271801/1‑139 (MQ=255)                  ttGCTCAAAGAATCATTTTATGAATTACAAAGCCTTCACCCAGATCGCTATTGACCTGTTGTCCGCGAAATTGTGCAACTGCACGCAGGCGATCATGACGCATATAATCGCGTCGTTCTTAGCTTTTATGTTTTTTTaa                                                                                                                     <  2:90524/139‑1 (MQ=255)                    gCTCAAAGAATCATTTTATGAATTACAAAGCCTTCACCCAGATCGCTATTGACCTGTTGTCCGCGAAATTGTGCAACTGCACGCAGGCGATCATGACGCATATAATCGCGTCGTTCTTAGCTTTTATGTTTTTTTAATg                                                                                                                   <  1:192430/139‑1 (MQ=255)                      tCAAAGAATCATTTTATGAATTACAAAGCCTTCACCCAGATCGCTATTGACCTGTTGTCCGCGAAATTGTGCAACTGCACGCAGGCGATCATGACGCATATAATCGCGTCGTTCTTAGCTTTTATGTTTTTTTAATgca                                                                                                                 <  1:271801/139‑1 (MQ=255)                                  tttATGAATTACAAAGCCTTCACCCAGATCGCTATTGACCTGTTGTCCGCGAAATTGTGCAACTGCACGCAGGCGATCATGACGCATATAATCGCGTCGTTCTTAGCTTTTATGTTTTTTTAATGCAGCAAGTTTGACt                                                                                                     <  2:416892/139‑1 (MQ=255)                                        aaTTACAAAGCCTTCACCCAGATCGCTATTGACCTGTTGTCCGCGAAATTGTGCAACTGCACGCAGGCGATCATGACGCATATAATCGCGTCGTTCTTAGCTTTTATGTTTTTTTAATgcagc                                                                                                               >  1:615897/1‑123 (MQ=255)                                        aaTTACAAAGCCTTCACCCAGATCGCTATTGACCTGTTGTCCGCGAAATTGTGCAACTGCACGCAGGCGATCATGACGCATATAATCGCGTCGTTCTTAGCTTTTATGTTTTTTTAATgcagc                                                                                                               <  2:615897/123‑1 (MQ=255)                                            aCAAAGCCTTCACCCAGATCGCTATTGACCTGTTGTCCGCGAAATTGTGCAACTGCACGCAGGCGATCATGACGCATATAATCGCGTCGTTCTTAGCTTTTATGTTTTTTTAATGCAGCAAGTTTGACTGTGAAATAtt                                                                                           <  2:678897/139‑1 (MQ=255)                                                  ccTTCACCCAGATCGCTATTGACCTGTTGTCCGCGAAATTGTGCAACTGCACGCAGGCGATCATGACGCATATAATCGCGTCGTTCTTAGCTTTTATg                                                                                                                              <  1:294745/98‑1 (MQ=255)                                                  ccTTCACCCAGATCGCTATTGACCTGTTGTCCGCGAAATTGTGCAACTGCACGCAGGCGATCATGACGCATATAATCGCGTCGTTCTTAGCTTTTATg                                                                                                                              >  2:294745/1‑98 (MQ=255)                                                     tCACCCAGATCGCTATTGACCTGTTGTCCGCGAAATTGTGCAACTGCACGCAGGCGATCATGACGCATATAATCGCGTCGTTCTTAGCTTTTATGTTTTTTTAATGCAGCAAGTTTGACTGTGAAATAttc                                                                                          >  1:443673/1‑131 (MQ=255)                                                     tCACCCAGATCGCTATTGACCTGTTGTCCGCGAAATTGTGCAACTGCACGCAGGCGATCATGACGCATATAATCGCGTCGTTCTTAGCTTTTATGTTTTTTTAATGCAGCAAGTTTGACTGTGAAATAttc                                                                                          <  2:443673/131‑1 (MQ=255)                                                      cACCCAGATCGCTATTGACCTGTTGTCCGCGAAATTGTGCAACTGCACGCAGGCGATCATGACGCATATAATCGCGTCGTTCTTAGCTTTTATGTTTTTTTAATGCAGCAAGTTTGACTGTGAAATATTCTTCTTTAAc                                                                                 >  2:582141/1‑139 (MQ=255)                                                            gATCGCTATTGACCTGTTGTCCGCGAAATTGTGCAACTGCACGCAGGCGATCATGACGCATATAATCGCGTCGTTCTTAGCTTTTATGTTTTTTTAATGCAGCAAGTTTGACTGTGAAATATTCTTCTTTAACGGATTc                                                                           >  1:632355/1‑139 (MQ=255)                                                            gATCGCTATTGACCTGTTGTCCGCGAAATTGTGCAACTGCACGCAGGCGATCATGACGCATATAATCGCGTCGTTCTTAGCTTTTATGTTTTTTTAATGCAGCAAGTTTGACTGTGAAATATTCTTCTTTAACGGATTc                                                                           >  2:10073/1‑139 (MQ=255)                                                                       aCCTGTTGTCCGCGAAATTGTGCAACTGCACGCAGGCGATCATGACGCATATAATCGCGTCGTTCTTAGCTTTTATGTTTTTTTAATGCAGCAAGTTTGACTGTGAAATAttc                                                                                          >  1:398967/1‑113 (MQ=255)                                                                       aCCTGTTGTCCGCGAAATTGTGCAACTGCACGCAGGCGATCATGACGCATATAATCGCGTCGTTCTTAGCTTTTATGTTTTTTTAATGCAGCAAGTTTGACTGTGAAATAttc                                                                                          <  2:398967/113‑1 (MQ=255)                                                                            ttgtCCGCGAAATTGTGCAACTGCACGCAGGCGATCATGACGCATATAATCGCGTCGTTCTTAGCTTTTATGTTTTTTTAATGCAGCAAGTTTGACTGTGAAATATTCTTCTTTAACgg                                                                               <  1:285404/119‑1 (MQ=255)                                                                            ttgtCCGCGAAATTGTGCAACTGCACGCAGGCGATCATGACGCATATAATCGCGTCGTTCTTAGCTTTTATGTTTTTTTAATGCAGCAAGTTTGACTGTGAAATATTCTTCTTTAACgg                                                                               >  2:285404/1‑119 (MQ=255)                                                                            ttgtCCGCGAAATTGTGCAACTGCACGCAGGCGATCATGACGCATATAATCGCGTCGTTCTTAGCTTTTATGTTTTTTTAATGCAGCAAGTTTGACTGTGAAATATTCTTCTTTAACGGATTCAAAAACCTGAGGCATa                                                           >  2:627134/1‑139 (MQ=255)                                                                              gtCCGCGAAATTGTGCAACTGCACGCAGGCGATCATGACGCATATAATCGCGTCGTTCTTAGCTTTTATGTTTTTTTAATGCAGCAAGTTTGACTGTGAAATATTCTTCTTTAACGGATTCAAAAACCTGAGGCATaaa                                                         >  1:737873/1‑139 (MQ=255)                                                                                           tgCAACTGCACGCAGGCGATCATGACGCATATAATCGCGTCGTTCTTAGCTTTTATGTTTTTTTAATGCAGCAAGTTTGACTGTGAAATATTCTTCTTTaa                                                                                  >  2:160055/1‑101 (MQ=255)                                                                                           tgCAACTGCACGCAGGCGATCATGACGCATATAATCGCGTCGTTCTTAGCTTTTATGTTTTTTTAATGCAGCAAGTTTGACTGTGAAATATTCTTCTTTaa                                                                                  <  1:160055/101‑1 (MQ=255)                                                                                                   cACGCAGGCGATCATGACGCATATAATCGCGTCGTTCTTAGCTTTTATGTTTTTTTAATGCAGCAAGTTTGACTGTGAAATATTCTTCTTTAACGGATTCAAAAACCTGAGGCATAAATGAAAGCCAGGTACTTGGGGt                                    >  1:709295/1‑139 (MQ=255)                                                                                                     cGCAGGCGATCATGACGCATATAATCGCGTCGTTCTTAGCTTTTATGTTTTTTTAATGCAGCAAGTTTGACTGTGAAATATTCTTCTTTAACGGATTCAAAAACCTGAGGCATAAATGAAAGCCAGGTACTTGGGGttt                                  >  2:75503/1‑139 (MQ=255)                                                                                                        aGGCGATCATGACGCATATAATCGCGTCGTTCTTAGCTTTTATGTTTTTTTAATGCAGCAAGTTTGACTGTGAAATATTCTTCTTTAACGGATTCAAAAACCTGAGGCATAAATGAAAGCCAGGTACTTGGGGTTTCGt                               >  1:29201/1‑139 (MQ=255)                                                                                                                            aTCGCGTCGTTCTTAGCTTTTATGTTTTTTTAATGCAGCAAGTTTGACTGTGAAATATTCTTCTTTAACGGATTCAAAAACCTGAGGCATAAATGAAAGCCAGGTACTTGGGGTTTCGTAGCCGAGAATTTGCGGAATa           >  1:633794/1‑139 (MQ=255)                                                                                                                                cgtcgtTCTTAGCTTTTATGTTTTTTTAATGCAGCAAGTTTGACTGTGAAATATTCTTCTTTAACGGATTCAAAAACCTGAGGCATAAATGAAAGCCAGGTACTTGGGGTTTCGTAGCCGAGAATTTGCGGAATAGTTc       >  1:444468/1‑139 (MQ=255)                                                                                                                                     ttCTTAGCTTTTATGTTTTTTTAATGCAGCAAGTTTGACTGTGAAATATTCTTCTTTAACGGATTCAAAAACCTGAGGCATAAATGAAAGCCAGGTACTTGGGGTTTCGTAGCCGAGAATTTGCGGAATAGTTCggcag  >  2:266426/1‑139 (MQ=255)                                                                                                                                        |                                                                                                                                        CCAGGTTTTCATGCAATTTGCTCAAAGAATCATTTTATGAATTACAAAGCCTTCACCCAGATCGCTATTGACCTGTTGTCCGCGAAATTGTGCAACTGCACGCAGGCGATCATGACGCATATAATCGCGTCGTTCCTAGCTTTTATGTTTTTTTAATGCAGCAAGTTTGACTGTGAAATATTCTTCTTTAACGGATTCAAAAACCTGAGGCATAAATGAAAGCCAGGTACTTGGGGTTTCGTAGCCGAGAATTTGCGGAATAGTTCGGCAG  >  NC_000913/383924‑384194

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