Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	338,125	A→C	I450L (ATC→CTC)	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,125	0	A	C	100.0%	102.7 / NA	29	I450L (ATC→CTC)	yahF	putative acyl‑CoA synthetase YahF
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base C (9/20);  total (9/20)

CCCGATGATTGACCCCACCAACCGCATCAGTCGCTTGATCGAAGAGGCGCGCGATCCAGAAGTGGCGGTGATCGTGA‑TGGATTTTGTGCTCGGATTTGGATCGCATGAAGATCCGGTCGGCTCCACCATCGAGACGAT‑CAAAGAAGCGAAAGCGATCGCCGC‑TGCCGAAGGACGCGAGTTGATCATTCTCGCCTATGTGCTGGGTACCGATCTTGATACGCCATCGTTAGAACAACAAAGCCAGATGCTGCTTGATGCCGGAGTGATTC  >  NC_000913/337989‑338257                                                                                                                                          |                                                                                                                                       cccGATGATTGACCCCACCAACCGCATCAGTCGCTTGATCGAAGAGGCGCGCGATCCAGAAGTGGCGGTGATCGTGA‑TGGATTTTGTGCTCGGATTTGGATCGCATGAAGATCCGGTCGGCTCCACCATCGAGACGCT‑Ca                                                                                                                                    <  1:852333/140‑1 (MQ=255)                caccaACCGCATCAGTCGCTTGATCGAAGAGGCGCGCGATCCAGAAGTGGCGGTGATCGTGA‑TGGATTTTGTGCTCGGATTTGGATCGCATGAAGATCCGGTCGGCTCCACCATCGAGACGCT‑CAAAGaagcgaaagcga                                                                                                                     <  2:302569/140‑1 (MQ=255)                            cAGTCGCTTGATCGAAGAGGCGCGCGATCCAGAAGTGGCGGTGATCGTGA‑TGGATTTTGTGCTCGGATTTGGATCGCATGAAGATCCGGTCGGCTCCACCATCGAGACGCT‑CAAAGAAGCGAAAGCGATCGCCGC‑TGCCg                                                                                                        <  1:671202/140‑1 (MQ=255)                                 gtttGATCGAAGAGGCGCGCGATCCAGAAGTGGCGGTGATCGTGA‑TGGATTTTGTGCTCGGATTTGGATCGCATGAAGATCCGGTCTGCTCCACCATCGAGACGCT‑CAAAGAAGCGAAAGCGATCGCCGC‑TGCCGAAGGa                                                                                                   <  1:722963/138‑1 (MQ=255)                                 gCTTGATCGAAGAGGCGCGCGATCCAGAAGTGGCGGTGATCGTGA‑TGGGTTTTTTGGTCGGATTTTGATTGCATGAAAGTCCGGGCGGCCCCCCCATCCAGACGGT‑CAAAGAAACGAAAAACATCGCCCC‑TGCCggaaga                                                                                                   >  2:72024/1‑135 (MQ=255)                                 gCTTGATCGAAGAGGCGCGCGATCCAGAAGTGGCGGAGTTCGGGATTGGATTTT‑TGCTCGGGGTTGGATTGCAAGCAGATCCGGCCGGCCTCACCCTTGCGGCGCG‑CAAAAAAACGAAAGGGAACGCCCCGTCCCGAAgg                                                                                                    >  2:342962/1‑140 (MQ=255)                                      aTCGAAGAGGCGCGCGATCCAGAAGTGGCGGTGATCGTGA‑TGGATTTTGTGCTCGGATTTGGATCGCATGAAGATCCGGTCGGCTCCACCATCGAGACGCT‑CAAAGAAGCGAAAGCGATCGCCGC‑TGCCGAAGGACGCGa                                                                                              <  1:975139/140‑1 (MQ=255)                                           agagGCGCGCGATCCAGAAGTGGCGGTGATCGTGA‑TGGATTTTGTGCTCGGATTTGGATCGCATGAAGATCCGGTCGGCTCCACCATCGAGACGCT‑CAAAGAAGCGAAAGCGATCGCCGC‑TGCCGAAGGACGCGAGTTGa                                                                                         <  2:404958/140‑1 (MQ=255)                                               gcgcgcGATCCAGAAGTGGCGGTGATCGTGA‑TGGGGTTTTTGCTCGGGTTTGGGTTGCATGAAGAGCCGGGCGGGTCCCCCACCCAGACACTCCAAAAAAGGGAAA‑CGAACGCCCC‑CGCCCCAGGACGGGGGGTGGTCAt                                                                                     >  2:970989/1‑140 (MQ=255)                                                cgcgcgATCCAGAAGTGGCGGTGATCGTGA‑TGGATTTTGTGCTCGGATTTGGATCGCATGAAGATCCGGTCGGCTCCACCATCGAGACGCT‑CAAAGAAGCGAAAGCGATCGCCGC‑TGCCGAAGGACGCGAGTTGATCAtt                                                                                    <  1:409193/140‑1 (MQ=255)                                                      aTCCAGAAGTGGCGGTGATCGTGA‑TGGATTTTGTGCTCGGATTTGGATCGCATGAAGATCCGGTCGGCTCCACCATCGAGACGCT‑CAAAGAAGCGAAAGCGATCGCCGC‑TGCCGAAGGACGCGAGTTGATCATTCTCGcc                                                                              <  1:90893/140‑1 (MQ=255)                                                       tCCAGAAGTGGCGGTGATCGTGA‑TGGATTTTGTGCTCGGATTTGGATCGCATGAAGATCCGGTCGGCTCCACCATCGAGACGCT‑CAAAGAAGCGAAAGCGATCGCCGC‑TGCCGAAGGACGCGAGTTGATCATTCTCGCCt                                                                             <  1:639753/140‑1 (MQ=255)                                                                 gCGGTGATCGTGA‑TGGGTTTTTTGCTCGGGTTTTGATCGCATGAAGATCCGGGCGGCTCCACCCTCGAGACGCT‑CAAAAAAGCGAAAGCGATCGCCGC‑TGCCGAAAGACGCGAGGTGATCATTCTCGCCTATGTGCTggg                                                                   >  2:424134/1‑140 (MQ=255)                                                                   ggTGATCGTGA‑TGGATTTTGTGCTCGGATTTGGATCGCATGAAGATCCGGTCGGCTCCACCATCGAGACGCT‑CAAAGAAGCGAAAGCGATCGCCGC‑TGCCGAAGGACGCGAGTTGATCATTCTCGCCCATGTGCTGGGTa                                                                 <  2:56718/140‑1 (MQ=255)                                                                    gTGATCGTGA‑TGGATTTTGTGCTCGGATTTGGATCGCATGAAGATCCGGTCGGCTCCGCCATCGAGACGCT‑CAAAGAAGCGAAAGCGATCGCCGC‑TGCCGAAGGACGCGAGTTGATCATTCTCGCCTATGTGCTGGGTAc                                                                <  1:791520/140‑1 (MQ=255)                                                                        tCGTGA‑TGGATTTTGTGCTCGGATTTGGATCGCATGAAGATCCGGTCGGCTCCACCATCGAGACGCT‑CAAAGAAGCGAAAGCGATCGCCGC‑TGCCGAAGGACGCGAGTTGATCATTCTCGCCTATGTGCTGGGTACCGAt                                                            <  1:179744/140‑1 (MQ=255)                                                                                    tttGTGCTCGGATTTGGATCGCATGAAGATCCGGTCGGCTCCACCATCGAGACGCT‑CAAAGAAGCGAAAGCGATCGCCGC‑TGCCGAAGGACGCGAGTTGATCATTCTCGCCTATGTGCTGGGTACCGATCTTGATACGcg                                                 <  1:970989/140‑2 (MQ=255)                                                                                    tttGTGCTCGGATTTGGATCGCATGAAGATCCGGTCGGCTCCACCATCGAGACGCT‑CAAAGAAGCGAAAGCGATCGCCGC‑TGCCGAAGGACGCGAGTTGATCATTCTCGCCTATGTGCTGGGTACCGATCTTGATACGcc                                                 <  1:727525/140‑1 (MQ=255)                                                                                    tttGTGCTCGGATTTGGATCGCATGAAGATCCGGTCGGCTCCACCATCGAGACGCA‑CAAAGAAGCGAAAGCGATCGCCGC‑TGCCGAAGGACGCGAGTTGATCATTCTCGCCTATGTGCTGGGTACCGATCTTGATACGcc                                                 <  1:732352/140‑1 (MQ=255)                                                                                          cTCGGATTTGGATCGCATGAAGATCCGGTCGGCTCCACCATCGAGACGCT‑CAAAGAAGCGAAAGCGATCGCCGC‑TGCCGAAGGACGCGAGTTGATCATTCTCGCCTATGTGCTGGGTACCGATCTTGATACGCCATCGtt                                           >  2:709652/1‑140 (MQ=255)                                                                                          cTCGGATTTGGATCGCATGAAGATCCGGTCGGCTCCACCATCGAGACGCT‑CAAAGAAGCGAAAGCGATCGCCGC‑TGCCGAAGGACGCGAGTTGATCATTCTCGCCTATGTGCTGGGTACCGATCTTGATACGCCATCGtt                                           >  2:951584/1‑140 (MQ=255)                                                                                                  tGGATCGCATGAAGATCCGGTCGGCTCCACCATCGAGACGCT‑CAAAGAAGCGAAAGCGATCGCCGC‑TGCCGAAGGACGCGAGTTGATCATTCTCGCCTATGTGCTGGGTACCGATCTTGATACGCCATCGTTAGaacaac                                   <  1:612738/140‑1 (MQ=255)                                                                                                     aTCGCATGAAGATCCGGTCGGCTCCACCATCGAGACGCT‑CAAAGAAGCGAAAGCGATCGCCGC‑TGCCGAAGGACGCGAGTTGATCATTCTCGCCTATGTGCTGGGTACCGATCTTGATACGCCATCGTTAGAACAACaaa                                <  1:96180/140‑1 (MQ=255)                                                                                                         cATGAAGATCCGGTCGGCTCCACCATCGAGACGCT‑CAAAGAAGCGAAAGCGATCGCCGC‑TGCCGAAGGACGCGAGTTGATCATTCTCGCCTATGTGCTGGGTACCGATCTTGATACGCCATCGTTAGAACAACAAAGCCa                            >  2:582242/1‑140 (MQ=255)                                                                                                         cATGAAGATCCGGTCGGCTCCACCATCGAGACGCT‑CAAAGAAGCGAAAGCGATCGCCGC‑TGCCGAAGGACGCGAGTTGATCATTCTCGCCTATGTGCTGGGTACCGATCTTGATACGCCATCGTTAGAACAACAAAGCCa                            >  2:694111/1‑140 (MQ=255)                                                                                                         cATGAAGATCCGGTCGGCTCCACCATCGAGACGCT‑CAAAGAAGCGAAAGCGATCGCCGC‑TGCCGAAGGACGCGAGTTGATCATTCTCGCCTATGTGCTGGGTACCGATCTTGATACGCCATCGTTAGAACAACAAAGCCa                            >  2:861561/1‑140 (MQ=255)                                                                                                         cATGAAGATCCGGTCGGCTCCACCATCGAGACGCT‑CAAAGAAGCGAAAGCGATCGCCGC‑TGCCGAAGGACGCGAGTTGATCATTCTCGCCTATGTGCTGGGTACCGATCTTGATACGCCATCGTTAGAACAACAAAGCCa                            >  1:185363/1‑140 (MQ=255)                                                                                                         cATGAAGATCCGGTCGGCTCCACCATCGAGACGCT‑CAAAGAAGCGAAAGCGATCGCCGC‑TGCCGAAGGACGCGAGTTGATCATTCACGCCTATGTGCTGGGTACCGATCTTGATACGCCATCGTTAGAacaaca                                  <  1:654909/134‑1 (MQ=255)                                                                                                         cATGAAGATCCGGTCGGCTCCACCATCGAGACGCT‑CAAAGAAGCGAAAGCGATCGCCGC‑TGCCGAAGGACGCGAGTTGATCATTCACGCCTATGTGCTGGGTACCGATCTTGATACGCCATCGTTAGAacaaca                                  >  2:654909/1‑134 (MQ=255)                                                                                                           tGAAGATCCGGTCGGCTCCACCATCGAGACGCT‑CAAAGAAGCGAAAGCGATCGCCGC‑TGCCGAAGGACGCGAGTTGATCATTCTCGCCTATGTGCTGGGTACCGATCTTGATACGCCATCGTTAGAACAACAAAGCCAGa                          >  2:128682/1‑140 (MQ=255)                                                                                                              aGATCCGGTCGGCTCCACCATCGAGACGCT‑CAAAGAAGCGAAAGCGATCGCCGC‑TGCCGAAGGACGCGAGTTGATCATTCTCGCCTATGTGCTGGGTACCGATCTTGATACGCCATCGTTAGAACAACAAAGCCAGAtgc                       <  1:582242/140‑1 (MQ=255)                                                                                                                             caccaTCGAGACGCT‑CAAAGAAGCGAAAGCGATCGCCGC‑TGCCGAAGGACGCGAGTTGATCATTCTCGCCTATGTGCTGGGTACCGATCTTGATACGCCATCGTTAGAACAACAAAGCCAGATGCTGCTTGATGCCGGAg        <  1:951584/140‑1 (MQ=255)                                                                                                                                   cGAGACGCT‑CAAAGAAGCGAAAGCGATCGCCGC‑TGCCGAAGGACGCGAGTTGATCATTCTCGCCTATGTGCTGGGTACCGATCTTGATACGCCATCGTTAGAACAACAAAGCCAGATGCTGCTTGATGCCGGAGTGATTc  >  2:711728/1‑140 (MQ=255)                                                                                                                                          |                                                                                                                                       CCCGATGATTGACCCCACCAACCGCATCAGTCGCTTGATCGAAGAGGCGCGCGATCCAGAAGTGGCGGTGATCGTGA‑TGGATTTTGTGCTCGGATTTGGATCGCATGAAGATCCGGTCGGCTCCACCATCGAGACGAT‑CAAAGAAGCGAAAGCGATCGCCGC‑TGCCGAAGGACGCGAGTTGATCATTCTCGCCTATGTGCTGGGTACCGATCTTGATACGCCATCGTTAGAACAACAAAGCCAGATGCTGCTTGATGCCGGAGTGATTC  >  NC_000913/337989‑338257

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