Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	338,109	C→T	G444G (GGC→GGT)	yahF →	putative acyl‑CoA synthetase YahF

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	338,109	0	C	T	100.0%	87.6 / NA	26	G444G (GGC→GGT)	yahF	putative acyl‑CoA synthetase YahF
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (7/19);  total (7/19)

AGCCGCACCCGATGATTGACCCCACCAACCGCATCAGTCGCTTGATCGAAGAGGCGCGCGATCCAGAAGTGGCGGTGATCGTGATGGATTTTGTGCTCGGATTTGGATCGCATGAAGATCCGGTCGGCTCCACCATCGAGACGATCAAAGAAGCGAAAGCGATCGCCGCTGCCGAAGGACGCGAGTTGATCATTCTCGCCTATGTGCTGGGTACCGATCTTGATACGCCATCGTTAGAACAACAAAGCCAGATGCTGCTTG  >  NC_000913/337982‑338242                                                                                                                                |                                                                                                                                      aGCCGCACCCGATGATTGACCCCACCAACCGCATCAGTCGCTTGATCGAAGAGGCGCGCGATCCAGAAGTGGCGGTGATCGTGATGGGTTTTTTGCTCGGGTTTTGATCGCATGAAGAACCGGTCGGGTCCCCCATCgag                                                                                                                           >  2:416819/1‑140 (MQ=255) aGCCGCACCCGATGATTGACCCCACCAACCGCATCAGTCGCTTGATCGAAGAGGCGCGCGATCCAGAAGTGGCGGTGATCGGGGTGGGTTTTGTGCTCGGATTTGGGTTGCCTCTAAATCCGGTGGGTTTCTCCATcccg                                                                                                                           >  2:65429/1‑137 (MQ=255)                   aCCCCACCAACCGCATCAGTCGCTTGATCGAAGAGGCGCGCGATCCAGAAGTGGCGGTGATCGGGGTGGGTTTTTTGCTCGGGTTTGGATCGCATGAAGAACCGGGCGGTTCCCCCATCGAGACGATCAAAAaagcgaaa                                                                                                         >  2:325317/1‑140 (MQ=255)                          caACCGCATCAGTCGCTTGATCGAAGAGGCGCGCGATCCAGAAGTGGCGGTGATCGTGATGGATTTTGTGCTCGGATTTGGATCGCATGAAGATCCGGTCGGTTCCACCATCGAGACGATCAAAGAAGCGAAAGCGATcg                                                                                                  <  1:172203/140‑1 (MQ=255)                             ccGCATCAGTCGCTTGATCGAAGAGGCGCGCGATCCAGAAGTGGCGGTGATCGTGATGGATTTTGTGCTCGGATTTGGATCGCATGAAGATCCGGTCGGTTCCACCATCGAGACGATCAAAGAAGCGAAAGCGATcgccg                                                                                               <  2:523112/140‑1 (MQ=255)                                    aGTCGCTTGATCGAAGAGGCGCGCGATCCAGAAGTGGCGGTGATCGTGATGGATTTTGTGCTCGGATTTGGATCGCATGAAGATCCGGTCGGTTCCACCATCGAGACGATCAAAGAAGCGAAAGCGATCGCCGCTGCCGa                                                                                        <  2:370859/140‑1 (MQ=255)                                         cTTGATCGAAGAGGCGCGCGATCCAGAAGTGGCGGTGATCGTGATGGATTTTGTGCTCGGATTTGGATCGCATGAAGATCCGGTCGGTTCCACCATCGAGACGATCAAAGAAGCGAAAGCGATCGCCGCTGCCGAAGGAc                                                                                   <  1:157598/140‑1 (MQ=255)                                          ttGATCGAAGAGGCGCGCGATCCAGAAGTGGCGGTGATCGTGATGGATTTTGTGCTCGGATTTGGATCGCATGAAGATCCGGTCGGTTCCACCATCGAGACGATCAAAGAAGCGAAAGCGATCGCCGCTGCCGAAGGAcg                                                                                  <  1:738196/140‑1 (MQ=255)                                               cGAAGAGGCGCGCGATCCAGAAGTGGCGGTGATCGTGATGGATTTTGTGCTCGGATTTGGATCGCATGAAGATCCGGTCGGTTCCACCATCGAGACGATCAAAGAAGCGAAAGCGATCGCCGCTGCCGAAGGACGCGAGt                                                                             <  1:200807/140‑1 (MQ=255)                                                  agagGCGCGCGATCCAGAAGTGGCGGTGATCGTGATGGATTTTGTGCTCGGATTTGGATCGCATGAAGATCCGGTCGGTTCCACCATCGAGACGATCAAAGAAGCGAAAGCGATCGCCGCTGCCGAAGGACGCGAGTTGa                                                                          <  1:65429/140‑1 (MQ=255)                                                      gcgcgcGATCCAGAAGTGGCGGTGATCGTGATGGATTTTGTGCTCGGATTTGGATCGCATGAAGATCCGGTCGGTTCCACCATCGAGACGATCAAAGAAGCGAAAGCGATCGCCGCTGCCGAAGGACGCGAGTTGATCAt                                                                      <  2:8251/140‑1 (MQ=255)                                                       cgcgcgATCCAGAAGTGGCGGTGATCGTGATGGGTTTTTTGCTCGGGTTTGGATCGCATGAAGAACCGGGCGGGTCCACCATCCAGACGATCAAAAAAGCGAAAGCGAACGCCGCCGCCGGAGGACGCGGGGTGGTCAtt                                                                     >  2:524496/1‑140 (MQ=255)                                                           agATCCAGAAGTGGCGGTGATCGTGATGGATTTTGTGCTCGGATTTGGATCGCATGAAGATCCGGTCGGTTCCACCATCGAGACGATCAAAGAAGCGAAAGCGATCGCCGCTGCCGAAGGACGCGAGTTGATCATTCTCg                                                                 <  1:150424/139‑1 (MQ=255)                                                                cAGAAGTGGCGGTGATCGTGATGGATTTTGTGCTCGGATTTGGATCGCATGAAGATCCGGTCGGTTCCACCATCGAGACGATCAAAGAAGCGAAAGCGATCGCCGCTGCCGAAGGACGCGAGTTGATCATTCTCGCCTAt                                                            <  1:436841/140‑1 (MQ=255)                                                                 aGAAGTGGCGGTGATCGTGATGGATTTTGTGCTCGGATTTGGATCGCATGAAGATCCGGTCGGTTCCACCATCGAGACGATCAAAGAAGCGAAAGCGATCGCCGCTGCCGAAGGACGCGAGTTGATCATTCTCGCCTAtg                                                           <  1:289174/140‑1 (MQ=255)                                                                   aaGTGGCGGTGATCGTGATGGATTTTGTGCTCGGATTTGGATCGCATGAAGATCCGGTCGGTTCCACCATCGAGACGATCAAAGAAGCGAAAGCGATCGCCGCTGCCGAAGGACGCGAGTTGATCATTCTCGCCTAtgtg                                                         <  1:325317/140‑1 (MQ=255)                                                                        gCGGTGATCGTGGTGGGTTTTTTGCTCGGGGTTGGGTCGCATGAAGAACCCGGCGGGTCCACCATCGAGACGATTAAAAAAACGAAAGCGATCGCCGCCGCCGAAGGACGCGGGGTGATCATTCTCGCCTATGTGCTggg                                                    >  2:588179/1‑140 (MQ=255)                                                                        gCGGTGATCGTGATGGGTTTTTTGCTCGGGTTTTGATCGCATGAAGATCCGGTCGGTTCCACCATCGAGACGATCAAAAAAGCGAAAACGATCGCCGCCGCCGAAGGACGCGAGGTGATCATTCTCGCCCATTTGCTggg                                                    >  2:284365/1‑140 (MQ=255)                                                                        gCGGTGATCGTGATGGGTTTTTTGCTCGGGTTTTGATCGCATGAAGATCCGGGCGGGTCCACCATCGAGACGATCAAAAAAACGAAAACGAACGCCGCCGCCGAAGGACGCGGGTTGATCATTCTCGCCTATTTGCTggg                                                    >  2:551925/1‑140 (MQ=255)                                                                        gCGGTGATCGTGATGGGTTTTGTGGTCGGGTTTTGGTCGCATGAAGATCCGGGCGGGTCCACCCTCGAGACGATCAAAAAAGCGAAAACGATCGCCGCTGCCGAAGGACGCGGGGTGATCATTTTCGCCTATGTGCTggg                                                    >  2:514654/1‑140 (MQ=255)                                                                        gCGGTGATCGTGATGGATTTTTTGCTCGGGTTTGGATCGCATGAAGATCCGGTCGGTTCCACCATCGAGACGATCAAAAAAACGAAAACGAACGCCGCCGCCGAAGGACGCGGGGTGATCATTCTCTCCTATGTGCCggg                                                    >  2:289174/1‑140 (MQ=255)                                                                                    aTGGATTTTGTGCTCGGATTTGGATCGCATGAAGATCCGGTCGGTTCCACCATCGAGACGATCAAAGAAGCGAAAGCGATCGCCGCTGCCGAAGGACGCGAGTTGATCATTCTCGCCTATGTGCTGGGTACCGATCTTGa                                        <  1:551925/140‑1 (MQ=255)                                                                                        aTTTTGTGCTCGGATTTGGATCGCATGAAGATCCGGTCGGTTCCACCATCGAGACGATCAAAGAAGCGAAAGCGATCGCCGCTGCCGAAGGACGCGAGTTGATCATTCTCGCCTATGTGCTGGGTACCGATCTTGATACg                                    <  1:524496/140‑1 (MQ=255)                                                                                        aTTTTGTGCTCGGATTTGGATCGCATGAAGATCCGGTCGGTTCCACCATCGAGACGATCAAAGAAGCGAAAGCGATCGCCGCTGCCGAAGGACGCGAGTTGATCATTCTCGCCTATGTGCTGGGTACCGATCTTGATACg                                    <  1:284365/140‑1 (MQ=255)                                                                                          tttGTGCTCGGATTTGGATCGCATGAAGATCCGGTCGGTTCCACCATCGAGACGATCAAAGAAGCGAAAGCGATCGCCGCTGCCGAAGGACGCGAGTTGATCATTCTCGCCTATGTGCTGGGTACCGATCTTGATACGcc                                  <  1:416819/140‑1 (MQ=255)                                                                                             gtgCTCGGATTTGGATCGCATGAAGATCCGGTCGGTTCCACCATCGAGACGATCAAAGAAGCGAAAGCGATCGCCGCTGCCGAAGGACGCGAGTTGATCATTCTCGCCTATGTGCTGGGTACCGATCTTGATACGCCATc                               <  1:108779/140‑1 (MQ=255)                                                                                                   ggATTTGGATCGCATGAAGATCCGGTCGGTTCCACCATCGAGACGATCAAAGAAGCGAAAGCGATCGCCGCTGCCGAAGGACGCGAGTTGATCATTCTCGCCTATGTGCTGGGTACCGATCTTGATACGcc                                  >  2:511339/1‑131 (MQ=255)                                                                                                   ggATTTGGATCGCATGAAGATCCGGTCGGTTCCACCATCGAGACGATCAAAGAAGCGAAAGCGATCGCCGCTGCCGAAGGACGCGAGTTGATCATTCTCGCCTATGTGCTGGGTACCGATCTTGATACGcc                                  <  1:511339/131‑1 (MQ=255)                                                                                                    gATTTGGATCGCATGAAGATCCGGTCGGTTCCACCATCGAGACGATCAAAGAAGCGAAAGCGATCGCCGCTGCCGAAGGACGCGAGTTGATCATTCTCGCCTATGTGCTGGGTACCGATCTTGATACGCCATCGTTAGaa                        >  1:753296/1‑140 (MQ=255)                                                                                                         ggATCGCATGAAGATCCGGTCGGTTCCACCATCGAGACGATCAAAGAAGCGAAAGCGATCGCCGCTGCCGAAGGACGCGAGTTGATCATTCTCGCCTATGTGCTGGGTACCGATCTTGATACGCCATCGTTAGAacaaca                   <  2:627312/140‑1 (MQ=255)                                                                                                           aTCGCATGAAGATCCGGTCGGTTCCACCATCGAGACGATCAAAGAAGCGAAAGCGATCGCCGCTGCCGAAGGACGCGAGTTGATCATTCTCGCCTATGTGCTGGGTACCGATCTTGATACGCCATCGTTAGAACAACaaa                 >  2:397116/1‑140 (MQ=255)                                                                                                               cATGAAGATCCGGTCGGTTCCACCATCGAGACGATCAAAGAAGCGAAAGCGATCGCCGCTGCCGAAGGACGCGAGTTGATCATTCTCGCCTATGTGCTGGGTACCGATCTTGATACGCCATCGTTAGAACAACAAAGCCa             >  1:24767/1‑140 (MQ=255)                                                                                                               cATGAAGATCCGGTCGGTTCCACCATCGAGACGATCAAAGAAGCGAAAGCGATCGCCGCTGCCGAAGGACGCGAGTTGATCATTCTCGCCTATGTGCTGGGTACCGATCTTGATACGCCATCGTTAGAACAACAAAGCCa             >  2:430463/1‑140 (MQ=255)                                                                                                               cATGAAGATCCGGTCGGTTCCACCATCGAGACGATCAAAGAAGCGAAAGCGATCGCCGCTGCCGAAGGACGCGAGTTGATCATTCTCGCCTATGTGCTGGGTACCGATCTTGATACGCCATCGTTAGAACAACAAAGCCa             >  1:363063/1‑140 (MQ=255)                                                                                                                          ggTCGGTTCCACCATCGAGACGATCAAAGAAGCGAAAGCGATCGCCGCTGCCGAAGGACGCGAGTTGATCATTCTCGCCTATGTGCTGGGTACCGATCTTGATACGCCATCGTTAGAACAACAAAGCCAGATGCTGCTTg  >  1:454888/1‑140 (MQ=255)                                                                                                                                |                                                                                                                                      AGCCGCACCCGATGATTGACCCCACCAACCGCATCAGTCGCTTGATCGAAGAGGCGCGCGATCCAGAAGTGGCGGTGATCGTGATGGATTTTGTGCTCGGATTTGGATCGCATGAAGATCCGGTCGGCTCCACCATCGAGACGATCAAAGAAGCGAAAGCGATCGCCGCTGCCGAAGGACGCGAGTTGATCATTCTCGCCTATGTGCTGGGTACCGATCTTGATACGCCATCGTTAGAACAACAAAGCCAGATGCTGCTTG  >  NC_000913/337982‑338242

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