Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,897,866	T→C	N2S (AAC→AGC)	ygcS ←	putative transporter YgcS

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,897,866	0	T	C	100.0%	50.1 / NA	15	N2S (AAC→AGC)	ygcS	putative transporter YgcS
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (8/7);  total (8/7)

TGGCGTAACCAATGACGCCGAGAACATAACCGTCGGTCAGGTGTGCGCCGAAAGTGAGCGCAGCAATGCGGCAGTGAAAACGGTTAAGCGGTAAATCATCCATTCGCACCGGTGAAGTGTTCATCTTTACCCCTTCATTGTCCGGTAATAAGAACGGCTTCTCGCCTGAGAAGCCGAAAAGGCAACGCCCCGTCATCGCGGAACGCTGCCTGACACATTATTTTTCAATCGGATAGA  >  NC_000913/2897747‑2897983                                                                                                                        |                                                                                                                      tGGCGTAACCAATGACGCCGAGAACATAACCGTCGGTCAGGTGTGCGCCGAAAGTGAGCGCAGCAATGCGGCAGTGAAAACGGTTAAGCGGTAAATCATCCATTCGCACCGGTGAAGTGCTCATCTTTAcccc                                                                                                          >  1:611813/1‑133 (MQ=255)      tAACCAATGACGCCGAGAACATAACCGTCGGTCAGGTGTGCGCCGAAAGTGAGCGCAGCAATGCGGCAGTGAAAACGGTTAAGCGGTAAATTATCCATTCCCACCGGTCAAGTGCTCATCTTTACCCCTTCAt                                                                                                     >  2:247373/1‑133 (MQ=255)            aTGACGCCGAGAACATAAACGTCGGTCAGGTGTGCGCCGAAAGTGAGCGCAGCAATGCGGCAGTGAAAACGGTTAAGCGGTAAATCATCCATTCGCACCGGTGAAGTGCTCATCTTTACCCCTTCATTGTCCg                                                                                               <  1:247373/133‑1 (MQ=255)                 gCCGAGAACAGAACCTTCGGTCAGCTGTTCGCAGAAAGTGACCGCAGCAATGCGGCAGTGAAAACGGTTAAGCGGTAAATCATCCATTCGCACCGGTGAAGTGCTCATCTTTACCCCTTCATTGTCCGGtaat                                                                                          <  2:201282/133‑1 (MQ=255)                   cGAGAACATAACCGTCGGTCAGGTGTGCGCCGAAAGTGAGCGCAGCAATGCGGCAGTGAAAACGGTTAAGCGGTAAATCATCCATTCGCACCGGTGAAGTGCTCATCTTTACCCCTTCATTGTCCGGtaataa                                                                                        >  1:546536/1‑133 (MQ=255)                            aaCCGTCGGTCAGGTGTGCGCCGAAAGTGAGCGCAGCAATGCGGCAGTGAAAACGGTTAAGCGGTAAATCATCCATTCGCACCGGTGAAGTGCTCATCTTTACCCCTTCATTGTCCGGTAATAAGAACGGCtt                                                                               <  1:588366/133‑1 (MQ=255)                                                         aGCGCAGCAATGCGGCAGTGAAAACGGTTAAGCGGTAAATCATCCATTCGCACCGGTGAAGTGCTCATCTTTACCCCTTCATTGTCCGGTAATAAGAACGTCTTTTCGGCTGAGAAGGCGAAAAGGCAACGcc                                                  >  2:235278/1‑133 (MQ=255)                                                           cgcAGCAATGCGGCAGTGAAAACGGTTAAGCGGTAAATCATCCATTCGCACCGGTGAAGTGCTCATCTTTACCCCTTCATTGTCCGGTAATAAGAACGGCTTCTCGCCTGAGAAGCCGAAAAGCCAACGcctc                                                >  2:295075/1‑131 (MQ=255)                                                                           tGAAAACGGTTAAGCGGTAAATCATCCATTCGCACCGGTGAAGTGCTCATCTTTACCCCTTCATTGTCCGGTAATAAGAACGGCTTCTCGCCTGAGAAGCCGAAAAGGCAACGCCCCGTCATCGCGGAACGCt                                >  2:521225/1‑133 (MQ=255)                                                                                   gTTAAGCGGTAAATCCTCAATTCGCACCGGTGAAGTGCTCATCTTTACCCCTTCATTGTCCGGTAATAAGAACGGCTTCTCGCCTGAGAAGCCGAAAAGGCAACGCCCCGTCATCGCGGAACGCTGCCTGaca                        <  1:235278/133‑1 (MQ=255)                                                                                               atcatcCATTCGCACCGGTGAAGTGCTCATCTTTACCCCTTCATTGTCCGGTAATAAGAACGGCTTCTCGCCTGAGAAGCCGAAAAGGCAACGCCCCGTCATCGCGGAACGCTGCCTGACACATTATTTTTCa            <  2:546536/133‑1 (MQ=255)                                                                                                      aTTCGCACCGGTGAAGTGCTCATCTTTACCCCTTCATTGTCCGGTAATAAGAACGGCTTCTCGCCTGAGAAGCCGAAAAGGCAACGCCCCGTCATCGCGGAACGCTGCCTGACACATTATTTTTCAATCGgat     >  2:366737/1‑133 (MQ=255)                                                                                                        tCGCACCGGTGAAGTGCTCATCTTTACCCCTTCATTGTCCGGTAATAAGAACGGCTTCTCGCCTGAGAAGCCGAAAAGGCAACGCCCCGTCATCGCGGAACGCTGCCTGACACATTATTTTTCAATCGga      >  1:528972/1‑130 (MQ=255)                                                                                                        tCGCACCGGTGAAGTGCTCATCTTTACCCCTTCATTGTCCGGTAATAAGAACGGCTTCTCGCCTGAGAAGCCGAAAAGGCAACGCCCCGTCATCGCGGAACGCTGCCTGACACATTATTTTTCAATCGga      <  2:528972/130‑1 (MQ=255)                                                                                                         cGCACCGGTGAAGTGCTCATCTTTACCCCTTCATTGTCCGGTAATAAGAACGGCTTCTCGCCTGAGAAGCCGAAAAGGCAACGCCCCGTCATCGCGGAACGCTGCCTGACACATTATTTTTCAATCGgataga  <  2:62992/133‑1 (MQ=255)                                                                                                                        |                                                                                                                      TGGCGTAACCAATGACGCCGAGAACATAACCGTCGGTCAGGTGTGCGCCGAAAGTGAGCGCAGCAATGCGGCAGTGAAAACGGTTAAGCGGTAAATCATCCATTCGCACCGGTGAAGTGTTCATCTTTACCCCTTCATTGTCCGGTAATAAGAACGGCTTCTCGCCTGAGAAGCCGAAAAGGCAACGCCCCGTCATCGCGGAACGCTGCCTGACACATTATTTTTCAATCGGATAGA  >  NC_000913/2897747‑2897983

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