Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,661,961	G→T	S57R (AGC→AGA)	iscR ←	DNA‑binding transcriptional dual regulator IscR

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,661,961	0	G	T	100.0%	74.9 / NA	24	S57R (AGC→AGA)	iscR	DNA‑binding transcriptional dual regulator IscR
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base T (10/14);  total (10/14)

GACACATCCAGCACTTCCTGGTTATTAACCAGTTCGCCTAAAGTAATGTTGTTGAGAAAACCGGTGAGACGGTCGCTCAAATCACGCCACAGCGCGTGGGTCAGGCATTTATCGCCGCCCTGGCAGCCGCCTTTACCCTGACAACGGGTGGCATCTACAGATTCGTCAACGGCGCTAATTAC‑TTCGCCAACGGCGATGCTGCTGGCATCTTTGCCTAACAGATAACCACCGCCTGGTCCACGTACGCTGGAAACCAGACCATTTTTACGCAGACGGGAAAACAGTTGTTCCAGATAAGAAAGGGAAATTCCCTGACGTTCGGAAATATCAGCCAACGGTACCGGGCCCGCTTCAGAGTTGAGCGCAACGTCAAGCATTGCGGTCACGGCATAGCGCCCTTTAGATGTCAGTCTCATGTCTTACTTCACCTCAAACTCGCCCCTGCCCGGGGTTTTTTATTGTAAAAGTGGGGGTATTGCATAGC  >  NC_000913/2661716‑2662199                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                               gACACATCCAGCACTTCCTGGTTATTAACCAGTTCGCCTACAGTAATGTTGTTGAGAAAACCGGTGAGACGGTCGCTCAAATCACGCCACAGCGCGTGGGTCAGGCATTTATCGCCGCCCTGGCAGCCGCCTTTACCCTGACAACGGGTGGCATCTACAGATTCGTCAACGGCGCTAATTAC‑TTCGCCAACGGCGATGCTGCAGGCATCTTTGCCTAACAGATAACCCCCGCCTGGTCCACGTACTCTgg                                                                                                                                                                                                                                            >  2:97471/1‑250 (MQ=255) gACACATCCAGCACTTCCTGGTTATTAACCAGTTCGCCTAAAGTAATGTTGTTGAGAAAACCGGTGAGACGGTCGCTCAAATCACGCCACAGCGCGTGGGTCAGGCATTTATCGCCGCCCTGGCAGCCGCCTTTACCCTGACAACGGGTGGCATCTACAGATTCGTCAACGGCGCTAATTAC‑TTCGCCAACGGCGATGCTGCTGGCATCTTTGCCTAACAGATAACCACCGCCTGGTCCACGTACTCTGGa                                                                                                                                                                                                                                           <  1:45871/251‑1 (MQ=255)      aTCCAGCACTTCCTGGTTATTAACCAGTTCGCCTAAAGTAATGTTGTTGAGAAAACCGGTGAGACGGTCGCTCAAATCACGCCACAGCGCGTGGGTCAGGCATTTATCGCCGCCCTGGCAGCCGCCTTTACCCTGACAACGGGTGGCATCTACAGATTCGTCAACGGCGCTAATTAC‑TTCGCCAACGGCGATGCTGCTGGCATCTTTGCCTAACAGATAACCACCGCCTGGTCCACGTACTCTGGAAAcc                                                                                                                                                                                                                                       >  2:296788/1‑250 (MQ=255)                tctttgTTATTAACCAGTTCGCCTAAAGTAATTTTGTTGAGAAAACCGGGGAGCAGGTCGCTCAAATCACGCCACAGCGCTGGGGTCAGGCGTTTATCCCCGCCCGGGCGGCCGCCTTTCCCCTGACAGCGGGTGGCATCTCCAGATTCGTCAACGGCGCTAATTAC‑TTCGCCAACGGCGATGCTGCTGGCATCTTTGCCTAACAGATAACCACCGCCTGGTCCACGTACTCTGGAAACCATACCAtttt                                                                                                                                                                                                                             <  2:266047/245‑1 (MQ=255)                        atgaaCCAGTTCGCCAAAATAAATGTTGTTGAGAAAACCGGTGAGACGGTCGCTCAAATCACGCCACAGCGCGTGGGTCAGGCATTTATCGCCGCCCTGGCAGCCGCCTTTACCCTGACAACGGGTGGCATCTACAGATTCGTCAACGGCGCTAATTAC‑TTCGCCAACGGCGATGCTGCTGGCATCTTTGCCTAACAGATAACCACCGCCTGGTCCACGTACTCTGGAAACCAGACCATTTTTACGCAGa                                                                                                                                                                                                                     <  1:343632/247‑1 (MQ=255)                             ccAGTTCGCCTAAAGTAATGTTGTTGAGAAAACCGGTGAGACGGTCGCTCAAATCACGCCACAGCGCGTGGGTCAGGCATTTATCGCCGCCCTGGCAGCCGCCTTTACCCTGACAACGGGTGGCATCTCCAGATTCGTCAACGGCGCTAATTAC‑TTCGCCAACGGCGATGCTGCTGGCATCTTTGCCTAACAGATAACCACCGCCTGGTCCACGTACTCTGGAAACCAGACCATTTTTACGCAGCCGGGa                                                                                                                                                                                                                >  2:36870/1‑250 (MQ=255)                              cAGTTCGCCAAAAGTAATGTTTTTGAGAAAACCGGTGAGACGGTCGCCCAAAACACGCCACAGCGCGGGGGTCAGGCATTTATCGCCGCCCTGGCAGCCGCCTTTACCCTGACAACGGGTGGCATCTACAGATTCGTCAACGGCGCTAATTAC‑TTCGCCAACGGCGATGCTGCTGGCATCTTTGCCTAACAGATAACCACCGCCTGGTCCACGTACTCTGGAAACCAGACCATTTTTACGCAGACGGGaaa                                                                                                                                                                                                              <  2:360782/251‑1 (MQ=255)                                    gCCTAAAGTAATGTTGTTGAGAAAACCGGTGAGACGGTCGCTCAAATCACGCCACAGCGCGTGGGTCAGGCATTTATCGCCGCCCTGGCAGCCGCCTTTACCCTGACAACGGGTGGCATCTACAGATTCGTCAACGGCGCTAATTAC‑TTCGCCAACGGCGATGCTGCTGGCATCTTTGCCTAACAGATAACCACCGCCTGGTCCACGTACTCTGGAAACCAGACCATTTTTACGCAGACGGGAAAACAgtt                                                                                                                                                                                                        <  1:97471/251‑1 (MQ=255)                                                tgttgggagaACACCGGGGAGACGGTCGCCCAAATCACGCCACAGCGCGTGGGTCAGGCATTTATCGCCGCCCTGGCAGCCGCCTTTACCCTGACAACGGGTGGCATCTACAGAGTCGTCATCGCCGCTAATTAC‑TTCGCCAACGGCGATGCTGCTGGCATCGTGGCCTAACAGATAACCACCGCCTGTTCCACGTACTCTGGAAACCAGACCATTTTTACGCAGACGGGAAAACAGTTGTTCCAGATaa                                                                                                                                                                                             <  2:69960/244‑1 (MQ=255)                                                   gttgAGAAAACCGGTGAGACGGTCGCTCAAATCACGCCACAGCGCGTGGGTCAGGCATTTATCGCCGCCCTGGCAGCCGCCTTTACCCTGACAACGGGTGGCATCTACAGATTCGTCAACGGCGCTAATTAC‑TTCGCCAACGGCGATGCTGCTGGCATCTTTGCCTAACAGATAACCACCGCCTGGTCCACGTACTCTGGAAACCAGACCATTTTTACGCAGACGGGAAAACAGTTGTTCCAGATAAGaa                                                                                                                                                                                          <  1:282843/250‑1 (MQ=255)                                                           aaCCGGTGAGACGGTCGCCCAAATCACGCCACAGCGCGTGGGTCAGGCATTTATCGCCGCCCTGGCAGCCGCCTTTACCCTGACAACGGGTGGCATCTACAGATTCGTCAACGGCGCTAATTAC‑TTCGCCAACGGCGATGCTGCTGGCATCTTTGCCTAACAAATAACCACCGCCTGGTCCACGTACTCTGGAAACCAGACCATTTTTACGCAGACGGGAAAACAGTTGTTCCAGATAAGAAAGGGAAAtt                                                                                                                                                                                 <  1:253146/251‑1 (MQ=255)                                                                                               cgTGGGTCAGGCATTTATCGCCGCCCTGGCAGCCGCCTTTACCCTGACAACGGGTGGCATCTACAGATTCGTCAACGGCGCTAATTAC‑TTCGCCAACGGCGATGCTGCTGGCATCTTTGCCTAACAGATAACCACCGCCTGGTCCACGTACTCTGGAAACCAGACCATTTTTACGCAGACGGGAAAACAGTTGTTCCAGATAAGAAAGGGAAATTCCCTGACGTTCGGAAATATCAGCCAACGGTACCggg                                                                                                                                             >  1:73684/1‑251 (MQ=255)                                                                                                 ggtgtcagtcaTTGAGCGCCGGCCTGGCAGCCGCCTTTACCCTGACAACGGGTGGCACTTACAGTTTCGTCAACGGCGCGAATTTC‑TTCGCCAACGGTGATGCTGCTGGCATTTTTGCCATACAGATAACCACCGCCTGGTCCACGTACTCTGGAAACCAGACCATTTTTACGCAGACTGGAAAACAGTTGTTCCAGAGAAGAAAGGGAAATTTCCTGACGTTCGGAAATATCAGCCAACGGTACCGGGcc                                                                                                                                           <  1:136836/248‑1 (MQ=255)                                                                                                    gTCAGGCATTTATCGCCGCCCTGGCAGCCGCCTTTACCCTGACAACGGGTGGCATCTACAGATTCGTCAACGGCGCTAATTAC‑TTCGCCAACGGCGATGCTGCTGGCATCTTTGCCTAACAGATAACCACCGCCTGGTCCACGTACTCTGGAAACCAGACCATTTTTACGCAGACGGGAAAACAGTTGTTCCAGATAAGAAAGGGAAATTCCCTGACGTTCGGAAATATCAGCCAACGGTACCGGGCCCGc                                                                                                                                        >  2:145963/1‑251 (MQ=255)                                                                                                            tttATCGCCGCCCTGGCAGCCGCCTTTACCCTGACAACGGGTGGCATCTACAGATTCGTCAACGGCGCTAATTAC‑TTCGCCAACGGCGATGCTGCTGGCATCTTTGCCTAACAGATAACCACCGCCTGGTCCACGTACTCTGGAAACCAGACCATTTTTACGCAGACGGGAAAACAGTTGTTCCAGATAAGAAAGGGAAATTCCCTGACGTTCGGAAATATCAGCCAACGGTACCGGGCCCGCTTCAGAGt                                                                                                                                >  2:220777/1‑251 (MQ=255)                                                                                                                            cAGCCGCCTTTACCCTGACAACGGGTGGCATCTACAGATTCGTCAACGGCGCTAATTAC‑TTCGCCAACGGCGATGCTGCTGGCATCTTTGCCTAACAGATAACCACCGCCTGGTCCACGTACTCTGGAAACCAGACCATTTTTACGCAGACGGGAAAACAGTTGTTCCAGATAAGAAAGGGAAATTCCCTGACGTTCGGAAATATCAGCCAACGGTACCGGGCCCGCTTCAGAGTTGAGCGCAACGTCa                                                                                                                  >  2:9969/1‑249 (MQ=255)                                                                                                                            cAGCCGCCTTTACCCTGACAACGGGTGGCATCTACAGATTCGTCAACGGCGCTAATTAC‑TTCGCCAACGGCGATGCTGCTGGCATCTTTGCCTAACAGATAACCACCGCCTGGTCCACGTACTCTGGAAACCAGACCATTTTTACGCAGACGGGAAAACAGTTGTTCCAGATAAGAAAGGGAAATTCCCTGACGTTCGGAAATATCAGCCAACGGTACCGGGCCCGCTTCAGAGTTGAGCGCAACGTCAAg                                                                                                                <  1:272162/251‑1 (MQ=255)                                                                                                                                                                    cGTCAACGGCTCTAATTCC‑TTCGCCAACGGCGATGCTGCTGGCATCTTTGCCTAACAGATAACCACCGCCTGGTCCACGTACTCTGGAAACCAGACCATTTGTACGCAGACGGGAAAACAGTTGTTCCAGATAAGAAAGGGAAATTCCCTGACGTTCGGAAATATCAGCCAACGGTACCGGGCCCGCTTCAGAGTTGAGCGCAACGTCAAGCATTGCGGTCACGGCATAGCGCCCTTTAGATGTCAGtctc                                                                        <  1:64024/251‑1 (MQ=255)                                                                                                                                                                    cGTCAACGGCG‑TAAATACCATCGCCAACGGCGATGCTGCTGGCATCTTTGCCTAACAGATAACCACCGCCTGGTCCACGTACTCTGGAAACCAGACCATTTTTACGCAGACGGGAAAACAGTTGTTCCAGATAAGAAAGGGAAATTCCCTGACGTTCGGAAATATCAGCCAACGGTACCGGGCCCGCTTCAGAGTTGAGCGCAACGTCAAGCATTGCGGTCACGGCATAGCGCCCTTTAGATGTCAGtct                                                                         <  2:92278/250‑1 (MQ=255)                                                                                                                                                                      tCAACGGCGCTAATTAC‑TTCGCCAACGGCGATGCTGCTGGCATCTTTGCCTAACAGATAACCACCGCCTGGTCCACGTACTCTGGAAACCAGACCATTTTTACGCAGACGGGAAAACAGTTGTTCCAGATAAGAAAGGGAAATTCCCTGACGTTCGGAAATATCAGCCAACGGTACCGGGCCCGCTTCAGAGTTGAGCGCAACGTCAAGCATTGCGGTCACGGCATAGCGCCCTTTAGATGTCAGTCTCAt                                                                      >  1:263348/1‑251 (MQ=255)                                                                                                                                                                              ccTAATTAC‑TTCGCCAACGGCGATGCTGCTGGCATCTTTGCCTAACAGATAACCACCGCCTGGTCCACGTACTCTGGAAACCAGACCATTTTTACGCAGACGGGAAAACAGTTGTTCCAGATAAGAAAGGGAAATTCCCTGACGTTCGGAAATATCAGCCAACGGTACCGGGCCCGCTTCAGAGTTGAGCGCAACGTCAAGCATTGCGGTCACGGCATAGCGCCCTTTAGATGTCAGTCTCATGTCTTACt                                                              <  1:94924/250‑1 (MQ=255)                                                                                                                                                                                                   cGATGCTGCTGGCATCTTTGCCTAACAGATAACCACCGCCTGGTCCACGTACTCTGGAAACCAGACCATTTTTACGCAGACGGGAAAACAGTTGTTCCAGATAAGAAAGGGAAATTCCCTGACGTTCGGAAATATCAGCCAACGGTACCGGGCCCGCTTCAGAGTTGAGCGCAACGTCAAGCATTGCGGTCACGGCATAGCGCCCTTTAGATGTCAGTCTCATGTCTTACTTCACCTCAAACTCGCCCCTg                                          <  1:296788/251‑1 (MQ=255)                                                                                                                                                                                                                          aaCAGATACCCACCGCCTGATCCACGTACTCTGGAAACCAGACCATTTTTACGCAGACGGGAAAACCGTTGTTCCAGATAAGAAAGGGAAATTCCCTGACGTTCGGAAATATAAGCCAACGGTACCGGGCCCGCTTCAGAGTTGAGCGCAACGTCCAGCATTGCGGTCACGGCATAGCGCCCTTTAGCTGTCAGTCTCATGTCTTACTTCACCTCAAACTCGCCCCTGCCCGGGGTTTTTTATTGTaaa                     >  2:227139/1‑249 (MQ=255)                                                                                                                                                                                                                                            ggTCCACGTACTCTGGAAACCAGACCATTTTTACGCAGACGGGAAAACAGTTGTACCAGATAAGAAAGGGAAATTCCCTGACGTTCGGAAATATCAGCCAACGGTACCGGGCCCGCTTCAGAGTTGAGCGCAACGTCAAGCATTGCGGTCACGGCATAGCGCCCTTTAGATGTCAGTCTCATGTCTTACTTCACCTCAAACTCGCCCCTGCCCGGGGTTTTTTATTGTAAAAGTGGGGGTATTGCATAGc  >  1:349073/1‑250 (MQ=255)                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                               GACACATCCAGCACTTCCTGGTTATTAACCAGTTCGCCTAAAGTAATGTTGTTGAGAAAACCGGTGAGACGGTCGCTCAAATCACGCCACAGCGCGTGGGTCAGGCATTTATCGCCGCCCTGGCAGCCGCCTTTACCCTGACAACGGGTGGCATCTACAGATTCGTCAACGGCGCTAATTAC‑TTCGCCAACGGCGATGCTGCTGGCATCTTTGCCTAACAGATAACCACCGCCTGGTCCACGTACGCTGGAAACCAGACCATTTTTACGCAGACGGGAAAACAGTTGTTCCAGATAAGAAAGGGAAATTCCCTGACGTTCGGAAATATCAGCCAACGGTACCGGGCCCGCTTCAGAGTTGAGCGCAACGTCAAGCATTGCGGTCACGGCATAGCGCCCTTTAGATGTCAGTCTCATGTCTTACTTCACCTCAAACTCGCCCCTGCCCGGGGTTTTTTATTGTAAAAGTGGGGGTATTGCATAGC  >  NC_000913/2661716‑2662199

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