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%	108.8 / NA	35	*123* (TAG→TAA)	sufA	Fe‑S cluster assembly protein
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (18/17);  total (18/17)

TCTTCGTTTATCCCCTTTGCCAGCTCATCGGTGGCTAACTGGGTGAAGAATCCTTCTTTATAATTCAGCGGGCCGCCGGTCCAGGTTTTGACATCGTCAGTTGCTTCAGTATTACGAGACATAGTACCGCCTATACCCCAAAGCTTTCGCCACAGCCACATTCATTCTGGGCTTTAGGGTTGTGAAATTTGAATATCTGATTAAGTCCTTCACGAACGAAATCGACTTCCGTGCCATCAATAAACGGCATCGCTTGCAGCGGG  >  NC_000913/1763888‑1764150                                                                                                                                   |                                                                                                                                     tctTCGTTTATCCCCTTTGCCAGCTCATCGGTGGCTAACTGGGTGAAGAATCCTTCTTTATAATTCAGCGGGCCGCCGGTCCAGGTTTTGACATCGTCA‑TTGCTTCAGTATTACGAGACATAGTACCGCTTATACCCCa                                                                                                                             <  2:76298/139‑1 (MQ=255)       tttATCCCCTTTGCCAGCTCATCGGTGGCTAACTGGGTGAAGAATCCTTCTTTATAATTCAGCGGGCCGCCGGTCCAGGTTTTGACATCGTCA‑TTGCTTCAGTATTACGAGACATAGTACCGCTTATACCCCAAAGCtt                                                                                                                       >  1:591348/1‑139 (MQ=255)          aTCCCCTTTGCCAGCTCATCGGTGGCTAACTGGGTGAAGAATCCTTCTTTATAATTCAGCGGGCCGCCGGTCCAGGTTTTGACATCGTCA‑TTGCTTCAGTATTACGAGACATAGTACCGCTTATACCCCAAAGCTTTCg                                                                                                                    >  2:239147/1‑139 (MQ=255)               cTTTGCCAGCTCATCGGTGGCTAACTGGGTGAAGAATCCTTCTTTATAATTCAGCGGGCCGCCGGTCCAGGTTTTGACATCGTCA‑TTGCTTCAGTATTACGAGACATAGTACCGCTTATACCCCAAAGCTTTCgccaca                                                                                                               <  1:239147/139‑1 (MQ=255)                tttGCCAGCTCATCGGTGGCTAACTGGGTGAAGAATCCTTCTTTATAATTCAGCGGGCCGCCGGTCCAGGTTTTGACATCGTCA‑TTGCTTCAGTATTACGAGACATAGTACCGCTTATACCCCAAAGCTTTCgcc                                                                                                                  >  2:20368/1‑135 (MQ=255)                tttGCCAGCTCATCGGTGGCTAACTGGGTGAAGAATCCTTCTTTATAATTCAGCGGGCCGCCGGTCCAGGTTTTGACATCGTCA‑TTGCTTCAGTATTACGAGACATAGTACCGCTTATACCCCAAAGCTTTCgcc                                                                                                                  <  1:20368/135‑1 (MQ=255)                     cAGCTCATCGGTGGCTAACTGGGTGAAGAATCCTTCTTTATAATTCAGCGGGCCGCCGGTCCAGGTTTTGACATCGTCA‑TTGCTTCAGTATTACGAGACATAGTACCGCTTATACCCCAAAGCTTTCgccacagccaca                                                                                                         >  1:508920/1‑139 (MQ=255)                        cTCATCGGTGGCTAACTGGGTGAAGAATCCTTCTTTATAATTCAGCGGGCCGCCGGTCCAGGTTTTGACATCGTCA‑TTGCTTCAGTATTACGAGACATAGTACCGCTTATACCCCAAAGCTTTCGCCACAGCCAcattc                                                                                                      >  2:250299/1‑139 (MQ=255)                                         gggTGAAGAATCCTTCTTTATAATTCAGCGGGCCGCCGGTCCAGGTTTTGACATCGTCA‑TTGCTTCAGTATTACGAGACATAGTACCGCTTATACCCCAAAGCTTTCGCCACAGCCACATTCATTCTGGGCTTTAGGGt                                                                                     <  2:591348/139‑1 (MQ=255)                                               agaaTCCTTCTTTATAATTCAGCGGGCCGCCGGTCCAGGTTTTGACATCGTCA‑TTGCTTCAGTATTACGAGACATAGTACCGCTTATACCCCAAAGCTTTCGCCACAGCCACATTCATTCTGGGCTTTAGGGTTGTGaa                                                                               <  2:478972/139‑1 (MQ=255)                                                        cttTATAATTCAGCGGGCCGCCGGTCCAGGTTTTGACATCGTCA‑TTGCTTCAGTATTACGAGACATAGTACCGCTTATACCCCAAAGCTTTCgccacagc                                                                                                             >  1:256291/1‑100 (MQ=255)                                                        cttTATAATTCAGCGGGCCGCCGGTCCAGGTTTTGACATCGTCA‑TTGCTTCAGTATTACGAGACATAGTACCGCTTATACCCCAAAGCTTTCgccacagc                                                                                                             <  2:256291/100‑1 (MQ=255)                                                               aTTCAGCGGGCCGCCGGTCCAGGTTTTGACATCGTCA‑TTGCTTCAGTATTACGAGACATAGTACCGCTTATACCCCAAAGCTTTCGCCACAGCCACattcattc                                                                                                  >  2:707142/1‑104 (MQ=255)                                                               aTTCAGCGGGCCGCCGGTCCAGGTTTTGACATCGTCA‑TTGCTTCAGTATTACGAGACATAGTACCGCTTATACCCCAAAGCTTTCGCCACAGCCACattcattc                                                                                                  <  1:707142/104‑1 (MQ=255)                                                               aTTCAGCGGGCCGCCGGTCCAGGTTTTGACATCGTCA‑TTGCTTCAGTATTACGAGACATAGTACCGCTTATACCCCAAAGCTTTCGCCACAGCCACATTCATTCTGGGCTTTAGGGTTGTGAAATTTGAATATCTGAtt                                                               <  1:250299/139‑1 (MQ=255)                                                                     cGGGCCGCCGGTCCAGGTTTTGACATCGTCA‑TTGCTTCAGTATTACGAGACATAGTACCGCTTATACCCCAAAGCTTTCgccacagc                                                                                                             <  2:894982/87‑1 (MQ=255)                                                                     cGGGCCGCCGGTCCAGGTTTTGACATCGTCA‑TTGCTTCAGTATTACGAGACATAGTACCGCTTATACCCCAAAGCTTTCgccacagc                                                                                                             >  1:894982/1‑87 (MQ=255)                                                                     cGGGCCGCCGGTCCAGGTTTTGACATCGTCA‑TTGCTTCAGTATTACGAGACATAGTACCGCTTATACCCCAAAGCTTTCGCCACAGCCACATTCATTCTg                                                                                                <  2:9965/100‑1 (MQ=255)                                                                     cGGGCCGCCGGTCCAGGTTTTGACATCGTCA‑TTGCTTCAGTATTACGAGACATAGTACCGCTTATACCCCAAAGCTTTCGCCACAGCCACATTCATTCTg                                                                                                >  1:9965/1‑100 (MQ=255)                                                                     cGGGCCGCCGGTCCAGGTTTTGACATCGTCA‑TTGCTTCAGTATTACGAGACATAGTACCGCTTATACCCCAAAGCTTTCGCCACAGCCACATTCATTCTGGGCTTTAGGGTTGTGAAATTTGAATATCTGATTAAGTcc                                                         >  1:270088/1‑139 (MQ=255)                                                                            ccgGTCCAGGTTTTGACATCGTCA‑TTGCTTCAGTATTACGAGACATAGTACCGCTTATACCCCAAAGCTTTCGCCACAGCCACATTCATTCTGGGCTTTAGGGTTGTGAAATTTGAATATCTGATTAAGTCCTTCacga                                                  <  2:270088/139‑1 (MQ=255)                                                                                                           cAGTATTACGAGACATAGTACCGCTTATACCCCAAAGCTTTCGCCACAGCCACATTCATTCTGGGCTTTAGGGTTGTGAAATTTGAATATCTGAtt                                                               >  2:877580/1‑96 (MQ=255)                                                                                                           cAGTATTACGAGACATAGTACCGCTTATACCCCAAAGCTTTCGCCACAGCCACATTCATTCTGGGCTTTAGGGTTGTGAAATTTGAATATCTGAtt                                                               <  1:877580/96‑1 (MQ=255)                                                                                                           cAGTATTACGAGACATAGTACCGCTTATACCCCAAAGCTTTCGCCACAGCCACATTCATTCTGGGCTTTAGGGTTGTGAAATTAGAATATCTGATTAAGTCCTTCACGAACGAAATCGACTTCCGTGCCATCAATAAAc                    >  1:590042/1‑139 (MQ=255)                                                                                                                  aCGAGACATAGTACCGCTTATACCCCAAAGCTTTCGCCACAGCCACATTCATTCTGGGCTTTAGGGTTGTGAAATTTGAATATCTGATTAAGTCCTTCACGAACGAAATCGACTTCCGTGCCATCAATAAACGGCATCg             >  2:232741/1‑139 (MQ=255)                                                                                                                   cGAGACATAGTACCGCTTATACCCCAAAGCTTTCGCCACAGCCACATTCATTCTGGGCTTTAGGGTTGTGAAATTTGAATATCTGATTAAGTCCTTCACGAACGAAATCGACTTCCGTGCCATCAATAAACGGCATCGc            >  2:550514/1‑139 (MQ=255)                                                                                                                      gaCATAGTACCGCTTATACCGCAAAGCTTTCGCCACAGCCACATTCATTCTGTGCTTTAGGGTTGTGAAATTTGAATATCTGATTAAGTCCTTCACGAACGAAATCGACTTCCGTGCCATCAATAAACGGCATCGCTTg         <  1:735568/139‑1 (MQ=255)                                                                                                                        cATAGTACCGCTTATACCCCAAAGCTTTCGCCACAGCCACATTCATTCTGGGCTTTAGGGTTGTGAAATTTGAATATCTGATTAAGTCCTTTACGAACGAAATCGACTTCCGTGCCATCAATAAAGGGCATCGCTTACa       >  2:735568/1‑139 (MQ=255)                                                                                                                        cATAGTACCGCTTATACCCCAAAGCTTTCGCCACAGCCACATTCATTCTGGGCTTTAGGGTTGTGAAATTTGAATATCTGATTAAGTCCTTCACGAACGAAATCGACTTCCGt                                 <  2:213366/113‑1 (MQ=255)                                                                                                                        cATAGTACCGCTTATACCCCAAAGCTTTCGCCACAGCCACATTCATTCTGGGCTTTAGGGTTGTGAAATTTGAATATCTGATTAAGTCCTTCACGAACGAAATCGACTTCCGt                                 >  1:213366/1‑113 (MQ=255)                                                                                                                        cATAGTACCGCTTATACCCCAAAGCTTTCGCCACAGCCACATTCATTCTGGGCTTTAGGGTTGTGAAATTTGAATATCTGATTAAGTCCTTCACGAACGAAATCGACTTCCGTGCCATCAATAAACGGCATCGCTTGCa       >  1:993764/1‑139 (MQ=255)                                                                                                                        cATAGTACAGCTTATACCCCAAAGCTTTCGCCACAGCCACATTCATTCTGGGCTTTAGGGTTGTGAAATTTGAATATCTGATTAAGTCCTTCACGAACGAAATCGACTTCCGTGCCATCAATAAACGTCATCGCTTGCa       >  1:722751/1‑139 (MQ=255)                                                                                                                           aGTACCGCTTATACCCCAAAGCTTTCGCCACAGCCACATTCATTCTGGGCTTTAGGGTTGTGAAATTTGAATATCTGATTAAGTCCTTCACGAACGAAATCGACTTCCGTGCCATCAATAAACGGCATCGCTTGCAGCg    <  2:508920/139‑1 (MQ=255)                                                                                                                           aGTACCGCTTATACCCCAAAGCTTTCGCCACAGCCACATTCATTCTGGGCTTTAGGGTTGTGAAATTTGAATATCTGATTAAGTCCTTCACGAACGAAATCGACTTCCGTGCCATCAATAAACGGCATCGCTTGCAGCg    <  1:232741/139‑1 (MQ=255)                                                                                                                             tACCGCTTATACCCCAAAGCTTTCGCCACAGCCACATTCATTCTGGGCTTTAGGGTTGTGAAATTTGAATATCTGATTAAGTCCTTCACGAACGAAATCGACTTCCGTGCCATCAATAAACGGCATCGCTTGCAGCggg  <  2:590042/139‑1 (MQ=255)                                                                                                                                   |                                                                                                                                     TCTTCGTTTATCCCCTTTGCCAGCTCATCGGTGGCTAACTGGGTGAAGAATCCTTCTTTATAATTCAGCGGGCCGCCGGTCCAGGTTTTGACATCGTCAGTTGCTTCAGTATTACGAGACATAGTACCGCCTATACCCCAAAGCTTTCGCCACAGCCACATTCATTCTGGGCTTTAGGGTTGTGAAATTTGAATATCTGATTAAGTCCTTCACGAACGAAATCGACTTCCGTGCCATCAATAAACGGCATCGCTTGCAGCGGG  >  NC_000913/1763888‑1764150

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