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%	133.6 / NA	42	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 (21/21);  total (21/21)

TAACGGGCGTTCTGAGGTGCTGGCCTCTGAATTTCGGGATGTGCTGCGCTGTATTCGCTGCGGGGCTTGTATGAATACTTGTCCGGCATATCGCCATATTGGCGGTCATGGATATGGCTCTATTTATCCAGGGCCAATTGGTGCGGTGATTTCTCCGCTACTTGGCGGCTATAAAGATTTTAAAGATTTACCCTACGCCTGCTCTTTATGCACAGCTTGTGACAACGTGTGTCCGGTGCGTATTCCGCTGTCAAAACTGATTTTGC  >  NC_000913/323222‑323487                                                                                                                                      |                                                                                                                                     tAACGGGCGTTCTGAGGTGCTGGCCTCTGAATTTCGGGATGTGCTGCGCTGTATTCGCTGCGGGGCTTGTATGAATACTTGTCCGGCATATCGCCATATTGGCGGTCATGGATATGGCTCTATTTATCCAGGGTCAAtt                                                                                                                                 >  1:1193115/1‑139 (MQ=255)            cTGAGGTGCTGGCCTCTGAATTTCGGGATGTGCTGCGCTGTATTCGCTGCGGGGCTTGTATGAATACTTGTCCGGCATATCGCCATATTGGCGGTCATGGATATGGCTCTATTTATCCAGGGTCAATTGGTGCGGTGAt                                                                                                                      <  1:579893/139‑1 (MQ=255)                                     ggATGTGCTGCGCTGTATTCGCTGCGGGGCTTGTATGAATACTTGTCCGGCATATCGCCATATTGGCGGTCATGGATATGGCTCTATTTATCCAGGGTCAATTGGTGCGGTGATTTCTCCGCTACTTGGCGGCTAtaaa                                                                                             <  2:1168184/139‑1 (MQ=255)                                                  tGTATTCGCTGCGGGGCTTGTATGAATACTTGTCCGGCATATCGCCATATTGGCGGTCATGGATATGGCTCTATTTATCCAGGATCAATTGGTGCGGTGATTTCTCCGCTACTTGGCGGCTAtaaagattttaaagatt                                                                                <  1:341818/139‑1 (MQ=255)                                                    tATTCGCTGCGGGGCTTGTATGAATACTTGTCCGGCATATCGCCATATTGGCGGTCATGGATATGGCTCTATTTATCCAGGGTCAATTGGTGCGGTGATTTCTCCGCTACTTGGCGGCTATAAAGATTTTAAAGATTTa                                                                              >  1:413460/1‑139 (MQ=255)                                                       tCGCTGCGGGGCTTGTATGAATACTTGTCCGGCATATCGCCATATTGGCGGTCATGGATATGGCTCTATTTATCCAGGGTCAATTGGTGCGGTGATTTCTCCGCTACTTGGCGGCTATAAAGATTTTAAAGATTTAccc                                                                           <  1:7311/139‑1 (MQ=255)                                                           tGCGGGGCTTGTATGAATACTTGTCCGGCATATCGCCATATTGGCGGTCATGGATATGGCTCTATTTATCCAGGGTCAATTGGTGCGGTGATTTCTCCGCTACTTGGCGGCTATAAAGATTTTAAAGATTTACCCTACg                                                                       <  2:413460/139‑1 (MQ=255)                                                              ggggCTTGTATGAATACTTGTCCGGCATATCGCCATATTGGCGGTCATGGATATGGCTCTATTTATCCAGGGTCAATTGGTGCGGTGATTTCTCCGCTACTTGGCGGCTAtaaagattttaaagatt                                                                                <  2:157600/127‑1 (MQ=255)                                                              ggggCTTGTATGAATACTTGTCCGGCATATCGCCATATTGGCGGTCATGGATATGGCTCTATTTATCCAGGGTCAATTGGTGCGGTGATTTCTCCGCTACTTGGCGGCTAtaaagattttaaagatt                                                                                >  1:157600/1‑127 (MQ=255)                                                                     gTATGAATACTTGTCCGGCATATCGCCATATTGGCGGTCATGGATATGGCTCTATTTATCCAGGGTCAATTGGTGCGGTGATTTCTCCGCTACTTggc                                                                                                      <  2:1041621/98‑1 (MQ=255)                                                                     gTATGAATACTTGTCCGGCATATCGCCATATTGGCGGTCATGGATATGGCTCTATTTATCCAGGGTCAATTGGTGCGGTGATTTCTCCGCTACTTggc                                                                                                      >  1:1041621/1‑98 (MQ=255)                                                                        tGAATACTTGTCCGGCATATCGCCATATTGGCGGTCATGGATATGGCTCTATTTATCCAGGGTCAATTGGTGCGGTGATTTCTCCGCTACTTGGCGGCTATAAAGATTTTAAAGATTTACCCTACGCCTGCTCTTTATg                                                          <  2:1193115/139‑1 (MQ=255)                                                                               ttGTCCGGCATATCGCCATATTGGCGGTCATGGATATGGCTCTATTTATCCAGGGTCAATTGGTGCGGTGATTTCTCCGCTACTTGGCGGCTATAAAGATTTTAAAGATTTACCCTACGCCTGCTCTTTATGCACAGCt                                                   >  1:350584/1‑139 (MQ=255)                                                                                 gTCCGGCATATCGCCATATTGGCGGTCATGGATATGGCTCTATTTATCCAGGGTCAATTGGTGCGGTGATTTCTCCGCTACTTGGCGGCTAtaaagatttta                                                                                      >  2:508191/1‑102 (MQ=255)                                                                                 gTCCGGCATATCGCCATATTGGCGGTCATGGATATGGCTCTATTTATCCAGGGTCAATTGGTGCGGTGATTTCTCCGCTACTTGGCGGCTAtaaagatttta                                                                                      <  1:508191/102‑1 (MQ=255)                                                                                 gTCCGGCATATCGCCATATTGGCGGTCATGGATATGGCTCTATTTATCCAGGGTCAATTGGTGCGGTGATTTCTCCGCTACTTGGCGGCTATAAAGATTTTAAAGATTTACCCTACGCCTGCTCTTTATGCACAGCTtg                                                 >  2:524418/1‑139 (MQ=255)                                                                                      gCATATCGCCATATTGGCGGTCATGGATATGGCTCTATTTATCCAGGGTCAATTGGTGCGGTGATTTCTCCGCTACTTggc                                                                                                      <  2:1121495/81‑1 (MQ=255)                                                                                      gCATATCGCCATATTGGCGGTCATGGATATGGCTCTATTTATCCAGGGTCAATTGGTGCGGTGATTTCTCCGCTACTTggc                                                                                                      >  1:1121495/1‑81 (MQ=255)                                                                                      gCATATCGCCATATTGGCGGTCATGGATATGGCTCTATTTATCCAGGGTCAATTGGTGCGGTGATTTCTCCGCTACTTGGCGGCTATAAAGATTTTAAAGATTTACCCTACGCCTGCTCTTTATGCACAGCTTGTGACa                                            >  2:1014534/1‑139 (MQ=255)                                                                                            cGCCATATTGGCGGTCATGGATATGGCTCTATTTATCCAGGGTCAATTGGTGCGGTGATTTCTCCg                                                                                                               >  2:336003/1‑66 (MQ=255)                                                                                            cGCCATATTGGCGGTCATGGATATGGCTCTATTTATCCAGGGTCAATTGGTGCGGTGATTTCTCCg                                                                                                               <  1:336003/66‑1 (MQ=255)                                                                                                  atTGGCGGTCATGGATATGGCTCTATTTATCCAGGGTCAATTGGTGCGGTGATTTCTCCGCTACTTGGCGGCTATAAAGATTTTAAAGATTTACCCTACGCCTGCTCTTTATGCACAGCTTGTGACAACGTGTGTCCgg                                <  2:350584/139‑1 (MQ=255)                                                                                                      gCGGTCATGGATATGGCTCTATTTATCCAGGGTCAATTGGTGCGGTGATTTCTCCGCTACTTGGCGGCTATAAAGATTTTAAAGATTTAccc                                                                           <  1:1079653/92‑1 (MQ=255)                                                                                                      gCGGTCATGGATATGGCTCTATTTATCCAGGGTCAATTGGTGCGGTGATTTCTCCGCTACTTGGCGGCTATAAAGATTTTAAAGATTTAccc                                                                           >  2:1079653/1‑92 (MQ=255)                                                                                                                    ggCTCTATTTATCCAGGGTCAATTGGTGCGGTGATTTCTCCGCTACTTGGCGGCtata                                                                                               <  1:420959/58‑1 (MQ=255)                                                                                                                    ggCTCTATTTATCCAGGGTCAATTGGTGCGGTGATTTCTCCGCTACTTGGCGGCtata                                                                                               >  2:420959/1‑58 (MQ=255)                                                                                                                    ggCTCTATTTATCCAGGGTCAATTGGTGCGGTGATTTCTCCGCTACTTGGCGGCTAtaaagattttaaagatt                                                                                <  1:857336/73‑1 (MQ=255)                                                                                                                    ggCTCTATTTATCCAGGGTCAATTGGTGCGGTGATTTCTCCGCTACTTGGCGGCTAtaaagattttaaagatt                                                                                >  2:857336/1‑73 (MQ=255)                                                                                                                    ggCTCTATTTATCCAGGGTCAATTGGTGCGGTGATTTCTCCGCTACTTGGCGGCTAtaaagatttta                                                                                      >  2:280715/1‑67 (MQ=255)                                                                                                                    ggCTCTATTTATCCAGGGTCAATTGGTGCGGTGATTTCTCCGCTACTTGGCGGCTAtaaagatttta                                                                                      <  1:280715/67‑1 (MQ=255)                                                                                                                    ggCTCTATTTATCCAGGGTCAATTGGTGCGGTGATTTCTCCGCTACTTGGCGGCTATAAAGATTTTAAAGATTTACCCTACGCCTGc                                                                  >  1:373353/1‑87 (MQ=255)                                                                                                                    ggCTCTATTTATCCAGGGTCAATTGGTGCGGTGATTTCTCCGCTACTTGGCGGCTATAAAGATTTTAAAGATTTACCCTACGCCTGc                                                                  <  2:373353/87‑1 (MQ=255)                                                                                                                    ggCTCTATTTATCCAGGGTCAATTGGTGCGGTGATTTCTCCGCTACTTGGCGGCTATAAAGATTTTAAAGATTTACCCTACGCCTGCTCTTTATGCACAGCTTGTGACAACGTGTGTCCGGTGCGTATTCCGCTGTCaa              >  2:142845/1‑139 (MQ=255)                                                                                                                      ctctATTTATCCAGGGTCAATTGGTGCGGTGATTtctc                                                                                                                 <  2:1118247/38‑1 (MQ=255)                                                                                                                      ctctATTTATCCAGGGTCAATTGGTGCGGTGATTtctc                                                                                                                 >  1:1118247/1‑38 (MQ=255)                                                                                                                         tATTTATCCAGGGTCAATTGGTGCGGTGATTTCTCCGCTACTTGGCGGCTATAAAGATTTTAAAGATTTACCCTACGCCTGCTCTTTATGCACAGCTTGTGACAACGTGTGTCCGGTGCGTATTCCGCTGTCAAAACTg         >  1:654027/1‑139 (MQ=255)                                                                                                                               tCCAGGGTCAATTGGTGCGGTGATTTCTCCGCTACTTgg                                                                                                       <  2:322394/39‑1 (MQ=255)                                                                                                                               tCCAGGGTCAATTGGTGCGGTGATTTCTCCGCTACTTgg                                                                                                       >  1:322394/1‑39 (MQ=255)                                                                                                                                ccAGGGTCAATTGGTGCGGTGATTTCTCCGCTACTTGGCGGCTATAAAGATTTTAAAGATTTACCCTACGCCTGCTCTTTATGCACAGc                                                    <  1:628417/89‑1 (MQ=255)                                                                                                                                ccAGGGTCAATTGGTGCGGTGATTTCTCCGCTACTTGGCGGCTATAAAGATTTTAAAGATTTACCCTACGCCTGCTCTTTATGCACAGc                                                    >  2:628417/1‑89 (MQ=255)                                                                                                                                ccAGGGTCAATTGGTGCGGTGATTTCTCCGCTACTTGGCGGCTATAAAGATTTTAAAGATTTACCCTACGCCTGCTCTTTATGCACAGCTTGTGACAACGTGTGTCCGGTGCGTATTCCGCTGTCAAAACTGATTTTGc  >  2:233539/1‑139 (MQ=255)                                                                                                                                   gggTCAATTGGTGCGGTGATTTCTCCGCTACTTGGCGGCTAtaaagatttta                                                                                      <  2:943895/52‑1 (MQ=255)                                                                                                                                      |                                                                                                                                     TAACGGGCGTTCTGAGGTGCTGGCCTCTGAATTTCGGGATGTGCTGCGCTGTATTCGCTGCGGGGCTTGTATGAATACTTGTCCGGCATATCGCCATATTGGCGGTCATGGATATGGCTCTATTTATCCAGGGCCAATTGGTGCGGTGATTTCTCCGCTACTTGGCGGCTATAAAGATTTTAAAGATTTACCCTACGCCTGCTCTTTATGCACAGCTTGTGACAACGTGTGTCCGGTGCGTATTCCGCTGTCAAAACTGATTTTGC  >  NC_000913/323222‑323487

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