Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,657,434	G→A	S501F (TCT→TTT)	hscA ←	iron‑sulfur cluster biosynthesis chaperone HscA

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,657,434	0	G	A	100.0%	94.2 / NA	28	S501F (TCT→TTT)	hscA	iron‑sulfur cluster biosynthesis chaperone HscA
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (18/10);  total (18/10)

CCAGCACACGCGCCGCTTCTACTTTTTGTTCTGCCAGCATTCGGGCTTTTACGTCCTGCTCGGCATAGCTCATTGAGTCTTTGATCATCGAAGCGATTTCGCTATCGGTCAGACCGTAAGACGGTTTGACCTGAATAGACGCCTCAACGCCGGTGGATTTCTCCATCGCCGTCACGCTCAAAAGACCGTCGGCATCGACCTGGAACGTCACGCGAATATGCGCACCGCCAGCCGGTAG  >  NC_000913/2657315‑2657552                                                                                                                        |                                                                                                                       ccAGCACACGCGCCGCTTCTACTTTTTGTTCTGCCAGCATTCGGGCTTTTACGTCCTGCTCGGCATAGCTCATTGAGTCTTTGATCATCGAAGCGATTTCGCTATCGGTCAGACCGTAAAACGGTTTGACCTGAATAGAc                                                                                                    >  1:36888/1‑140 (MQ=255) ccAGCACACGCGCCGCTTCTACTTTTTGTTCTGCCAGCATTCGGGCTTTTACGTCCTGCTCGGCATAGCTCATTGAGTCTTTGATCATCGAAGCGATTTCGCTATCGGTCAGACCGTAAAACGGTTTGACCTGAATAGAc                                                                                                    >  1:635190/1‑140 (MQ=255)      acacGCGCCGCTTCTACTTTTTGTTCTGCCAGCATTCGGGCTTTTACGTCCTGCTCGGCATAGCTCATTGAGTCTTTGATCATCGAAGCGATTTCGCTATCGGTCAGACCGTAAAACGGTTTGACCTGAATAGACGCCTc                                                                                               >  1:648184/1‑140 (MQ=255)               gcTTCTACTTTTTGTTCTGCCAGCATTCGGGCTTTTACGTCCTGCTCGGCATAGCTCATTGAGTCTTTGATCATCGAAGCGATTTCGCTATCGGTCAGACCGTAAAACGGTTTGACCTGAATAGACGCCTCAACGCCGGt                                                                                      >  1:21897/1‑140 (MQ=255)                                tGCCAGCATTCGGGCTTTTACGTCCTGCTCGGCATAGCTCATTGAGTCTTTGATCATCGAAGCGATTTCGCTATCGGTCAGACCGTAAAACGGTTTGACCTGAATAGACGCCTCAACGCCGGTGGATTTCTCCATCGCCg                                                                     >  2:649999/1‑140 (MQ=255)                                  ccAGCATTCGGGCTTTTACGTCCTGCTCGGCATAGCTCATTGAGTCTTTGATCATCGAAGCGATTTCGCTATCGGTCAGACCGTAAAACGGTTTGACCTGAATAGACGCCTCAACGCCGGTGGATTTCTCCATCGCCGTc                                                                   <  2:753548/140‑1 (MQ=255)                                       aTTCGGGCTTTTACGTCCTGCTCGGCATAGCTCATTGAGTCTTTGATCATCGAAGCGATTTCGCTATCGGTCAGACCGTAAAACGGTTTGACCTGAATAGACGCCTCAACGCCGGTGGATTTCTCCATCGCCGTCACGCt                                                              >  1:522025/1‑140 (MQ=255)                                           gggCTTTTACGTCCTGCTCGGCATAGCTCATTGAGTCTTTGATCATCGAAGCGATTTCGCTATCGGTCAGACCGTAAAACGGTTTGACCTGAATAGACGCCTCAACGCCGGTGGATTTCTCCATCGCCGTCACGCTCaaa                                                          >  2:626363/1‑140 (MQ=255)                                                       ccTGCTCGGCATAGCTCATTGAGTCTTTGATCATCGAAGCGATTTCGCTATCGGTCAGACCGTAAAACGGTTTGACCTGAATAGACGCCTCAACGCCGGTGGATTTCTCCATCGCCGTCACGCTCAAAAGACCGt                                                   >  2:83732/1‑135 (MQ=255)                                                       ccTGCTCGGCATAGCTCATTGAGTCTTTGATCATCGAAGCGATTTCGCTATCGGTCAGACCGTAAAACGGTTTGACCTGAATAGACGCCTCAACGCCGGTGGATTTCTCCATCGCCGTCACGCTCAAAAGACCGt                                                   <  1:83732/135‑1 (MQ=255)                                                        cTGCTCGGCATAGCTCATTGAGTCTTTGATCATCGAAGCGATTTCGCTATCGGTCAGACCGTAAAACGGTTTGACCTGAATAGACGCCTCAACGCCGGTGGATTTCTCCATCGCCGTCACGCTCAAAAGACCGTCGGCAt                                             <  1:649999/140‑1 (MQ=255)                                                                     cTCATTGAGTCTTTGATCATCGAAGCGATTTCGCTATCGGTCAGACCGTAAAACGGTTTGACCTGAATAGACGCCTCAACGCCGGTGGATTTCTCCATCGCCGTCACGCTCAAAAGACCGTCGGCATCGACCTGGAACGt                                >  2:472906/1‑140 (MQ=255)                                                                          tGAGTCTTTGATCATCGAAGCGATTTCGCTATCGGTCAGACCGTAAAACGGTTTGACCTGAATAGACGCCTCAACGCCGGTGGATTTCTCCATCGCCGTCACGCTCAAAAGACCGTCGGCATCGACCTGGAACGTCAcgc                           >  1:90781/1‑140 (MQ=255)                                                                          tGAGTCTTTGATCATCGAAGCGATTTCGCTATCGGTCAGACCGTAAAACGGTTTGACCTGAATAGACGCCTCAACGCCGGTGGATTTCTCCATCGCCGTCACGCTCAAAAGACCGTCGGCATCGACCTGGAACGTCAcgc                           >  2:37629/1‑140 (MQ=255)                                                                           gAGTCTTTGATCATCGAAGCGATTTCGCTATCGGTCAGACCGTAAAACGGTTTGACCTGAATAGACGCCTCAACGCCGGTGGATTTCTCCATCGCCGTCACGCTCAAAAGACCGTCGGCATCGACCTGGAACGTCAcgcg                          <  2:90781/140‑1 (MQ=255)                                                                               cTTTGATCATCGAAGCGATTTCGCTATCGGTCAGACCGTAAAACGGTTTGACCTGAATAGACGCCTCAACGCCGGTGGATTTCTCCATCGCCGTCACGCTCAAAAGACCGTCGGCATCGACCTGGAACGTCACGCGAata                      >  1:721021/1‑140 (MQ=255)                                                                                  tGATCATCGAAGCGATTTCGCTATCGGTCAGACCGTAAAACGGTTTGACCTGAATAGACGCCTCAACGCCGGTGGATTTCTCCATCGCCGTCACGCTCAAAAGACCGTCGGCATCGACCTGGAACGTCACGCGAATATgc                   >  2:239100/1‑140 (MQ=255)                                                                                   gATCATCGAAGCGATTTCGCTATCGGTCAGACCGTAAAACGGTTTGACCTGAATAGACGCCTCAACGCCGGTGGATTTCTCCATCGCCGTCACGCTCAAAAGACCGTCGGCATCGACCTGGAACGTCACGCGAATATgcg                  <  2:620692/140‑1 (MQ=255)                                                                                   gATCATCGAAGCGATTTCGCTATCGGTCAGACCGTAAAACGGTTTGACCTGAATAGACGCCTCAACGCCGGTGGATTTCTCCATCGCCGTCACGCTCAAAAGACCGTCGGCATCGACCTGGAACGTCACGCGAATATgcg                  <  1:472906/140‑1 (MQ=255)                                                                                   gATCATCGAAGCGATTTCGCTATCGGTCAGACCGTAAAACGGTTTGACCTGAATAGACGCCTCAACGCCGGTGGATTTCTCCATCGCCGTCACGCTCAAAAGACCGTCGGCATCGACCTGGAACGTCACGCGAATATgcg                  <  1:37629/140‑1 (MQ=255)                                                                                      catcGAAGCGATTTCGCTATCGGTCAGACCGTAAAACGGTTTGACCTGAATAGACGCCTCAACGCCGGTGGATTTCTCCATCGCCGTCACGCTCAAAAGACCGTCGGCATCGACCTGGAACGTCACGCGAATATGCGCAc               >  2:28716/1‑140 (MQ=255)                                                                                        tcGAAGCGATTTCGCTATCGGTCAGACCGTAAACCGGTTTGACCTGAATAGACGCCTCAACGCCGGTGGATTTCTCCATCGCCGTCACGCTCAAAAGACCGTCGGCATCGACCTGGAACGTCACGCGAATATGCGCACCg             >  1:531365/1‑140 (MQ=255)                                                                                        tcGAAGCGATTTCGCTATCGGTCAGACCGTAAAACGGTTTGACCTGAATAGACGCCTCAACGCCGGTGGATTTCTCCATCGCCGTCACGCTCAAAAGACCGTCGGCATCGACCTGGAACGTCACGCGAATATGCGCACCg             >  2:300475/1‑140 (MQ=255)                                                                                           aaGCGATTTCGCTATCGGTCAGACCGTAAAACGGTTTGACCTGAATAGACGCCTCAACGCCGGTGGATTTCTCCATCGCCGTCACGCTCAAAAGACCGTCGGCATCGACCTGGAACGTCACGCGAATATGCGCACCGCCa          <  1:774199/140‑1 (MQ=255)                                                                                                 tttCGCTATCGGTCAGACCGTAAAACGGTTTGACCTGAATAGACGCCTCAACGCCGGTGGATTTCTCCATCGCCGTCACGCTCAAAAGACCGTCGGCATCGACCTGGAACGTCACGCGAATATGCGCACCGCCAGCCGGt    <  2:807060/140‑1 (MQ=255)                                                                                                   tCGCTATCGGTCAGACCGTAAAACGGTTTGACCTGAATAGACGCCTCAACGCCGGTGGATTTCTCCATCGCCGTCACGCTCAAAAGACCGTCGGCATCGACCTGGAACGTCACGCGAATATGCGCACCGCCAGCCGGTAg  >  2:386124/1‑140 (MQ=255)                                                                                                              cAGACCGTAAAACGGTTTGACCTGAATAGACGCCTCAACGCCGGTGGATTTCTCCATCGCCGTCACGCTCAAAAGACCGTCGGCATCGACCTGGAACGTCACGCGAATATGCGCACCGCCAGCCGGt    >  2:228784/1‑127 (MQ=255)                                                                                                              cAGACCGTAAAACGGTTTGACCTGAATAGACGCCTCAACGCCGGTGGATTTCTCCATCGCCGTCACGCTCAAAAGACCGTCGGCATCGACCTGGAACGTCACGCGAATATGCGCACCGCCAGCCGGt    <  1:228784/127‑1 (MQ=255)                                                                                                                        |                                                                                                                       CCAGCACACGCGCCGCTTCTACTTTTTGTTCTGCCAGCATTCGGGCTTTTACGTCCTGCTCGGCATAGCTCATTGAGTCTTTGATCATCGAAGCGATTTCGCTATCGGTCAGACCGTAAGACGGTTTGACCTGAATAGACGCCTCAACGCCGGTGGATTTCTCCATCGCCGTCACGCTCAAAAGACCGTCGGCATCGACCTGGAACGTCACGCGAATATGCGCACCGCCAGCCGGTAG  >  NC_000913/2657315‑2657552

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