Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,764,018	C→T	*123* (TAG→TAA)	sufA ←	Fe‑S cluster assembly protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,764,018	0	C	T	100.0%	88.3 / NA	28	*123* (TAG→TAA)	sufA	Fe‑S cluster assembly protein
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (18/10);  total (18/10)

CTTTGCCAGCTCATCGGTGGCTAACTGGGTGAAGAATCCTTCTTTATAATTCAGCGGGCCGCCGGTCCAGGTTTTGACATCGTCAGTTGCTTCAGTATTACGAGACATAGTACCGCCTATACCCCAAAGCTTTCGCCACAGCCACATTCATTCTGGGCTTTAGGGTTGTGAAATTTGAATATCTGATTAAGTCCTTCACGAACGAAATCGACTTCCGTGCCATCAATAAACGGCATCGCTTGCA  >  NC_000913/1763902‑1764145                                                                                                                     |                                                                                                                                cTTTGCCAGCTCATCGGTGGCTAACTGGGTGAAGAATCCTTCTTTATAATTCAGCGGGCCGCCGGTCCAGGTTTTGACATCGTCA‑TTGCTTCAGTATTACGAGACATAGTACCGCTTATACCCCAAAGCTTTCgccaca                                                                                                          >  2:493011/1‑139 (MQ=255)     gCCAGCTCATCGGTGGCTAACTGGGTGAAGAATCCTTCTTTATAATTCAGCGGGCCGCCGGTCCAGGTTTTGACATCGTCA‑TTGCTTCAGTATTACGAGACATAGTACCGCTTATACCCCAAAGCTTTCgccacagcca                                                                                                      <  1:92803/139‑1 (MQ=255)               cGGTGGCTAACTGGGTGAAGAATCCTTCTTTATAATTCAGCGGGCCGCCGGTCCAGGTTTTGACATCGTCA‑TTGCTTCAGTATTACGAGACATAGTACCGCTTATACCCCAAAGCTTTCGCCACAGCCACATTCATTCt                                                                                            >  1:703582/1‑139 (MQ=255)               cGGTGGCTAACTGGGTGAAGAATCCTTCTTTATAATTCAGCGGGCCGCCGGTCCAGGTTTTGACATCCTCA‑TTGCTTCAGTATTACGAGACATAGTACCGCTTATACCCCAAAGCTTTCGCCACAGCCACATTCATTCt                                                                                            <  1:458191/139‑1 (MQ=255)                            ggTGAAGAATCCTTCTTTATAATTCAGCGGGCCGCCAGTCCAGGTTTCGACATCGTCA‑TTGCTGCAGTATTACGAGACATAGTACCGATTATACCCCAAAGCTTTCGCCACAGCCACATTCATTCTGGGCTTTAGGGtt                                                                               >  2:509500/1‑139 (MQ=255)                                            ttATAATTCAGCGGGCCGCCGGTCCAGGTTTTGACATCGTCA‑TTGCTTCAGTATTACGAGACATAGTACCGCTTATACCCCAAAGCTTTCGCCACAGCCACATTCATTCTGGGCTTTAGGGTTGTGAAAttt                                                                      <  2:742368/132‑1 (MQ=255)                                            ttATAATTCAGCGGGCCGCCGGTCCAGGTTTTGACATCGTCA‑TTGCTTCAGTATTACGAGACATAGTACCGCTTATACCCCAAAGCTTTCGCCACAGCCACATTCATTCTGGGCTTTAGGGTTGTGAAAttt                                                                      >  1:742368/1‑132 (MQ=255)                                            ttATAATTCAGCGGGCCGCCGGTCCAGGTTTTGACATCGTCA‑TTGCTTCAGTATTACGAGACATAGTACCGCTTATACCCCAAAGCTTTCGCCACAGCCACATTCATTCTGGGCTTTAGGGTTGTGAAATTTGAATATc                                                               <  2:698838/139‑1 (MQ=255)                                                  ttCAGCGGGCCGCCGGTCCAGGTTTTGACATCGTCA‑TTGCTTCAGTATTACGAGACATAGTACCGCTTATACCCCAAAGCTTTCGCCACAGCCACATTCATTCTGGGCTTTAGGGTTGTGAAATTTGAATATCTGATTa                                                         >  2:22454/1‑139 (MQ=255)                                                       cGGGCCGCCGGTCCAGGTTTTGACATCGTCA‑TTGCTTCAGTATTACGAGACATAGTACCGCTTATACCCCAAAGCTTTCGCCACAGCCACATTCATTCTGGGCTTTa                                                                                    <  2:597893/107‑1 (MQ=255)                                                       cGGGCCGCCGGTCCAGGTTTTGACATCGTCA‑TTGCTTCAGTATTACGAGACATAGTACCGCTTATACCCCAAAGCTTTCGCCACAGCCACATTCATTCTGGGCTTTa                                                                                    >  1:597893/1‑107 (MQ=255)                                                       cGGGCCGCCGGTCCAGGTTTTGACATCGTCA‑TTGCTTCAGTATTACGAGACATAGTACCGCTTATACCCCAAAGCTTTCGCCACAGCCACATTCATTCTGGGCTTTAGGGTTGTGAAATTTGAATATCTGATTAAGTcc                                                    >  1:150420/1‑139 (MQ=255)                                                       cGGGCCGCCGGTCCAGGTTTTGACATCGTCA‑TTGCTTCAGTATTACGAGACATAGTACCGCTTATACCCCAAAGCTTTCGCCACAGCCACATTCATTCTGGGCTTTAGGGTTGTGAAATTTGAATATCTGATTAAGTcc                                                    >  1:415324/1‑139 (MQ=255)                                                       cGGGCCGCCGGTCCAGGTTTTGACATCGTCA‑TTGCTTCAGTATTACGAGACATAGTACCGCTTATACCCCAAAGCTTTCGCCACAGCCACATTCATTCTGGGCTTTAGGGTTGTGAAATTTGAATATCTGATTAAGTcc                                                    >  1:222350/1‑139 (MQ=255)                                                            cgccgGTCCAGGTTTTGACATCGTCA‑TTGCTTCAGTATTACGAGACATAGTACCGCTTATACCCCAAAGCTTTCGCCACAGCCACATTCATTCTGGGCTTTAGGGTTGTGAAATTTGAATATCTGATTAAGTCCTTCac                                               >  2:461477/1‑139 (MQ=255)                                                                               aTCGTCA‑TTGCTTCAGTATTACGAGACATAGTACCGCTTATACCCCAAAGCTTTCGCCACAGCCACATTCATTCTGGGCTTTAGGGTTGTGAAATTTGAATATCTGATTAAGTCCTTCACGAACGAAATCGACTTCCGt                            <  1:493011/139‑1 (MQ=255)                                                                                     cattGCTTCAGTATTACGAGACATAGTACCGCTTATAGCCCAAAGCTTTCGCCACAGCCACATTCATTCTGGGCTTTAGGGTTGTGAAATTTGAATATCTGATTAAGTCCTTCACGAACGAAATCGACTTCCGTGCCAt                       <  1:509500/137‑1 (MQ=255)                                                                                            tCAGTATTACGAGACATAGTACCGCTTATACCCCAAAGCTTTCGCCACAGCCACATTCATTCTGGGCTTTAGGGTTGTGAAATTTGAATATCTGATTAAGTCCTTCACGAACGAAATCGACTTCCGTGCCATCAATaaa                >  1:711288/1‑139 (MQ=255)                                                                                              aGTATTACGAGACATAGTACCGCTTATACCCCAAAGCTTTCGCCACAGCCACATTCATTCTGGGCTTTAGGGTTGTGAAATTTGAATATCTGATTAAGTCCTTCACGAACGAAATCGACTTCCGTGCCATCAATAAACg              >  2:68412/1‑139 (MQ=255)                                                                                              aGTATTACGAGACATAGTACCGCTTATACCCCAAAGCTTTCGCCACAGCCACATTCATTCTGGGCTTTAGGGTTGTGAAATTTGAATATCTGATTAAGTCCTTCACGAACGAAATCGACTTCCGTGCCATCAATAAACg              >  2:606234/1‑139 (MQ=255)                                                                                                tATTACGAGACATAGTACCGCTTATACCCCAAAGCTTTCGCCACAGCCACATTCATTCTGGGCTTTAGGGTTGTGAAATTTGAATATCTGATTAAGTCCTTCACGAACGAAATCGACTTCCGTGCCATCAATAAACGGc            >  2:645877/1‑139 (MQ=255)                                                                                                tATTACGAGACATAGTACCGCTTATACCCCAAAGCTTTCGCCACAGCCACATTCATTCTGGGCTTTAGGGTTGTGAAATTTGAATATCTGATTAAGTCCTTCACGAACGAAATCGACTTCCGTGCCATCAATAAACGGc            <  2:150420/139‑1 (MQ=255)                                                                                                tATTACGAGACATAGTACCGCTTATACCCCAAAGCTTTCGCCACAGCCACATTCATTCTGGGCTTTAGGGTTGTGAAATTTGAATATCTGATTAAGTCCTTCACGAACGAAATCGACTTCCGTGCCATCAATAAACGGc            <  2:222350/139‑1 (MQ=255)                                                                                                  ttACGAGACATAGTACCGCTTATACCCCAAAGCTTTCGCCACAGCCACATTCATTCTggg                                                                                         <  1:136743/60‑1 (MQ=255)                                                                                                  ttACGAGACATAGTACCGCTTATACCCCAAAGCTTTCGCCACAGCCACATTCATTCTggg                                                                                         >  2:136743/1‑60 (MQ=255)                                                                                                     cGAGACATAGTACCGCTTATACCCCAAAGCTTTCGCCACAGCCACATTCATTCTGGGCTTTAGGGTTGTGAAATTTGAATATCTGATTAAGTCCTTCACGAACGAAATCGACTTCCGTGCCATCAATAAACGGCATCGc       >  2:124573/1‑139 (MQ=255)                                                                                                     cGAGACATAGTACCGCTTATACCCCAAAGCTTTCGCCACAGCCACATTCATTCTGGGCTTTAGGGTTGTGAAATTTGAATATCTGATTAAGTCCTTCACGAACGAAATCGACTTCCGTGCCATCAATAAACGGCATCGc       >  2:755427/1‑139 (MQ=255)                                                                                                          cATAGTACCGCTTATACCCCAAAGCTTTCGCCACAGCCACATTCATTCTGGGCTTTAGGGTTGTGAAATTTGAATATCTGATTAAGTCCTTCACGAACGAAATCGACTTCCGTGCCATCAATAAACGGCATCGCTTGCa  >  2:744418/1‑139 (MQ=255)                                                                                                                     |                                                                                                                                CTTTGCCAGCTCATCGGTGGCTAACTGGGTGAAGAATCCTTCTTTATAATTCAGCGGGCCGCCGGTCCAGGTTTTGACATCGTCAGTTGCTTCAGTATTACGAGACATAGTACCGCCTATACCCCAAAGCTTTCGCCACAGCCACATTCATTCTGGGCTTTAGGGTTGTGAAATTTGAATATCTGATTAAGTCCTTCACGAACGAAATCGACTTCCGTGCCATCAATAAACGGCATCGCTTGCA  >  NC_000913/1763902‑1764145

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