Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	323,355	C→T	P340S (CCA→TCA)	ykgF →	ferridoxin‑like LutB family protein, putative electron transport chain YkgEFG component

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	323,355	0	C	T	100.0%	125.3 / NA	40	P340S (CCA→TCA)	ykgF	ferridoxin‑like LutB family protein, putative electron transport chain YkgEFG component
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (19/21);  total (19/21)

CGGGCGTTCTGAGGTGCTGGCCTCTGAATTTCGGGATGTGCTGCGCTGTATTCGCTGCGGGGCTTGTATGAATACTTGTCCGGCATATCGCCATATTGGCGGTCATGGATATGGCTCTATTTATCCAGGGCCAATTGGTGCGGTGATTTCTCCGCTACTTGGCGGCTATAAAGATTTTAAAGATTTACCCTACGCCTGCTCTTTATGCACAGCTTGTGACAACGTGTGTCCGGTGCG  >  NC_000913/323225‑323461                                                                                                                                   |                                                                                                           cGGGCGTTCTGAGGTGCTGGCCTCTGAATTTCGGGATGTGCTGCGCTGTATTCGCTGCGGGGCTTGTATGAATACTTGTCCGACATATCGCCATATTGGCGGTCATGGATATGGCTCTATTTATCCAGGGTCAATTGGt                                                                                                    <  1:587529/139‑1 (MQ=255)           gAGGTGCTGGCCTCTGAATTTCGGGATGTGCTGCGCTGTATTCGCTGCGGGGCTTGTATGAATACTTGTCCGGCATATCGCCATATTGGCGGTCATGGATATGGCTCTATTTATCCAGGGTCAATTGGTGCGGTGAttt                                                                                          <  2:361997/139‑1 (MQ=255)                 cTGGCCTCTGAATTTCGGGATGTGCTGCGCTGTATTAGCTGCGGGGCGTGTATGAATACTTGTACGGCATATCGCCATATTGGCGGTGATGGATATGGCTCTATTTATCCAGGGTCAATTGGTGCGGTGATTTCTCCGc                                                                                    >  1:857774/1‑139 (MQ=255)                  tGGCCTCTGAATTTCGGGATGTGCTGCGCTGTATTCGCTGCGGGGCTTGTATGAATACTTGTCCGGCATATCGCCATATTGGCGGTCATGGATATGGCTCTATTTATCCAGGGTCAATTGGTGCGGTGATTTCTCCGCt                                                                                   >  1:1121592/1‑139 (MQ=255)                      ctctGAATTTCGGGATGTGCTGCGCTGTATTCGCTGCGGGGCTTGTATGAATACTTGTCCGGCATATCGCCATATTGGCGGTCATGGATATGGCTCTATTTATCCAGGGTCAAt                                                                                                        >  2:203730/1‑114 (MQ=255)                      ctctGAATTTCGGGATGTGCTGCGCTGTATTCGCTGCGGGGCTTGTATGAATACTTGTCCGGCATATCGCCATATTGGCGGTCATGGATATGGCTCTATTTATCCAGGGTCAAt                                                                                                        <  1:203730/114‑1 (MQ=255)                      ctctGAATTTCGGGATGTGCTGCGCTGTATTCGCTGCGGGGCTTGTATGAATACTTGTCCGGCATATCGCCATATTGGCGGTCATGGATATGGCTCTATTTATCCAGGGTCAATTGGTGCGGTGATTTCTCCGCTACtt                                                                               >  2:1150184/1‑139 (MQ=255)                          gAATTTCGGGATGTGCTGCGCTGTATTCGCTGCGGGGCTTGTATGAATACTTGTCCGGCATATCGCCATATTGGCGGTCATGGATATGGCTCTATTTATCCAGGGTCAATTGGTGCGGTGATTTCTCCGCTACTTggcg                                                                           <  1:1150184/139‑1 (MQ=255)                                                      gCTGCGGGGCTTGTATGAATACTTGTCCGGCATATCGCCATATTGGCGGTCATGGATATGGCTCTATTTATCCAGGGTCAATTGGTGCGGTGATTTCTCCGCTACTTGGCGGCTATAAAGATTTTAAAGATTTACCCTa                                               <  2:1121592/139‑1 (MQ=255)                                                           ggggCTTGTATGAATACTTGTCCGGCTTATCGCCATAATGGCGGTCATGGATATGGCTCTATTTATCCAGGGTCAATTGGTGCGGTGATTTCTCCGCTACTTGGCGGCTATAAAGATTTTAAAGATTTACCCTACGCCt                                          <  2:857774/139‑1 (MQ=255)                                                           ggggCTTGTATGAATACTTGTCCGGCATATCGCCATATTGGCGGTCATGGATATGGCTCTATTTATCCAGGGTCAATTGGTGCGGTGATTTCTCCGCTACTTGGCGGCTAtaaagatttta                                                            >  1:370189/1‑121 (MQ=255)                                                           ggggCTTGTATGAATACTTGTCCGGCATATCGCCATATTGGCGGTCATGGATATGGCTCTATTTATCCAGGGTCAATTGGTGCGGTGATTTCTCCGCTACTTGGCGGCTAtaaagatttta                                                            <  2:370189/121‑1 (MQ=255)                                                               cTTGTATGAATACTTGTCCGGCATATCGCCATATTGGCGGTCATGGATATGGCTCTATTTATCCAGGGTCAATTGGTGCGGTGATTTCTCCGCTACTTGGCGGCTATAAAGATTTTAAAGATTTACCCTACGCCTGctc                                      <  1:485898/139‑1 (MQ=255)                                                                  gTATGAATACTTGTCCGGCATATCGCCATATTGGCGGTCATGGATATGGCTCTATTTATCCAGGGTCAATTGGTGCGGTGATTTCTCCGCTACTTGGCGGCTATAAAGATTTTAAAGATTTACCCTACGCCTGCTCttt                                   >  2:110072/1‑139 (MQ=255)                                                                             tGTCCGGCATATCGCCATATTGGCGGTCATGGATATGGCTCTATTTATCCAGGGTCAATTGGTGCGGTGAttt                                                                                          >  2:775635/1‑73 (MQ=255)                                                                             tGTCCGGCATATCGCCATATTGGCGGTCATGGATATGGCTCTATTTATCCAGGGTCAATTGGTGCGGTGAttt                                                                                          <  1:775635/73‑1 (MQ=255)                                                                             tGTCCGGCATATCGCCATATTGGCGGTCATGGATATGGCTCTATTTATCCAGGGTCAATTGGTGCGGTGATTTCTCCGCTACTTGGCGGCTATAAAGATTTTAAAGATTTACCCTACGCCTGCTCTTTATGca                              >  2:1298742/1‑133 (MQ=255)                                                                             tGTCCGGCATATCGCCATATTGGCGGTCATGGATATGGCTCTATTTATCCAGGGTCAATTGGTGCGGTGATTCCTCCGCTACTTGGCGGCTATAAAGATTTTAAAGATTTACCCTACGCCTGCTCTTTATGca                              <  1:1298742/133‑1 (MQ=255)                                                                              gTCCGGCATATCGCCATATTGGCGGTCATGGATATGGCTCTATTTATCCAGGGTCAATTGGTGCGGTGATTTCTCCGCTACTTGGCGGCTATAAAGATTTTAAAGATTTACCCTACGCCTGCTCTTTATGca                              <  1:1064724/132‑1 (MQ=255)                                                                              gTCCGGCATATCGCCATATTGGCGGTCATGGATATGGCTCTATTTATCCAGGGTCAATTGGTGCGGTGATTTCTCCGCTACTTGGCGGCTATAAAGATTTTAAAGATTTACCCTACGCCTGCTCTTTATGca                              >  2:1064724/1‑132 (MQ=255)                                                                               tCCGGCATATCGCCATATTGGCGGTCATGGATATGGCTCTATTTATCCAGGGTCAATTGGTGCGGTGATTTCTCCGCTACTTGGCGGCTAtaaagattttaaagatt                                                      <  1:403755/107‑1 (MQ=255)                                                                               tCCGGCATATCGCCATATTGGCGGTCATGGATATGGCTCTATTTATCCAGGGTCAATTGGTGCGGTGATTTCTCCGCTACTTGGCGGCTAtaaagattttaaagatt                                                      >  2:403755/1‑107 (MQ=255)                                                                                              tatTGGCGGTCATGGATATGGCTCTATTTATCCAGGGTCAATTGGTGCGGTGATTTCTCCGCTACTTGGCGGCTATAAAGATTTTAAAGATTTACCCTACGCCTGc                                        <  1:610509/106‑1 (MQ=255)                                                                                              tatTGGCGGTCATGGATATGGCTCTATTTATCCAGGGTCAATTGGTGCGGTGATTTCTCCGCTACTTGGCGGCTATAAAGATTTTAAAGATTTACCCTACGCCTGc                                        >  2:610509/1‑106 (MQ=255)                                                                                                   gCGGTCATGGATATGGCTCTATTTATCCAGGGTCAATTGGTGCGGTGATTTCTCCGCTACTTGGCGGCTATAAAGATTTTAAAGATTTACCCTACGCCTGCTCTTTATGCACAGCTTGTGACAACGTGTGTCCGGTGCg  >  1:103483/1‑139 (MQ=255)                                                                                                         aTGGATATGGCTCTATTTATCCAGGGTCAATTGGTGCGGTGATTTCTCCGCTACTTGGCGGCTATAAAGATTTTAAAGATTTACCCTACGcc                                           >  2:282831/1‑92 (MQ=255)                                                                                                         aTGGATATGGCTCTATTTATCCAGGGTCAATTGGTGCGGTGATTTCTCCGCTACTTGGCGGCTATAAAGATTTTAAAGATTTACCCTACGcc                                           <  1:282831/92‑1 (MQ=255)                                                                                                                 ggCTCTATTTATCCAGGGTCAATTGGTGCGGTGATTTCTCCGCTACTTGGCGGCtata                                                                     <  2:113244/58‑1 (MQ=255)                                                                                                                 ggCTCTATTTATCCAGGGTCAATTGGTGCGGTGATTTCTCCGCTACTTGGCGGCtata                                                                     >  1:113244/1‑58 (MQ=255)                                                                                                                 ggCTCTATTTATCCAGGGTCAATTGGTGCGGTGATTTCTCCGCTACTTGGCGGCTATAAAGATTTTAAAGATTTACCCTACGCCt                                          <  2:427654/85‑1 (MQ=255)                                                                                                                 ggCTCTATTTATCCAGGGTCAATTGGTGCGGTGATTTCTCCGCTACTTGGCGGCTATAAAGATTTTAAAGATTTACCCTACGCCt                                          >  1:427654/1‑85 (MQ=255)                                                                                                                 ggCTCTATTTATCCAGGGTCAATTGGTGCGGTGATTTCTCCGCTACTTGGCGGCTATAAAGATTTTAAAGATTTACCCTACGCCTGCTCTTTATGCACAGCTTGTGACAACGTGTGTCCGGt     >  1:654468/1‑122 (MQ=255)                                                                                                                 ggCTCTATTTATCCAGGGTCAATTGGTGCGGTGATTTCTCCGCTACTTGGCGGCTATAAAGATTTTAAAGATTTACCCTACGCCTGCTCTTTATGCACAGCTTGTGACAACGTGTGTCCGGt     <  2:654468/122‑1 (MQ=255)                                                                                                                      tATTTATCCAGGGTCAATTGGTGCGGTGATTTCTCCGCTACTTGGCGGCTAtaaagatttta                                                            <  1:64900/62‑1 (MQ=255)                                                                                                                      tATTTATCCAGGGTCAATTGGTGCGGTGATTTCTCCGCTACTTGGCGGCTAtaaagatttta                                                            >  2:64900/1‑62 (MQ=255)                                                                                                                             ccAGGGTCAATTGGTGCGGTGATTTCTCCGCTACTTgg                                                                             >  1:653401/1‑38 (MQ=255)                                                                                                                             ccAGGGTCAATTGGTGCGGTGATTTCTCCGCTACTTgg                                                                             <  2:613429/38‑1 (MQ=255)                                                                                                                             ccAGGGTCAATTGGTGCGGTGATTTCTCCGCTACTTgg                                                                             <  2:653401/38‑1 (MQ=255)                                                                                                                             ccAGGGTCAATTGGTGCGGTGATTTCTCCGCTACTTgg                                                                             >  1:613429/1‑38 (MQ=255)                                                                                                                                gggTCAATTGGTGCGGTGATTTCTCCGCTACTTGGCGGCTATAAAGATTTTAAAGATTTAccc                                                 <  1:274919/63‑1 (MQ=255)                                                                                                                                   |                                                                                                           CGGGCGTTCTGAGGTGCTGGCCTCTGAATTTCGGGATGTGCTGCGCTGTATTCGCTGCGGGGCTTGTATGAATACTTGTCCGGCATATCGCCATATTGGCGGTCATGGATATGGCTCTATTTATCCAGGGCCAATTGGTGCGGTGATTTCTCCGCTACTTGGCGGCTATAAAGATTTTAAAGATTTACCCTACGCCTGCTCTTTATGCACAGCTTGTGACAACGTGTGTCCGGTGCG  >  NC_000913/323225‑323461

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