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%	149.4 / NA	48	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 (26/22);  total (26/22)

ATAACGGGCGTTCTGAGGTGCTGGCCTCTGAATTTCGGGATGTGCTGCGCTGTATTCGCTGCGGGGCTTGTATGAATACTTGTCCGGCATATCGCCATATTGGCGGTCATGGATATGGCTCTATTTATCCAGGGCCAATTGGTGCGGTGATTTCTCCGCTACTTGGCGGCTATAAAGATTTTAAAGATTTACCCTACGCCTGCTCTTTATGCACAGCTTGTGACAACGTGTGTCCGGTGCGTATTCCGCTGTCAAAACTGATTTTGCGTC  >  NC_000913/323221‑323490                                                                                                                                       |                                                                                                                                        aTAACGGGCGTTCTGAGGTGCTGGCCTCTGAATTTCGGGATGTGCTGCGCTGTATTCGCTGCGGGGCTTGTATGAATACTTGTCCGGCATATCGCCATATTGGCGGTCATGGATATGGCTCTATTTATCCAGGGTCAAt                                                                                                                                     >  2:304761/1‑139 (MQ=255) aTAACGGGCGTTCTGAGGTGCTGGCCTCTGAATTTCGGGATGTGCTGCGCTGTATTCGCTGCGGGGCTTGTATGAATACTTGTCCGGCATATCGCCATATTGGCGGTCATGGATATGGCTCTATTTATCCAGGGTCAAt                                                                                                                                     <  2:1557818/139‑1 (MQ=255)     cGGGCGTTCTGAGGTGCTGGCCTCTGAATTTCGGGATGTGCTGCGCTGTATTCGCTGCGGGGCTTGTATGAATACTTGTCCGGCATATCGCCATATTGGCGGTCATGGATATGGCTCTATTTATCCAGGGTCAATTGGt                                                                                                                                 <  2:904543/139‑1 (MQ=255)                     cTGGCCTCTGAATTTCGGGATGTGCTGCGCTGTATTCGCTGCGGGGCTTGTATGAATACTTGTCCGGCATATCGCCATATTGGCGGTCATGGATATGGCTCTATTTATCCAGGGTCAATTGGTGCGGTGATTTCTCCGc                                                                                                                 >  1:511428/1‑139 (MQ=255)                            ctGAATTTCGGGATGTGCTGCGCTGTATTCGCTGCGGGGCTTGTATGAATACTTGTCCGGCATATCGCCATATTGGCGGTCATGGATATGGCTCTATTTATCCAGGGTCAATTGGTGCGGTGATTTCTCCGCTACTTgg                                                                                                          <  2:511428/139‑1 (MQ=255)                             tGAATTTCGGGATGTGCTGCGCTGTATTCGCTGCGGGGCTTGTATGAATACTTGTCCGGCATATCGCCATATTGGCGGTCATGGATATGGCTCTATTTATCCAGGGTCAATTGGTGCGGTGATTTCTCCGCTACTTggc                                                                                                         <  1:304761/139‑1 (MQ=255)                             tGAATTTCGGGATGTGCTGCGCTGTATTCGCTGCGGGGCTTGTATGAATACTTGTCCAGCATATCGCCATATTGGCGGTCATGGATATGGCTCTATTTATCCAGGGTCAATTGGTGCGGTGATTTCTCCGCTACTTggc                                                                                                         <  1:179200/139‑1 (MQ=255)                                            gctgcGCTGTATTCGCTGCGGGGCTTGTATGAATACTTGTCCGGCATATCGCCATATTGGCGGTCATGGATATGGCTCTATTTATCCAGGGTCAATTGGt                                                                                                                                 >  2:1082885/1‑100 (MQ=255)                                            gctgcGCTGTATTCGCTGCGGGGCTTGTATGAATACTTGTCCGGCATATCGCCATATTGGCGGTCATGGATATGGCTCTATTTATCCAGGGTCAATTGGt                                                                                                                                 <  1:1082885/100‑1 (MQ=255)                                             ctgcGCTGTATTCGCTGCGGGGCTTGTATGAATACTTGTCCGGCATATCGCCATATTGGCGGTCATGGATATGGCTCTATTTATCCAGGGTCAATTGGTGCGGTGATTTCTCCGCTACTTGGCGGCTAtaaagatttta                                                                                         >  2:33340/1‑139 (MQ=255)                                             ctgcGCTGTATTCGCTGCGGGGCTTGTATGAATACTTGTCCGGCATATCGCCATATTGGCGGTCATGGATATGGCTCTATTTATCCAGGGTCAATTGGTGCGGTGATTTCTCCGCTACTTGGCGGCTAtaaagatttta                                                                                         <  1:33340/139‑1 (MQ=255)                                                     tATTCGCTGCGGGGCTTGTATGAATACTTGTCCGGCATATCGCCATATTGGCGGTCATGGATATGGCTCTATTTATCCAGGGTCAATTGGTGCGGTGATTTCTCCGCTACTTGGCGGCTATAAAGATTTTAAAGATTTa                                                                                 >  1:150582/1‑139 (MQ=255)                                                               ggggTTTGTATGAATACTTGTCCGGCATATCGCCATATTGGCGGTCATGGATATGGCTCTATTTATCCACGGTCAATTGGTGCGGTGATTTCTCCGCTACTTGGCGGCTATAAAGATTTTAAAGATTTACCCTACGCCt                                                                       >  2:1009427/1‑139 (MQ=255)                                                               ggggCTTGTATGAATACTTGTCCGGCATATCGCCATATTGGCGGTCATGGATATGGCTCTATTTATACAGGGTCAATTGGTGCGGTGATTTCTCCGCTACTTGGCGGCTATAAAGATTTTAAAGATTTACCCTACGCCt                                                                       >  2:368916/1‑139 (MQ=255)                                                                             tACTTGTCCGGCATATCGCCATATTGGCGGTCATGGATATGGCTCTATTTATCCAGGGTCAATTGGt                                                                                                                                 >  1:1002026/1‑67 (MQ=255)                                                                             tACTTGTCCGGCATATCGCCATATTGGCGGTCATGGATATGGCTCTATTTATCCAGGGTCAATTGGt                                                                                                                                 <  2:1002026/67‑1 (MQ=255)                                                                                 tGTCCGGCATATCGCCATATTGGCGGTCATGGATATGGCTCTATTTATCCAGGGTCAATTGGTGCGGTGATTTCTCCGCTACTTGGCGGCTATAAAGATTTTAAAGATTTACCCTACGCCTGCTCTTTATGCACAGCtt                                                     >  2:124395/1‑139 (MQ=255)                                                                                  gTCCGGCATATCGCCATATTGGCGGTCATGGATATGGCTCTATTTATCCAGGGTCAATTGGTGCGGTGATTTCTCCGCTACTTGGCGGCTATAAAGATTTTAAAGATTTACCCTACGCCTGCTCTTTATGca                                                           >  1:1324954/1‑132 (MQ=255)                                                                                  gTCCGGCATATCGCCATATTGGCGGTCATGGATATGGCTCTATTTATCCAGGGTCAATTGGTGCGGTGATTTCTCCGCTACTTGGCGGCTATAAAGATTTTAAAGATTTACCCTACGCCTGCTCTTTATGCACAGCTtg                                                    >  1:1457150/1‑139 (MQ=255)                                                                                  gTCCGGCATATCGCCATATTGGCGGTCATGGATATGGCTCTATTTATCCAGGGTCAATTGGTGCGGTGATTTCTCCGCTACTTGGCGGCTATAAAGATTTTAAAGATTTACCCTACGCCTGCTCTTTATGCACAGCTtg                                                    >  1:1461665/1‑139 (MQ=255)                                                                                  gTCCGGCATATCGCCATATTGGCGGTCATGGATATGGCTCTAGTTATCCAGGGTCAATTGGTGCGGTGATTTCTCCGCTACTTGGCGGCTATAAAGATTTTAAAGATTTACCCTACGCCTGCTCTTTATGca                                                           <  2:1324954/132‑1 (MQ=255)                                                                                     cGGCATATCGCCATATTGGCGGTCATGGATATGGCTCTATTTATCCAGGGTCAATTGGTGCGGTGATTTCTCCGCTACTTgg                                                                                                          >  2:327581/1‑82 (MQ=255)                                                                                     cGGCATATCGCCATATTGGCGGTCATGGATATGGCTCTATTTATCCAGGGTCAATTGGTGCGGTGATTTCTCCGCTACTTgg                                                                                                          <  1:327581/82‑1 (MQ=255)                                                                                                cATATTGGCGGTCATGGATATGGCTCTATTTATCCAGGGTCAATTGGTGCGGTGATTTCTCCGCTACTTGGCGGCTATAAAGATTTTAAAGATTTACCCTACGCCTGCTCTTTATGCACAGCTTGTGACAACGTGTGTc                                      <  2:1635254/139‑1 (MQ=255)                                                                                                 atgttGGCGTTCATGGATATGGCTTTCTTTATACAGGGTCAATTGGTGCGGTGATTTATCCGCTACTTGGCGGCTATAAAGATTTTAAAGATTTACCCTACGCCTGCTCTTTATGCACAGCTTGTGACAACGTGTGTcc                                     <  2:150582/136‑1 (MQ=255)                                                                                                 atatTGGCGGTCATGGATATGGCTCTATTTATCCAGGGTCAATTGGTGCGGTGATTTCTCCGCTACTTGGCGGCTATAAAGATTTTAAAGATTTACCCTACGCCTGCTCTTTATGCACAGCTTGTGACAACGTGTGTcc                                     <  2:117407/139‑1 (MQ=255)                                                                                                       gCGGTCATGGATATGGCTCTATTTATCCAGGGTCAATTGGTGCGGTGATTTCTCCGCTACTTGGCGGCTATAAAGATTTTAAAGATTTGCCCTACGCCTGCTCTTTATGCACAGCTTGTGACAACGTGTGTCCGGTGCg                               >  2:596272/1‑139 (MQ=255)                                                                                                            cATGGATATGGCTCTATTTATCCAGGGTCAATTGGTGCGGTGATTTCTCCGCTACTTggc                                                                                                         <  1:1394819/60‑1 (MQ=255)                                                                                                            cATGGATATGGCTCTATTTATCCAGGGTCAATTGGTGCGGTGATTTCTCCGCTACTTggc                                                                                                         >  2:1394819/1‑60 (MQ=255)                                                                                                            cATGGATATGGCTCTATTTATCCAGGGTCAATTGGTGCGGTGATTTCTCCGCTACTTGGCGGCTATAAAGATTTTAAAGATTTACCCTACGCCTGCTCTTTATGCACAGCTTGTGACAACGTGTGTCCGGTGCGTATTc                          >  2:530966/1‑139 (MQ=255)                                                                                                             aTGGATATGGCTCTATTTATCCAGGGTCAATTGGTGCGGTGATTTCTCCGCTACTTGGCGGCTATAAAGATTTTAAAGATTTACCCTACGCCTGCTCTTTATGCACAGCTTGTGACAACGTGTGTCCGGTGCGTATTcc                         <  2:1461665/139‑1 (MQ=255)                                                                                                              tGGATATGGCTCTATTTATCCAGGGTCAATTGGTGCGGt                                                                                                                            <  2:1213595/39‑1 (MQ=255)                                                                                                              tGGATATGGCTCTATTTATCCAGGGTCAATTGGTGCGGt                                                                                                                            >  1:1213595/1‑39 (MQ=255)                                                                                                              tGGATATGGCTCTATTTATCCAGGGTCAATTGGTGCGGTGATTTCTCCGCTACTTgg                                                                                                          >  2:1285341/1‑57 (MQ=255)                                                                                                              tGGATATGGCTCTATTTATCCAGGGTCAATTGGTGCGGTGATTTCTCCGCTACTTgg                                                                                                          <  1:1285341/57‑1 (MQ=255)                                                                                                              tGGATATGGCTCTATTTATCCAGGGTCAATTGGTGCGGTGATTTCTCCGCTACTTGGCGGCTAtaaagatttta                                                                                         <  2:1097684/74‑1 (MQ=255)                                                                                                              tGGATATGGCTCTATTTATCCAGGGTCAATTGGTGCGGTGATTTCTCCGCTACTTGGCGGCTAtaaagatttta                                                                                         >  1:1097684/1‑74 (MQ=255)                                                                                                              tGGATATGGCTCTATTTATCCAGGGTCAATTGGTGCGGTGATTTCTCCGCTACTTGGCGGCTAtaaagatttta                                                                                         <  2:534977/74‑1 (MQ=255)                                                                                                              tGGATATGGCTCTATTTATCCAGGGTCAATTGGTGCGGTGATTTCTCCGCTACTTGGCGGCTAtaaagatttta                                                                                         >  1:534977/1‑74 (MQ=255)                                                                                                                     ggCTCTATTTATCCAGGGTCAATTGGTGCGGTGATTTCTCCGCTACTTGGCGGCTAtaaa                                                                                                <  2:1387622/60‑1 (MQ=255)                                                                                                                     ggCTCTATTTATCCAGGGTCAATTGGTGCGGTGATTTCTCCGCTACTTGGCGGCTAtaaa                                                                                                >  1:1387622/1‑60 (MQ=255)                                                                                                                     ggCTCTATTTATCCAGGGTCAATTGGTGCGGTGATTTCTCCGCTACTTGGCGGCTATAAAGATTTTAAAGATTTACCCTACGCCTGCTCTTTATGCACAGCTTGTGACAACGTGTGTCCGGTGCGTATTCCGCTGTCaa                 >  2:1656466/1‑139 (MQ=255)                                                                                                                            tttATCCAGGGTCAATTGGTGCGGTGATTTCTCCGCTACTTGGCGGCTAtaaagatttta                                                                                         >  1:301720/1‑60 (MQ=255)                                                                                                                            tttATCCAGGGTCAATTGGTGCGGTGATTTCTCCGCTACTTGGCGGCTAtaaagatttta                                                                                         <  2:301720/60‑1 (MQ=255)                                                                                                                            tttATCCAGGGTCAATTGGTGCGGTGATTTCTCCGCTACTTGGCGGCTATAAAGATTTTAAAGATTTACCCTACGCCTGCTCTTTATGCACAGCTTGTGACAACGTGTGTCCGGTGCGTATTCCGCTGTCAAAACTGAt          >  1:1035668/1‑139 (MQ=255)                                                                                                                                 ccAGGGTCAATTGGTGCGGTGATTTCTCCGCTACTTGGCGGCTATAAAGATTTTAAAGATTTACCCTACGCCTGCTCTTTATGCACAGCTTGTGACAACGTGTGTCCGGTGCGTATTCCGCTGTCAAAACTGATTTTGc     >  2:1010981/1‑139 (MQ=255)                                                                                                                                    gggTCAATTGGTGCGGTGATTTCTCCGCTACTTGGCGGCTATAAAGATTTTAAAGATTTACCCTACGCCTGCTCTTTATGCACAGCTTGTGACAACGTGTGTCCGGTGCGTATTCCGCTGTc                   <  1:474381/122‑1 (MQ=255)                                                                                                                                    gggTCAATTGGTGCGGTGATTTCTCCGCTACTTGGCGGCTATAAAGATTTTAAAGATTTACCCTACGCCTGCTCTTTATGCACAGCTTGTGACAACGTGTGTCCGGTGCGTATTCCGCTGTCAAAACTGATTTTGCGTc  >  1:969138/1‑139 (MQ=255)                                                                                                                                       |                                                                                                                                        ATAACGGGCGTTCTGAGGTGCTGGCCTCTGAATTTCGGGATGTGCTGCGCTGTATTCGCTGCGGGGCTTGTATGAATACTTGTCCGGCATATCGCCATATTGGCGGTCATGGATATGGCTCTATTTATCCAGGGCCAATTGGTGCGGTGATTTCTCCGCTACTTGGCGGCTATAAAGATTTTAAAGATTTACCCTACGCCTGCTCTTTATGCACAGCTTGTGACAACGTGTGTCCGGTGCGTATTCCGCTGTCAAAACTGATTTTGCGTC  >  NC_000913/323221‑323490

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