Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,413,180	(A)5→4	coding (153/1164 nt)	yjfC →	putative acid‑‑amine ligase YjfC

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,413,176	0	A	.	100.0%	182.3 / NA	42	coding (149/1164 nt)	yjfC	putative acid‑‑amine ligase YjfC
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base . (28/14);  total (28/14)

CATATCATCGACAATGAAATCTATTGGGATGAGAGTCGGGCTTACCGTTTTACTCTGCGCCAGATTGAAGAGCAGATCGAAAAACCGACTGCGGAACTGCATCAGATGTGCCTTGAGGTGGTGGATCGCGCGGTAAAAGATGAAGAGATCCT  >  NC_000913/4413097‑4413248                                                                                |                                                                         cATATCATCGACAATGAAATCTATTGGGATGAGAGTCGGGCTTACCGTTTTACTCTGCGCCAGATTGAAGAGCAGATCGa                                                                          >  2:889190/1‑80 (MQ=255)     tcatcGACAATGAAATCTATTGGGATGAGAGTCGGGCTTACCGTTTTACTCTGCGCCAGATTGAAGAGCAGATCGaaaacc                                                                     >  1:889189/1‑79 (MQ=255)     tcatcGACAATGAAATCTATTGGGATGAGAGTCGGGCTTACCGTTTTACTCTGCGCCAGATTGAAGAGCAGATCGaaaac                                                                      <  1:889128/80‑2 (MQ=255)     tcatcGACAATGAAATCTATTGGGATGAGAGTCGGGCTTACCGTTTTACTCTGCGCCAGATTGAAGAGCAGATCGaaaac                                                                      <  1:889127/80‑2 (MQ=255)       atcGACAATGAAATCTATTGGGATGAGAGTCGGGCTTACCGTTTTACTCTGCGCCAGATTGAAGAGCAGATCGaaaaccg                                                                    <  1:889129/80‑4 (MQ=255)        tcGACAATGAAATCTATTGGGATGAGAGTCGGGCTTACCGTTTTACTCTGCGCCAGATTGAAGAGCAGATCG‑AAAACCGAc                                                                 <  2:889130/81‑1 (MQ=255)            cAATGAAATCTATTGGGATGAGAGTCGGGCTTACCGTTTTACTCTGCGCCAGATTGAAGAGCAGATCG‑AAAACCGACTGCg                                                             <  1:889131/81‑1 (MQ=255)              aTGAAATCTATTGGGATGAGAGTCGGGCTTACCGTTTTACTCTGCGCCAGATTGAAGAGCAGATCG‑AAAACCGACTGCGGa                                                           <  1:889132/81‑1 (MQ=255)                gAAATCTATTGGGATGAGAGTCGGGCTTACCGTTTTACTCTGCGCCAGATTGAAGAGCAGATCG‑AAAACCGACTGCGGAAc                                                         >  2:889199/1‑81 (MQ=255)                gAAATCTATTGGGATGAGAGTCGGGCTTACCGTTTTACTCTGCGCCAGATTGAAGAGCAGATCG‑AAAACCGACTGCGGAAc                                                         >  1:889184/1‑81 (MQ=255)                     cTATTGGGATGAGAGTCGGGCTTACCGTTTTACTCTGCGCCAGATTGAAGAGCAGATCG‑AAAACCGACTGCGGAACTGCAt                                                    >  1:889188/1‑81 (MQ=255)                         tGGGATGAGAGTCGGGCTTACCGTTTTACTCTGCGCCAGATTGAAGAGCAGATCG‑AAAACCGACTGCGGAACTGCAt                                                    >  2:889207/1‑77 (MQ=255)                            gATGAGAGTCGGGCTTACCGTTTTACTCTGCGCCAGCTTGAAGAGCAGATCG‑AAAACCGACTGCGGAACTGCATCAGAtgt                                             >  2:889267/1‑81 (MQ=255)                               gagaGTCGGGCTTACCGTTTTACTCTGCGCCAGATTGAAGAGCAGATCG‑AAAACCGACTGCGGAACTGCATCAGATGTGcc                                          >  2:889183/1‑81 (MQ=255)                                 gagTCGGGCTTACCGTTTTACTCTGCGCCAGATTGAAGAGCAGATCG‑AAAACCGACTGCGGAACTGCATCAGATGTGCCtt                                        >  1:889182/1‑81 (MQ=255)                                    tCGGTCTTACCGTTTTACTCTGCGCCAGATTGAAGAGCAGATCG‑AAAACCGACTGCGGAACTGCATCAGATGTGCCTTGAg                                     >  1:889195/1‑81 (MQ=255)                                    tCGGGCTTACCGTTTTACTCTGCGCCAGATTGAAGAGCAGATCG‑AAAACCGACTGCGGAACTGCATCAGATGTGCCTTGAg                                     >  2:889224/1‑81 (MQ=255)                                     cGGGCTTACCGTTTTACTCTGCGCCAGATTGAAGAGCAGATCG‑AAAACCGACTGCGGAACTGCATCAGATGTGCCTTGAgg                                    >  1:889208/1‑81 (MQ=255)                                       ggCTTACCGTTTTACTCTGCGCCAGATTGAAGAGCAGATCG‑AAAACCGACTGCGGAACTGCATCAGATGTGCCTTGAggtg                                  >  2:889209/1‑81 (MQ=255)                                         cTTACCGTTTTACTCTGCGCCAGATTGAAGAGCAGATCG‑AAAACCGACTGCGGAACTGCATCAGATGTGCCTTGAggtggt                                <  1:889133/81‑1 (MQ=255)                                         cTTACCGTTTTACTCTGCGCCAGATTGAAGAGCAGATCG‑AAAACCGACTGCGGAACTGCATCAGATGTGCCTTGAggtggt                                >  2:889194/1‑81 (MQ=255)                                          ttACCGTTTTACTCTGCGCCAGATTGAAGAGCAGATCG‑AAAACCGACTGCGGAACTGCATCAGATGTGCCTTGAGgtggtg                               >  2:889212/1‑81 (MQ=255)                                             ccGTTTTACTCTGCGCCAGATTGAAGAGCAGATCG‑AAAACCGACTGCGGAACTGCATCAGATGTGCCTTGAGGTGGTGGAt                            <  1:889134/81‑1 (MQ=255)                                              cGTTTTACTCTGCGCCAGATTGAAGAGCAGATCG‑AAAACCGACTGCGGAACTGCATCAGATGTGCCTTGAGGTGGTGGATc                           <  1:889137/81‑1 (MQ=255)                                              cGTTTTACTCTGCGCCAGATTGAAGAGCAGATCG‑AAAACCGACTGCGGAACTGCATCAGATGTGCCTTGAGGTGGTGGATc                           <  1:889135/81‑1 (MQ=255)                                              cGTTTTACTCTGCGCCAGATTGAAGAGCAGATCG‑AAAACCGACTGCGGAACTGCATCAGATGTGCCTTGAGGTGGTGGATc                           <  1:889136/81‑1 (MQ=255)                                               gTTTTACTCTGCGCCAGATTGAAGAGCAGATCG‑AAAACCGACTGCGGAACTGCATCAGATGTGCCTTGAGGTGGTGGATcg                          >  2:889243/1‑81 (MQ=255)                                               gTTTTACTCTGCGCCAGATTGAAGAGCAGATCG‑AAAACCGACTGCGGAACTGCATCAGATGTGCCTTGAGGTGGTGGATcg                          >  1:889201/1‑81 (MQ=255)                                                     ctctGCGCCAGATTGAAGAGCAGATCG‑AAAACCGACTGCGGAACTGCATCAGATGTGCCTTGAGGTGGTGGATCGCGCgg                     <  1:889139/80‑1 (MQ=255)                                                     ctctGCGCCAGATTGAAGAGCAGATCG‑AAAACCGACTGCGGAACTGCATCAGATGTGCCTTGAGGTGGTGGATCGCGCGGt                    <  1:889138/81‑1 (MQ=255)                                                        tGCGCCAGATTGAAGAGCAGATCG‑AAAACCGACTGCGGAACTGCATCAGATGTGCCTTGAGGTGGTGGATCGCGCGGTa                   <  1:889140/79‑1 (MQ=255)                                                          cgcCAGATTGAAGAGCAGATCG‑AAAACCGACTGCGGAACTGCATCAGATGTGCCTTGAGGTGGTGGATCGCGCGGTaa                  <  2:889142/78‑1 (MQ=255)                                                          cgcCAGATTGAAGAGCAGATCG‑AAAACCGACTGCGGAACTGCATCAGATGTGCCTTGAGGTGGTGGATCGCGCGGTAAAAg               <  2:889141/81‑1 (MQ=255)                                                          cgcCAGATTGAAGAGCAGATCG‑AAAACCGACTGCGGAACTGCATCAGATGTGCCTTGAGGTGGTGGAACGCGCGGTaa                  >  2:889178/1‑78 (MQ=255)                                                             cAGATTGAAGAGCAGATCG‑AAAACCGACTGCGGAACTGCATCAGATGTGCCTTGAGGTGGTGGATCGCGCGGTAAAAGAt             >  2:889203/1‑80 (MQ=255)                                                             cAGATTGAAGAGCAGATCG‑AAAACCGACTGCGGAACTGCATCAGATGTGCCTTGAGGTGGTGGATCGCGCGGTAAAAGAt             >  1:889216/1‑80 (MQ=255)                                                             cAGATTGAAGAGCAGATCG‑AAAACCGACTGCAGAACTGCATCAGATGTGCCTTGAGGTGGTGGATCGCGCGGTAAAAGATg            <  2:889143/81‑1 (MQ=255)                                                              aGATTGAAGAGCAGATCG‑AAAACCGACTGCGGAACTGCATCAGATGTGCCTTGAGGTGGTGGATCGCGCGGTAAAAGATg            >  1:889186/1‑80 (MQ=255)                                                                  tGAAGAGCAGATCG‑AAAACCGACTGCGGAACTGCATCAGATGTGCCTTGAGGTGGTGGATCGCGCGGTAAAAGATGaa          >  2:889221/1‑78 (MQ=255)                                                                  tGAAGAGCAGATCG‑AAAACCGACTGCGGAACTGCATCAGATGTGCCTTGAGGTGGTGGATCGCGCGGTAAAAGATGaa          >  1:889196/1‑78 (MQ=255)                                                                   gAAGAGCAGATCG‑AAAACCGACTGCGGAACTGCATCAGATGTGCCTTGAGGTGGTGGATCGCGCGGTAAAAGATGaa          >  2:889205/1‑77 (MQ=255)                                                                   gAAGAGCAGATCG‑AAAACCGACTGCGGAACTGCATCAGATGTGCCTTGAGGTGGTGGATCGCGCGGTAAAAGATGaa          >  1:889219/1‑77 (MQ=255)                                                                   gAAGAGCAGATCG‑AAAACCGACTGCGGAACTGCATCAGATGTGCCTTGAGGTGGTGGATCGCGCGGTAAAAGATGaa          >  2:889220/1‑77 (MQ=255)                                                                   gAAGAGCAGATCG‑AAAACCGACTGCGGAACTGCATCAGATGTGCCTTGAGGTGGTGGATCGCGCGGTAAAAGATGaa          >  1:889197/1‑77 (MQ=255)                                                                   gAAGAGCAGATCG‑AAAACCGACTGCGGAACTGCATCAGATGTGCCTTGAGGTGGTGGATCGCGCGGTAAAAGATGaa          >  2:889227/1‑77 (MQ=255)                                                                    aaGAGCAGATCG‑AAAACCGACTGCGGAACTGCATCAGATGTGCCTTGAGGTGGTGGATCGCGCGGTAAAAGATGaa          >  2:889230/1‑76 (MQ=255)                                                                       agCAGATCG‑AAAACCGACTGCGGAACTGCATCAGATGTGCCTTGAGGTGGTGGATCGCGCGGTAAAAGATGAAGAGATCCt  >  1:889222/1‑81 (MQ=255)                                                                                |                                                                         CATATCATCGACAATGAAATCTATTGGGATGAGAGTCGGGCTTACCGTTTTACTCTGCGCCAGATTGAAGAGCAGATCGAAAAACCGACTGCGGAACTGCATCAGATGTGCCTTGAGGTGGTGGATCGCGCGGTAAAAGATGAAGAGATCCT  >  NC_000913/4413097‑4413248

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