Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,661,821	C→G	C104S (TGC→TCC)	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,821	0	C	G	100.0%	63.3 / NA	26	C104S (TGC→TCC)	iscR	DNA‑binding transcriptional dual regulator IscR
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base G (13/13);  total (13/13)

TGAGTATGCTGACGACCAGACACATCCAGCACTTCCTGGTTATTAACCAGTTCGCCTAAAGTAATGTTGTTGAGAAAACCGGTGAGACGGTCGCTCAAATCACGCCACAGCGCGTGGGTCAGGCATTTATCGCCGCCCTGGCAGCCGCCTTTACCCTGACAACGGGTGGCATCTACAGATTCGTCAACGGCGCTAATTACTTCGCCAACGGCGATGCTGCTGGCATCTTTGCCTAACAGATAACCACC  >  NC_000913/2661698‑2661945                                                                                                                            |                                                                                                                             ggAGCATGCTGACGACCAGACACAGCCAGCACTTCCTGGTTATTAACCAGTTCGCCTAAAGTAATGTTGTTGAGAAAACCGGTGAGACTGTCGCTCAAATCACGCCACAGCGCGTGCGTCATGGATTTa                                                                                                                         <  2:306424/128‑1 (MQ=255)    gTATGCTGACGACCCGACACCTCCAGCACTTCCTGGTTATTAACCAGTTCGAATAAATTAATGTTGTTGAGAAAACCGGTGAGACGGTCGCTCAAATCACGCCACAGCGCGTGGGTCAGGGATTTATcg                                                                                                                      >  1:506932/1‑129 (MQ=255)         cTGACGACCAGCCACATCCAGCACTTCCTGGATATTAACAAGAACGCAAAAAGAAATGTAATAGAGAAAACCGGTGAGACGGTCGCTCAAATCACGCCACAGCGCGTGGGTCAGGGATTTATCgccgcc                                                                                                                 >  1:329898/1‑129 (MQ=255)            ccgacCAGACACATCCAGCACTTCCTGGTTATTAACCAGTTCGCCTAAAATAATGTTGTTGAGAAAACCAGTGAGCCCGTCGCTCAAATCACGCCCCAGCGCGTGCGTCCGAGATGTATCGCCGCCCTg                                                                                                              >  2:257917/2‑129 (MQ=255)                  aGACACATCCAGCACTTCCTGGTTATTAACCAGTTCGCCTAAAGTAATATTGTTGAGAAAACCGGTGAGACGGTCGCGCAAATCACGCCACAGCGCGTCGGTCAGGGATGCATCGCCGCCCTTGCAccc                                                                                                        >  2:459470/1‑129 (MQ=255)                        aTCCAACACTTCCTGGTTATAAACCAGTTCGCCTAAAGTAATGATGTTGAGAAAACCGGTGAGACGGTCGCTCAAATCACGCCACACCGCGAGGGCCAGCGATATATCGCCGACCCGGCCGACGCCatt                                                                                                  >  2:661288/1‑129 (MQ=255)                            aGCACTTCCTGGTTATTATCCAGTTCGCCTAAAGTAATGTTGTTGAGAAAACCGGTGAGACGGTCGCTCAAATCACGCCACAGCGCGTGGGTCAGGGATTTATCGCCGCCCTGGCAGCCGCCTTTAccc                                                                                              <  1:259037/129‑1 (MQ=255)                            aGCACTTCCTGGTTATTAACCAGTTCGACTAAAGTAATGTTGTTGAGAAAACCGGTGAGACGGTCGCTCAAATCACGCCACAGCGCGTGGGTCAGGGATGTATCGCCGCCCTGACAGCCGCCTTTAccc                                                                                              >  2:416236/1‑129 (MQ=255)                               aCTTCCTGGTTATTATCCATTTCGCCTAAAGTAATGTTGTTGAGAAAACCGGTGAGACGGTCGCTCAAATCACGCCACAGCGCGTGGGTCAGGGATTTATCGCCGCCCTGGCAGCCGCCTTTACCCTGa                                                                                           <  2:233767/129‑1 (MQ=255)                                                caggtggCCTAAAGTAAGGTTGTTGAGAAAACCGGTGAGACGGTCGCTCAAATCACGCCACAGCGCGTGGGTCAGGGATTTATTGCCGTTTTGGCAGCCGCCTTTACCCTGACAACGGGTGGCATCTAc                                                                          <  1:491375/123‑1 (MQ=255)                                                    tCGCCTAAACTAATGTTGTTGAGAAAACCGGTGAGACGGTCGCTCAAATCACGCCACAGCGCGTGGGTCAGGGATTTATCGTCGCCCTGGCAGCCGCCTTTACCCTGACAACGGGTGGCATCTACAGAt                                                                      <  2:105826/129‑1 (MQ=255)                                                       cctgtaGTAATGTTGTTGAGAAAACGCTTGTTTCGTACGAACAAGTCTGTGCACAGCGCGTGGGTCAGTGATTGATCTCTGCCCTGGCAGCCTCCTTTACCCTGACAACGGGTGGCATCTACAGATTCg                                                                   <  2:355594/124‑1 (MQ=255)                                                        cTAAAGTAATGTTGTTGAGAAACCCGGTGTGACGGTCGCTCAAATCACGCCACAGCGCGTGGGTCAGGGATTTATCGCCGCTCTGTTTGTTTTTTTTTTCCTGACAACGGGTGGCATCTACAGATTCGt                                                                  <  1:159004/129‑1 (MQ=255)                                                                       tgAGAAAACCGGTGAGACGGTCGCTCAAATCACACCACAACGCATAGAACAGAGATTTATCGACACCCTGGCAGCCGCCTTTACCCTGACAACGGGTGGCATCTACAGATTCGTCAACGGCGCTAATTa                                                   >  1:559372/1‑129 (MQ=255)                                                                                          gTCGCTCAAATCACGCCACAGCGCGTGGGACAGGGATTTATCGCCGCCCTGGCAGCCGCCTTTACCCTGACAACGGGTGGCATCTACAGATTCGTCAACGGCGCTAATTACTTCGCCAACGGCGAtgct                                >  1:32335/1‑129 (MQ=255)                                                                                               tCAAATCAGGCCACAGCGCGTGGGTCAGGGATTTATCGCCGCCCTGGCAGCCGCCTTTACCCTGACAACGTGTGGCATCTACAGATTTTTTAACGGTGTTAATTACTTCGCCAACGGCGATGCTGCTgg                           <  1:43025/129‑1 (MQ=255)                                                                                                  aaTCACGCCACAGCGCGTGGGTCAGGGATTTATAGCCTCGCGGGCAGCCGCCTTTACCCTGACAACGGGTGGCATCTACAGATTCGTCAACGGTGCTAATTACTTCGCCAACGGCGATGCTGCTGGCAt                        <  1:417689/129‑1 (MQ=255)                                                                                                    tCACGCCACAGCGCGTGGGTCAGGGATTTATCGCCGCCCTGGCAGCCGCCTTTACCCTGACAACGGGTGCCATCTACAGATTCGTCAACGGCGCTAATTACTTCGCCAACGGCGATGCGGATGGccgct                      >  2:92478/1‑125 (MQ=255)                                                                                                     gaCGCCACAGGGCGTGGGTCAGGGATTTATCGCCGCCCTGGCAGCCGCCTTTACCCTGACAACGGGTGGCATCTACAGATTCGTCAACGGCGCTAATTACTTCGCCAACGGCGATGCTGCTGGCATCtt                     <  2:590196/128‑1 (MQ=255)                                                                                                     cACGCCACAGCGCGTGGGTCAGGGATTTATCGCCGCCCTGGCAGCAGCCTTTACCCTGACAACGGGTGGCATCTACAGCTTCGTCAACGCCGCTCATTACTTCGCCACACGCCATGCTGCTGGCATcgt                     >  2:387517/1‑127 (MQ=255)                                                                                                      aCGCCACAGCGCGTGGGTCAGGGATTTATCGCCGCCCTGGCAGCCGCCTTTACCCTGACAGCGGGTGGCATCTACAGATTCGTTTACGTCGCTATTTTTTTCGCTAACGGCGATGCTGCTGGCATCttt                    <  1:497355/129‑1 (MQ=255)                                                                                                       cGCCACAGCGCGTGGGTCACGGATTTATCGCCGCCCTGACAACCGCCTATAACCAGAAAACGGGTGGCATCTACAGATTCGTCAACGGCGCTAATTACTTCGCCAACGGCGATGCTGCTGGCATCTTTg                   >  1:604720/1‑129 (MQ=255)                                                                                                             aGCGCGTGGGTCCGGGCTTTATCGCCGCCCTGGCAGCCACCATTACACTGACAACGGGAGGCATCTACAGATTCGTCAACGGCGCTAATTACTTCGCCAACGGCGATGCTGCTGGCATCTTTGCCTagc             >  1:343219/1‑127 (MQ=255)                                                                                                             aGCGCGTGGGTCAGGGATTTATCGCCGCCCTGGCAGCCGCCTTTACCCTGACAACGGGTGGCATCTACAGATTCGTCAACGGCGCTAATTACTTCGCCAACGGCGATGCTGCTGGCATCTTTGCCTAAc             >  2:492340/1‑129 (MQ=255)                                                                                                               cgcgTGGGTCAGGGATTTATCGCCGCCCTGGCAGCCGCCTTTACCCTGACAACGGTTTGCATCTACATATTCGTCAACGGCGCTAATTTCTTCGCCAACGGCGATGCTGCTGGCATCTTTGCCTAACAg           <  1:459470/129‑1 (MQ=255)                                                                                                                        cAGGGATTTATCGCCGCCCTGGCAGCCGCCTTTACCCTGACAACGGGTGGCATCTACAGATTCGTCAACGGCGCTAATTACTTCGCCAACGGCGATGCTGCTGGCATCTTTGCCTAACAGATAaccacc  <  2:32335/129‑1 (MQ=255)                                                                                                                            |                                                                                                                             TGAGTATGCTGACGACCAGACACATCCAGCACTTCCTGGTTATTAACCAGTTCGCCTAAAGTAATGTTGTTGAGAAAACCGGTGAGACGGTCGCTCAAATCACGCCACAGCGCGTGGGTCAGGCATTTATCGCCGCCCTGGCAGCCGCCTTTACCCTGACAACGGGTGGCATCTACAGATTCGTCAACGGCGCTAATTACTTCGCCAACGGCGATGCTGCTGGCATCTTTGCCTAACAGATAACCACC  >  NC_000913/2661698‑2661945

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