Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,370,044	(T)9→8	intergenic (+41/‑172)	pspE → / → ycjM	thiosulfate sulfurtransferase PspE/glucosylglycerate phosphorylase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,370,036	0	T	.	100.0%	148.9 / NA	35	intergenic (+33/‑180)	pspE/ycjM	thiosulfate sulfurtransferase PspE/glucosylglycerate phosphorylase
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base . (18/17);  total (18/17)

ACGCCAGTCAGGGCAAGCAAAAGAGATCCTTAGCGAGATGGGATATACCCACGTTGAGAACGCCGGTGGCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCGTTTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTGGTAAAAAATGATTTCCTGCGTGACTAAAACCCTTGTGCTCAATTGACAGTTTATT  >  NC_000913/1369898‑1370164                                                                                                                                           |                                                                                                                                 aCGCCAGTCAGGGCAAGCAAAAGAGATCCTTAGCGAGATGGGATATACCCACGTTGAGAACGCCGGTGGCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCGtt                                                                                                                                 <  2:1168355/140‑1 (MQ=255)     cAGTCAGGGCAAGCAAAAGAGATCCTTAGCGAGATGGGATATACCCACGTTGAGAACGCCGGTGGCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCGtttttt                                                                                                                             <  2:1107609/140‑1 (MQ=255)       gTCAGGGCAAGCAAAAGAGACCCTTAGCGAGATGGGATATACCCACGTTGAGAACGCCGGTGGCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCGtttttttt                                                                                                                           <  1:157321/140‑1 (MQ=255)            ggCAAGCAAAAGAGATCCTTAGCGAGATGGGATATACCCACGTTGAGAACGCCGGTGGCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCG‑TTTTTTTTCTGTc                                                                                                                     <  1:424566/140‑1 (MQ=255)              caagcaaAAGAGATCCTTAGCGAGATGGGATATACCCACGTTGAGAACGCCGGTGGCCTGAAAGACATTGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCG‑TTTTTTTTCTGTcac                                                                                                                   <  1:181617/140‑1 (MQ=255)               aagcaaAAGAGATCCTTAGCGAGATGGGATATACCCACGTTGAGAACGCCGGTGGCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCG‑TTTTTTTTCTGTcaca                                                                                                                  >  2:698567/1‑140 (MQ=255)                 gcaaAAGAGATCCTTAGCGAGATGGGATATACCCACGTTGAGAACGCCGGTGGCCTGAAAGACATCGCAATGCGGAAGGTCAAAGGTTAAAAGAATTCGCCATGAGCGGCGCTTATGCGCCG‑TTTTTTTTCTGTCACAcc                                                                                                                <  2:922301/140‑1 (MQ=255)                  caaAAGAGATCCTTAGCGAGATGGGATATACCCACGTTGAGAACGCCGGTGGCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCG‑TTTTTTTTCTGTCACACCt                                                                                                               <  1:327214/140‑1 (MQ=255)                    aaaGAGATCCTTAGCGAGATGGGATATACCCACGTTGAGAACGCCGGTGGCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCG‑TTTTTTTTCTGTCACACCTTa                                                                                                             >  1:416035/1‑140 (MQ=255)                      agagATCCTTAGCGAGATGGGATATACCCACGTTGAGAACGCCGGTGGCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCG‑TTTTTTTTCTGTCACACCTTAtt                                                                                                           <  2:773484/140‑1 (MQ=255)                          aTCCTTAGCGAGATGGGATATACCCACGTTGAGAACGCCGGTGGCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCG‑TTTTTTTTCTGTCACACCTTATTTACa                                                                                                       <  1:303035/140‑1 (MQ=255)                          aTCCTTAGCGAGATGGGATATACCCACGTTGAGAACGCCGGTGGCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCG‑TTTTTTTTCTGTCACACCTTATTTACa                                                                                                       <  1:698567/140‑1 (MQ=255)                           tCCTTAGCGAGATGGGATATACCCACGTTGAGAACGCCGGTGGCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCG‑TTTTTTTTCTGTCACACCTTATTTACAt                                                                                                      >  2:81125/1‑140 (MQ=255)                                  cGAGATGGGATATACCCACGTTGAGAACGCCGGTGGCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCG‑TTTTTTTTCTGTCACACCTTATTTACATCCCCATa                                                                                               <  1:81125/140‑1 (MQ=255)                                       tGGGATATACCCACGTTGAGAACGCCGGTGGCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCG‑TTTTTTTTCTGTCACACCTTATTTACATCCCCATAGa                                                                                             <  2:566222/137‑1 (MQ=255)                                       tGGGATATACCCACGTTGAGAACGCCGGTGGCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCG‑TTTTTTTTCTGTCACACCTTATTTACATCCCCATAGa                                                                                             >  1:566222/1‑137 (MQ=255)                                              taCCCACGTTGAGAACGCCGGTGGCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCG‑TTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGt                                                                                   <  1:370683/140‑1 (MQ=255)                                              taCCCACGTTGAGAACGCCGGTGGCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCG‑TTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGt                                                                                   <  2:126450/140‑1 (MQ=255)                                                       tGAGAACGCCGGTGGCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCG‑TTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTcaca                                                                          <  2:1216875/140‑1 (MQ=255)                                                         agaACGCCGGTGGCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGTGGCGCTTATGCGCCG‑TTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACaaa                                                                        >  2:1234521/1‑140 (MQ=255)                                                               ccGGTGGCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCG‑TTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCt                                                                  >  1:137581/1‑140 (MQ=255)                                                                          aaaGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCG‑TTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTGGt                                                       >  2:521010/1‑140 (MQ=255)                                                                                      aTGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCG‑TTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTGGTAAAAAATGAttt                                           <  2:416035/140‑1 (MQ=255)                                                                                       tGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCG‑TTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTGGTAAAAAATGATTTc                                          >  1:517361/1‑140 (MQ=255)                                                                                        gCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCG‑TTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTGGTAAAAAATGATTTcc                                         >  1:976135/1‑140 (MQ=255)                                                                                             aGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCG‑TTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTGGTAAAAAATGATTTCCTGCGt                                    <  2:137581/140‑1 (MQ=255)                                                                                                tCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCG‑TTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTGGTAAAAAATGATTTCCTGCGTGAc                                 <  1:521010/140‑1 (MQ=255)                                                                                                 cAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCG‑TTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTGGTAAAAAATGATTTCCTGCGTGACt                                >  2:780391/1‑140 (MQ=255)                                                                                                         aaaaGAATTCACCATGAGCGGCGCTTATGCGCCG‑TTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTGGTAAAAAATGATTTCCTGCGTGACTAAAACCCt                        >  2:666564/1‑140 (MQ=255)                                                                                                         aaaaGAATTCACCATGAGCGGCGCTTATGCGCCG‑TTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTGGTAAAAAATGATTTCCTGCGTGACTAAAACCCt                        >  1:398591/1‑140 (MQ=255)                                                                                                         aaaaGAATTCACCATGAGCGGCGCTTATGCGCCG‑TTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTGGTAAAAAATCATTTCCTGCGTGACTAAAACCCt                        >  2:231815/1‑140 (MQ=255)                                                                                                          aaaGAATTCACCATGAGCGGCGCTTATGCGCCG‑TTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTGGTAAAAAATGATTTCCTGCGTGACTAAAACCCtt                       >  1:220923/1‑140 (MQ=255)                                                                                                              aaTTCACCATGAGCGGCGCTTATGCGCCG‑TTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTGGTAAAAAATGATTTCCTGCGTGACTAAAACCCTTGTGc                   >  2:970393/1‑140 (MQ=255)                                                                                                                       tGAGTGGCGCTTATGCGCCG‑TTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTGGTAAAAAATGATTTCCTGCGTGACTAAAACCCTTGTGCTCAATTGAc          <  1:1234521/140‑1 (MQ=255)                                                                                                                       tGAGCGGCGCTTATGCGCCG‑TTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTGGTAAAAAATGATTTCCTGCGTGACTAAAACCCTTGTGCTCAATTGAc          <  1:1067971/140‑1 (MQ=255)                                                                                                                        gAGCGGCGCTTATGCGCCG‑TTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTGGTAAAAAATGATTTCCTGCGTGACTAAAACCCTTGTGCTCAATTGACa         >  2:685370/1‑140 (MQ=255)                                                                                                                               gcTTATGCGCCG‑TTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTGGTAAAAAATGATTTCCTGCGTGACTAAAACCCTTGTGCTCAATTGACAGTTTAtt  >  2:56964/1‑140 (MQ=255)                                                                                                                               gcTTATGCGCCG‑TTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTGGTAAAAAATGATTTCCTGCGTGACTAAAACCCTTGTGCTCAATTGACAGTTTAtt  >  1:758517/1‑140 (MQ=255)                                                                                                                                           |                                                                                                                                 ACGCCAGTCAGGGCAAGCAAAAGAGATCCTTAGCGAGATGGGATATACCCACGTTGAGAACGCCGGTGGCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCGTTTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTGGTAAAAAATGATTTCCTGCGTGACTAAAACCCTTGTGCTCAATTGACAGTTTATT  >  NC_000913/1369898‑1370164

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: ‑