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

CCCTTTGCCAGCTCATCGGTGGCTAACTGGGTGAAGAATCCTTCTTTATAATTCAGCGGGCCGCCGGTCCAGGTTTTGACATCGTCAGTTGCTTCAGTATTACGAGACATAGTACCGCCTATACCCCAAAGCTTTCGCCACAGCCACATTCATTCTGGGCTTTAGGGTTGTGAAATTTGAATATCTGATTAAGTCCTTCACGAACGAAATCGACTTCCGTGCCATCAATAAACGGCATCGCTTGCAGCG  >  NC_000913/1763900‑1764148                                                                                                                       |                                                                                                                                   cccTTTGCCAGCTCGTCGGTGGCTAACTGGGTGAAGAATCCTTCTTTATAAGTCAGCGGGCCGCCGGTCCAGGTTTTGACATCGTCA‑TTGCTTCAGTATTACGAGACATAGTACCGCTTATACCCCAAAGCTTTCgcca                                                                                                               <  2:290365/139‑1 (MQ=255) cccTTTGCCAGCTCATCGGTGGCTAACTGGGTGAAGAATCCTTCTTTATAATTCAGCGGGCCGCCGGTCCAGGTTTTGACATCGTCA‑TTGCTTCAGTATTACGAGACATAGTACCGCTTATACCCCAAAGCTTTCgcca                                                                                                               <  2:645068/139‑1 (MQ=255)  ccTTTGCCAGCTCATCGGTGGCTAACTGGGTGAAGAATCCTTCTTTATAATTCAGCGGGCCGCCGGTCCAGGTTTTGACATCGTCA‑TTGCTTCAGTATTACGAGACATAGTACCGCTTATACCCCAAAGCTTTCgccac                                                                                                              <  2:1644412/139‑1 (MQ=255)     ttGCCAGCTCATCGGTGGCTAACTGGGTGAAGAATCCTTCTTTATAATTCAGCGGGCCGCCGGTCCAGGTTTTGACATCGTCA‑TTGCTTCAGTATTACGAGACATAGTACCGCTTATACCCCAAAGCTTTCgccacagc                                                                                                           <  2:17287/139‑1 (MQ=255)     ttGCCAGCTCATCGGTGGCTAACTGGGTGAAGAATCCTTCTTTATAATTCAGCGGGCCGCCGGTCCAGGTTTTGACATCGTCA‑TTGCTTCAGTATTACGAGACATAGTACCGCTTATACCCCAAAGCTTTCgccacagc                                                                                                           <  1:520368/139‑1 (MQ=255)                 cGGTGGCTAACTGGGTGAAGAATCCTTCTTTATAATTCAGCGGGCCGCCGGTCCAGGTTTTGACATCGTCA‑TTGCTTCAGTATTACGAGACATAGTACCGCTTATACCCCAAAGCTTTCGCCACAGCCACATTCATTCt                                                                                               >  2:1595581/1‑139 (MQ=255)                 cGGTGGCTAACTGGGTGAAGAATCCTTCTTTATAATTCAGCGGGCCGCCGGTCCAGGTTTTGACATCGTCA‑TTGCTTCAGTATTACGAGACATAGTACCGCTTATACCCCAAAGCTTTCGCCACAGCCACATTCATTCt                                                                                               >  2:1283616/1‑139 (MQ=255)                  ggTGGCTAACTGGGTGAAGAATCCTTCTTTATAATTCAGCGGGCCGCCGGTCCAGGTTTTGACATCGTCA‑TTGCTTCAGTATTACGAGACATAGTACCGCTTATACCCCAAAGCTTTCGCCACAGCCACATTCATTCTg                                                                                              <  2:739919/139‑1 (MQ=255)                  ggTGGCTAACTGGGTGAAGAATCCTTCTTTATAATTCAGCGGGCCGCCGGTCCAGGTTTTGACATCGTCA‑TTGCTTCAGTATTACGAGACATAGTACCGCTTATACCCCAAAGCTTTCGCCACAGCCACATTCATTCTg                                                                                              <  2:96098/139‑1 (MQ=255)                  ggTGGCTAACTGGGTGAAGAATCCTTCTTTATAATTCAGCGGGCCGCCGGTCCAGGTTTTGACATCGTCA‑TTGCTTCAGTATTACGAGACATAGTACCGCTTATACCCCAAAGCTTTCGCCACAGCCACATTCATTCTg                                                                                              <  2:1300934/139‑1 (MQ=255)                          aCTGGGTGAAGAATCCTTCTTTATAATTCAGCGGGCCGCCGGTCCAGGTTTTGACATCGTCA‑TTGCTTCAGTATTACGAGACATAGTACCGCTTATACCCCAAAGCTTTCGCCACAGCCACATTCATTCTGGGCTTTAg                                                                                      >  2:135961/1‑139 (MQ=255)                              ggTGAAGAATCCTTCTTTATAATTCAGCGGGCCGCCGGTCCAGGTTTTGACATCGTCA‑TTGCTTCAGTATTACGAGACATAGTACCGCTTATACCCCa                                                                                                                           >  1:1229826/1‑98 (MQ=255)                              ggTGAAGAATCCTTCTTTATAATTCAGCGGGCCGCCGGTCCAGGTTTTGACATCGTCA‑TTGCTTCAGTATTACGAGACATAGTACCGCTTATACCCCa                                                                                                                           <  2:1229826/98‑1 (MQ=255)                                tGAAGAATCCTTCTTTATAATTCAGCGGGCCGCCGGTCCAGGTTTTGACATCGTCA‑TTGCTTCAGTATTACGAGACATAGTACCGCTTATACCCCAAAGCTTTCGCCACAGCCACATTCATTCTGGGCTTTAGGGTtgt                                                                                >  1:1109128/1‑139 (MQ=255)                                         cttcttTATAATTCAGCGGGCCGCCGGTCCAGGTTTTGACATCGTCA‑TTGCTTCAGTATTACGAGACATAGTACCGCTCATACCCCAAAGCTTTCGCCACAGCCAcat                                                                                                      >  1:1249751/1‑108 (MQ=255)                                         cttcttTATAATTCAGCGGGCCGCCGGTCCAGGTTTTGACATCGTCA‑TTGCTTCAGTATTACGAGACATAGTACCGCTCATACCCCAAAGCTTTCGCCACAGCCAcat                                                                                                      <  2:1249751/108‑1 (MQ=255)                                          ttcttTATAATTCAGCGGGCCGCCGGTCCAGGTTTTGACATCGTCA‑TTGCTTCAGTATTACGAGACATAGTACCGCTTATACCCCAAAGCTTTCGCCACAGCCACATTCATTCTGGGCTTTAGGGTTGTGAAATTTGaa                                                                      >  1:1175367/1‑139 (MQ=255)                                            cttTATAATTCAGCGGGCCGCCGGTCCAGGTTTTGACATCGTCA‑TTGCTTCAGTATTACGAGACATAGTACCGCTTATACCCCAAAGCTTTCGCCACAGCCACATTCATTCTGGGCTTTAGGGTTGTGAAATTTGAata                                                                    <  2:796642/139‑1 (MQ=255)                                             tttATAATTCAGCGGGCCGCCGGTCCAGGTTTTGACATCGTCA‑TTGCTTCAGTATTACGAGACATAGTACCGCTTATACCCCAAAGCTTTCGCCACAGCCACATTCATTCTGGGCTTTAGGGTTGTGAAATTTGAatat                                                                   <  2:1175367/139‑1 (MQ=255)                                              ttATAATTCAGCGGGCCGCCGGTCCAGGTTTTGACATCGTCA‑TTGCTTCAGTATTACGAGACATAGTACCGCTTATACCCCAAAGCTTTCGCCACAGCCACATTCATTCTGGGCTTTAGGGTTGTGAAATTTGAATATc                                                                  <  1:135961/139‑1 (MQ=255)                                                         cGGGCCGCCGGTCCAGGTTTTGACATCGTCA‑TTGCTTCAGTATTACGAGACATAGTACCGCTTATACCCCAGAGCTTTCGCCACAGCCACATTCATTCTGGGCTTTAGGGTTGTGAAATTTGAATATCTGATTAAGTcc                                                       >  2:948535/1‑139 (MQ=255)                                                         cGGGCCGCCGGTCCAGGTTTTGACATCGTCA‑TTGCTTCAGTATTACGAGACATAGTACCGCTTATACCCCAAAGCTTTCgccacagc                                                                                                           <  1:1235120/87‑1 (MQ=255)                                                         cGGGCCGCCGGTCCAGGTTTTGACATCGTCA‑TTGCTTCAGTATTACGAGACATAGTACCGCTTATACCCCAAAGCTTTCgccacagc                                                                                                           >  2:1235120/1‑87 (MQ=255)                                                         cGGGCCGCCGGTCCAGGTTTTGACATCGTCA‑TTGCTTCAGTATTACGAGACATAGTACCGCTTATACCCCAAAGCTTTCGCCACAGCCACATTCATTCTGGGCTTTAGGGTTGTGAAATTTGAATATCTGATTAAGTcc                                                       >  2:657165/1‑139 (MQ=255)                                                         cGGGCCGCCGGTCCAGGTTTTGACATCGTCA‑TTGCTTCAGTATTACGAGACATAGTACCGCTTATACCCCAAAGCTTTCGCAACAGCCACATTCATTCTGGGCtt                                                                                         >  2:1532070/1‑105 (MQ=255)                                                         cGGGCCGCCGGTCCAGGTTTTGACATCGTCA‑TTGCTTCAGTATTACGAGACATAGTACCGCTTATACCCCAAAGCTTTCGCAACAGCCACATTCATTCTGGGCtt                                                                                         <  1:1532070/105‑1 (MQ=255)                                                             ccgccgGTCCAGGTTTTGACATCGTCA‑TTGCTTCAGTATTACGAGACATAGTACCGCTTATACCCCAAAGCTTTCGCCACAGCCACATTCATTCTGGGCTTTAGGGTTGTGAAATTTGAATATCTGATTAAGTCCTTCa                                                   <  1:657165/139‑1 (MQ=255)                                                              cgccgGTCCAGGTTTTGACATCGTCA‑TTGCTTCAGTATTACGAGACATAGTACCGCTTATACCCCAAAGCTTTCGCCACAGCCACATTCATTCTGGGCTTTAGGGTTGTGAAATTTGAATATCTGATTAAGTCCTTCac                                                  >  1:478361/1‑139 (MQ=255)                                                              cgccgGTCCAGGTTTTGACATCGTCA‑TTGCTTCAGTATTACGAGACATAGTACCGCTTATACCCCAAAGCTTTCGCCACAGCCACATTCATTCTGGGCTTTAGGGTTGTGAAATTTGAATATCTGATTAAGTCCTTCac                                                  >  2:1481033/1‑139 (MQ=255)                                                                     ccAGGTTTTGACATCGTCA‑TTGCTTCAGTATTACGAGACATAGTACCGCTTATACCCCAAAGCTTTCGCCACAGCCACATTCATTCTGGGCTTTAGGGTTGTGAAATTTGAATATCTGATTAAGTCCTTCAcgaacgaa                                           >  1:1551910/1‑139 (MQ=255)                                                                             tGACATCGTCA‑TTGCTTCAGTATTACGAGACATAGTACCGCTTATACCCCAAAGCTTTCGCCACAGCCACATTCATTCTGGGCTTTAGGGTTGTGAAATTTGAat                                                                     >  1:253775/1‑105 (MQ=255)                                                                             tGACATCGTCA‑TTGCTTCAGTATTACGAGACATAGTACCGCTTATACCCCAAAGCTTTCGCCACAGCCACATTCATTCTGGGCTTTAGGGTTGTGAAATTTGAat                                                                     <  2:253775/105‑1 (MQ=255)                                                                                  tCGTCA‑TTGCTTCAGTATTACGAGACATAGTACCGCTTATACCCCAAAGCTTTCGCCACAGCCACATTCATTCTGGGCTTTAGGGTTGTGAAATTTGAATATCTGATTAAGTCCTTCACGAACGAAATCGACTTCCGTg                              >  1:73579/1‑139 (MQ=255)                                                                                    cgtcattGCTTCAGTATTACGAGACATAGTACCGCTTATACCCCAAAGCTTTCgccacagc                                                                                                           <  2:1453097/56‑1 (MQ=255)                                                                                    cgtcattGCTTCAGTATTACGAGACATAGTACCGCTTATACCCCAAAGCTTTCgccacagc                                                                                                           >  1:1453097/6‑61 (MQ=255)                                                                                    cgtcattGCTTCAGTATTACGAGACATAGTACCGCTTATACCCCAAAGCTTTCGCCACAGCCACATTCATTCTGGGCTTTAGGGTTGTGAAATTTGAATATCTGATTAAGTCCTTCACGAACGAAATCGACTTCCGTGc                             <  1:1283616/134‑1 (MQ=255)                                                                                    cgtcattGCTTCAGTATTACGAGACATAGTACCGCTTATACCCCAAAGCTTTCGCCACAGCCACATTCATTCTGGGCTTTAGGGTTGTGAAATTTGAATATCTGATTAAGTCCTTCACGAACGAAATCGACTTCCGTGc                             >  1:1465027/6‑139 (MQ=255)                                                                                      tcattGCTTCAGTATTACGAGACATAGTACCGCTTATACCCCAAAGCTTTCGCCACAGCCACATTCATTCTGGGCTTTAGGGTTGTGAAATTTGAATATCTGATTAAGTCCTTCACGAACGAAATCGACTTCCGTGCCa                           <  1:1481033/136‑1 (MQ=255)                                                                                        attGCTTCAGTATTACGAGACATAGTACCGCTTATACCCCAAAGCTTTCGCCACAGCCACATTCATTCTGGGCTTTAGGGTTGTGAAATTTGAATATCTGATTAAGTCCTTCacgaacg                                             >  2:331586/2‑119 (MQ=255)                                                                                        attGCTTCAGTATTACGAGACATAGTACCGCTTATACCCCAAAGCTTTCGCCACAGCCACATTCATTCTGGGCTTTAGGGTTGTGAAATTTGAATATCTGATTAAGTCCTTCacgaacg                                             <  1:331586/118‑1 (MQ=255)                                                                                               cAGTATTACGAGACATAGTACCGCTTATACCCCAAAGCTTTCgccacagc                                                                                                           <  2:1647118/50‑1 (MQ=255)                                                                                               cAGTATTACGAGACATAGTACCGCTTATACCCCAAAGCTTTCgccacagc                                                                                                           >  1:1647118/1‑50 (MQ=255)                                                                                                  tATTACGAGACATAGTACCGCTTATACCCCAAAGCTTTCGCCACAGCCACATTCATTCTGGGCTTTAGGGTTGTGAAATTTGAATATCTGATTAAGTCCTTCACGAACGAAATCGACTTCCGTGCCATCAATAAACGGc               <  1:1595581/139‑1 (MQ=255)                                                                                                  tATTACGAGACATAGTACCGCTTATACCCCAAAGCTTTCGCCACAGCCACATTCATTCTGGGCTTTAGGGTTGTGAAATTTGAATATCTGATTAAGTCCTTCACGAACGAAATCGACTTCCGTGCCATCAATAAACGGc               <  2:530545/139‑1 (MQ=255)                                                                                                    ttACGAGACATAGTACCGCTTATACCCCAAAGCTTTCGCCACAGCCACATTCATTCTGGGCTTTAGGGTTGTGAAATTTGAATATCTGATTAAGTCCTTCACGAACGAAATCGACTTCCGTGCCATCAATAAACGGCAt             >  1:1275661/1‑139 (MQ=255)                                                                                                    ttACGAGACATAGTACCGCTTATACCCCAAAGCTTTCGCCACAGCCACATTCATTCTGGGCTTTAGGGTTGTGAAATTTGAATATCTGATTAAGTCCTTCACGAACGAAATCGACTTCCGTGCCATCAATAAACGGCAt             <  1:1196147/139‑1 (MQ=255)                                                                                                            cATAGTACCGCTTATACCCCAAAGCTTTCTCCACAGCCACCTTCACTCTGGGTTTTAGTGTTGTGAAATTTGAATATCTGATTAAGTCCTTCACGAACGAAATCGACTTCCGTGCCATCAATAAACGGCATCGCTTGCa     >  2:1396101/1‑139 (MQ=255)                                                                                                            cATAGTACCGCTTATACCCCAAAGCTTTCGCCACAGCCACATTCATTCTGGGCTTTAGGGTTGTGAAATTTGAat                                                                     >  2:718721/1‑75 (MQ=255)                                                                                                            cATAGTACCGCTTATACCCCAAAGCTTTCGCCACAGCCACATTCATTCTGGGCTTTAGGGTTGTGAAATTTGAat                                                                     <  1:718721/75‑1 (MQ=255)                                                                                                            cATAGTACCGCTTATACCCCAAAGCTTTCGCCACAGCCACATTCATTCTGGGCTTTAGGGTTGTGAAATTTGAATATCTGa                                                               >  1:893789/1‑81 (MQ=255)                                                                                                            cATAGTACCGCTTATACCCCAAAGCTTTCGCCACAGCCACATTCATTCTGGGCTTTAGGGTTGTGAAATTTGAATATCTGa                                                               <  2:893789/81‑1 (MQ=255)                                                                                                            cATAGTACCGCTTATACCCCAAAGCTTTCGCCACAGCCACATTCATTCTGGGCTTTAGGGTTGTGAAATTTGAATATCTGATTAAGTCCTTCACGAACGAAATCGACTTCCGTGCCATCAATAAACGGCATCGCTTGCa     <  2:367287/139‑1 (MQ=255)                                                                                                            cATAGTACCGCTTATACCCCAAAGCTTTCGCCACAGCCACATTCATTCTGGGCTTTAGGGTTGTGAAATTTGAATATCTGATTAAGTCCTTCACGAACGAAATCGACTTCCGTGCCATCAATAAACGGCATCGCTTGCa     >  1:367287/1‑139 (MQ=255)                                                                                                              tAGTACCGCTTATACCCCAAAGCTTTCGCCACAGCCACATTCATTCTGGGCTTTAGGGTTGTGAAATTTGAATATCTGATTAAGTCCTTCACGAACGAAATCGACTTCCGTGCCATCAATa                     >  2:1599291/1‑121 (MQ=255)                                                                                                              tAGTACCGCTTATACCCCAAAGCTTTCGCCACAGCCACATTCATTCTGGGCTTTAGGGTTGTGAAATTTGAATATCTGATTAAGTCCTTCACGAACGAAATCGACTTCCGTGCCATCAATa                     <  1:1599291/121‑1 (MQ=255)                                                                                                               aGTACCGCTTATACCCCAAAGCTTTCGCCACAGCCACATTCATTCTGGGCTTTAGGGTTGTGAAATTTGAATATCTGATTAAGTCCTTCACGAACGAAATCGACTTCCGTGCCATCAATAAACGGCATCGCTTGCAGCg  <  2:630660/139‑1 (MQ=255)                                                                                                                       |                                                                                                                                   CCCTTTGCCAGCTCATCGGTGGCTAACTGGGTGAAGAATCCTTCTTTATAATTCAGCGGGCCGCCGGTCCAGGTTTTGACATCGTCAGTTGCTTCAGTATTACGAGACATAGTACCGCCTATACCCCAAAGCTTTCGCCACAGCCACATTCATTCTGGGCTTTAGGGTTGTGAAATTTGAATATCTGATTAAGTCCTTCACGAACGAAATCGACTTCCGTGCCATCAATAAACGGCATCGCTTGCAGCG  >  NC_000913/1763900‑1764148

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: ‑