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

TCTTCGTTTATCCCCTTTGCCAGCTCATCGGTGGCTAACTGGGTGAAGAATCCTTCTTTATAATTCAGCGGGCCGCCGGTCCAGGTTTTGACATCGTCAGTTGCTTCAGTATTACGAGACATAGTACCGCCTATACCCCAAAGCTTTCGCCACAGCCACATTCATTCTGGGCTTTAGGGTTGTGAAATTTGAATATCTGATTAAGTCCTTCACGAACGAAATCGACTTCCGTGCCATCAATAAACGGCATCGCTTGCAGCGGGACA  >  NC_000913/1763888‑1764153                                                                                                                                   |                                                                                                                                        tctTCGTTTATCCCCTTTGCCAGCTCATCGGTGGCTAACTGGGTGAAGAATCCTTCTTTATAATTCAGCGGGCCGCCGGTCCAGGTTTTGACATCGTCA‑TTGCTTCAGTATTACGAGACATAGTACCGCTTATACCCCa                                                                                                                                <  2:1674712/139‑1 (MQ=255)     cGTTTATCCCCTTTGCCAGCTCATCGGTGGCTAACTGGGTGAAGAATCCTTCTTTATAATTCAGCGGGCCGCCGGTCCAGGTTTTGACATCGTCA‑TTGCTTCAGTATTACGAGACATAGTACCGCTTATACCCCAAAGc                                                                                                                            >  1:590268/1‑139 (MQ=255)            ccccTTTGCCAGCTCATCGGTGGCTAACTGGGTGAAGAATCCTTCTTTATAATTCAGCGGGCCGCCGGTCCAGGTTTTGACATCGTCA‑TTGCTTCAGTATTACGAGACATAGTACCGCTTATACCCCAAAGCTTTCgcc                                                                                                                     <  2:72962/139‑1 (MQ=255)            ccccTTTGCCAGCTCATCGGTGGCTAACTGGGTGAAGAATCCTTCTTTATAATTCAGCGGGCCGCCGGTCCAGGTTTTGACATCGTCA‑TTGCTTCAGTATTACGAGACATAGTACCGCTTATACCCCAAAGCTTTCgcc                                                                                                                     >  1:1363475/1‑139 (MQ=255)                          cATCGGTGGCTAACTGGGTGAAGAATCCTTCTTTATAATTCAGCGGGCCGCCGGTCCAGGTTTTGACATCGTCA‑TTGCTTCAGTATTACGAGACATAGTACCGCTTATACCCCAAAGCTTTCGCCACAGCCAcattcat                                                                                                       >  2:930337/1‑139 (MQ=255)                          cATCGGTGGCTAACTGGGTGAAGAATCCTTCTTTATAATTCAGCGGGCCGCCGGTCCAGGTTTTGACATCGTCA‑TTGCTTCAGTATTACGAGACATAGTACCGCTTATACCCCAAAGCTTTCGCCACAGCCAcattcat                                                                                                       >  2:850277/1‑139 (MQ=255)                           aTCGGTGGCTAACTGGGTGAAGAATCCTTCTTTATAATTCAGCGGGCCGCCGGTCCAGGTTTTGACATCGTCA‑TTGCTTCAGTATTACGAGACATAGTACCGCTTATACCCCAAAGCTTTCGCCACAGCCAcattcatt                                                                                                      <  1:850277/139‑1 (MQ=255)                           aTCGGTGGCTAACTGGGTGAAGAATCCTTCTTTATAATTCAGCGGGCCGCCGGTCCAGGTTTTGACATCGTCA‑TTGCTTCAGTATTACGAGACATAGTACCGCTTATACCCCAAAGCTTTCGCCACAGCCAcattcatt                                                                                                      <  1:930337/139‑1 (MQ=255)                            tCGGTGGCTAACTGGGTGAAGAATCCTTCTTTATAATTCAGCGGGCCGCCGGTCCAGGTTTTGACATCGTCA‑TTGCTTCAGTATTACGAGACATAGTACCGCTTATACCCCAAAGCTTTCGCCACAGCCACattcattc                                                                                                     <  2:1363475/139‑1 (MQ=255)                            tCGGTGGCTAACTGGGTGAAGAATCCTTCTTTATAATTCAGCGGGCCGCCGGTCCAGGTTTTGACATCGTCA‑TTGCTTCAGTATTACGAGACATAGTACCGCTTATACCCCAAAGCTTTCGCCACAGCCACattcattc                                                                                                     >  2:1365146/1‑139 (MQ=255)                                            tGAAGAATCCTTCTTTATAATTCAGCGGGCCGCCGGTCCAGGTTTTGACATCGTCA‑TTGCTTCAGTATTACGAGACATAGTACCGCTTATACCCCAAAGCTTTCGCCACAGCCACATTCATTCTg                                                                                                   <  2:1259847/125‑1 (MQ=255)                                            tGAAGAATCCTTCTTTATAATTCAGCGGGCCGCCGGTCCAGGTTTTGACATCGTCA‑TTGCTTCAGTATTACGAGACATAGTACCGCTTATACCCCAAAGCTTTCGCCACAGCCACATTCATTCTg                                                                                                   >  1:1259847/1‑125 (MQ=255)                                                     cttcttTATAATTCAGCGGGCCGCCGGTCCAGGTTTTGACATCGTCA‑TTGCTTCAGTATTACGAGACATAGTACCGCTTATACCCCAAAGCTTTCGCCACAGCCACATTCATTCTGGGCTTTAGGGTTGTGAAATTTGa                                                                            <  2:590268/139‑1 (MQ=255)                                                          ttATAATTCAGCGGGCCGCCGGTCCAGGTTTTGACATCGTCA‑TTGCTTCAGTATTACGAGACATAGTACCGCTTATACCCCAAAGCTTTCGCCACAGCCACATTCATTCTGGGCTTTAGGGTTGTGAAATTTGAATATc                                                                       <  1:199597/139‑1 (MQ=255)                                                                  cAGCGGGCCGCCGGTCCAGGTTTTGACATCGTCA‑TTGCTTCAGTATTACGAGACATAGTACCGCTTATACCCCAAAGCTTTCGCCACAGCCACATTCATTCTGGGCTTTAGGGTTGTGaa                                                                                  <  1:842115/120‑1 (MQ=255)                                                                  cAGCGGGCCGCCGGTCCAGGTTTTGACATCGTCA‑TTGCTTCAGTATTACGAGACATAGTACCGCTTATACCCCAAAGCTTTCGCCACAGCCACATTCATTCTGGGCTTTAGGGTTGTGaa                                                                                  >  2:842115/1‑120 (MQ=255)                                                                   aGCGGGCCGCCGGTCCAGGTTTTGACATCGTCA‑TTGCTTCAGTATTACGAGACATAGTACCGCTTATACCCCAAAGCTTTCgcca                                                                                                                    <  1:1190112/85‑1 (MQ=255)                                                                   aGCGGGCCGCCGGTCCAGGTTTTGACATCGTCA‑TTGCTTCAGTATTACGAGACATAGTACCGCTTATACCCCAAAGCTTTCgcca                                                                                                                    >  2:1190112/1‑85 (MQ=255)                                                                   aGCGGGCCGCCGGTCCAGGTTTTGACATCGTCA‑TTGCTTCAGTATTACGAGACATAGTACCGCTTATACCCCAAAGCTTTCGCCACAGCCACATTCATTCTGGGCTTTAGGGTTGTGAAATTTg                                                                             <  1:1323959/124‑1 (MQ=255)                                                                   aGCGGGCCGCCGGTCCAGGTTTTGACATCGTCA‑TTGCTTCAGTATTACGAGACATAGTACCGCTTATACCCCAAAGCTTTCGCCACAGCCACATTCATTCTGGGCTTTAGGGTTGTGAAATTTg                                                                             >  2:1323959/1‑124 (MQ=255)                                                                     cGGGCCGCCGGTCCAGGTTTTGACATCGTCA‑TTGCTTCAGTATTACGAGACATAGTACCGCTTATACCGCAAAGCTTTCgccacagccaca                                                                                                            <  2:1280738/91‑1 (MQ=255)                                                                     cGGGCCGCCGGTCCAGGTTTTGACATCGTCA‑TTGCTTCAGTATTACGAGACATAGTACCGCTTATACCGCAAAGCTTTCgccacagccaca                                                                                                            >  1:1280738/1‑91 (MQ=255)                                                                     cGGGCCGCCGGTCCAGGTTTTGACATCGTCA‑TTGCTTCAGTATTACGAGACATAGTACCGCTTATACCCCAAAGCttt                                                                                                                         >  2:1349229/1‑78 (MQ=255)                                                                     cGGGCCGCCGGTCCAGGTTTTGACATCGTCA‑TTGCTTCAGTATTACGAGACATAGTACCGCTTATACCCCAAAGCttt                                                                                                                         <  1:1349229/78‑1 (MQ=255)                                                                     cGGGCCGCCGGTCCAGGTTTTGACATCGTCA‑TTGCTTCAGTATTACGAGACATAGTACCGCTTATACCCCAAAGCTTTCGCCACAGCCACATTCATTCTGGGCttt                                                                                             <  2:1007905/106‑1 (MQ=255)                                                                     cGGGCCGCCGGTCCAGGTTTTGACATCGTCA‑TTGCTTCAGTATTACGAGACATAGTACCGCTTATACCCCAAAGCTTTCGCCACAGCCACATTCATTCTGGGCttt                                                                                             >  1:1007905/1‑106 (MQ=255)                                                                     cGGGCCGCCGGTCCAGGTTTTGACATCGTCA‑TTGCTTCAGTATTACGAGACATAGTACCGCTTATACCCCAAAGCTTTCGCCACAGCAACATTCATTCTGGGCTTTAGGGTTGTGAAATTTGa                                                                            >  2:1129363/1‑123 (MQ=255)                                                                     cGGGCCGCCGGTCCAGGTTTTGACATCGTCA‑TTGCTTCAGTATTACGACACATAGTACCGCTTATACCCCAAAGCTTTCGCCACAGCCACATTCATTCTGGGCTTTAGGGTTGTGAAATTTGa                                                                            <  1:1129363/123‑1 (MQ=255)                                                                     cGGGCCGCCGGTCCAGGTTTTGACATCGTCA‑TTGCATCAGTATTACGAGACATAGTACCGCTTATACCCCAAAGCTTTCGCCACAGCCACATTCATACTGGGATTTAGGGTAGAGAAATTTGAATATCTCATTAAGTcc                                                            >  2:469197/1‑139 (MQ=255)                                                                          cgccgGTCCAGGTTTTGACATCGTCA‑TTGCTTCAGTATTACGAGACATAGTACCGCTTATACCCCAAAGCTTTCGCCACAGCCACATTCATTCTGGGCTTTAGGGTTGTGAAATTTGAATATc                                                                       >  1:661137/1‑123 (MQ=255)                                                                          cgccgGTCCAGGTTTTGACATCGTCA‑TTGCTTCAGTATTACGAGACATAGTACCGCTTATACCCCAAAGCTTTCGCCACAGCCACATTCATTCTGGGCTTTAGGGTTGTGAAATTTGAATATc                                                                       <  2:661137/123‑1 (MQ=255)                                                                           gccgGTCCAGGTTTTGACATCGTCA‑TTGCTTCAGTATTACGAGACATAGTACCGCTTATACCCCAAAGCTTTCgcca                                                                                                                    >  1:366535/1‑77 (MQ=255)                                                                           gccgGTCCAGGTTTTGACATCGTCA‑TTGCTTCAGTATTACGAGACATAGTACCGCTTATACCCCAAAGCTTTCgcca                                                                                                                    <  2:366535/77‑1 (MQ=255)                                                                               gTCCAGGTTTTGACATCGTCA‑TTGCTTCAGTATTACGAGACATAGTACCGCTTATACCCCAAAGCTTTCGCCACAGCCACATTCATTCTGGGCTTTAGGGTTGTGAAATTTGAATATCTGATTAAGTCCTTCacgaacg                                                  <  1:1365146/139‑1 (MQ=255)                                                                                    ggTTTTGACATCGTCA‑TTGCTTCAGTATTACGAGACATAGTACCGCTTATACCCCAAAGCTTTCGCCACAGCCAcattcatt                                                                                                      >  2:737367/1‑82 (MQ=255)                                                                                    ggTTTTGACATCGTCA‑TTGCTTCAGTATTACGAGACATAGTACCGCTTATACCCCAAAGCTTTCGCCACAGCCAcattcatt                                                                                                      <  1:737367/82‑1 (MQ=255)                                                                                     gTTTTGACATCGTCA‑TTGCTTCAGTATTACGAGACATAGTACCGCTTATACCCCAAAGCTTTCGCCACAGCCACATTCATTCTGGGCTTTAGGGTtgt                                                                                     <  1:476203/98‑1 (MQ=255)                                                                                     gTTTTGACATCGTCA‑TTGCTTCAGTATTACGAGACATAGTACCGCTTATACCCCAAAGCTTTCGCCACAGCCACATTCATTCTGGGCTTTAGGGTtgt                                                                                     >  2:476203/1‑98 (MQ=255)                                                                                          gACATCGTCA‑TTGCTTCAGTATTACGAGACATAGTACCGCTTATACCCCAAAGCTTTCGCCACAGCCACATTCATTCTg                                                                                                   <  1:405343/79‑1 (MQ=255)                                                                                          gACATCGTCA‑TTGCTTCAGTATTACGAGACATAGTACCGCTTATACCCCAAAGCTTTCGCCACAGCCACATTCATTCTg                                                                                                   >  2:405343/1‑79 (MQ=255)                                                                                          gACATCGTCA‑TTGCTTCAGTATTACGAGACATAGTACCGCTTATACCCCAAAGCTTTCGCCACAGCCACATTCATTCTGGGCTTTa                                                                                            >  2:149988/1‑86 (MQ=255)                                                                                          gACATCGTCA‑TTGCTTCAGTATTACGAGACATAGTACCGCTTATACCCCAAAGCTTTCGCCACAGCCACATTCATTCTGGGCTTTa                                                                                            <  1:149988/86‑1 (MQ=255)                                                                                                cgtcattGCTTCAGTATTACGAGACATAGTACCGCTTATACCCCAAAGCTTTCgcca                                                                                                                    >  2:372056/6‑57 (MQ=255)                                                                                                cgtcattGCTTCAGTATTACGAGACATAGTACCGCTTATACCACAAAGCTTTCgcca                                                                                                                    <  1:372056/52‑1 (MQ=38)                                                                                                            aGTATTACGAGACATAGTACCGCTTATACCCCAAAGCTTTCGCCACAGCCACATTCATTCTGGGCTTTAGGGTTGTGAAATTTGAATATCTGATTAAGTCCTTCACGAACGAAATCGACTTCCGTGCCATCAATAAACg                      >  1:646225/1‑139 (MQ=255)                                                                                                                   cGAGACATAGTACCGCTTATACCCCAAAGCTTTCGCCACAGCCACATTCATTCTGGGCTTTAGGGTTGTGAAATTTGAATATCTGATTAAGTCCTTCACGAACGAAATCGACTTCCGTGCCATCAATAAACGGCATCGc               >  1:802581/1‑139 (MQ=255)                                                                                                                           aGTACCGCTTATACCCCAAAGCTTTCGCCACAGCCACATTCATTCTGGGCTTTa                                                                                            >  1:1194915/1‑54 (MQ=255)                                                                                                                           aGTACCGCTTATACCCCAAAGCTTTCGCCACAGCCACATTCATTCTGGGCTTTa                                                                                            <  2:1194915/54‑1 (MQ=255)                                                                                                                           aGTACCGCTTATACCCCAAAGCTTTCGCCACAGCCACATTCATTCTGGGCTTTAGGGTTGTGAAAt                                                                                >  2:1253786/1‑66 (MQ=255)                                                                                                                           aGTACCGCTTATACCCCAAAGCTTTCGCCACAGCCACATTCATTCTGGGCTTTAGGGTTGTGAAAt                                                                                <  1:1253786/66‑1 (MQ=255)                                                                                                                             tACCGCTTATACCCCAAAGCTTTCGCCACAGCCACATTCATTCTGGGCTTTAGGGTTGTGAAATTTGAATATCTGATTAAGTCCTTCACGAACGAAATCGACTTCCGTGCCATCAATAAACGGCATCGCTTGCAGCggg     <  2:802581/139‑1 (MQ=255)                                                                                                                                cGCTTATACCCCAAAGCTTTCGCCACAGCCACATTCATTCTGGGCTTTAGGGTTGTGAAATTTGAATATCTGATTAAGTCCTTCACGAACGAAATCGACTTCCGTGCCATCAATAAACGGCATCGCTTGCAGCGGGACa  <  2:646225/139‑1 (MQ=255)                                                                                                                                   |                                                                                                                                        TCTTCGTTTATCCCCTTTGCCAGCTCATCGGTGGCTAACTGGGTGAAGAATCCTTCTTTATAATTCAGCGGGCCGCCGGTCCAGGTTTTGACATCGTCAGTTGCTTCAGTATTACGAGACATAGTACCGCCTATACCCCAAAGCTTTCGCCACAGCCACATTCATTCTGGGCTTTAGGGTTGTGAAATTTGAATATCTGATTAAGTCCTTCACGAACGAAATCGACTTCCGTGCCATCAATAAACGGCATCGCTTGCAGCGGGACA  >  NC_000913/1763888‑1764153

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