Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,797,916	T→C	Y133H (TAT→CAT)	ygaP →	thiosulfate sulfurtransferase YgaP

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,797,916	0	T	C	100.0%	185.3 / NA	51	Y133H (TAT→CAT)	ygaP	thiosulfate sulfurtransferase YgaP
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (26/25);  total (26/25)

CGAAGATGGGATTGATGGCTGGAAAAAAGCAGGATTGCCAGTAGCGGTAAATAAATCTCAACCCTTGCCGTTAATGCGCCAGGTGCAGATCGCTGCGGGTGGTTTAATATTAATCGGCGTTGTACTGGGTTATACCGTAAATAGCGGTTTCTTCTTATTAAGTGGCTTTGTTGGTGCCGGGTTACTGTTTGCAGGAATCAGCGGTTTTTGTGGAATGGCAAGGTTGTTAGATAAGATGCCGTGGAACCAACGAGCTTGAGAAGCGAC  >  NC_000913/2797786‑2798052                                                                                                                                   |                                                                                                                                         cGAAGATGGGATTGATGGCTGGAAAAAAGCAGGATTGCCAGTAGCGGTAAATAAATCTCAACCCTTGCCGTTAATGCGCCAGGTGCAGATCGCTGCGGGTGGTTTAATATTAATCGGCGTTGTACTGGGTCATACCGTaa                                                                                                                                 >  1:823023/1‑140 (MQ=255)          gATTGATGGCTGGAAAAAAGCAGGATTGCCAGTAGCGGTAAATAAATCTCAACCCTTGCCGTTAATGCGCCAGGTGCAGATCGCTGCGGGTGGTTTAATATTAATCGGCGTTGTACTGGGTCATACCGTAAATAGCGGtt                                                                                                                        >  2:259870/1‑140 (MQ=255)          gATTGATGGCTGGAAAAAAGCAGGATTGCCAGTAGCGGTAAATAAATCTCAACCCTTGCCGTTAATGCGCCAGGTGCAGATCGCTGCGGGTGGTTTAATATTAATCGGCGTTGTACTGGGTCATACCGTAAATAGCGGtt                                                                                                                        <  2:717339/140‑1 (MQ=255)             tGATGGCTGGAAAAAAGCAGGATTGCCAGTAGCGGTAAATAAATCTCAACCCTTGCCGTTAATGCGCCAGGTGCAGATCGCTGCGGGTGGTTTAATATTAATCGGCGTTGTACTGGGTCATACCGTAAATAGCGGTttct                                                                                                                     <  2:867869/140‑1 (MQ=255)              gATGGCTGGAAAAAAGCAGGATTGCCAGTAGCGGTAAATAAATCTCAACCCTTGCCGTTAATGCGCCAGGTGCAGATCGCTGCGGGTGGTTTAATATTAATCGGCGTTGTACTGGGTCATACCGTAAATAGCGGTttctt                                                                                                                    <  2:158373/140‑1 (MQ=255)              gATGGCTGGAAAAAAGCAGGATTGCCAGTAGCGGTAAATAAATCTCAACCCTTGCCGTTAATGCGCCAGGTGCAGATCGCTGCGGGTGGTTTAATATTAATCGGCGTTGTACTGGGTCATACCGTAAATAGCGGTttctt                                                                                                                    >  1:653087/1‑140 (MQ=255)              gATGGCTGGAAAAAAGCAGGATTGCCAGTAGCGGTAAATAAATCTCAACCCTTGCCGTTAATGCGCCAGGTGCAGATCGCTGCGGGTGGTTTAATATTAATCGGCGTTGTACTGGGTCATACCGTAAATAGCGGTttctt                                                                                                                    >  2:397151/1‑140 (MQ=255)               aTGGCTGGAAAAAAGCAGGATTGCCAGTAGCGGTAAATAAATCTCAACCCTTGCCGTTAATGCGCCAGGTGCAGATCGCTGCGGGTGGTTTAATATTAATCGGCGTTGTACTGGGTCATACCGTAAATAGCGGTttcttc                                                                                                                   >  1:38473/1‑140 (MQ=255)                      gAAAAAAGCAGGATTGCCAGTAGCGGTAAATAAATCTCAACCCTTGCCGTTAATGCGCCAGGTGCAGATCGCTGCGGGTGGTTTAATATTAATCGGCGTTGTACTGGGTCATACCGTAAATAGCGGTTTCTTCTTATTaa                                                                                                            <  2:891172/140‑1 (MQ=255)                                ggATTGCCAGTAGCGGTAAATAAATCTCAACCCTTGCCGTTAATGCGCCAGGTGCAGATCGCTGCGGGTGGTTTAATATTAATCGGCGTTGTACTGGGTCATACCGTAAATAGCGGTTTCTTCTTATTAAGTGGCTttgt                                                                                                  >  1:134633/1‑140 (MQ=255)                                ggATTGCCAGTAGCGGTAAATAAATCTCAACCCTTGCCGTTAATGCGCCAGGTGCAGATCGCTGCGGGTGGTTTAATATTAATCGGCGTTGTACTGGGTCATACCGTAAATAGCGGTTTCTTCTTATTAAGTGGCTttgt                                                                                                  <  2:363493/140‑1 (MQ=255)                                ggATTGCCAGTAGCGGTAAATAAATCTCAACCCTTGCCGTTAATGCGCCAGGTGCAGATCGCTGCGGGTGGTTTAATATTAATCGGCGTTGTACTGGGTCATACCGTAAATAGCGGTTTCTTCTTATTAAGTGGCTttgt                                                                                                  <  2:823023/140‑1 (MQ=255)                                   ttGCCAGTAGCGGTAAATAAATCTCAACCCTTGCCGTTAATGCGCCAGGTGCAGATCGCTGCGGGTGGTTTAATATTAATCGGCGTTGTACTGGGTCATACCGTAAATAGCGGTTTCTTCTTATTAAGTGGCTTTGTTgg                                                                                               <  1:397151/140‑1 (MQ=255)                                   ttGCCAGTAGCGGTAAATAAATCTCAACCCTTGCCGTTAATGCGCCAGGTGCAGATCGCTGCGGGTGGTTTAATATTAATCGGCGTTGTACTGGGTCATACCGTAAATAGCGGTTTCTTCTTATTAAGTGGCTTTGTTgg                                                                                               <  2:1052281/140‑1 (MQ=255)                                     gCCAGTAGCGGTAAATAAATCTCAACCCTTGCCGTTAATGCGCCAGGTGCAGATCGCTGCGGGTGGTTTAATATTAATCGGCGTTGTACTGGGTCATACCGTAAATAGCGGTTTCTTCTTATTAAGTGGCTTTGTTGGTg                                                                                             <  2:717129/140‑1 (MQ=255)                                        aGTAGCGGTAAATAAATCTCAACCCTTGCCGTTAATGCGCCAGGTGCAGATCGCTGCGGGTGGTTTAATATTAATCGGCGTTGTACTGGGTCATACCGTAAATAGCGGTTTCTTCTTATTAAGTGGCTTTGTTGGTGCCg                                                                                          >  2:115595/1‑140 (MQ=255)                                             cGGTAAATAAATCTCAACCCTTGCCGTTAATGCGCCAGGTGCAGATCGCTGCGGGTGGTTTAATATTAATCGGCGTTGTACTGGGTCATACCGTAAATAGCGGTTTCTTCTTATTAAGTGGCTTTGTTGGTGCCGGGTTa                                                                                     >  1:987380/1‑140 (MQ=255)                                               gTAAATAAATCTCAACCCTTGCCGTTAATGCGCCAGGTGCAGATCGCTGCGGGTGGTTTAATATTAATCGGCGTTGTACTGGGTCATACCGTAAATAGCGGTTTCTTCTTATTAAGTGGCTTTGTTGGTGCCGGGTTACt                                                                                   >  1:708801/1‑140 (MQ=255)                                               gTAAATAAACCTCAACCCTTGCCGTTAATGCGCCAGGTGCAGATCGCTGCGGGTGGTTTAATATTAATCGGCGTTGTACTGGGTCATACCGTAAATAGCGGTTTCTTCTTATTAAGTGGCTTTGTTGGTGCCGGGTTACt                                                                                   >  2:985885/1‑140 (MQ=255)                                                taaataaaTCTCAACCCTTGCCGTTAATGCGCCAGGTGCAGATCGCTGCGGGTGGTTTAATATTAATCGGCGTTGTACTGGGTCATACCGTAAATAGCGGTTTCTTCTTATTAAGTGGCTTTGTTGGTGcc                                                                                           >  1:983829/1‑131 (MQ=255)                                                taaataaaTCTCAACCCTTGCCGTTAATGCGCCAGGTGCAGATCGCTGCGGGTGGTTTAATATTAATCGGCGTTGTACTGGGTCATACCGTAAATAGCGGTTTCTTCTTATTAAGTGGCTTTGTTGGTGcc                                                                                           <  2:983829/131‑1 (MQ=255)                                                taaataaaTCTCAACCCTTGCCGTTAATGCGCCAGGTGCAGATCGCTGCGGGTGGTTTAATATTAATCGGCGTTGTACTGGGTCATACCGTAAATAGCGGTTTCTTCTTATTAAGTGGCTTTGTTGGTGCCGGGTTACTg                                                                                  >  2:342122/1‑140 (MQ=255)                                                 aaataaatCTCAACCCTTGCCGTTAATGCGCCAGGTGCAGATCGCTGCGGGTGGTTTAATATTAATCGGCGTTGTACTGGGTCATACCGTAAATAGCGGTTTCTTCTTATTAAGTGGCTTTGTTGGTGCCGGGTTACTGt                                                                                 <  2:1007825/140‑1 (MQ=255)                                                          tcAACCCTTGCCGTTAATGCGCCAGGTGCAGATCGCTGCGGGTGGTTTAATATTAATCGGCGTTGTACTGGGTCATACCGTAAATAGCGGTTTCTTCTTATTAAGTGGCTTTGTTGGTGCCGGGTTACTGTTTGCAGGaa                                                                        >  2:1045318/1‑140 (MQ=255)                                                            aaCCCTTGCCGTTAATGCGCCAGGTGCAGATCGCTGCGGGTGGTTTAATATTAATCGGCGTTGTACTGGGTCATACCGTAAATAGCGGTTTCTTCTTATTAAGTGGCTTTGTTGGTGCCGGGTTACTGTTTGCAGGAATc                                                                      <  1:115595/140‑1 (MQ=255)                                                              cccTTGCCGTTAATGCGCCAGGTGCAGATCGCTGCGGGTGGTTTAATATTAATCGGCGTTGTACTGGGTCATACCGTAAATAGCGGTTTCTTCTTATTAAGTGGCTTTGTTGGTGCCGGGTTACt                                                                                   <  1:823935/125‑1 (MQ=255)                                                              cccTTGCCGTTAATGCGCCAGGTGCAGATCGCTGCGGGTGGTTTAATATTAATCGGCGTTGTACTGGGTCATACCGTAAATAGCGGTTTCTTCTTATTAAGTGGCTTTGTTGGTGCCGGGTTACt                                                                                   >  2:823935/1‑125 (MQ=255)                                                                         aaTGCGCCAGGTGCAGATCGCTGCGGGTGGTTTAATATTAATCGGCGTTGTACTGGGTCATACCGTAAATAGCGGTTTCTTCTTATTAAGTGGCTTTGTTGGTGCCGGGTTACTGTTTGCAGGAATCAGCGGTTTTtgtg                                                         <  2:38473/140‑1 (MQ=255)                                                                          aTGCGCCAGGTGCAGATCGCTGCGGGTGGTTTAATATTAATCGGCGTTGTACTGGGTCATACCGTAAATAGCGGTTTCTTCTTATTAAGTGGCTTTGTTGGTGCCGGGTTACTGTTTGCAGGAATCAGCGGTTTTTGTgg                                                        <  1:259870/140‑1 (MQ=255)                                                                           tGCGCCAGGTGCAGATCGCTGCGGGTGGTTTAATATTAATCGGCGTTGTACTGGGTCATACCGTAAATAGCGGTTTCTTCTTATTAAGTGGCTTTGTTGGTGCCGGGTTACTGTTTGCAGGAATCAGCGGTTTTTGTGGa                                                       >  2:1070138/1‑140 (MQ=255)                                                                               ccAGGTGCAGATCGCTGCGGGTGGTTTAATATTAATCGGCGTTGTACTGGGTCATACCGTAAATAGCGGTTTCTTCTTATTAAGTGGCTTTGTTGGTGCCGGGTTACTGTTTGCAGGAATCAGCGGTTTTTGTGGAATgg                                                   >  1:669305/1‑140 (MQ=255)                                                                                 aGGTGCAGATCGCTGCGGGTGGTTTAATATTAATCGGCGTTGTACTGGGTCATACCGTAAATAGCGGTTTCTTCTTATTAAGTGGCTTTGTTGGTGCCGGGTTACTGTTTGCAGGAATCAGCGGTTTTTGTGGAATGGCa                                                 <  1:1045318/140‑1 (MQ=255)                                                                                  ggTGCAGATCGCTGCGGGTGGTTTAATATTAATCGGCGTTGTACTGGGTCATACCGTAAATAGCGGTTTCTTCTTATTAAGTGGCTTTGTTGGTGCCGGGTTACTGTTTGCAGGAATCAGCGGTTTTTGTGGAATGGCaa                                                <  2:371967/140‑1 (MQ=255)                                                                                    tGCAGATCGCTGCGGGTGGTTTAATATTAATCGGCGTTGTACTGGGTCATACCGTAAATAGCGGTTTCTTCTTATTAAGTGTCTTTGTTGGTGCCGGGTTACTGTTTGCAGGAATCAGCGGTTTTTGTGGAATGGCAAgg                                              <  2:266372/140‑1 (MQ=255)                                                                                     gCAGATCGCTCCGGCTGGTTTAATATTAATCGGCGTTGTACTGGGTCATACCGTAAATAGCGGTTTCTTCTTATTAAGTGGCTTTGTTGGTGCCGGGTTACTGTTTGCAGGAATCAGCGGTTTTTGTGGAATGGCAAGgt                                             <  2:708801/140‑1 (MQ=255)                                                                                      cAGATCGCTGCGGGTGGTTTAATATTAATCGGCGTTGTACTGGGTCATACCGTAAATAGCGGTTTCTTCTTATTAAGTGGCTTTGTTGGTGCCGGGTTACTGTTTGCAGGAATCAGCGGTTTTTGTGGAATGGCAAGgtt                                            <  1:985885/140‑1 (MQ=255)                                                                                      cAGATCGCTGCGGGTGGTTTAATATTAATCGGCGTTGTACTGGGTCATACCGTAAATAGCGGTTTCTTCTTATTAAGTGGCTTTGTTGGTGCCGGGTTACTGTTTGCAGGAATCAGCGGTTTTTGTGGAATGGCAAGgtt                                            <  1:342122/140‑1 (MQ=255)                                                                                         aTCGCTGCGGGTGGTTTAATATTAATCGGCGTTGTACTGGGTCATACCGTAAATAGCGGTTTCTTCTTATTAAGTGGCTTTGTTGGTGCCGGGTTACTGTTTGCAGGAATCAGCGGTTTTTGTGGAATGGCAAGgttgtt                                         <  2:987380/140‑1 (MQ=255)                                                                                            gCTGCGGGTGGTTTAATATTAATCGGCGTTGTACTGGGTCATACCGTAAATAGCGGTTTCTTCTTATTAAGTGGCTTTGTTGGTGCCGGGTTACTGTTTGCAGGAATCAGCGGTTTTTGTGGAATGGCAAGGTTGTTAGa                                      >  1:1103949/1‑140 (MQ=255)                                                                                             cTGCGGGTGGTTTAATATTAATCGGCGTTGTACTGGGTCATACCGTAAATAGCGGTTTCTTCTTATTAAGTGGCTTTGTTGGTGCCGGGTTACTGTTTGCAGGAATCAGCGGTTTTTGTGGAATGGCAAGGTTGTTAGAt                                     >  2:923745/1‑140 (MQ=255)                                                                                              tGCGGGTGGTTTAATATTAATCGGCGTTGTACTGGGTCATACCGTAAATAGCGGTTTCTTCTTATTAAGTGGCTTTGTTGGTGCCGGGTTACTGTTTGCAGGAATCAGCGGTTTTTGTGGAATGGCAAGGTTGTTAGATa                                    >  2:204003/1‑140 (MQ=255)                                                                                              tGCGGGTGGTTTAATATTAATCGGCGTTGTACTGGGTCATACCGTAAATAGCGGTTTCTTCTTATTAAGTGGCTTTGTTGGTGCCGGGTTACTGTTTGCAGGAATCAGCGGTTTTTGTGGAATGGCAAGGTTGTTAGATa                                    >  1:154131/1‑140 (MQ=255)                                                                                                     ggtTTAATATTAATCGGCGTTGTACTGGGTCATACCGTAAATAGCGGTTTCTTCTTATTAAGTGGCTTTGTTGGTGCCGGGTTACTGTTTGCAGGAATCAGCGGTTTTTGTGGAATGGCAAGGTTGTTAGATAAGATGcc                             >  1:1034992/1‑140 (MQ=255)                                                                                                        ttAATATTAATCGGCGTTGTACTGGGTCATACCGTAAATAGCGGTTTCTTCTTATTAAGTGGCTTTGTTGGTGCCGGGTTACTGTTTGCAGGAATCAGCGGTTTTTGTGGAATGGCAAGGTTGTTAGATAAGATGCCGTg                          <  2:653087/140‑1 (MQ=255)                                                                                                              ttAATCGGCGTTGTACTGGGTCATACCGTAAATAGCGGTTTCTTCTTATTAAGTGGCTTTGTTGGTGCCGGGTTACTGTTTGCAGGAATCAGCGGTTTTTGTGGAATGGCAAGGTTGTTAGATAAGATGCCGTGGAACCa                    <  2:1034992/140‑1 (MQ=255)                                                                                                              ttAATCGGCGTTGTACTGGGTCATACCGTAAATAGCGGTTTCTTCTTATTAAGTGGCTTTGTTGGTGCCGGGTTACTGTTTGCAGGAATCAGCGGTTTTTGTGGAATGGCAAGGTTGTTAGATAAGATGCCGTGGAACCa                    <  2:134633/140‑1 (MQ=255)                                                                                                                 aTCGGCGTTGTACTGGGTCATACCGTAAATAGCGGTTTCTTCTTATTAAGTGGCTTTGTTGGTGCCGGGTTACTGTTTGCAGGAATCAGCGGTTTTTGTGGAATGGCAAGGTTGTTAGATAAGATGCCGTGGAACCAACg                 >  2:1027672/1‑140 (MQ=255)                                                                                                                 aTCGGCGTTGTACTGGGTCATACCGTAAATAGCGGTTTCTTCTTATTAAGTGGCTTTGTTGGTGCCGGGTTACTGTTTGCAGGAATCAGCGGTTTTTGTGGAATGGCAAGGTTGTTAGATAAGATGCCGTGGAACCAACg                 >  2:915901/1‑140 (MQ=255)                                                                                                                   cGGCGTTGTACTGGGTCATACCGTAAATAGCGGTTTCTTCTTATTAAGTGGCTTTGTTGGTGCCGGGTTACTGTTTGCAGGAATCAGCGGTTTTTGTGGAATGGCAAGGTTGTTAGATAAGATGCCGTGGAACCAACGAg               >  2:750583/1‑140 (MQ=255)                                                                                                                      cGTTGTACTGGGTCATACCGTAAATAGCGGTTTCTTCTTATTAAGTGGCTTTGTTGGTGCCGGGTTACTGTTTGCAGGAATCAGCGGTTTTTGTGGAATGGCAAGGTTGTTAGATAAGATGCCGTGGAACCAACGAGCtt            >  1:538500/1‑140 (MQ=255)                                                                                                                                ggTCATACCGTAAATAGCGGTTTCTTCTTATTAAGTGGCTTTGTTGGTGCCGGGTTACTGTTTGCAGGAATCAGCGGTTTTTGTGGAATGGCAAGGTTGTTAGATAAGATGCCGTGGAACCAACGAGCTTGAGAAGCGAc  >  1:901256/1‑140 (MQ=255)                                                                                                                                   |                                                                                                                                         CGAAGATGGGATTGATGGCTGGAAAAAAGCAGGATTGCCAGTAGCGGTAAATAAATCTCAACCCTTGCCGTTAATGCGCCAGGTGCAGATCGCTGCGGGTGGTTTAATATTAATCGGCGTTGTACTGGGTTATACCGTAAATAGCGGTTTCTTCTTATTAAGTGGCTTTGTTGGTGCCGGGTTACTGTTTGCAGGAATCAGCGGTTTTTGTGGAATGGCAAGGTTGTTAGATAAGATGCCGTGGAACCAACGAGCTTGAGAAGCGAC  >  NC_000913/2797786‑2798052

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