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

TTCGTTTATCCCCTTTGCCAGCTCATCGGTGGCTAACTGGGTGAAGAATCCTTCTTTATAATTCAGCGGGCCGCCGGTCCAGGTTTTGACATCGTCAGTTGCTTCAGTATTACGAGACATAGTACCGCCTATACCCCAAAGCTTTCGCCACAGCCACATTCATTCTGGGCTTTAGGGTTGTGAAATTTGAATATCTGATTAAGTCCTTCACGAACGAAATCGACTTCCGTGCCATCAATAAACGGCATCGCTTGCA  >  NC_000913/1763890‑1764145                                                                                                                                 |                                                                                                                                ttCGTTTATCCCCTTTGCCAGCTCATCGGTGGCTAACTGGGTGAAGAATCCTTCTTTATAATTCAGCGGGCCGCCGGTCCAGGTTTTGACATCGTCA‑TTGCTTCAGTATTACGAGACATAGTACCGCTTATACCCCaaa                                                                                                                      >  1:487330/1‑139 (MQ=255)   cGTTTATCCCCTTTGCCAGCTCATCGGTGGCTAACTGGGTGAAGAATCCTTCTTTATAATTCAGCGGGCCGCCGGTCCAGGTTTTGACATCGTCA‑TTGCTTCAGTATTACGAGACATAGTACCGCTTATACCCCAAAGc                                                                                                                    <  2:254948/139‑1 (MQ=255)      ttATTCCCTTTGCCAGGTCATCGGTGGCTAACTTGGTGAAGAATCCTTCTTTATAATTTCGCGGGCCGGCCGTCCAAGTTTTGACATCGTCA‑TTTCTTCAGTATTACGAAACATAGTACCGCTTTAACCCCAAAGCttt                                                                                                                 <  2:96546/139‑1 (MQ=255)       tATCCCCTTTGCCAGCTCATCGGTGGCTAACTGGGTGAAGAATCCTTCTTTATAATTCAGCGGGCCGCCGGTCCAGGTTTTGACATCGTCA‑TTGCTTCAGTATTACGAGACATAGTACCGCTTATACCCCAAAGCTTTc                                                                                                                >  1:96546/1‑139 (MQ=255)          ccccTTTGCCAGCTCATCGGTGGCTAACTGGGTGAAGAATCCTTCTTTATAATTCAGCGGGCCGCCGGTCCAGGTTTTGACATCGTCA‑TTGCTTCAGTATTACGAGACATAGTACCGCTTATACCCCAAAGCTTTCgcc                                                                                                             <  2:584533/139‑1 (MQ=255)               ttGCCAGCTCATCGGTGGCTAACTGGGTGAAGAATCCTTCTTTATAATTCAGCGGGCCGCCGGTCCAGGTTTTGACATCGTCA‑TTGCTTCAGTATTACGAGACATAGTACCGCTTATACCCCAAAGCTTTCgccacagc                                                                                                        >  2:269065/1‑139 (MQ=255)               ttGCCAGCTCATCGGTGGCTAACTGGGTGAAGAATCCTTCTTTATAATTCAGCGGGCCGCCGGTCCAGGTTTTGACATCGTCA‑TTGCTTCAGTATTACGAGACATAGTACCGCTTATACCCCAAAGCTTTCgccacagc                                                                                                        <  1:269065/139‑1 (MQ=255)                           cGGTGGCTAACTGGGTGAAGAATCCTTCTTTATAATTCAGCGGGCCGCCGGTCCAGGTTTTGACATCGTCA‑TTGCTTCAGTATTACGAGACATAGTACCGCTTATACCCCAAAGCTTTCGCCACAGCCACATTCATTCt                                                                                            <  2:487330/139‑1 (MQ=255)                           cGGTGGCTAACTGGGTGAAGAATCCTTCTTTATAATTCAGCGGGCCGCCGGTCCAGGTTTTGACATCGTCA‑TTGCTTCAGTATTACGAGACATAGTACCGCTTATACCCCAAAGCTTTCGCCACAGCCACATTCATTCt                                                                                            <  1:719420/139‑1 (MQ=255)                            ggTGGCTAACTGGGTGAAGAATCCTTCTTTATAATTCAGCGGGCCGCCGGTCCAGGTTTTGACATCGTCA‑TTGCTTCAGTATTACGAGACATAGTACCGCTTATACCCCAAAGCTTTCGCCACAGCCACATTCATTCTg                                                                                           <  1:685121/139‑1 (MQ=255)                                     cTGGGTGAAGAATCCTTCTTTATAATTCAGCGGGCCGCCGGTCCAGGTTTTGACATCGTCA‑TTGCTTCAGTATTACGAGACATAGTACCGCTTATACCCCAAAGCTTTCGCCACAGCCACATTCATTCTGGGCTTTAgg                                                                                  <  2:647151/139‑1 (MQ=255)                                                              ttCAGCGGGCCGCCGGTCCAGGTTTTGACATCGTCA‑TTGCTTCAGTATTACGAGACATAGTACCGCTTATACCCCAAAGCTTTCGCCACAGCCACATTCATTCTGGGCTTTAGGGTTGTGAAATTTGAATATCTGATTa                                                         >  2:816707/1‑139 (MQ=255)                                                                 aGCGGGCCGCCGGTCCAGGTTTTGACATCGTCA‑TTGCTTCAGTATTACGAGACATAGTACCGCTTATACCCCAAAGCTTTCGCCACAGCCACATTCATTCTGGGCTTTAGGGTTGTGAAATTTGAATATCTGATTAAGt                                                      >  1:67662/1‑139 (MQ=255)                                                                 aGCGGGCCGCCGGTCCAGGTTTTGACATCGTCA‑TTGCTTCAGTATTACGAGACATAGTACCGCTTATACCCCAAAGCTTTCGCCACAGCCACATTCATTCTGGGCTTTAGGGTTGTGAAATTTGAATATCTGATTAAGt                                                      >  2:185121/1‑139 (MQ=255)                                                                     ggCCGCCGGTCCAGGTTTTGACATCGTCA‑TTGCTTCAGTATTACGAGACATAGTACCGCTTATACCCCAAAGCTTTCGCCACAGCCACATTCATTCTGGGCTTTAGGGTTGTGAAATTTGAATATCTGATTAAGTCCtt                                                  >  1:1016528/1‑139 (MQ=255)                                                                                     tttGACATCGTCA‑TTGCTTCAGTATTACGAGACATAGTACCGCTTATACCCCAAAGCTTTCGCCACAGCCACATTCATTCTGGGCTTTAGGGTTGTGAAATTTGAATATCTGATTAAGTCCTTCACGAACGAAATCGAc                                  >  1:1079879/1‑139 (MQ=255)                                                                                        gACATCGTCA‑TTGCTTCAGTATTACGAGACATAGTACCGCTTATACCCCAAAGCTTTCGCCACAGCCACATTCATTCTGGGCTTTAGGGTTGTGAAATTTGAATATc                                                               <  1:571236/107‑1 (MQ=255)                                                                                        gACATCGTCA‑TTGCTTCAGTATTACGAGACATAGTACCGCTTATACCCCAAAGCTTTCGCCACAGCCACATTCATTCTGGGCTTTAGGGTTGTGAAATTTGAATATc                                                               >  2:571236/1‑107 (MQ=255)                                                                                         aCATCGTCA‑TTGCTTCAGTATTACGAGACATAGTACCGCTTATACCCCAAAGCTTTCGCCACAGCCACATTCATTCTGGGCTTTAGGGTTGTGAAATTTGAATATCTGATTAAGTCCTTCACGAACGAAATCGACTTcc                              <  2:67662/139‑1 (MQ=255)                                                                                            tCGTCA‑TTGCTTCAGTATTACGAGACATAGTACCGCTTATACCCCAAAGCTTTCgcc                                                                                                             >  2:1134063/1‑57 (MQ=255)                                                                                            tCGTCA‑TTGCTTCAGTATTACGAGACATAGTACCGCTTATACCCCAAAGCTTTCgcc                                                                                                             <  1:1134063/57‑1 (MQ=255)                                                                                                 cattGCTTCAGTATTACGAGACATAGTACCGCTTATACCCCAAAGCTTTCGCCACAGCCACATTCATTCTGGGCTTTAGGGTTGTGAAATTTGAATATCTGATTAAGTCCTTCACGAACGAAATCGACTTCCGTGCCAt                       <  1:816707/137‑1 (MQ=255)                                                                                                         cAGTATTACGAGACATAGTACCGCTTATACCCCAAAGCTTTCGCCACAGCCACATTCATTCTGGGCTTTAGGGTTGTGAAATTTGAATATCTGATTAAGTCCTTCACGAACGAAATCGACTTCCGTGCCATCAATAAAc               >  2:411065/1‑139 (MQ=255)                                                                                                                 cGAGACATAGTACCGCTTATACCCCAAAGCTTTCGCCACAGCCACATTCATTCTGGGCTTTAGGGTTGTGAAATTTGAATATCTGATTAAGTCCTTCACGAACGAAATCGACTTCCGTGCCATCAATAAACGGCATCGc       >  1:475737/1‑139 (MQ=255)                                                                                                                  gagaCATAGTACCGCTTATACCCCAAAGCTTTCGCCACAGCCACATTCATTCTGGGCTTTAGGGTTGTGAAATTTGAATATCTGATTAAGTCCTTCACGAACGAAATCGACTTCCGTGCCATCAATAAACGGCATCGCt      <  2:475737/139‑1 (MQ=255)                                                                                                                    gaCATAGTACCGCTTATACCCCAAAGCTTTCGCCACAGCCACATTCATTCTGGGCTTTAGGGTTGTGaa                                                                          >  1:1229775/1‑69 (MQ=255)                                                                                                                    gaCATAGTACCGCTTATACCCCAAAGCTTTCGCCACAGCCACATTCATTCTGGGCTTTAGGGTTGTGaa                                                                          <  2:1229775/69‑1 (MQ=255)                                                                                                                      cATAGTACCGCTTATACCCCAAAGCTTTCGCCACAGCCACATTCATTCTGGGCTTTAGGGTTGTGAAATTTGAATATCTGATTAAGTCCTTCACGAACGAAATCGACTTCCGTGCCATCAATAAACGGCATCGCTTGCa  <  2:1016528/139‑1 (MQ=255)                                                                                                                        tAGTACCGCTTATACCCCAAAGCTTTCGCCACAGCCAcat                                                                                                   >  1:706583/1‑40 (MQ=255)                                                                                                                        tAGTACCGCTTATACCCCAAAGCTTTCGCCACAGCCAcat                                                                                                   <  2:706583/40‑1 (MQ=255)                                                                                                                        tAGTACCGCTTATACCCCAAAGCTTTCGCCACAGCCACATTCATTCTGGGCTTTAGGGTTGTGAAATTTGAATATCTGATTAAGTCCTTCACGAACGAAATCGACTTcc                              <  1:1270403/109‑1 (MQ=255)                                                                                                                        tAGTACCGCTTATACCCCAAAGCTTTCGCCACAGCCACATTCATTCTGGGCTTTAGGGTTGTGAAATTTGAATATCTGATTAAGTCCTTCACGAACGAAATCGACTTcc                              >  2:1270403/1‑109 (MQ=255)                                                                                                                         aGTACCGCTTATACCCCAAAGCTTTCGCCACAGCCACattcattc                                                                                             >  2:1022969/1‑45 (MQ=255)                                                                                                                         aGTACCGCTTATACCCCAAAGCTTTCGCCACAGCCACattcattc                                                                                             <  1:1022969/45‑1 (MQ=255)                                                                                                                                 |                                                                                                                                TTCGTTTATCCCCTTTGCCAGCTCATCGGTGGCTAACTGGGTGAAGAATCCTTCTTTATAATTCAGCGGGCCGCCGGTCCAGGTTTTGACATCGTCAGTTGCTTCAGTATTACGAGACATAGTACCGCCTATACCCCAAAGCTTTCGCCACAGCCACATTCATTCTGGGCTTTAGGGTTGTGAAATTTGAATATCTGATTAAGTCCTTCACGAACGAAATCGACTTCCGTGCCATCAATAAACGGCATCGCTTGCA  >  NC_000913/1763890‑1764145

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