Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,532,981	C→T	G191G (GGC→GGT)	cysK →	O‑acetylserine sulfhydrylase A

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,532,981	0	C	T	100.0%	107.7 / NA	31	G191G (GGC→GGT)	cysK	O‑acetylserine sulfhydrylase A
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (18/13);  total (18/13)

CGACGGTCAGGTTGATGTATTTATTGCTGGCGTTGGGACTGGCGGTACGCTGACTGGCGTCAGCCGCTACATTAAAGGCACCAAAGGCAAGACCGATCTTATCTCTGTCGCCGTTGAGCCAACCGATTCTCCAGTTATCGCCCAGGCGCT  >  NC_000913/2532903‑2533052                                                                               |                                                                        cGACGGTCAGGTTGATGTATTTATTGCTGGCGTTGGGACTGGCGGTACGCTGACTGGCGTCAGCCGCTACATTAAAGGTAc                                                                       <  2:401163/81‑1 (MQ=255) cGACGGTCAGGTTGATGTATTTATTGCTGGCGTTGGGACTGGCGGTACGCTGACTGGCGTCAGCCGCTACATTAAAGGTAc                                                                       <  2:401162/81‑1 (MQ=255)  gACGGTCAGGTTGATGTATTTATTGCTGGCGTTGGGACTGGCGGTACGCTGACTGGCGTCAGCCGCTACATTAAAGGTAc                                                                       >  2:401204/1‑80 (MQ=255)    cGGTCAGGTTGATGTATTTATTGCTGGCGTTGGGACTGGCGGTACGCTGACTGGCGTCAGCCGCTACATTAAAGGTACCa                                                                     >  1:401192/1‑80 (MQ=255)      gTCAGGTTGATGTATTTATTGCTGGCGTTGGGACTGGCGGTACGCTGACTGGCGTCAGCCGCTACATTAAAGGTACCAAAg                                                                  >  1:401201/1‑81 (MQ=255)       tCAGGTTGATGTATTTATTGCTGGCGTTGGGACTGGCGGTACGCTGACTGGCGTCAGCCGCTACATTAAAGGTACCAAAgg                                                                 >  2:401211/1‑81 (MQ=255)        cAGGTTGATGTATTTATTGCTGGCGTTGGGACTGGCGGTACGCTGACTGGCGTCAGCCGCTACATTAAAGGTACCAAAgg                                                                 >  2:401196/1‑80 (MQ=255)         aGGTTGATGTATTTATGGCTGGCGTTGGGACTGGCGGTACGCTGACTGGCGTCAGCCGCTACATTAAAGGTACCAAAGGc                                                                <  1:401164/80‑1 (MQ=255)             tGATGTATTTATTGCTGGCGTTGGGACTGGCGGTACGCTGACTGGCGTCAGCCGCTACATTAAAGGTACCAAAGGCAAGAc                                                           >  2:401208/1‑81 (MQ=255)               aTGTATTTATTGCTGGCGTTGGGACTGGCGGTACGCTGACTGGCGTCAGCCGCTACATTAAAGGTACCAAAGGCAAGACCg                                                         <  1:401165/81‑1 (MQ=255)                tGTATTTATTGCTGGCGTTGGGACTGGCGGTACGCTGACTGGCGTCAGCCGCTACATTAAAGGTACCAGAGGCAAGACCGa                                                        <  2:401166/81‑1 (MQ=255)                    tttattGCTGGCGTTGGGACTGGCGGTACGCTGACTGGCGTCAGCCGCTACATTAAAGGTACCAAAGGCAAGACCGATCt                                                     <  1:401167/80‑1 (MQ=255)                          gCTGGCGTTGGGACTGGCGGTACGCTGACTGGCGTCAGCCGCTACATTAAAGGTACCAAAGGCAAGACCGATCTTAtctc                                               <  2:401168/80‑1 (MQ=255)                             ggCGTTGGGACTGGCGGTACGCTGACTGGCGTCAGCCGCTACATTAAAGGTACCAAAGGCAAGACCGATCTTATCTCTGTc                                           >  1:401226/1‑81 (MQ=255)                                 ttGGGACTGGCGGTACGCTGACTGGCGTCAGCCGCTACATTAAAGGTACCAAAGGCAAGACCGATCTTATCTCTGTCGCCg                                       >  2:401203/1‑81 (MQ=255)                                   gggAGTGGGGGTGCGGTGACTGGCGTCAGCCGGTACATGAGGGGTACCAAAGGCAGGACCGATCTTATCTCTGTCGCCGtt                                     <  2:401173/81‑1 (MQ=255)                                   gggACTGGCGGTACGCTGACTGGCGTCAGCCGCTACATTAAAGGTACCAAAGGCAAGACCGATCTTATCTCTGTCGCCGt                                      <  2:401169/80‑1 (MQ=255)                                       cTGGCGGTACGCTGACTGGCGTCAGCCGCTACATTAAAGGTACCAAAGGCAAGACCGATCTTATCTCTGTCGCCGTTGAGc                                 <  1:401170/81‑1 (MQ=255)                                        tGGCGGTACGCTGACTGGCGTCAGCCGCTACATTAAAGGTACCAAAGGCAAGACCGATCTTATCTCTGTCGCCGTTGAGc                                 >  1:401221/1‑80 (MQ=255)                                            ggTACGCTGACTGGCGTCAGCCGCTACATTAAAGGTACCAAAGGCAAGACCGATCTTATCTCTGTCGCCGTTGAGCCAAcc                            >  1:401217/1‑81 (MQ=255)                                                cGCTGACTGGCGTCAGCCGCTACATTAAAGGTACCAAAGGCAAGACCGATCTTATCTCTGTCGCCGTTGAGCCAACCGAtt                        >  2:401200/1‑81 (MQ=255)                                                cGCTGACTGGCGTCAGCCGCTACATTAAAGGTACCAAAGGCAAGACCGATCTTATCTCTGTCGCCGTTGAGCCAACCGAt                         <  2:401171/80‑1 (MQ=255)                                                 gCTGACTGGCGTCAGCCGCTACATTAAAGGTACCAAAGGCAAGACCGATCTTATCTCTGTCGCCGTTGAGCCAACCGATtc                       <  1:401172/81‑1 (MQ=255)                                                  cTGACTGGCGTCAGCCGCTACATTAAAGGTACCAAAGGCAAGACCGATCTTATCTCTGTCGCCGTTGAGCCAACCGATtct                      >  1:401198/1‑81 (MQ=255)                                                    gACTGGCGTCAGCCGCTACATTAAAGGTACCAAAGGCAAGACCGATCTTATCTCTGTCGCCGTTGAGCCAACCGATtctc                     >  1:401214/1‑80 (MQ=255)                                                          cGTCAGCCGCTACATTAAAGGTACCAAAGGCAAGACCGATCTTATCTCTGTCGCCGTTGAGCCAACCGATTCTCCAGTTAt              >  1:401210/1‑81 (MQ=255)                                                          cGTCAGCCGCTACATTAAAGGTACCAAAGGCAAGACCGATCTTATCTCTGTCGCCGTTGAGCCAACCGATTCTCCAGTTAt              >  2:401228/1‑81 (MQ=255)                                                             cAGCCGCTACATTAAAGGTACCAAAGGCAAGACCGATCTTATCTCTGTCGCCGTTGAGCCAACCGATTCTCCAGTTATCGc           >  1:401225/1‑81 (MQ=255)                                                              aGCCGCTACATTAAAGGTACCAAAGGCAAGACCGATCTTATCTCTGTCGCCGTTGAGCCAACCGATTCTCCAGTTATCGcc          >  2:401209/1‑81 (MQ=255)                                                              aGCCGCTACATTAAAGGTACCAAAGGCAAGACCGATCTTATCTCTGTCGCCGTTGAGCCAACCGATTCTCCAGTTATCGcc          >  2:401218/1‑81 (MQ=255)                                                                      cATTAAAGGTACCAAAGGCAAGACCGATCTTATCTCTGTCGCCGTTGAGCCAACCGATTCTCCAGTTATCGCCCAGGCGCt  <  2:401174/81‑1 (MQ=255)                                                                               |                                                                        CGACGGTCAGGTTGATGTATTTATTGCTGGCGTTGGGACTGGCGGTACGCTGACTGGCGTCAGCCGCTACATTAAAGGCACCAAAGGCAAGACCGATCTTATCTCTGTCGCCGTTGAGCCAACCGATTCTCCAGTTATCGCCCAGGCGCT  >  NC_000913/2532903‑2533052

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