Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,797,709	T→C	C64R (TGC→CGC)	ygaP →	thiosulfate sulfurtransferase YgaP

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,797,709	0	T	C	100.0%	112.5 / NA	33	C64R (TGC→CGC)	ygaP	thiosulfate sulfurtransferase YgaP
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (12/21);  total (12/21)

CAGATCTGGCTCCATTATCCGTGCTGGAACAGTCAGGTCTTCCGGCTAAATTACGTCACGAGCAAATTATATTCCACTGCCAGGCAGGTAAACGCACCAGTAATAATGCCGATAAATTAGCAGCGATTGCCGCCCCCGCAGAAATCTTTTTACTC  >  NC_000913/2797632‑2797786                                                                              |                                                                              cAGATCTGGCTCCATTATCCGTGCTGGAACAGTCAGGTCTTCCGGCTAAATTATGTCACGAGCAAATTATATTCCACCGc                                                                             >  1:450963/1‑80 (MQ=255) cAGATCTGGCTCCATTATCCGTGCTGGAACAGTCAGGTCTTCCGGCTAAATTACGTCACGAGCAAATTATATTCCACCGcc                                                                            <  1:450912/81‑1 (MQ=255) cAGATCTGGCTCCATTATCCGTGCTGGAACAGTCAGGTCTTCCGGCTAAATTACGTCACGAGCAAATTATATTCCACCGc                                                                             >  1:450951/1‑80 (MQ=255)   gATCTGGCTCCATTATCCGTGCTGGAACAGTCAGGTCTTCCGGCTAAATTACGTCACGAGCAAATTATATTCCACCGCcag                                                                          >  1:450952/1‑81 (MQ=255)       tGGCTCCATTATCCGTGCTGGAACAGTCAGGTCTTCCGGCTAAATTACGTCACGAGCAAATTATATTCCACCGCcaggcag                                                                      >  2:450956/1‑81 (MQ=255)       tGGCTCCATTATCCGTGCTGGAACAGTCAGGTCTTCCGGCTAAATTACGTCACGAGCAAATTATATTCCACCGCcaggcag                                                                      >  1:450961/1‑81 (MQ=255)         gCTCCATTATCCGTGCTGGAACAGTCAGGTCTTCCGGCTAAATTACGTCACGAGCAAATTATATTCCACCGCcaggcag                                                                      <  2:450913/79‑1 (MQ=255)          cTCCATTATCCGTGCTGGAACAGTCAGGTCTTCCGGCTAAATTACGTCACGAGCAAATTATATTCCACCGCCAGGCAGGt                                                                    >  2:450954/1‑80 (MQ=255)          cTCCATTATCCGTGCTGGAACAGTCAGGTCTTCCGGCTAAATTACGTCACGAGCAAATTATATTCCACCGCCAGGCAGGt                                                                    >  2:450965/1‑80 (MQ=255)             cATTATCCGTGCTGGAACAGTCAGGTCTTCCGGCTAAATTACGTCACGAGCAAATTATATTCCACCGCCAGGCAGGTaa                                                                  <  1:450914/79‑1 (MQ=255)                    cGTGCTGGAACAGTCAGGTCTTCCGGCTAAATTACGTCACGAGCAAATTATATTCCACCGCCAGGCAGGTAAACGCACCa                                                          <  1:450915/80‑1 (MQ=255)                        cTGGAACAGTCAGGTCTTCCGGCTAAATTACGTCACGAGCAAATTATATTCCACCGCCAGGCAGGTAAACGCACCAGtaat                                                     <  1:450919/81‑1 (MQ=255)                        cTGGAACAGTCAGGTCTTCCGGCTAAATTACGTCACGAGCAAATTATATTCCACCGCCAGGCAGGTAAACGCACCAGtaat                                                     <  2:450916/81‑1 (MQ=255)                        cTGGAACAGTCAGGTCTTCCGGCTAAATTACGTCACGAGCAAATTATATTCCACCGCCAGGCAGGTAAACGCACCAGtaat                                                     <  2:450917/81‑1 (MQ=255)                        cTGGAACAGTCAGGTCTTCCGGCTAAATTACGTCACGAGCAAATTATATTCCACCGCCAGGCAGGTAAACGCACCAGtaat                                                     <  2:450918/81‑1 (MQ=255)                                gTCAGGTCTTCCGGCTAAATTACGTCACGAGCAAATTATATTCCACCGCCAGGCAGGTAAACGCACCAGTAATAATGCCGa                                             <  1:450920/81‑1 (MQ=255)                                gTCAGGTCTTCCGGCTAAATTACGTCACGAGCAAATTATATTCCACCGCCAGGCAGGTAAACGCACCAGTAATAATGCCGa                                             <  2:450921/81‑1 (MQ=255)                                 tCAGGTCTTCCGGCTAAATTACGTCACGAGCAAATTATATTCCACCGCCAGGCAGGTAAACGCACCAGTAATAATGCCGAt                                            <  2:450922/81‑1 (MQ=255)                                   aGGTCTTCCGGCTAAATTACGTCACGAGCAAATTATATTCCACCGCCAGGCAGGTAAACGCACCAGTAATAATGCCGATaa                                          <  2:450923/81‑1 (MQ=255)                                   aGGTCTTCCGGCTAAATTACGTCACGAGCAAATTATATTCCACCGCCAGGCAGGTAAACGCACCAGTAATAATGCCGATa                                           <  1:450924/80‑1 (MQ=255)                                       cTTCCGGCTAAATTACGTCACGAGCAAATTATATTCCACCGCCAGGCAGGTAAACGCACCAGTAATAATGCCGATAAAt                                        >  1:450959/1‑79 (MQ=255)                                        ttCCGGCTAAATTACGTCACGAGCAAATTATATTCCACCGCCAGGCAGGTAAACGCACCAGTAATAATGCCGATAAATTag                                     >  1:450958/1‑81 (MQ=255)                                          ccGGCTAAATTACGTCACGAGCAAATTATATTCCACCGCCAGGCAGGTAAACGCACCAGTAATAATGCCGATAAATTagc                                    <  2:450925/80‑1 (MQ=255)                                            ggCTAAATTACGTCACGAGCAAATTATATTCCACCGCCAGGCAGGTAAACGCACCAGTAATAATGCCGATAAATTagcagc                                 <  1:450926/81‑1 (MQ=255)                                                 aaTTACGTCACGAGCAAATTATATTCCACCGCCAGGCAGGTAAACGCACCAGTAATAATGCCGATAAATTAGCAGCGAtt                             <  2:450927/80‑1 (MQ=255)                                                 aaTTACGTCACGAGCAAATTATATTCCACCGCCAGGCAGGTAAACGCACCAGTAATAATGCCGATAAATTAGCAGCGAtt                             <  1:450928/80‑1 (MQ=255)                                                   ttACGTCACGAGCAAATTATATTCCACCGCCAGGCAGGTAAACGCACCAGTAATAATGCCGATAAATTAGCAGCGATTgcc                          <  1:450929/81‑1 (MQ=255)                                                     aCGTCACGAGCAAATTATATTCCACCGCCAGGCAGGTAAACGCACCAGTAATAATGCCGATAAATTAGCAGCGATTgccgc                        >  2:450979/1‑81 (MQ=255)                                                     aCGTCACGAGCAAATTATATTCCACCGCCAGGCAGGTAAACGCACCAGTAATAATGCCGATAAATTAGCAGCGATTgccgc                        >  1:450972/1‑81 (MQ=255)                                                      cGTCACGAGCAAATTATATTCCACCGCCAGGCAGGTAAACGCACCAGTAATAATGCCGATAAATTAGCAGCGATTgccgcc                       <  2:450930/81‑1 (MQ=255)                                                           cGAGCAAATTATATTCCACCGCCAGGCAGGTAAACGCACCAGTAATAATGCCGATAAATTAGCAGCGATTGCCGCCCCCGc                  <  2:450931/81‑1 (MQ=255)                                                                      tatTCCACCGCCAGGCAGGTAAACGCACCAGTAATAATGCCGATAAATTAGCAGCGATTGCCGCCCCCGCAGAAATCtttt       >  1:450966/1‑81 (MQ=255)                                                                           cACCGCCAGGCAGGTAAACGCACCAGTAATAATGCCGATAAATTAGCAGCGATTGCCGCCCCCGCAGAAATCTTTTTACTc  <  1:450932/81‑1 (MQ=255)                                                                              |                                                                              CAGATCTGGCTCCATTATCCGTGCTGGAACAGTCAGGTCTTCCGGCTAAATTACGTCACGAGCAAATTATATTCCACTGCCAGGCAGGTAAACGCACCAGTAATAATGCCGATAAATTAGCAGCGATTGCCGCCCCCGCAGAAATCTTTTTACTC  >  NC_000913/2797632‑2797786

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