Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,370,155	A→G	intergenic (+152/‑61)	pspE → / → ycjM	thiosulfate:cyanide sulfurtransferase (rhodanese)/alpha amylase catalytic domain family protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,370,155	0	A	G	100.0%	142.4 / NA	38	intergenic (+152/‑61)	pspE/ycjM	thiosulfate:cyanide sulfurtransferase (rhodanese)/alpha amylase catalytic domain family protein
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (24/14);  total (24/14)

AATTCACCATGAGCGGCGCTTATGCGCCGTTTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTGGTAAAAAATGATTTCCTGCGTGACTAAAACCCTTGTGCTCAATTGACAGTTTATTTTCTGCGGAGTAGTCTCTCGTTTCATGGGACCGCTACCACGGAAAGGCAACATGAAACAGAAAATTACGGATTACCTGGACGAAATCTACGGTGGAACATTTACCGCAACTCATTTACAGAAACTTGTAACGCGTCT  >  NC_000913/1370007‑1370301                                                                                                                                                     |                                                                                                                                                   aaTTCACCATGAGCGGCGCTTATGCGCCGTTTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTGGTAAAAAATGATTTCCTGCGTGACTAAAACCCTTGTGCTCAATTGGCa                                                                                                                                                  <  1:639479/151‑1 (MQ=255)     cACCATGAGCGGCGCTTATGCGCCGTTTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTGGTAAAAAATGATTTCCTGCGTGACTAAAACCCTTGTGCTCAATTGGCAGttt                                                                                                                                              >  2:140604/1‑151 (MQ=255)           gAGCGCCGCTTATGCGCCGTGTTTTTTTCTGGCAAACCTTATTTACATCCACATAGATTACTTGCGTCAGCTCACAAATACGGTTTTTCCCTGGTAAAAAATGATTTCCTGCGTGACTAAAACCCTTGTGCTCAATTGGCAGTTTAtttt                                                                                                                                         >  1:33107/1‑150 (MQ=255)                gcgcTTATGCGCCGTTTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTGGTAAAAAATGATTTCCTGCGTGACTAAAACCCTTGTGCTCAATTGGCAGTTTATTTTCTGCgg                                                                                                                                   <  2:373184/151‑1 (MQ=255)                  gcTTATGCGCCGTTTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTGGTAAAAAATGATTTCCTGCGTGACTAAAACCCTTGTGCTCAATTGGCAGTTTATTTTCTGCgg                                                                                                                                   >  1:114942/1‑149 (MQ=255)                      aTGCGCCGTTTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTGGTAAAAAATGATTTCCTGCGTGACTAAAACCCTTGTGCTCAATTGGCAGTTTATTTTCTGCGGagtagt                                                                                                                             <  2:411938/151‑1 (MQ=255)                          gcCGTTTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTGGTAAAAAATGATTTCCTGCGTGACTAAAACCCTTGTGCTCAATTGGCAGTTTATTTTCTGCGGAGTAGtctct                                                                                                                         >  2:40189/1‑151 (MQ=255)                               ttttttttCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTGGTAAAAAATGATTTCCTGCGTGACTAAAACCCTTGTGCTCAATTGGCAGTTTATTTTCTGCGGAGTAGTCTCTCGttt                                                                                                                    >  2:448556/1‑151 (MQ=255)                               ttttttttCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTGGTAAAAAATGATTTCCTGCGTGACTAAAACCCTTGTGCTCAATTGGCAGTTTATTTTCTGCGGAGTAGTCTCTCGtt                                                                                                                     <  1:588531/150‑1 (MQ=255)                                 ttttttCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTGGTAAAAAATGATTTCCTGCGTGACTAAAACCCTTGTGCTCAATTGGCAGTTTATTTTCTGCGGAGTAGTCTCTCGTTTc                                                                                                                   >  1:327111/1‑150 (MQ=255)                                    tttCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTGGTAAAAAATGATTTCCTGCGTGACTAAAACCCTTGTGCTCAATTGGCAGTTTATTTTCTGCGGAGTAGTCTCTCGTTTCATg                                                                                                                <  2:226744/150‑1 (MQ=255)                                         gTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTGGTAAAAAATGATTTCCTGCGTGACTAAAACCCTTGTGCTCAATTGGCAGTTTATTTTCTGCGGAGTAGTCTCTCGTTTCATGGGAc                                                                                                            <  1:379297/149‑1 (MQ=255)                                                 ttATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTGGTAAAAAATGATTTCCTGCGTGACTAAAACCCTTGTGCTCAATTGGCAGTTTATTTTCTGCGGAGTAGTCTCTCGTTTCATGGGACCGCTACCACg                                                                                                  >  1:617173/1‑151 (MQ=255)                                                  tATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTGGTAAAAAATGATTTCCTGCGTGACTAAAACCCTTGTGCTCAATTGGCAGTTTATTTTCTGCGGAGTAGTCTCTCGTTTCATGGGACCGCTACCACgg                                                                                                 >  2:448579/1‑151 (MQ=255)                                                               aTAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTGGTAAAAAATGATTTCCTGCGTGACTAAAACCCTTGTGCTCAATTGGCAGTTTATTTTCTGCGGAGTAGTCTCTCGTTTCATGGGACCGCTACCACGGAAAGGCAACATg                                                                                     >  1:584642/1‑150 (MQ=255)                                                                 aGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTGGTAAAAAATGATTTCCTGCGTGACTAAAACCCTTGTGCTCAATTGGCAGTTTATTTTCTGCGGAGTAGTCTCTCGTTTCATGGGACCGCTACCACGGAAAGGCAACATGaa                                                                                   >  2:389739/1‑150 (MQ=255)                                                                          gCGTCAGCTCACAAATACGCTTTTTCCCTGGTAAAAAATGATTTCCTGCGTGACTAAAACCCTTGTGCTCAATTGGCAGTTTATTTTCTGCGGAGTAGTCTCTCGTTTCATGGGACCGCTACCACGGAAAGGCAACATGAAACAGAAAAtt                                                                         <  1:370418/151‑1 (MQ=255)                                                                              cAGCTCACAAATACGCTTTTTCCCTGGTAAAAAATGATTTCCTGCGTGACTAAAACCCTTGTGCTCAATTGGCAGTTTATGTTCTGCGGAGTAGTCTCTCGTTTCATGGGACCGCTACCACGGAAAGGCAACATGAAACAGAAAATTACg                                                                      <  1:341586/150‑1 (MQ=255)                                                                                gCTCACAAATACGCTTTTTCCCTGGTAAAAAATGATTTCCTGCGTGACTAAAACCCTTGTGCTCAATTGGCAGTTTATTTTCTGCGGAGTAGTCTCTCGTTTCATGGGACCGCTACCACGGAAAGGCAACATGAAACAGAAAATTACGGAt                                                                   >  2:497580/1‑151 (MQ=255)                                                                                  tCACAAATACGATTTTTCCCTGGTAAAAAATGATTTCCTGCTTGACTAAAACCCTTGTGCTCAATTGGCAGTTTATTTTCTGCGTAGTAGTCTCTCGTTTCATGGGACCGCTACCACGGAAAGGCAACATGAAACAGAAAATTACGGAtt                                                                  >  1:262422/1‑150 (MQ=255)                                                                                           cGCTTTTTCCCTGGTAAAAAATGATTTCCTGCGTGACTAAAACCCTTGTGCTCAATTGGCAGTTTATTTTCTGCGGAGTAGTCTCTCGTTTCATGGGACCGCTACCACGGAAAGGCAACATGAAACAGAAAATTACGGATTACCTGGACg                                                         >  1:561281/1‑150 (MQ=255)                                                                                            gCTTTTTCCCTGGTAAAAAATGATTTCCTGCGTGACTAAAACCCTTGTGCTCAATTGGCAGTTTATTTTCTGCGGAGTAGTCTCTCGTTTCATGGGACCGCTACCACGGAAAGGCAACATGAAACAGAAAATTACGGATTACCTGGACGa                                                        >  1:554695/1‑150 (MQ=255)                                                                                                    ccTGGTAAAAAATGATTTCCTGCGTGACTAAAACCCTTGTGCTCAATTGGCAGTTTATTTTCTGCGGAGTAGTCTCTCGTTTCATGGGACCGCTACCACGGAAAGGCAACATGAAACAGAAAATTACGGATTACCTGGACGAAATCTACgg                                               >  2:632282/1‑151 (MQ=255)                                                                                                     cTGGTAAAAAATGATTTCCTGCGTGACTAAAACCCTTGTGCTCAATTGGCAGTTTATTTTCTGCGGAGTAGTCTCTCGTTTCATGGGACCGCTACCACGGAAAGGCAACATGAAACAGAAAATTACGGATTACCTGGACGAAATCTACGGt                                              >  2:129178/1‑151 (MQ=255)                                                                                                           aaaaaTGATTTCCTGCGTGACTAAAACCCTTGTGCTCAATTGGCAGTTTATTTTCTGCGGAGTAGTCTCTCGTTTCATGGGACCGCTACCACGGAAAGGCAACATGAAACAGAAAATTACGGATTACCTGGACGAAATCTACGGTGGAACa                                        <  2:548647/151‑1 (MQ=255)                                                                                                            aaaaTGATTTCCTGCGTGACTAAAACCCTTGTGCTCAATTGGCAGTTTATTTTCTGCGGAGTAGTCTCTCGTTTCATGGGACCGCTACCACGGAAAGGCAACATGAAACAGAAAATTACGGATTACCTGGACGAAATCTACGGTGGAACAt                                       >  1:116758/1‑151 (MQ=255)                                                                                                             aaaTGATTTCCTGCGTGACTAAAACCCTTGTGCTCAATTGGCAGTTTATTTTCTGCGGAGTAGTCTCTCGTTTCATGGGACCGCTACCACGGAAAGGCAACATGAAACAGAAAATTACGGATTACCTGGACGAAATCTACGGTGGAACAtt                                      >  2:525276/1‑151 (MQ=255)                                                                                                                  aTTTCCTGCGTGACTAAAACCCTTGTGCTCAATTGGCAGTTTATTTTCTGCGGAGTAGTCTCTCGTTTCATGGGACCGCTACCACGGAAAGGCAACATGAAACAGAAAATTACGGATTACCTGGACGAAATCTACGGTGGAACATTTACCg                                 <  1:240937/151‑1 (MQ=255)                                                                                                                   tttCCTGCGTGACTAAAACCCTTGTGCTCAATTGGCAGTTTATTTTCTGCGGAGTAGTCTCTCGTTTCATGGGACCGCTACCACGGAAAGGCAACATGAAACAGAAAATTACGGATTACCTGGACGAAATCTACGGTGGAACATTTACCGc                                >  2:309666/1‑151 (MQ=255)                                                                                                                   tttCCTGCGTGACTAAAACCCTTGTGCTCAATTGGCAGTTTATTTTCTGCGGAGTAGTCTCTCGTTTCATGGGACCGCTACCACGGAAAGGCAACATGAAACAGAAAATTACGGATTACCTGGACGAAATCTACGGTGGAACATTTACCGc                                >  2:610554/1‑151 (MQ=255)                                                                                                                      ccTGCGTGACTAAAACCCTTGTGCTCAATTGGCAGTTTATTTTCTGCGGAGTAGTCTCTCGTTTCATGGGACCGCTACCACGGAAAGGCAACATGAAACAGAAAATTACGGATTACCTGGACGAAATCTACGGTGGAACATTTACCGCa                               <  1:521378/149‑1 (MQ=255)                                                                                                                             gACTAAAACCCTTGTGCTCAATTGGCAGTTTATTTTCTGCGGAGTAGTCTCTCGTTTCATGGGACCGCTACCACGGAAAGGCAACATGAAACAGAAAATTACGGATTACCTGGACGAAATCTACGGTGGAACATTTACCGCAACTCAttt                       <  1:389735/150‑1 (MQ=255)                                                                                                                             gACTAAAACCCTTGTGCTCAATTGGCAGTTTATTTTCTGCGGAGTAGTCTCTCGTTTCATGGGACCGCTACCACGGAAAGGCAACATGAAACAGAAAATTACGGATTACCTGGACGAAATCTACGGTGGAACATTTACCGCAACTCAttt                       <  1:129368/150‑1 (MQ=255)                                                                                                                              aCTAAAACCCTTGTGCTCAATTGGCAGTTTATTTTCTGCGGAGTAGTCTCTCGTTTCATGGGACCGCTACCACGGAAAGGCAACATGAAACAGAAAATTACGGATTACCTGGACGAAATCTACGGTGGAACATTTACCGCAACTCATTTAc                     >  2:142078/1‑151 (MQ=255)                                                                                                                                   aaCCCTTGTGCTCAATTGGCAGTTTATTTTCTGCGGAGTAGTCTCTCGTTTCATGGGACCGCTACCACGGAAAGGCAACATGAAACAGAAAATTACGGATTACCTGGACGAAATCTACGGTGGAACATTTACCGCAACTCATTTACAGaaa                <  2:620385/151‑1 (MQ=255)                                                                                                                                        ttGTGCTCAATTGGCAGTTTATTTTCTGCGGAGTAGTCTCTCGTTTCATGGGACCGCTACCACGGAAAGGCAACATGAAACAGAAAATTACGGATTACCTGGACGAAATCTACGGCGGAACATTTACCGCAACTCATTTACAGAAACTTGt           >  1:652236/1‑151 (MQ=255)                                                                                                                                            gCTCAATTGGCAGTTTATTTTCTGCGGAGTAGTCTCTCGTTTCATGGGACCGCTACCACGGAAAGGCAACATGAAACAGAAAATTACGGATTACCTGGACGAAATCTACGGTGGAACATTTACCGCAACTCATTTACAGAAACTTGTAAcg       >  1:539396/1‑151 (MQ=255)                                                                                                                                                 aTTGGCAGTTTATTTTCTGCGGAGTAGTCTCTCGTTTCATGGGACCGCTACCACGGAAAGGCAACATGAAACAGAAAATTACGGATTACCTGGACGAAATCTACGGTGGAACATTTACCGCAACTCATTTACAGAAACTTGTAACGCGTCt  >  2:385085/1‑151 (MQ=255)                                                                                                                                                     |                                                                                                                                                   AATTCACCATGAGCGGCGCTTATGCGCCGTTTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTGGTAAAAAATGATTTCCTGCGTGACTAAAACCCTTGTGCTCAATTGACAGTTTATTTTCTGCGGAGTAGTCTCTCGTTTCATGGGACCGCTACCACGGAAAGGCAACATGAAACAGAAAATTACGGATTACCTGGACGAAATCTACGGTGGAACATTTACCGCAACTCATTTACAGAAACTTGTAACGCGTCT  >  NC_000913/1370007‑1370301

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