Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,656,514	G→A	G111S (GGC→AGC)	speG →	spermidine N‑acetyltransferase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,656,514	0	G	A	100.0%	112.5 / NA	32	G111S (GGC→AGC)	speG	spermidine N‑acetyltransferase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (11/21);  total (11/21)

TCAGGGGAAAGGTCTGGCAACCCGTGCCGCCAAATTAGCAATGGACTATGGCTTTACCGTTCTCAATCTCTATAAGCTGTATCTGATCGTTGATA  >  NC_000913/1656465‑1656559                                                  |                                              tCAGGGGAAAGGTCTGGCAACCCGTGCCGCCAAATTAGCAATGGACTATg                                               <  1:1994975/50‑1 (MQ=255)   aGGGGAAAGGTCTGGCAACCCGTGCCGCCAAATTAGCAATGGACTATAGc                                             >  1:2073633/1‑50 (MQ=255)    ggggAAAGGTCTGGCAACCCGTGCCGCCAAATTAGCAATGGACTATAGCt                                            <  1:555851/50‑1 (MQ=255)      ggAAAGGTCTGGCAACCCGTGCCGCCAAATTAGCAATGGACTATAGCtt                                           >  1:1303054/1‑49 (MQ=255)         aaGGTCTGGCAACCCGTGCCGCCAAATTAGCAATGGACTATAGCTTTAcc                                       >  1:90192/1‑50 (MQ=255)         aaGGTCTGGCAACCCGTGCCGCCAAATTAGCAATGGACTATAGCTTTAcc                                       <  1:654692/50‑1 (MQ=255)           ggTCTGGCAACCCGTGCCGCCAAATTAGCAATGGACTATAGCTTTACCGt                                     <  1:1824201/50‑1 (MQ=255)           ggTCTGGCAACCCGTACCGCCAAATTAGCAATGGACTATAGCTTTACCGt                                     >  1:1782380/1‑50 (MQ=255)                ggCAACCCGTGCCGCCAAATTAGCAATGGACTATAGCTTTACCGTtctc                                 >  1:254886/1‑49 (MQ=255)                 gCAACCCGTGCCGCCAAATTAGCAATGGACTATAGCTTTACCGTtctc                                 <  1:687036/48‑1 (MQ=255)                 gCAACCCGTGCCGCCAAATTAGCAATGGACTATAGCTTTACCGTTCTCaa                               <  1:649671/50‑1 (MQ=255)                   aaCCCGTGCCGCCAAATTAGCAATGGACTATAGCTTTACCGTTCTCAAtc                             >  1:1157861/1‑50 (MQ=255)                    aCCCGTGCCGCCAAATTAGCAATGGACTATAGCTTTACCGTTCTCAAtct                            <  1:2087348/50‑1 (MQ=255)                      ccGTGCCGCCAAATTAGCAATGGACTATAGCTTTACCGTTCTCAATctct                          <  1:1127024/50‑1 (MQ=255)                       cGTGCCGCCAAATTAGCAATGGACTATAGCTTTACCGTTCTCAATCTCta                         <  1:1613089/50‑1 (MQ=255)                        gTGCCGCCAAATTAGCAATGGACTATAGCTTTACCGTTCTCAATCTCta                         >  1:665982/1‑49 (MQ=255)                         tGCCGCCAAATTAGCAATGGACTATAGCTTTACCGTTCTCAATCTCtata                       <  1:1346532/50‑1 (MQ=255)                         tGCCGCCAAATTAGCAATGGACTATAGCTTTACCGTTCTCAATCTCtata                       <  1:353938/50‑1 (MQ=255)                           ccgccAAATTAGCAATGGACTATAGCTTTACCGTTCTCAATCTCTATaa                      <  1:1730842/49‑1 (MQ=255)                              ccAAATTAGCAATGGACTATAGCTTTACCGTTCTCAATCTCTATAAGCTg                  <  1:1613768/50‑1 (MQ=255)                                 aaTTAGCAATGGACTATAGCTTTACCGTTCTCAATCTCTATAAGCTGTAt               >  1:1167373/1‑50 (MQ=255)                                    tAGCAATGGACTATAGCTTTACCGTTCTCAATCTCTATAAGCTGTATCt             >  1:1331344/1‑49 (MQ=255)                                      gCAATGGACTATAGCTTTACCGTTCTCAATCTCTATAAGCTGTATCTGa           <  1:194809/49‑1 (MQ=255)                                      gCAATGGACTATAGCTTTACCGTTCTCAATCTCTATAAGCTGTATCTGa           <  1:897577/49‑1 (MQ=255)                                      gCAATGGACTATAGCTTTACCGTTCTCAATCTCTATAAGCTGTATCTGAt          <  1:999976/50‑1 (MQ=255)                                       cAATGGACTATAGCTTTACCGTTCTCAATCTCTATAAGCTGTATCTGATc         >  1:2235393/1‑50 (MQ=255)                                          tGGACTATAGCTTTACCGTTCTCAATCTCTATAAGCTGTATCTGATCGtt      <  1:1514519/50‑1 (MQ=255)                                            gACTATAGCTTTACCGTTCTCAATCTCTATAAGCTGTATCTGATCGTTg     <  1:1364823/49‑1 (MQ=255)                                            gACTATAGCTTTACCGTTCTCAATCTCTATAAGCTGTATCTGATCGTTg     <  1:118307/49‑1 (MQ=255)                                            gACTATAGCTTTACCGTTCTCAATCTCTATAAGCTGTATCTGATCGTTg     <  1:2071795/49‑1 (MQ=255)                                             aCTATAGCTTTACCGTTCTCAATCTCTATAAGCTGTATCTGATCGTTGa    >  1:2170766/1‑49 (MQ=255)                                             aCTATAGCTTTACCGTTCTCAATCTCTATAAGCTGTATCTGATCGTTGAt   <  1:78141/50‑1 (MQ=255)                                             aCTATAGCTTTACCGTTCTCAATCTCTATAAGCTGTATCTGATCGTTGAt   <  1:2143909/50‑1 (MQ=255)                                               tataGCTTTACCGTTCTCAATCTCTATAAGCTGTATCTGATCGTTGATa  >  1:291/1‑49 (MQ=255)                                                  |                                              TCAGGGGAAAGGTCTGGCAACCCGTGCCGCCAAATTAGCAATGGACTATGGCTTTACCGTTCTCAATCTCTATAAGCTGTATCTGATCGTTGATA  >  NC_000913/1656465‑1656559

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